CN103637815B - 乳腺自动曝光参考区域的确定方法及系统 - Google Patents
乳腺自动曝光参考区域的确定方法及系统 Download PDFInfo
- Publication number
- CN103637815B CN103637815B CN201310699971.5A CN201310699971A CN103637815B CN 103637815 B CN103637815 B CN 103637815B CN 201310699971 A CN201310699971 A CN 201310699971A CN 103637815 B CN103637815 B CN 103637815B
- Authority
- CN
- China
- Prior art keywords
- breast
- row
- breast area
- area
- pixel quantity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 210000005075 mammary gland Anatomy 0.000 title claims abstract description 93
- 238000000034 method Methods 0.000 title claims abstract description 29
- 210000000481 breast Anatomy 0.000 claims abstract description 323
- 210000002445 nipple Anatomy 0.000 claims abstract description 86
- 238000012545 processing Methods 0.000 claims abstract description 5
- 239000000284 extract Substances 0.000 claims description 10
- 238000000605 extraction Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 11
- 238000006243 chemical reaction Methods 0.000 description 9
- 210000000038 chest Anatomy 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 206010006187 Breast cancer Diseases 0.000 description 2
- 208000026310 Breast neoplasm Diseases 0.000 description 2
- 210000000746 body region Anatomy 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 210000003484 anatomy Anatomy 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及图像处理技术领域，尤其涉及一种乳腺自动曝光参考区域的确定方法及系统。与现有技术相比，本发明提供的方法及系统可以从乳腺曝光图像中检测出乳房方向及乳头的大致位置，并根据乳房方向及乳头的大致位置确定出最适合乳腺影像检测的自动曝光参考区域。根据该自动曝光参考区域在预曝光时的亮度值与曝光剂量的关系可推算出该区域要获得预期的亮度值所需的曝光剂量，从而使该区域在正式曝光时获得最接近预期的亮度值，为得出准确的乳腺影像诊断结果提供了技术保证。
Description
技术领域
本发明涉及图像处理技术领域，尤其涉及一种乳腺自动曝光参考区域的确定方法及系统。
背景技术
在医学影像曝光领域中，现有X射线自动曝光系统的原理是：先用小剂量X射线拍摄小尺寸的影像，在小尺寸影像上选取固定的几块区域作为参考，以第一张影像上该区域的亮度值与目标亮度值的关系，用当前剂量推算要获取目标亮度值时所需的剂量大小。用于推算剂量的区域的位置和大小目前是固定的，适合于拍胸、腰等部位，如果拍其它部位，如：四肢、乳腺，固定地选取参考区域，很可能存在所选的参考区域处于空气穿透区域，而不是人体组织部分。而选取空气区域作为参考区域，没有任何意义，不能正确的推算出实际所需的曝光剂量。如果能根据所拍的体位，检测出人体有效组织区域，在组织区域选参数区域，才有意义。
发明内容
本发明所要解决的技术问题是：提供一种乳腺自动曝光参考区域的确定方法及系统，用于从乳腺曝光图像中检测出乳头附近区域作为自动曝光参考区域，以解决传统乳腺影像诊断中固定选取自动曝光参考区域的方式不能使乳腺影像检测中最重要的参考区域(即乳头附近区域)准确地获得预期亮度值的问题。
本发明是通过如下方式实现的：
一种乳腺自动曝光参考区域的确定方法，包括如下步骤：
步骤S1：将乳腺曝光图像二值化，并从获得的二值图像中提取乳房区域；
步骤S2：提取所述乳房区域的外边缘轮廓图像并据此判断乳房方向；
步骤S3：根据乳房方向确定乳头的大致坐标；
步骤S4：根据乳房方向及乳头的大致坐标对乳房区域进行截取；
步骤S5：求取被截取的乳房区域的质心，确定所述被截取的乳房区域的以所述质心为中心的最大内接矩形区域，并从所述最大内接矩形区域的靠近乳头一侧起截取设定比例区域作为自动曝光参考区域；
所述步骤S2包括如下步骤：
计算所述外边缘轮廓图像的各行及各列中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最大的行与最大的列；
比较所述最大的行与最大的列中属于乳房区域的像素数量的大小，并根据比较出的大小关系判断乳房方向；当最大的行中属于乳房区域的像素数量大于最大的列中属于乳房区域的像素数量时，判断乳房在所述乳腺曝光图像中为水平方向，当最大的列中属于乳房区域的像素数量大于最大的行中属于乳房区域的像素数量时，判断乳房在所述乳腺曝光图像中为垂直方向。
进一步地，所述步骤S3包括如下步骤：
当乳房在所述乳腺曝光图像中为垂直方向时：
计算乳房区域各行中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最大的行，并确定该行为乳头大致所在的行；
计算乳房区域各列中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最小的列，并确定该列为乳头大致所在的列；
当乳房在所述乳腺曝光图像中为水平方向时：
计算乳房区域各列中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最大的列，并确定该列为乳头大致所在的列；
计算乳房区域各行中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最小的行，并确定该行为乳头大致所在的行。
进一步地，所述步骤S4包括如下步骤：
如果乳房在所述乳腺曝光图像中为垂直方向，以乳头大致所在的行为对称轴，从所述乳房区域中截取设定宽度；
如果乳房在所述乳腺曝光图像中为水平方向，以乳头大致所在的列为对称轴，从所述乳房区域中截取设定宽度。
进一步地，所述步骤S1与所述步骤S2之间还包括如下步骤：
步骤S1-2：对所述乳房区域进行腐蚀及膨胀处理。
一种乳腺自动曝光参考区域的确定系统，包括：
乳房区域提取模块，用于将乳腺曝光图像二值化，并从获得的二值图像中提取乳房区域；
乳房方向确定模块，用于提取所述乳房区域的外边缘轮廓图像并据此判断乳房方向；
乳头坐标确定模块，用于根据乳房方向确定乳头的大致坐标；
乳房区域截取模块，用于根据乳房方向及乳头的大致坐标对乳房区域进行截取；
自动曝光参考区域确定模块，用于求取被截取的乳房区域的质心，确定所述被截取的乳房区域的以所述质心为中心的最大内接矩形区域，并从所述最大内接矩形区域的靠近乳头一侧起截取设定比例区域作为自动曝光参考区域；
所述乳房方向确定模块包括：
像素计算模块，用于计算所述外边缘轮廓图像的各行及各列中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最大的行与最大的列；
方向判断模块，用于当最大的行中属于乳房区域的像素数量大于最大的列中属于乳房区域的像素数量时，判断乳房在所述乳腺曝光图像中为水平方向，当最大的列中属于乳房区域的像素数量大于最大的行中属于乳房区域的像素数量时，判断乳房在所述乳腺曝光图像中为垂直方向。
进一步地，所述乳头坐标确定模块包括：
乳头所在行判断模块，用于当乳房在所述乳腺曝光图像中为垂直方向时计算乳房区域各行中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最大的行，并确定该行为乳头大致所在的行，以及当乳房在所述乳腺曝光图像中为水平方向时计算乳房区域各行中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最小的行，并确定该行为乳头大致所在的行；
乳头所在列判断模块，用于当乳房在所述乳腺曝光图像中为垂直方向时计算乳房区域各列中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最小的列，并确定该列为乳头大致所在的列，以及当乳房在所述乳腺曝光图像中为水平方向时，计算乳房区域各列中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最大的列，并确定该列为乳头大致所在的列。
进一步地，所述乳房区域截取模块包括：
水平截取模块，用于当乳房在所述乳腺曝光图像中为垂直方向时以乳头大致所在的行为对称轴，从所述乳房区域中截取设定宽度；
垂直截取模块，用于当乳房在所述乳腺曝光图像中为水平方向时以乳头大致所在的列为对称轴，从所述乳房区域中截取设定宽度。
进一步地，所述系统还包括：乳房区域处理模块，用于对从所述二值图像中提取的乳房区域进行腐蚀及膨胀处理。
与现有技术相比，本发明提供的方法及系统可以从乳腺曝光图像中检测出乳房方向及乳头的大致位置，并根据乳房方向及乳头的大致位置确定出最适合乳腺影像检测的自动曝光参考区域。根据该自动曝光参考区域在预曝光时的亮度值与曝光剂量的关系可推算出该区域要获得预期的亮度值所需的曝光剂量，从而使该区域在正式曝光时获得最接近预期的亮度值，为得出准确的乳腺影像诊断结果提供了技术保证。
附图说明
图1：本发明实施例提供的乳腺自动曝光参考区域的确定方法流程示意图；
图2：一幅乳房方向为垂直方向的乳腺曝光图像示意图；
图3：对上述乳腺曝光图像进行二值化处理后的效果示意图；
图4：对上述二值化处理获得的二值图像进行灰度值翻转后的效果示意图；
图5：对上述经过灰度翻转后的二值图像提取外边缘轮廓后的效果示意图；
图6：乳腺X光机工作示意图；
图7：乳头大致坐标检测示意图；
图8：根据乳头位置及乳房方向对乳房区域进行截取示意图；
图9：乳房区域被截取部分的示意图；
图10：根据乳房区域被截取部分的质心确定该被截取部分的以该质心为中心的最大内接矩形区域在乳腺曝光图像中的位置示意图；
图11：本发明实施例提供的乳腺自动曝光参考区域的确定系统示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用于解释本发明，并不用于限定本发明。
图1示出了本发明所提供的乳腺自动曝光参考区域确定方法的流程示意图，根据图1所示，本方法包括如下步骤：
步骤S1：将乳腺曝光图像二值化，并从获得的二值图像中提取乳房区域；
步骤S2：提取所述乳房区域的外边缘轮廓图像并据此判断乳房方向；
步骤S3：根据乳房方向确定乳头的大致坐标；
步骤S4：根据乳房方向及乳头的大致坐标对乳房区域进行截取；
步骤S5：求取被截取的乳房区域的质心，确定所述被截取的乳房区域的以所述质心为中心的最大内接矩形区域，并从所述最大内接矩形区域的靠近乳头一侧起截取设定比例区域作为自动曝光参考区域。
下面分别就每一步骤进行详细揭示。
步骤S1：将乳腺曝光图像二值化，并从获得的二值图像中提取乳房区域。图2示出了一幅乳腺曝光图像，参考图2可以看出，未经处理的乳腺曝光图像中的各像素具有不同的灰度等级。而本发明所提供的方法主要是基于对乳房区域的分析而确定自动曝光参考区域，因此，乳腺曝光图像中存在多个灰度等级会对分析乳房区域产生干扰，有必要对乳腺曝光图像进行二值化处理，使二值化处理后的图像中的各像素只存在两个灰度等级，其中一个灰度等级的像素表示乳房区域，另一个则表示背景区域。对乳腺曝光图像进行二值化处理可以区分出乳房区域与背景区域，同时，由于多个灰度等级的像素被处理成了只有两个灰度等级的像素，因此，二值化处理还可以在不失去原有图像整体和局部特征的前提下将原来的乳腺曝光图像数据简单化，减少运算量，提高运算速度。基于以上考虑，在做二值化处理时，其阈值的设定以能够最准确地区分出乳房区域与背景区域为原则。图3示出了乳腺曝光图像经过二值化处理后的效果图，可以看出，该图中只有黑白两种颜色，实际上，该图中只有两个灰度等级，即255和0，因为原乳腺曝光图像中乳房区域的灰度值小于背景区域的灰度值，因此，二值化处理后的图像中，灰度值为255的像素表示背景区域，灰度值为0的像素表示乳房区域。由于根据一般习惯，在二值图像中，灰度值为255的像素用于表示主体区域，灰度值为0像素用于表示背景区域，因此，对上述二值图像进行灰度翻转处理，灰度翻转处理后的图像如图4所示，灰度值为255的像素用于表示主体区域，定义为属于乳房区域的像素，灰度值为0像素用于的表示背景区域，定义为属于背景区域的像素。
步骤S2进一步包括如下步骤：
计算所述外边缘轮廓图像的各行及各列中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最大的行与最大的列；
比较所述最大的行与最大的列中属于乳房区域的像素数量的大小，并根据比较出的大小关系判断乳房方向；当最大的行中属于乳房区域的像素数量大于最大的列中属于乳房区域的像素数量时，判断乳房在所述乳腺曝光图像中为水平方向，当最大的列中属于乳房区域的像素数量大于最大的行中属于乳房区域的像素数量时，判断乳房在所述乳腺曝光图像中为垂直方向。
如图5所示，乳房区域的外边缘轮廓图像是指该乳房区域最外围的那一层像素的集合，它勾勒出了乳房的外形轮廓。在对该步骤进行解释说明前，有必要对乳腺拍摄方法做一个简单说明，因为本发明是基于当前对乳房进行拍摄的一般方法而得出的。如图6所示，一般而言，在拍摄乳腺X射线图像时，都是用乳腺压迫板3将乳房2压扁后再由X射线球管1透过乳腺压迫板3从乳房2的一面对乳房2进行整体拍摄，再由乳房2另一面的平板探测器9进行X射线感光成像。一般而言，平板探测器9的曝光区域为矩形，因此，在拍摄出的乳腺X射线图像中，乳房靠近胸部的一侧所表现出来的轮廓理论上应为一条直线。但由于曝光区域边缘可能未被完全曝光的原因，因此在实际情况下，乳房靠近胸部的一侧所表现出来的轮廓为一条近似直线。当乳房在乳腺曝光图像中为垂直方向时(如图2所示)，该条近似直线在乳房区域的外边缘轮廓图像中也处于垂直方向，如果此时计算所述外边缘轮廓图像的各行及各列中属于乳房区域的像素数量，则属于乳房区域的像素数量最大的列中属于乳房区域的像素数量将会大于属于乳房区域的像素数量最大的行中属于乳房区域的像素数量。当乳房在乳腺曝光图像中为水平方向时，该条近似直线在该外边缘轮廓图像种也处于水平方向，如果此时计算所述外边缘轮廓图像的各行及各列中属于乳房区域的像素数量，则属于乳房区域的像素数量最大的行中属于乳房区域的像素数量将会大于属于乳房区域的像素数量最大的列中属于乳房区域的像素数量。因此，可以通过这种方式确定乳房在乳腺曝光图像中的方向为水平方向还是垂直方向。
为了进一步确保乳房方向检测的准确性，在步骤S1与步骤S2之间插入步骤S1-2:对乳房区域进行先腐蚀后膨胀的处理，使得乳房靠近胸部的一侧与该图像的边缘完全融合，即将乳房靠近胸部的一侧区域中处于曝光区域边缘、本应被曝光而未被曝光的部分用属于乳房区域的像素填充。进行腐蚀膨胀处理后的乳房区域中，乳房靠近胸部一侧的轮廓在实际就完全成为了一条直线，而非近似直线。对腐蚀膨胀处理后的乳房区域提取外边缘轮廓图像，再进行乳房方向判断时，可使得计算出的属于乳房区域的像素数量最大的行中属于乳房区域的像素数量与属于乳房区域的像素数量最大的列中属于乳房区域的像素数量差异更为明显，提高了乳房方向检测的精确性。
步骤S3进一步包括如下步骤：
当乳房在所述乳腺曝光图像中为垂直方向时：
计算乳房区域各行中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最大的行，并确定该行为乳头大致所在的行；
计算乳房区域各列中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最小的列，并确定该列为乳头大致所在的列；
当乳房在所述乳腺曝光图像中为水平方向时：
计算乳房区域各列中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最大的列，并确定该列为乳头大致所在的列；
计算乳房区域各行中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最小的行，并确定该行为乳头大致所在的行。
乳腺组织在乳头附近分布密度最大，乳头附近也是在乳腺检测中最为重要的参考区域，因此，乳头的位置检测意义重大。乳头的位置可以通过确定乳头在乳房区域中所在的行和列来确定。根据乳房的形状可知，如果乳房在乳腺曝光图像中为垂直方向(如图7所示)，则乳头所在行(图7横轴)的属于乳房区域的像素数量大于其他行中属于乳房区域的像素数量，而乳头所在列(图7纵轴)的属于乳房区域的像素数量小于其他列中属于乳房区域的像素数量，如果乳房在乳腺曝光图像中为水平方向，则乳头所在列中属于乳房区域的像素数量大于其他列中属于乳房区域的像素数量，乳头所在行中属于乳房区域的像素数量小于其他行中属于乳房区域的像素数量。因此可以根据这个原理检测乳头的大致坐标。
步骤S4进一步包括如下步骤：
如果乳房在所述乳腺曝光图像中为垂直方向，以乳头大致所在的行为对称轴，从所述乳房区域中截取设定宽度；
如果乳房在所述乳腺曝光图像中为水平方向，以乳头大致所在的列为对称轴，从所述乳房区域中截取设定宽度。
在上述步骤中，截取的设定宽度以能包含乳腺的主要分布区域为原则，一般而言，乳房的最下边缘到对称轴之间的乳房区域可以完全截取(即图8中底部直线与中部直线之间的部分)，再在对称轴上部截取与其下部相等宽度的乳房区域(及图8中上部直线与中部直线之间的部分)，这部分的乳腺组织密度较大，最能体现出乳腺分布，图9所示为乳房区域被截取后的示意图。
在步骤S5中，质心坐标为对所述乳房区域中被截取部分的所有象素点的横坐标及纵坐标求均值而成，求取的质心如图10所示，所述最大内接矩形区域如图10中矩形所示，其可在最大程度上体现乳腺的分布，由于乳腺密度离乳头越近，密度越大，因此，可以根据实际需要从最大内接矩形区域的靠近乳头一侧起截取设定比例区域作为自动曝光参考区域。
自动曝光参考区域确定后，根据在预曝光时该区域的亮度值与所用曝光剂量的对应关系可以推算出想要该区域获得某个亮度值时所应该使用的曝光剂量，从而在正式曝光时使该曝光区域获得最接近预期的亮度值，为得出最准确的乳腺检测结果提供了技术保证。
如图11所示，本发明还提供了一种乳腺自动曝光参考区域的确定系统。该系统包括：
乳房区域提取模块4，用于将乳腺曝光图像二值化，并从获得的二值图像中提取乳房区域；
乳房方向确定模块5，用于提取所述乳房区域的外边缘轮廓图像并据此判断乳房方向；
乳头坐标确定模块6，用于根据乳房方向确定乳头的大致坐标；
乳房区域截取模块7，用于根据乳房方向及乳头的大致坐标对乳房区域进行截取；
自动曝光参考区域确定模块8，用于求取被截取的乳房区域的质心，确定所述被截取的乳房区域的以所述质心为中心的最大内接矩形区域，并从所述最大内接矩形区域的靠近乳头一侧起截取设定比例区域作为自动曝光参考区域；
乳房区域处理模块9，用于对所述乳房区域进行腐蚀及膨胀处理。
其中，乳房方向确定模块5包括：
像素计算模块51，用于计算所述外边缘轮廓图像的各行及各列中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最大的行与最大的列；
方向判断模块52，用于当最大的行中属于乳房区域的像素数量大于最大的列中属于乳房区域的像素数量时，判断乳房在所述乳腺曝光图像中为水平方向，当最大的列中属于乳房区域的像素数量大于最大的行中属于乳房区域的像素数量时，判断乳房在所述乳腺曝光图像中为垂直方向；
乳头坐标确定模块6包括：
乳头所在行判断模块61，用于当乳房在所述乳腺曝光图像中为垂直方向时计算乳房区域各行中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最大的行，并确定该行为乳头大致所在的行，以及当乳房在所述乳腺曝光图像中为水平方向时计算乳房区域各行中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最小的行，并确定该行为乳头大致所在的行；
乳头所在列判断模块62，用于当乳房在所述乳腺曝光图像中为垂直方向时计算乳房区域各列中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最小的列，并确定该列为乳头大致所在的列，以及当乳房在所述乳腺曝光图像中为水平方向时，计算乳房区域各列中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最大的列，并确定该列为乳头大致所在的列；
乳房区域截取模块7包括：
水平截取模块71，用于当乳房在所述乳腺曝光图像中为垂直方向时以乳头大致所在的行为对称轴，从所述乳房区域中截取设定宽度；
垂直截取模块72，用于当乳房在所述乳腺曝光图像中为水平方向时以乳头大致所在的列为对称轴，从所述乳房区域中截取设定宽度。
上述系统各模块的工作过程与上述方法中各步骤相对应，在此不再赘述。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (8)
1.一种乳腺自动曝光参考区域的确定方法，其特征在于，包括如下步骤：
步骤S1：将乳腺曝光图像二值化，并从获得的二值图像中提取乳房区域；
步骤S2：提取所述乳房区域的外边缘轮廓图像并据此判断乳房方向；
步骤S3：根据乳房方向确定乳头的大致坐标；
步骤S4：根据乳房方向及乳头的大致坐标对乳房区域进行截取；
步骤S5：求取被截取的乳房区域的质心，确定所述被截取的乳房区域的以所述质心为中心的最大内接矩形区域，并从所述最大内接矩形区域的靠近乳头一侧起截取设定比例区域作为自动曝光参考区域；
所述步骤S2包括如下步骤：
计算所述外边缘轮廓图像的各行及各列中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最大的行与最大的列；
比较所述最大的行与最大的列中属于乳房区域的像素数量的大小，并根据比较出的大小关系判断乳房方向；当最大的行中属于乳房区域的像素数量大于最大的列中属于乳房区域的像素数量时，判断乳房在所述乳腺曝光图像中为水平方向，当最大的列中属于乳房区域的像素数量大于最大的行中属于乳房区域的像素数量时，判断乳房在所述乳腺曝光图像中为垂直方向。
2.如权利要求1所述的方法，其特征在于，所述步骤S3包括如下步骤：
当乳房在所述乳腺曝光图像中为垂直方向时：
计算乳房区域各行中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最大的行，并确定该行为乳头大致所在的行；
计算乳房区域各列中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最小的列，并确定该列为乳头大致所在的列；
当乳房在所述乳腺曝光图像中为水平方向时：
计算乳房区域各列中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最大的列，并确定该列为乳头大致所在的列；
计算乳房区域各行中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最小的行，并确定该行为乳头大致所在的行。
3.如权利要求1所述的方法，其特征在于，所述步骤S4包括如下步骤：
如果乳房在所述乳腺曝光图像中为垂直方向，以乳头大致所在的行为对称轴，从所述乳房区域中截取设定宽度；
如果乳房在所述乳腺曝光图像中为水平方向，以乳头大致所在的列为对称轴，从所述乳房区域中截取设定宽度。
4.如权利要求1所述的方法，其特征在于，所述步骤S1与所述步骤S2之间还包括如下步骤：
步骤S1-2：对所述乳房区域进行腐蚀及膨胀处理。
5.一种乳腺自动曝光参考区域的确定系统，其特征在于，包括：
乳房区域提取模块，用于将乳腺曝光图像二值化，并从获得的二值图像中提取乳房区域；
乳房方向确定模块，用于提取所述乳房区域的外边缘轮廓图像并据此判断乳房方向；
乳头坐标确定模块，用于根据乳房方向确定乳头的大致坐标；
乳房区域截取模块，用于根据乳房方向及乳头的大致坐标对乳房区域进行截取；
自动曝光参考区域确定模块，用于求取被截取的乳房区域的质心，确定所述被截取的乳房区域的以所述质心为中心的最大内接矩形区域，并从所述最大内接矩形区域的靠近乳头一侧起截取设定比例区域作为自动曝光参考区域；
所述乳房方向确定模块包括：
像素计算模块，用于计算所述外边缘轮廓图像的各行及各列中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最大的行与最大的列；
方向判断模块，用于当最大的行中属于乳房区域的像素数量大于最大的列中属于乳房区域的像素数量时，判断乳房在所述乳腺曝光图像中为水平方向，当最大的列中属于乳房区域的像素数量大于最大的行中属于乳房区域的像素数量时，判断乳房在所述乳腺曝光图像中为垂直方向。
6.如权利要求5所述的乳腺自动曝光参考区域的确定系统，其特征在于，所述乳头坐标确定模块包括：
乳头所在行判断模块，用于当乳房在所述乳腺曝光图像中为垂直方向时计算乳房区域各行中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最大的行，并确定该行为乳头大致所在的行，以及当乳房在所述乳腺曝光图像中为水平方向时计算乳房区域各行中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最小的行，并确定该行为乳头大致所在的行；
乳头所在列判断模块，用于当乳房在所述乳腺曝光图像中为垂直方向时计算乳房区域各列中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最小的列，并确定该列为乳头大致所在的列，以及当乳房在所述乳腺曝光图像中为水平方向时，计算乳房区域各列中属于乳房区域的像素数量，确定其中属于乳房区域的像素数量最大的列，并确定该列为乳头大致所在的列。
7.如权利要求5所述的乳腺自动曝光参考区域的确定系统，其特征在于，所述乳房区域截取模块包括：
水平截取模块，用于当乳房在所述乳腺曝光图像中为垂直方向时以乳头大致所在的行为对称轴，从所述乳房区域中截取设定宽度；
垂直截取模块，用于当乳房在所述乳腺曝光图像中为水平方向时以乳头大致所在的列为对称轴，从所述乳房区域中截取设定宽度。
8.如权利要求5所述的乳腺自动曝光参考区域的确定系统，其特征在于，还包括：乳房区域处理模块，用于对从所述二值图像中提取的乳房区域进行腐蚀及膨胀处理。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310699971.5A CN103637815B (zh) | 2013-12-18 | 2013-12-18 | 乳腺自动曝光参考区域的确定方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310699971.5A CN103637815B (zh) | 2013-12-18 | 2013-12-18 | 乳腺自动曝光参考区域的确定方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103637815A CN103637815A (zh) | 2014-03-19 |
| CN103637815B true CN103637815B (zh) | 2015-08-19 |
Family
ID=50243209
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310699971.5A Active CN103637815B (zh) | 2013-12-18 | 2013-12-18 | 乳腺自动曝光参考区域的确定方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103637815B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2989342A1 (en) * | 2015-06-12 | 2016-12-15 | The Trustees Of Dartmouth College | Systems and methods for guiding tissue resection |
| CN105662445A (zh) * | 2016-02-25 | 2016-06-15 | 深圳柏斯生物科技有限公司 | 数字乳腺x射线成像系统的全自动曝光控制方法和装置 |
| CN107862691B (zh) * | 2017-11-24 | 2020-06-19 | 上海联影医疗科技有限公司 | 检测乳房图像中非乳房区域的方法、装置、计算机和介质 |
| CN108364275B (zh) * | 2018-03-02 | 2022-04-12 | 成都西纬科技有限公司 | 一种图像融合方法、装置、电子设备及介质 |
| CN111528919A (zh) * | 2020-05-08 | 2020-08-14 | 上海联影医疗科技有限公司 | 一种x射线与超声融合的乳腺成像方法、系统和装置 |
| CN114155179A (zh) * | 2020-08-17 | 2022-03-08 | 浙江宇视科技有限公司 | 光源缺陷检测方法、装置、设备和存储介质 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101472381A (zh) * | 2007-12-29 | 2009-07-01 | Ge医疗系统环球技术有限公司 | 控制x射线ct系统中的x射线曝光的方法 |
| CN101627916A (zh) * | 2009-08-10 | 2010-01-20 | 深圳圣诺医疗设备有限公司 | 乳腺x光机及其实现全自动曝光的方法 |
| CN102956035A (zh) * | 2011-08-25 | 2013-03-06 | 深圳市蓝韵实业有限公司 | 用于乳腺x线图像中提取乳腺区域的预处理方法及系统 |
| CN103034857A (zh) * | 2012-12-18 | 2013-04-10 | 深圳市安健科技有限公司 | 自动检测影像图片中曝光区域的方法及系统 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006068373A (ja) * | 2004-09-03 | 2006-03-16 | Fuji Photo Film Co Ltd | 乳頭検出装置およびそのプログラム |
| US8135199B2 (en) * | 2006-12-19 | 2012-03-13 | Fujifilm Corporation | Method and apparatus of using probabilistic atlas for feature removal/positioning |
-
2013
- 2013-12-18 CN CN201310699971.5A patent/CN103637815B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101472381A (zh) * | 2007-12-29 | 2009-07-01 | Ge医疗系统环球技术有限公司 | 控制x射线ct系统中的x射线曝光的方法 |
| CN101627916A (zh) * | 2009-08-10 | 2010-01-20 | 深圳圣诺医疗设备有限公司 | 乳腺x光机及其实现全自动曝光的方法 |
| CN102956035A (zh) * | 2011-08-25 | 2013-03-06 | 深圳市蓝韵实业有限公司 | 用于乳腺x线图像中提取乳腺区域的预处理方法及系统 |
| CN103034857A (zh) * | 2012-12-18 | 2013-04-10 | 深圳市安健科技有限公司 | 自动检测影像图片中曝光区域的方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103637815A (zh) | 2014-03-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103637815B (zh) | 乳腺自动曝光参考区域的确定方法及系统 | |
| US8855387B2 (en) | System for detecting bone cancer metastases | |
| JP2021509721A5 (zh) | ||
| WO2020182036A1 (zh) | 图像处理方法、装置、服务器及存储介质 | |
| CN110246580B (zh) | 基于神经网络和随机森林的颅侧面影像分析方法和系统 | |
| CN102920477A (zh) | 医学影像的目标区域边界确定装置和方法 | |
| CN110189307B (zh) | 一种基于多模型融合的肺结节检测方法及系统 | |
| US10568600B2 (en) | System and method for detecting anatomical regions | |
| CN113450399A (zh) | 一种正位胸片心胸比测量方法及装置 | |
| CN114022554A (zh) | 一种基于yolo的按摩机器人穴位检测与定位方法 | |
| JP5654561B2 (ja) | 輪郭形成のための自動コントラスト増強法 | |
| KR101128662B1 (ko) | 좌심실 내외벽 경계 검출장치 및 그 방법 | |
| US9320485B2 (en) | System and method for molecular breast imaging | |
| TWI675646B (zh) | 乳房影像分析方法、乳房影像分析系統及非暫態電腦可讀取媒體 | |
| CN105205784B (zh) | 标尺识别方法及影像拼接方法 | |
| US8831325B2 (en) | Radiographic image processing apparatus and radiographic image processing program | |
| Schreibmann et al. | Voxel clustering for quantifying PET‐based treatment response assessment | |
| CN106447664A (zh) | 确定匹配点对的方法及图像获取方法 | |
| CN115330896A (zh) | 金属物体识别方法、系统及装置 | |
| JP2001307064A (ja) | 照射領域抽出方法、装置および記録媒体 | |
| CN109602430B (zh) | 骨科射线成像机 | |
| CN118121164B (zh) | 基于风湿病变区域多维检测数据的骨质状态筛查方法和系统 | |
| US11948389B2 (en) | Systems and methods for automatic detection of anatomical sites from tomographic images | |
| Mandal et al. | Study of variation of pelvis positioning for patients suffering from rectal cancer using daily kilo-voltage cone beam ct images | |
| US20240029256A1 (en) | Determining target object type and position |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| C56 | Change in the name or address of the patentee | ||
| CP03 | Change of name, title or address |
Address after: 518000, Room 408, building A, Hua Han Innovation Park, Nanshan District, Shenzhen, Guangdong Patentee after: SHENZHEN ANGELL TECHNOLOGY CO., LTD. Address before: 518000 Hua Han Road, 16, Nanshan District, Guangdong, Shenzhen Province, China Han technology, 3A Patentee before: Shenzhen Angell Technology Co., Ltd. |