[发明专利]一种分割乳腺病灶的方法

权利要求书

1.一种分割乳腺病灶的方法，其特征在于，包括如下步骤：

(1)对注射造影剂前的磁共振图像基于椭圆模型进行胸部粗分割，获得第一图像；

(2)对所述注射造影剂前、后的磁共振图像进行剪影，在获得的剪影图像上基于所述第一图像提取乳腺部分的剪影图像；

(3)对每层乳腺部分的剪影图像进行病灶检测，获得第二图像；

(4)在所述第二图像上依次找到三维连通域，依次对每个连通域以重心为种子点在所述剪影图像上进行自适应的区域增长，获得分割后的乳腺可疑病灶的磁共振图像。

2.如权利要求1所述的分割乳腺病灶的方法，其特征在于，所述胸部粗分割前进行滤波器去噪处理和闭操作。

3.如权利要求2所述的分割乳腺病灶的方法，其特征在于，所述滤波器去噪处理利用中值滤波器、高斯滤波器、高通滤波器、均值滤波器或低通滤波器。

4.如权利要求1所述的分割乳腺病灶的方法，其特征在于，所述第一图像的取得包括如下步骤：

(1)对所述注射造影剂前的磁共振图像从上到下逐行进行检测，当检测到连通域的范围达到阈值M时，停止检测并记录行数N′；

(2)以行数N₁＝N′+n作为椭圆的起始行，以所述注射造影剂前的磁共振图像的最后一行作为椭圆的终止行为行数N₂，进行椭圆分割处理；

(3)对椭圆分割后的图像从上到下逐行检测，当测得的灰度值为非零值时，停止检测，记录行数N′₁作为乳腺部分的起始行，并记录所述非零值的像素点所在的列数L；

(4)对椭圆分割后的图像从下到上逐行检测第L列，当测得的灰度值为非零值时，停止检测，记录行数N′₂作为乳腺部分的终止行，以所述起始行N′₁和终止行N′₂提取所述乳腺部分，获得第一图像；

所述阈值M的取值范围为所述磁共振图像宽度的40～60％，所述行数n取自然数，取值范围为5～20行。

5.如权利要求1所述的分割乳腺病灶的方法，其特征在于，所述剪影步骤前对所述注射造影剂前、后的磁共振图像进行配准操作。

6.如权利要求1～5中任意一项所述的分割乳腺病灶的方法，其特征在于，所述注射造影剂后的磁共振图像为造影剂初始增强的磁共振图像。

7.如权利要求1所述的分割乳腺病灶的方法，其特征在于，所述病灶检测的方法包括如下步骤：

(1)对每层乳腺部分的剪影图像根据灰度直方图进行阈值分割；

(2)在阈值分割后的图像上进行检测，若三维连通区域的体积小于阈值b，则去除所述三维连通区域；

(3)在阈值分割后的图像上，依次得到三维连通域，对每个连通域进行二维投影，对所述二维投影进行检测，当满足时，则去除所述三维连通域，获得第二图像；

所述l_max为所述二维投影最长轴的长度，所述S为该连通域二维投影的面积，所述阈值b的取值范围为10～50立方毫米，所述c为倍数，其取值范围为0.9～1.2。

8.如权利要求7所述的分割乳腺病灶的方法，其特征在于，所述阈值分割方法包括如下步骤：

(1)对每层乳腺部分的剪影图像计算灰度直方图，根据所述灰度直方图分别定义第k层的高阈值T_k1和低阈值T_k2；

(2)若第k层图像的高阈值T_k1＜a×T_ave，则以高阈值T_k1对该层进行阈值分割；

(3)若第k层图像的高阈值T_k1≥a×T_ave，则以低阈值T_k2对该层进行阈值分割；

所述k取自然数，所述阈值T_ave为各层剪影图的高阈值平均值，所述a为倍数，取值范围为0.9～1.2。

9.如权利要求7或8所述的分割乳腺病灶的方法，其特征在于，对所述阈值分割后的图像进行滤波器去噪处理。

10.如权利要求7所述的分割乳腺病灶的方法，其特征在于，所述连通域检测前进行闭操作。