本文简要介绍尺度空间理论。如果你对这个理论是如何提出、发展的,可以详细地阅读相关文献,见参考文献。
简单地说,尺度空间理论(scale-space theory)的提出主要是为了解决如何识别不同大小物体的问题。通过将一系列不同“大小”的高斯滤波器 G ( x , y ; t ) G(x,y;t) G(x,y;t)与原始图像 I ( x , y ) I(x,y) I(x,y)进行卷积,获得不同分辨率的图像 L ( x , y ; t ) L(x,y;t) L(x,y;t)。
对于二维图像 I ( x , y ) I(x,y) I(x,y)在不同尺度 t t t下的图像, L ( x , y ; t ) = G ( x , y ; t ) ∗ I ( x , y ) L(x,y;t)=G(x,y;t)*I(x,y) L(x,y;t)=G(x,y;t)∗I(x,y)
其中, G ( x , y ; t ) = 1 2 π t e − ( x 2 + y 2 ) / 2 t G(x,y;t)=\frac{1}{2\pi t}e^{-(x^2+y^2)/2t} G(x,y;t)=2πt1e−(x2+y2)/2t
由卷积的微分性质,可以得到
L ( x , y ; t ) x x = ∂ 2 ( G ( x , y ; t ) ∗ I ( x , y ) ) ∂ 2 x = ∂ 2 G ( x , y ; t ) ∂ 2 x ∗ I ( x , y ) L(x,y;t)_{xx}=\frac{\partial^2{(G(x,y;t)*I(x,y))}}{\partial^2{x}}=\frac{\partial^2{G(x,y;t)}}{\partial^2{x}}*I(x,y) L(x,y;t)xx=∂2x∂2(G(x,y;t)∗I(x,y))=∂2x∂2G(x,y;t)∗I(x,y)
L ( x , y ; t ) y y = ∂ 2 ( G ( x , y ; t ) ∗ I ( x , y ) ) ∂ 2 y = ∂ 2 G ( x , y ; t ) ∂ 2 y ∗ I ( x , y ) L(x,y;t)_{yy}=\frac{\partial^2{(G(x,y;t)*I(x,y))}}{\partial^2{y}}=\frac{\partial^2{G(x,y;t)}}{\partial^2{y}}*I(x,y) L(x,y;t)yy=∂2y∂2(G(x,y;t)∗I(x,y))=∂2y∂2G(x,y;t)∗I(x,y)
L ( x , y ; t ) x y = ∂ 2 ( G ( x , y ; t ) ∗ I ( x , y ) ) ∂ x ∂ y = ∂ 2 G ( x , y ; t ) ∂ x ∂ y ∗ I ( x , y ) L(x,y;t)_{xy}=\frac{\partial^2{(G(x,y;t)*I(x,y))}}{\partial{x}\partial{y}}=\frac{\partial^2{G(x,y;t)}}{\partial{x}\partial{y}}*I(x,y) L(x,y;t)xy=∂x∂y∂2(G(x,y;t)∗I(x,y))=∂x∂y∂2G(x,y;t)∗I(x,y)
因此,图像的Hessian矩阵为
H
(
x
,
y
)
=
[
∂
2
G
(
x
,
y
;
t
)
∂
2
x
∂
2
G
(
x
,
y
;
t
)
∂
x
∂
y
∂
2
G
(
x
,
y
;
t
)
∂
x
∂
y
∂
2
G
(
x
,
y
;
t
)
∂
2
y
]
∗
I
(
x
,
y
)
H(x,y)=\begin{bmatrix} \frac{\partial^2{G(x,y;t)}}{\partial^2{x}} & \frac{\partial^2{G(x,y;t)}}{\partial{x}\partial{y}} \\ \frac{\partial^2{G(x,y;t)}}{\partial{x}\partial{y}} & \frac{\partial^2{G(x,y;t)}}{\partial^2{y}} \\ \end{bmatrix}*I(x,y)
H(x,y)=[∂2x∂2G(x,y;t)∂x∂y∂2G(x,y;t)∂x∂y∂2G(x,y;t)∂2y∂2G(x,y;t)]∗I(x,y)
参考文献
- 1996, Lindeberg Tony, Scale-space: A framework for handling image structures at multiple scales
- 2000, Krissian Karl, Model Based Detection of Tubular Structures in 3D Images
- 2009, Marcin Rudzki, Vessel Detection Method Based on Eigenvalues of the Hessian Matrix and its Applicability to Airway Tree Segmentation