DoG算子
D
o
G
(
D
i
f
f
e
r
e
n
c
e
o
f
G
a
u
s
s
i
a
n
)
DoG(Difference of Gaussian)
DoG(DifferenceofGaussian)是高斯函数的差分。它是可以通过将图像与高斯函数进行卷积得到一幅图像的低通滤波结果,即去噪过程。
高斯函数:
G
σ
(
x
,
y
)
=
1
2
π
σ
2
e
x
p
(
−
x
2
+
y
2
2
σ
2
)
G_{\sigma}(x,y)=\frac{1}{\sqrt{2\pi \sigma^{2}}}exp(-\frac{x^2+y^2}{2\sigma^2})
Gσ(x,y)=2πσ21exp(−2σ2x2+y2)
对图像高斯滤波:
g
(
x
,
y
)
=
G
σ
(
x
,
y
)
∗
f
(
x
,
y
)
g(x,y)=G_{\sigma}(x,y)*f(x,y)
g(x,y)=Gσ(x,y)∗f(x,y)
对同一图像进行不同
σ
\sigma
σ的高斯滤波结果相减得到
D
o
G
DoG
DoG图:
g
1
(
x
,
y
)
−
g
2
(
x
,
y
)
=
G
σ
1
(
x
,
y
)
∗
f
(
x
,
y
)
−
G
σ
2
(
x
,
y
)
∗
f
(
x
,
y
)
=
(
G
σ
1
−
G
σ
2
)
∗
f
(
x
,
y
)
=
D
o
G
∗
f
(
x
,
y
)
g_{1}(x,y)-g_{2}(x,y)=G_{\sigma_1}(x,y)*f(x,y)-G_{\sigma_2}(x,y)*f(x,y)=(G_{\sigma_1}-G_{\sigma_2})*f(x,y)=DoG*f(x,y)
g1(x,y)−g2(x,y)=Gσ1(x,y)∗f(x,y)−Gσ2(x,y)∗f(x,y)=(Gσ1−Gσ2)∗f(x,y)=DoG∗f(x,y)
D
o
G
DoG
DoG表示为:
D
o
G
≜
G
σ
1
−
G
σ
2
=
1
2
π
(
1
σ
1
e
−
(
x
2
+
y
2
)
/
(
2
σ
1
2
)
−
1
σ
2
e
−
(
x
2
+
y
2
)
/
(
2
σ
2
2
)
)
DoG\triangleq G_{\sigma1}-G_{\sigma2}=\\ \frac{1}{\sqrt{2\pi}}(\frac{1}{\sigma_{1}}e^{-(x^2+y^2)/(2\sigma_{1}^2)}-\\ \frac{1}{\sigma_{2}}e^{-(x^2+y^2)/(2\sigma_{2}^2)})
DoG≜Gσ1−Gσ2=2π1(σ11e−(x2+y2)/(2σ12)−σ21e−(x2+y2)/(2σ22))
LoG算子
拉普拉斯算子为
▽
2
f
=
∂
2
f
∂
x
2
+
∂
2
f
∂
y
2
\bigtriangledown ^2f=\frac{\partial ^2f}{\partial x^2}+\frac{\partial ^2f}{\partial y^2}
▽2f=∂x2∂2f+∂y2∂2f
对二维高斯函数应用拉普拉斯算子
▽
2
G
=
∂
2
G
∂
x
2
+
∂
2
G
∂
y
2
=
−
2
σ
2
+
x
2
+
y
2
2
π
σ
5
e
−
(
x
2
+
y
2
)
/
(
2
σ
2
)
\bigtriangledown ^2G=\frac{\partial ^2G}{\partial x^2}+\\ \frac{\partial ^2G}{\partial y^2}=\\ \frac{-2\sigma ^2+x^2+y^2}{2\pi\sigma ^5}e^{-(x^2+y^2)/(2\sigma ^2)}
▽2G=∂x2∂2G+∂y2∂2G=2πσ5−2σ2+x2+y2e−(x2+y2)/(2σ2)
L
o
G
(
L
a
p
l
a
c
i
a
n
o
f
G
a
u
s
s
i
a
n
)
LoG(Laplacian of Gaussian )
LoG(LaplacianofGaussian)算子定义为
L
o
G
=
σ
2
▽
2
G
LoG=\sigma ^2\bigtriangledown ^2G
LoG=σ2▽2G
DoG算子可以近似LoG算子
图像下采样(降采样)
下采样就是将图像缩小,对于一个图像I尺寸为MN,对其进行s倍下采样,即得到(M/s)(N/s)分辨率的图像,对于矩阵形式的图像,就是把原始图像s*s窗口内的像素变成一个像素,这个像素点的值是窗口内所有像素点的均值,对于第k个窗口:
p
k
=
∑
i
∈
w
i
n
(
k
)
I
(
i
)
/
s
2
p_k=\sum_{i\in win(k)}I_{(i)}/s^2
pk=i∈win(k)∑I(i)/s2
积分图
为了加速求局部矩阵求和运算,建立一个和原始图像相同的二维数组S,其中
S
(
x
,
y
)
=
∑
x
′
≤
x
,
y
′
≤
y
I
(
x
′
,
y
′
)
S(x,y)=\sum_{x'\leq x,y'\leq y}I(x',y')
S(x,y)=x′≤x,y′≤y∑I(x′,y′)
即
S
(
x
,
y
)
S(x,y)
S(x,y)相当于该点与图像原点构成的矩形的点的像素和
比如求区域D所有像素的和: S ( 4 ) − S ( 3 ) − S ( 2 ) + S ( 1 ) {S(4)-S(3)-S(2)+S(1)} S(4)−S(3)−S(2)+S(1)
三角剖分
对下面图像三角剖分
【定义】三角剖分:假设V是二维实数域上的有限点集,边e是由点集中的点作为端点构成的封闭线段,E为e的集合。那么该点集V的一个三角剖分T=(V,E)是一个平面图G,该平面图满足条件:
1、除了端点,平面图中的边不包含点集中的任何点。
2、没有相交边。//边和边没有交叉点
3、平面图中所有的面都是三角面,且所有三角面的合集是散点集V的凸包。
凸包的概念:用不严谨的话来讲,给定二维平面上的点集,凸包就是将最外层的点连接起来构成的凸多边型,它能包含点集中所有的点。
具体解释
https://blog.csdn.net/piaoxuezhong/article/details/68065170