在影像处理 及电脑视觉 领域中，非等向性扩散 （英语：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} }
∂
I
∂
t
=
d
i
v
(
c
(
x
,
y
,
t
)
∇
I
)
=
∇
c
⋅
∇
I
+
c
(
x
,
y
,
t
)
Δ
I
{\displaystyle {\frac {\partial I}{\partial t}}=\mathrm {div} \left(c(x,y,t)\nabla I\right)=\nabla c\cdot \nabla I+c(x,y,t)\Delta I}
Δ
{\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年最早提出非等向性扩散的概念，且提出了两种计算扩散系数的方程式：
c
(
‖
∇
I
‖
)
=
e
−
(
‖
∇
I
‖
/
K
)
2
{\displaystyle c\left(\|\nabla I\|\right)=e^{-\left(\|\nabla I\|/K\right)^{2}}}
和
c
(
‖
∇
I
‖
)
=
1
1
+
(
‖
∇
I
‖
K
)
2
{\displaystyle c\left(\|\nabla I\|\right)={\frac {1}{1+\left({\frac {\|\nabla I\|}{K}}\right)^{2}}}}
常数K控制方程式对于边界的敏感度，而其值通常是根据影像中的噪音所产生，或者根据实验所产生。
若
M
{\displaystyle M}
E
:
M
→
R
{\displaystyle E:M\rightarrow \mathbb {R} }
E
:
M
→
R
{\displaystyle E:M\rightarrow \mathbb {R} }
E
[
I
]
=
1
2
∫
Ω
g
(
‖
∇
I
(
x
)
‖
2
)
d
x
{\displaystyle E[I]={\frac {1}{2}}\int _{\Omega }g\left(\|\nabla I(x)\|^{2}\right)\,dx}
其中
g
:
R
→
R
{\displaystyle g:\mathbb {R} \rightarrow \mathbb {R} }
h
{\displaystyle h}
d
d
t
|
t
=
0
E
[
I
+
t
h
]
=
d
d
t
|
t
=
0
1
2
∫
Ω
g
(
‖
∇
(
I
+
t
h
)
(
x
)
‖
2
)
d
x
=
∫
Ω
g
′
(
‖
∇
I
(
x
)
‖
2
)
∇
I
⋅
∇
h
d
x
=
−
∫
Ω
d
i
v
(
g
′
(
‖
∇
I
(
x
)
‖
2
)
∇
I
)
h
d
x
{\displaystyle {\begin{aligned}\left.{\frac {d}{dt}}\right|_{t=0}E[I+th]&={\frac {d}{dt}}{\big |}_{t=0}{\frac {1}{2}}\int _{\Omega }g\left(\|\nabla (I+th)(x)\|^{2}\right)\,dx\\&=\int _{\Omega }g'\left(\|\nabla I(x)\|^{2}\right)\nabla I\cdot \nabla h\,dx\\&=-\int _{\Omega }\mathrm {div} (g'\left(\|\nabla I(x)\|^{2}\right)\nabla I)h\,dx\end{aligned}}}
假设
∇
E
I
{\displaystyle \nabla E_{I}}
L
2
(
Ω
,
R
)
{\displaystyle L^{2}(\Omega ,\mathbb {R} )}
∇
E
I
=
−
d
i
v
(
g
′
(
‖
∇
I
(
x
)
‖
2
)
∇
I
)
{\displaystyle \nabla E_{I}=-\mathrm {div} (g'\left(\|\nabla I(x)\|^{2}\right)\nabla I)}
因此，其梯度下降法的方程式可以表示成
∂
I
∂
t
=
−
∇
E
I
=
d
i
v
(
g
′
(
‖
∇
I
(
x
)
‖
2
)
∇
I
)
{\displaystyle {\frac {\partial I}{\partial t}}=-\nabla E_{I}=\mathrm {div} (g'\left(\|\nabla I(x)\|^{2}\right)\nabla I)}
我们假设
c
=
g
′
{\displaystyle c=g'}
修正后的Perona-Malik模型，又被成为正规化的P-M方程式，其未知部分在非线性部分与高斯函数进行卷积，得到
∂
I
∂
t
=
d
i
v
(
c
(
|
D
G
σ
∗
I
|
)
∇
I
)
{\displaystyle {\frac {\partial I}{\partial t}}=\mathrm {div} \left(c(|DG_{\sigma }*I|)\nabla I\right)}
其中
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 模型的话，其结果会在区域内进行扩散而且不会使其超过较强的边界，因此在移除噪声后，影像中的边界及结构仍可以保留下来。
除了移除噪声之外，非等向性扩散也可以用于边界侦测。只要根据边界侦测方程式来进行多次递回的非等向性扩散，其最终结果影像会趋向于剩下一个一个的色块，而相邻色块之间的区域则会被侦测为边界。
[ 1]
^ wiki_en . [2017-07-07 ] . （原始内容 存档于2021-01-19）.
![{\displaystyle E[I]={\frac {1}{2}}\int _{\Omega }g\left(\|\nabla I(x)\|^{2}\right)\,dx}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0827f250202d043383ae984378b55d2984020f85)
![{\displaystyle {\begin{aligned}\left.{\frac {d}{dt}}\right|_{t=0}E[I+th]&={\frac {d}{dt}}{\big |}_{t=0}{\frac {1}{2}}\int _{\Omega }g\left(\|\nabla (I+th)(x)\|^{2}\right)\,dx\\&=\int _{\Omega }g'\left(\|\nabla I(x)\|^{2}\right)\nabla I\cdot \nabla h\,dx\\&=-\int _{\Omega }\mathrm {div} (g'\left(\|\nabla I(x)\|^{2}\right)\nabla I)h\,dx\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/58c06d6bc4f85f8385223d0f228c107898ab1ab7)
