第二初始分割 点Os的横坐标,参见图5所示,图5示出了本实施例中所述第二初始分割点Os的横坐标和所述 第二分割点Ps的横坐标之间的关系。
[0105] 执行S1013,根据所述第二分割点Ps确定第一分割点Pf,本实施例中具体的,选择从 乳房最上方开始沿乳房边缘向下的乳房边缘的长度属于[0.2Lb,0.化b]的部分乳房边缘,其 中Lb为乳房边缘上的首末两点对应在竖直方向上的坐标之间的距离,参见图6,本实施例中 为Xf (首点横坐标)、Xe(末点横坐标)之间的距离,即Xe,图6中加粗显示的上侧乳房的乳房边 缘部分即为本实施例中选取的部分乳房边缘,连接所述第二分割点Ps与上侧乳房的乳房边 缘上的各个像素点,如图6中虚线所示,计算第二分割点Ps与上侧乳房的乳房边缘上的各个 像素点所在直线的斜率。取直线斜率最大时(图6中用实线表示)与所述第二分割点Ps连接 的位于上侧乳房的乳房边缘上的像素点为第一分割点Pf。
[0106] 执行S1014,连接所述第一分割点Pf和第二分割点Ps,二者之间的连线与所述乳房 图像的水平边界和竖直边界所组成的区域,如图7中所示的直角Ξ角形区域即为包含了胸 肌区域的初始区域,图7中直角Ξ角形区域的斜边为所述初始区域的胸肌壁分界线。
[0107] 需要说明的是,本实施例中是先在所述乳房图像的水平边界上选取第一初始分割 点,然后基于所述第一初始分割点和从乳房最下方开始(近胸墙侧)沿乳房边缘向上的部分 乳房边缘上各像素点之间连线的斜率的大小W确定第二初始分割点,根据所述第二初始分 割点的位置W确定所述第二分割点,进而根据所述第二分割点与上侧乳房的乳房边缘上的 像素点之间的连线斜率的大小确定第一分割点。而在其他实施例中,也可W先在所述乳房 图像的竖直边界上选取分割点,然后根据竖直边界上的分割点进而确定在所述乳房图像水 平边界上的分割点,连接竖直边界和水平边界上的分割点,二者的连线和所述乳房图像的 水平边界、竖直边界所组成的区域即为包含了胸肌区域的初始区域。因此,确定所述初始区 域过程中,先确定水平方向的分割点还是先确定竖直方向的分割点不应作为对本发明技术 方案的限定。
[0108] 此外,本实施例中根据第一初始分割点和部分乳房边缘上像素点连线的斜率的大 小W确定第二初始分割点,W及根据第二分割点和部分乳房边缘上像素点连线的斜率的大 小确定第一分割点,究竟取斜率大的点还是斜率小的点和胸肌区域的形态相关,而选取的 部分乳房边缘的长度,W及根据第二初始分割点的位置确定第二分割点的位置,则是根据 对实际采集到的多幅乳房图像样本进行统计W确定的。因此,在实际应用中,可W根据不同 地区的乳房形态的不同,对部分乳房边缘W及第二分割点的位置进行相应的调整。因此,本 实施例中选取的部分乳房边缘的长度,W及选取的第二分割点的位置,也不应作为对本发 明技术方案的限定。
[0109] 通过上述步骤在所述乳房图像中进行粗分割确定包括了胸肌区域的初始区域后, 执行S102,确定所述初始区域中胸肌区域和其他组织区域的分割阔值。本实施例中主要通 过对所述初始区域的灰度直方图进行形态学分析,如胸肌区域的灰度值相对于乳房区域的 灰度值而言,胸肌区域的灰度值低,从胸肌区域过渡到乳房区域时灰度值变化剧烈,灰度值 变化率大等等。参见图8,图8是本发明实施例一的确定胸肌区域和其他组织区域的分割阔 值的流程图,所述确定胸肌区域和其他组织区域的分割阔值包括:
[0110] S1020:对所述初始区域的灰度直方图归一化W获得第一灰度直方图;
[0111] S1021:确定所述第一灰度直方图中像素点个数最多对应的灰度值,在大于该灰度 值所对应的第一灰度直方图中捜索与该灰度值坐标的预设邻域内的任一灰度值坐标距离 最近的点,W该点对应的灰度值为所述乳房图像中胸肌区域和其他组织区域的分割阔值。
[0112] 本实施例中,在对所述初始区域的灰度直方图进行形态学分析前,可W先对所述 初始区域的灰度直方图进行平滑处理W消除灰度直方图的毛刺,如采用均值滤波器对所述 初始区域的灰度直方图进行平滑,滤波器核的大小可W通过初始区域的灰度直方图中像素 点的灰度分布范围及对应的像素点的个数来确定。举例来说:若初始区域的灰度直方图中 像素点的灰度分布范围为300~500,且在该灰度分布范围中包括了 10000个像素点,那么滤 波器核为取不小于
的奇数,即滤波器的核为9。在对所述初始区域的灰度直方图进 行均值滤波后,对其进形态学分析。W下结合图9和图10对本实施例中的对所述初始区域的 灰度直方图进行形态学分析进行相应的说明。
[0113] 参见图9,图9是本发明实施例一的初始区域的灰度直方图,其中横坐标为灰度值, 纵坐标为像素点个数,执行S1020,对图9所示的灰度直方图进行归一化,具体地就是对图9 中的纵坐标(像素点个数)进行归一化,使得图9中峰值处所对应的像素点的个数与图9中最 大灰度值与峰值处灰度值的差相等(图9中正方形的虚线框所示)。若图9中的最大灰度值为 1200,峰值处的灰度值为300且对应的像素点个数为4500,则就是对图9中纵坐标上的每一 个值除W5,W对所述初始区域的灰度直方图进行归一化获得第一灰度直方图。
[0114] 执行S1021,对归一化的初始区域的灰度直方图即第一灰度直方图进行形态学分 析W确定胸肌区域和其他组织区域的分割阔值,如图10所示,首先确定在所述第一灰度直 方图中像素点个数最多对应的灰度值Gm,在大于灰度值Gm所对应的第一灰度直方图中,即图 10中所示虚线的右半部分,捜索与该灰度值Gm坐标的预设邻域内的任一灰度值坐标距离最 近的点,本实施例中所述预设邻域的半径可W为[0,0.15Gl],其中Gl为第一灰度直方图中灰 度区间的长度,即第一灰度直方图的灰度分布范围,W图10为例,Gl为1000,若Gm的灰度值为 300,则可W在灰度值为300到450的范围内捜索图10中所示虚线的右半部分中与该范围内 灰度值坐标距离最近的点。图10中W预设邻域的半径为0进行了说明,也即在图10中所示虚 线的右半部分捜索与Gm坐标距离最近的点,也即WGm坐标为圆屯、做虚线右侧的灰度直方图 的内切圆,切点T即为与Gm坐标距离最近的点,切点T对应的灰度值则为所述乳房图像中胸 肌区域和其他组织区域的分割阔值。
[0115] 通过上述步骤确定了胸肌区域和其他组织区域的分割阔值后,大于等于该分割阔 值的像素点为其他组织区域的像素点,而小于该分割阔值的像素点则为胸肌区域的像素 点,进而可W在所述初始区域中分割出胸肌区域。
[0116] 至此,通过上述的对所述乳房图像先进行粗分割确定包括胸肌区域的呈Ξ角形的 初始区域,进而通过对该初始区域的灰度直方图进行形态学分析W确定胸肌区域和其他组 织区域的分割阔值,通过先粗分割后精确分割的方式提高了胸肌区域检测的准确度,同时 可W使得在对所述乳房图像进行后处理前将胸肌区域去除,避免胸肌区域对后续处理的干 扰,进而提高最终获得的乳房图像的质量。
[0117] 在另一实施例中,在对所述初始区域的灰度直方图进行归一化前,还可W先对其 进行预处理(预处理前也可W先对初始区域的灰度直方图进行平滑处理),W去除其他组织 区域对待检测的胸肌区域的影响获得第二灰度直方图;然后对所述第二灰度直方图归一化 W获得第Ξ灰度直方图;最后,确定所述第Ξ灰度直方图中像素点个数最多对应的灰度值, 在大于该灰度值所对应的第Ξ灰度直方图中捜索与该灰度值坐标的预设邻域内的任一灰 度值坐标距离最近的点,W该点对应的灰度值为所述乳房图像中胸肌区域和其他组织区域 的分割阔值。
[0118] W下结合图11对上述的预处理进行相应的说明。图11是对初始区域的灰度直方图 进行预处理的示意图,图11中横坐标为灰度值,纵坐标为像素点个数,1200为所述初始区域 的灰度直方图中的最大灰度值,而第一灰度值则为与最大灰度值距离最近的灰度值,也即 最大灰度值与所述第一灰度值之差最小,且第一灰度值对应的像素点个数不小于平均像素 点个数;或者说第一灰度值是W最大灰度值作为起点,在灰度值逐渐减小的过程中,第一个 灰度值对应的像素点个数不小于平均像素点个数的灰度值。图11中第一灰度值为Gf,去除 比第一灰度值Gf大的区域,也即去除图11中所示的阴影线所在的区域。Mp的值则通过初始区 域中像素点个数的总和与所述初始区域中像素点的灰度值范围的比值而定。举例来说:若 所述初始区域中像素点个数的总和为100000,像素点灰度值的范围为300~500,则平均像 素点个数Mp为500。
[0119] 在又一实施例中,还可W基于大津法在所述初始区域中确定所述胸肌区域和其他 组织区域的分割阔值,且采用大津法确定分割阔值时,也可W先对所述初始区域的灰度直 方图进行平滑处理,再对平滑处理后的初始区域的灰度直方图进行预处理,进而在进行了 平滑处理和预处理后的初始区域的灰度直方图中确定胸肌区域和其他组织区域的分割阔 值。
[0120] 基于上述的乳房图像中检测胸肌的方法,本发明实施例还提供一种乳房图像中检 测胸肌的装置,请参考图12,图12是本发明实施例一的乳房图像中检测胸肌的装置结构示 意图,如图12所示:所述乳房图像中检测胸肌的装置包括:
[0121] 第一区域确定单元101,用于在所述乳房图像中确定初始区域,所述初始区域至少 包括部分胸肌区域,所述乳房图像是指探测器采集到的原始乳房图像;
[0122] 第一确定单元102,用于在所述初始区域中确定所述胸肌区域和其他组织区域的 分割阔值W在所述初始区域中分割出胸肌区域。
[0123] 本实施例中,所述第一区域确定单元101包括:
[0124] 第一初始分割点确定单元,用于在所述乳房图像的边界上选取第一初始分割点;
[0125] 第二初始分割点确定单元,用于基于所述第一初始分割点与部分乳房边缘上各像 素点之间连线的斜率的大小确定第二初始分割点;
[0126] 第二分割点确定单元,用于根据所述第二初始分割点的位置确定第二分割点,所 述第二分割点坐标的绝对值小于所述第二初始分割点坐标的绝对值;
[0127] 第一分割点确定单元,用于基于所述第二分割点与部分乳房边缘上各像素点之间 连线的斜率的大小确定第一分割点;
[0128] 第四确定单元,用于确定所述第一分割点和第二分割点之间的连线与所述