ITK的血管提取的方法主要来自于文章:Three-dimensional multi-scale line filter for segmentation and visualization of curvilinear structures in medical images。
文章的大体思路是构造一个3D线性的血管滤波器,用于过滤噪声,保留血管结构。
首先介绍几个概念:一阶导和二阶导。
一阶导数描述的变换的情况,图像的一阶导数通常用于图像灰度梯度的变化情况,如图所示,一阶导那行,响应较大(绝对值较大)的地方为边缘。
二阶倒数还包含了变换的方向,且响应的位置比一阶导更为细致,如上第二行所示,响应的