非等向性扩散

✍ dations ◷ 2025-04-04 11:19:35 #Image processing,Image noise reduction techniques

在影像处理及电脑视觉领域中,Anistropic Diffusion(非等向性扩散)是一项用来减少影像噪声但却不会影响到影像中较重要成分的技术,像是边界、线条或者影像中较明显的细节。一般影像扩散处理是将原始影像与二维高斯滤波器进行卷积,这种扩散处理是线性且具有空间不变性的转换。而非等向性扩散处理则是会根据影像产生区域性的滤波器,再将原始影像与产生的滤波器进行卷积,所以非等向性扩散是一种非线性且不具有空间不变性的转换。

Perona和Malik在1987年提出不具有空间不变性的滤波器时,其原始的概念是等向性扩散但会根据影像内容产生不同的滤波器,这也使得在靠近边界的区域其产生的滤波器会很类似狄拉克δ函数,让边界及影像中较重要的结构能够在经过扩散处理后还能保留下来。而当初Perona和Malik称之为非等向性扩散,即使其产生的区域性滤波器是具有等向性的,而当时这种处理又被称为不均匀扩散、非线性扩散及Perona-Malik扩散。而实际上的非等向性扩散则是根据边界及结构的方向而产生非等向性的区域性滤波器,这种方法又被称为shape-adapted smoothing或coherence enhancing diffusion。其产生的影像可以同时进行平滑化并保留原本影像的结构,而这类方法所使用的扩散方程式通常是根据在原始影像中的位置及原始影像的像素值所产生。

虽然其结果是由原始影像及区域性滤波器卷积所产生,但实际应用上这样会需要大量的运算,所以通常会用近似法来进行加速,也就是说每一张新的影像是由上一张产生的影像套用非等向性扩散所产生。整体来说,非等向性扩散是一种迭代性的处理,其产生的结果会越来越平滑直到达到所需要的结果。

Ω R 2 {\displaystyle \Omega \subset \mathbb {R} ^{2}} 代表的是平面上的子集合,且 I ( , t ) : Ω R {\displaystyle I(\cdot ,t):\Omega \rightarrow \mathbb {R} } 是一组灰阶影像,则非等向性扩散可以定义为

Δ {\displaystyle \Delta } 代表的是拉普拉斯运算子, {\displaystyle \nabla } 代表的是梯度运算子, d i v ( ) {\displaystyle \mathrm {div} (\dots )} 则是散度运算子,而 c ( x , y , t ) {\displaystyle c(x,y,t)} 代表的是扩散系数. c ( x , y , t ) {\displaystyle c(x,y,t)} 控制扩散的程度,而且通常是根据影像梯度所产生的方程式,所以能够保存原本影像中的边界。 Pietro Perona 和 Jitendra Malik 在1990年最早提出非等向性扩散的概念,且提出了两种计算扩散系数的方程式:

常数K控制方程式对于边界的敏感度,而其值通常是根据影像中的噪音所产生,或者根据实验所产生。

M {\displaystyle M} 代表的是平滑的影像,则上面的扩散方程式就可以被转换成用梯度下降法寻找方程式 E : M R {\displaystyle E:M\rightarrow \mathbb {R} } 的最小能量,而 E : M R {\displaystyle E:M\rightarrow \mathbb {R} } 则定义为

其中 g : R R {\displaystyle g:\mathbb {R} \rightarrow \mathbb {R} } 是一个实数函数,其代表的是扩散系数之间的关系。对于可微函数 h {\displaystyle h}

假设 E I {\displaystyle \nabla E_{I}} 代表 E 对 L 2 ( Ω , R ) {\displaystyle L^{2}(\Omega ,\mathbb {R} )} 内积的梯度,则

因此,其梯度下降法的方程式可以表示成

我们假设 c = g {\displaystyle c=g'} 就可以得到非等向性方程式了。

修正后的Perona-Malik模型,又被成为正规化的P-M方程式,其未知部分在非线性部分与高斯函数进行卷积,得到

其中 G σ = C σ ( 1 / 2 ) e x p ( | x | 2 / 4 σ ) {\displaystyle G_{\sigma }=C{\sigma }^{-\left(1/2\right)}exp\left(-|x|^{2}/4{\sigma }\right)} .

正规化虽然可以增加其稳定性,但同时也会产生模糊效果,所以要在事前得知噪音的程度才能够决定正规化的所需要的常数。

非等向性扩散可以用来减少数位影像的噪声而不会模糊其边界。如果在固定的扩散系数下,非等向性扩散方程式所减少的heat equation与高斯模糊是相同的,但这样会在消除噪声时同时模糊边界。如果扩散系数是根据边界侦测方程式来决定,像是Perona Malik 模型的话,其结果会在区域内进行扩散而且不会使其超过较强的边界,因此在移除噪声后,影像中的边界及结构仍可以保留下来。

除了移除噪声之外,非等向性扩散也可以用于边界侦测。只要根据边界侦测方程式来进行多次递回的非等向性扩散,其最终结果影像会趋向于剩下一个一个的色块,而相邻色块之间的区域则会被侦测为边界。

相关

  • 情夫情夫是已婚女性的男性婚外情人,尤其指被亿万富婆包养的一类,中国古代称面首,又俗称小白脸,主人翁,牛郎, 现在又有二公、二爷、小王的俗称(相当于情妇被称为“二奶”或当今俗称的“
  • 阳明海洋文化艺术馆阳明海洋文化艺术馆位于台湾基隆市港西街,与基隆火车站及海港大楼相毗邻,是阳明海运创立的一所以海洋文化为主题的博物馆。阳明海洋文化艺术馆是一栋具有尖塔圆拱的建筑物,为基
  • 乌坦战争坦桑尼亚 乌干达全国解放阵线(UNLA)支援:坦桑尼亚:乌干达全国解放阵线:乌干达-坦桑尼亚战争,简称乌坦战争,在乌干达国内一般称为解放战争。这场战争起源于1978年乌干达侵略邻国坦桑
  • 约翰一世 (法兰西)(遗腹子)约翰一世(法语:Jean Ier de France,1316年11月15日-1316年11月20日),名义上的法兰西国王和纳瓦拉国王(称胡安一世)。约翰为卡佩王朝国王路易十世的遗腹子,母亲是匈牙利公主克莱
  • 联合抵制奥林匹克运动会列表奥林匹克运动会是国际性的综合运动会,在历史上它曾被抵制过六次,全都发生在夏季奥运。第一次抵制发生于1936年夏季奥林匹克运动会,最近一次则发生于1988年夏季奥林匹克运动会。
  • 陈冬龙陈冬龙(印尼语:Wilson Tjandinegara,1946年12月20日-),印尼文学社(印尼语:Komunitas Sastra Indonesia)的著名作家和翻译家,印华作协总会理事。1946年出生于望加锡,他的作品发表于《印
  • 浩吉铁路浩吉铁路,原称蒙西至华中地区铁路煤运通道或蒙华铁路,是中华人民共和国一条以煤炭运输为主的铁路。线路自内蒙古自治区鄂尔多斯市乌审旗浩勒报吉南站至江西省吉安市吉安站,全长
  • 阮文存阮文存(越南语:Nguyễn Văn Tồn/.mw-parser-output .han-nom{font-family:"Nom Na Tong","Han-Nom Gothic","Han-Nom Ming","HAN NOM A","HAN NOM B","Ming-Lt-HKSCS-UNI-H",
  • 冯淬冯淬(1940年-),女,汉族,北京人,中华人民共和国政治人物,曾任中华全国妇女联合会书记处书记,第九、十届全国政协委员。
  • 杨抡 (乾隆进士)杨抡,字方叔,号莲趺。清代常州府金匮县(今属无锡市)人。杨潮观之子。雍正十一年(1733年)邑庠生。乾隆三十九年(1774年)举人,乾隆四十三年(1778年)进士。官浙江天台县知县。无子,以杨芳灿