[发明专利]基于对称性匹配滤波器组和区域生长的三维血管分割方法

权利要求书

1.一种基于对称性匹配滤波器组和区域生长的三维血管分割方法，其特征在于，包括如下步骤：

步骤1：设计由三维高斯函数和复正弦函数实部乘积组成的符合血管管状结构的基本匹配滤波器模板；

步骤2：构建正交平面将基本匹配滤波器模板等分成四个象限的匹配滤波器模板之和；

步骤3：构造多尺度、多方向的对称性的匹配滤波器组，与待增强的图像进行空间卷积，对于尺度大于设定阈值的匹配滤波器，使用基本匹配滤波器减去四象限卷积两两做减法的绝对值，对于尺度不大于设定阈值的匹配滤波器，使用基本匹配滤波器减去对称象限的卷积响应的绝对值，同一尺度下的多个方向取最大值作为增强图像的卷积响应值，最后融合对称性匹配滤波器组的增强卷积响应，各个体素点选取与其尺度最相符的滤波器的空间卷积结果作为其增强结果值，得到增强图像；包括：

步骤3.1：选定m_j尺度，0<j<＝n，的l个方向的滤波器与CTA/MRA图像做空间卷积，根据图像实际血管尺度设定尺度阈值k，当m_j<＝k时，得到原图与对称性匹配滤波器的卷积响应值3DVSMFR_j为：

当m_j>k时，得到原图与对称性匹配滤波器的卷积响应值3DVSMFR_j为：

其中3DMFR_j表示基本匹配滤波器与原图的卷积响应，3DMFR_j^a、3DMFR_j^b、3DMFR_j^c、3DMFR_j^d分别表示四个象限的基本匹配滤波器与原图的卷积响应值，各个体素点选取与其方向最相符的滤波器的空间卷积结果作为该尺度下的最终响应值，最终形成m_j尺度下的响应图像；

步骤3.2：对m_j尺度下的响应图像，采取对于边缘高亮区域和响应异常点赋0值处理进行异常点的去除；

步骤3.3：对所有n个尺度的最终响应值进行归一化处理，各个体素点选取与其尺度最相符即卷积响应输出值最大的滤波器的空间卷积结果作为其增强结果值；

步骤4：利用固定阈值三维区域生长算法提取增强图像的血管，获取最终的分割结果。

2.根据权利要求1所述的基于对称性匹配滤波器组和区域生长的三维血管分割方法，其特征在于，步骤1中设计的在第j个尺度m_j下第i个方向的基本匹配滤波器模板h_ij(x_j,y_j,z_j)＝g_j(x_j,y_j,z_j)×sr_ij(x_j,y_j,z_j)；

其中(x_j,y_j,z_j)表示三维坐标，

表示尺度为m_j的3D高斯球的第j个点的值，其中σ_x表示g_j(x_j,y_j,z_j)在x轴上的标准差，σ_y表示g_j(x_j,y_j,z_j)在y轴上的标准差，σ_z表示g_j(x_j,y_j,z_j)在z轴上的标准差；

表示尺度为m_j、母线方向为(cosβ_icosα_i,cosβ_isinα_i，sinβ_i)的圆柱体的第j个点的值，其中F_j表示径向空间频率，β表示单位球面上某个离散点在球坐标系下的仰角，α是单位球面上某个离散点在球坐标系下的方位角。

3.根据权利要求1所述的基于对称性匹配滤波器组和区域生长的三维血管分割方法，其特征在于，所述步骤4中包括：

(4.1)在增强图像的血管主干入口处设置生长种子点；

(4.2)设置生长阈值范围；

(4.3)使用三维全局区域生长提取增强图像的血管部分，并将分割结果作为掩膜获取原图中的血管分割结果。