数字乳腺x射线成像系统曝光剂量的评价方法
【技术领域】
[0001] 本发明涉及一种医学影像处理技术，特别涉及一种医学影像曝光剂量的评价方 法。
【背景技术】
[0002] 乳腺作为女性健康的重要课题，已经越来越被人们所重视。近年来，乳腺癌的发病 率在我国逐年上升，在许多大城市中乳腺癌已经成为女性最常见的恶性肿瘤，严重威胁着 女性的身体健康，研究表明乳腺癌如果能早期发现，那么将大大提高其治愈率，所以乳腺癌 的早期诊断具有非常重要的临床价值。
[0003] 实际临床发现乳腺癌的早期表现为乳腺内存在微小的钙化点，这些微小的钙化点 在使用超声成像时很难发现，而使用数字乳腺X射线成像系统对乳腺进行检查是早期发现 钙化点的有效方式。对于X射线检查，辐射剂量是非常重要的，应在获得高质量图像的前提 下尽可能的降低曝光参数。因此需要一种评价乳腺图像曝光剂量的方法和相关标准。目前 能够搜索到的近似相关专利有：
[0004] 《一种乳腺机X射线剂量控制方法及系统》，申请人：北京国药恒瑞美联信息技术 有限公司。此方法通过针对不同性质的目标乳腺进行预曝光获得预曝光图像，并对预曝光 图像进行分析计算从而获得目标乳腺在不同投照角度的X射线剂量值，保证获取最佳的拍 摄图像。此方法只能控制当次曝光的曝光参数，而无法评价曝光剂量的高低。
[0005] 《乳腺X光机及其实现全自动曝光的方法》，申请人：深圳圣诺医疗有限公司。此方 法首先根据压迫厚度确定预曝光参数，然后利用预曝光图像的剂量信息确定第二次曝光 的剂量参数。此方法可以自动的计算合适的曝光参数，但是无法评价曝光剂量的高低。
【发明内容】
[0006] 为了解决上述的技术问题，本发明提出一种数字乳腺X射线成像系统曝光剂量的 评价方法，本发明采用如下技术方案来解决上述技术问题：
[0007] 一种数字乳腺X射线成像系统曝光剂量的评价方法，包括如下步骤：
[0008] a.制定目标曝光指数；
[0009] b.计算偏差指数DI。
[0010] 优选的，在上述的一种数字乳腺X射线成像系统曝光剂量的评价方法，中所述步 骤a包括如下步骤：
[0011] al.收集临床病例中对应体位的100幅原始图像；
[0012] a2?计算每幅图像曝光指数EI ;
[0013] a3.将所有图像的EI值按照从小到大排序得到EI值的序列Lei ;
[0014] a4.从序列Lei中取中间的10个EI值，求平均值，这个平均值即为目标曝光指数 EI t;
[0015] a5.将EIt写入配置文件或数据库，供计算偏差指数DI时调用。
[0016] 优选的，在上述一种数字乳腺X射线成像系统曝光剂量的评价方法中，所述步骤 a2包括如下步骤：
[0017] a21.校正曝光指数，获取图像斜率Slop和偏移量Offset ;
[0018] a22?读入原始图像；
[0019]a23.分割出乳腺组织区域ROI，计算ROI的像素均值Mean和像素点个数Count ;
[0020] a24.对乳腺组织区域进行二次分割，分割出组织区域中的致密区域secondROI， 计算secondROI的像素均值secondMean和像素点个数secondCount ;
[0021] a25.利用公式计算乳腺组织区域的感兴趣值Value，所述公式为：
[0022] Value = (secondMean*secondCount+Mean*(Count-secondCount))/Count ;
[0023] a26?利用公式EI = 10*(Slop*Value+0ffset)计算曝光指数EI，并将EI写入 Dicom文件。
[0024] 优选的，在上述的一种数字乳腺X射线成像系统曝光剂量的评价方法中，步骤a21 具体包括如下步骤：
[0025] A.将剂量仪放置于平板探测器的胸壁侧，确定RQM8辐射质量；
[0026] B.分别采集6不同mAs的图像，并记录每幅图像对应的剂量Dose ;
[0027] C.计算每幅图像的均值Mean ;
[0028] D?计算Mean与Dose间的线性关系，通过线性拟合公式Dose = Slop*Mean+0ffset 确定斜率Slop和偏移量Offset，将Slop和Offset写入配置文件或数据库，供计算曝光指 数EI时调用。
[0029] 优选的，在上述的一种数字乳腺X射线成像系统曝光剂量的评价方法中，
[0030] 所述步骤b为：使用公式计算偏差指数DI，
成的DI写入Dicom文件。
[0031] 不同的乳腺中致密腺体区域所占整个腺体区域的比例不同，本发明计算乳腺组织 区域的感兴趣值时考虑了致密腺体区域占整个腺体区域的比例，提高了相关数据的准确 性。使用本发明的技术方案，根据曝光指数与目标曝光指数的偏差来评价乳腺图像曝光剂 量高低的方法，从而方便医生后续调整曝光参数，获得曝光剂量适当的图像，算法简单，使 用方便。
【附图说明】
[0032] 图1为本发明流程图；
[0033] 图2为本发明曝光指数校正流程图；
[0034] 图3为计算曝光指数流程图；
[0035] 图4为制定目标曝光指数流程图。
【具体实施方式】
[0036] 下面结合附图对本发明作进一步详细说明，但是本发明不仅仅局限于以下实施 例。本发明所述的一种数字乳腺X射线成像系统曝光剂量的评价方法，具体包括以下步骤：[0037] 1.曝光指数校正：
[0038] 因不同平板探测器对剂量的响应不同，故使用不同类型的平板探测器前均需要首 先进行曝光指数的校正。
[0039] 将剂量仪放置于平板探测器的胸壁侧即患者拍摄时靠近身体的一侧，确定RQM8 辐射质量。RQM8辐射质量是本领域研究人员所熟知的相关标准，在此不做描述。
[0040] 分别采集6不同mAs的图像，并记录每幅图像对应的剂量Dose。
[0041] 计算每幅图像的均值Mean。采用OTSU分割算法分割掉剂量仪区域；为了克服边 界的影响，可适当裁剪掉图像的部分边界，本例中胸壁侧不裁剪，其余三个边界均裁剪掉20 个像素的宽度。OTSU分割算法是本领域研究人员所熟知的相关技术，故在此不做描述。
[0042] 计算Mean与Dose间的线性关系。通过线性拟合公式Dose = Slop*Mean+Offset 确定斜率Slop和偏移量Offset。可将Slop和Offset写入配置文件或数据库，供计算曝 光指数EI时调用。
[0043] 2.计算曝光指数EI :
[0044] 在得到Slop和Offset的前提下计算曝光指数EI。该数据表征拍摄图像时对患者 的曝光程度。
[0045] 读入原始图像0。图像允许相关的校正处理，但是不允许进行增强处理。
[0046] 利用OTSU分割算法分割出乳腺组织区域ROI。计算ROI的像素均值Mean和像素 点个数Count。OTSU分割算法是本领域研究人员所熟知的相关技术，故在此不做描述。
[0047] 利用OTSU分割算法分割对乳腺组织区域进行二次分割，分割出组织区域中的致 密区域secondROI，计算secondROI的像素均值secondMean和像素点个数secondCount。
[0048] 利用公式 Value = (secondMean*secondCount+Mean* (Count-secondCount)) / Count计算乳腺组织区域的感兴趣值Value。
[0049] 利用公式EI = 10* (Slop*Value+0ffset)计算曝光指数EI，并将EI写入Dicom文 件。
[0050] 3.目标曝光指数的制定：
[0051] 不同体位的目标曝光指数不同，本实施例以RCC位即右侧正位为例，介绍相关的 流程。
[0052] 收集100副RCC位的原始图像。
[0053] 根据步骤2计算每幅图像的曝光指数EI，利用快速排序法对所有图像的EI值由小 到大排序，得到EI值的序列Lei。快速排序法是本领域研究人员所熟知的相关技术，故在此 不做描述。
[0054] 从序列Lei中取中间的10个EI值，这10个EI值的平均值就作为RCC位的目标 曝光指数EI t。
[0055] 将EIt写入配置文件或数据库，供计算偏差指数DI时调用。
[0056] 4.计算偏差指数DI :
[0057] DI表示曝光剂量的高低，DI为正时表示曝光剂量偏高，DI为负时表示曝光剂量偏 低。由下式计算：
[0059] 将DI写入Dicom文件，方便以后查看图像曝光剂量的高低。
[0060] 以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精 神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
【主权项】
1. 一种数字乳腺X射线成像系统曝光剂量的评价方法，其特征在于，包括如下步骤： a. 制定目标曝光指数； b. 计算偏差指数DI。2. 根据权利要求1所述的一种数字乳腺X射线成像系统曝光剂量的评价方法，其特征 在于，所述步骤a包括如下步骤： al.收集临床病例对应体位的100幅原始图像；计算每幅图像曝光指数EI; a2.将所有图像的EI值按照从小到大排序得到EI值的序列Lu; a3.从序列Lu中取中间的10个EI值，求平均值，送个平均值即为目标曝光指数EIt; a4.将EIt写入配置文件或数据库，供计算偏差指数DI时调用。3. 根据权利要求2所述的一种数字乳腺X射线成像系统曝光剂量的评价方法，其特征 在于，所述步骤a2包括如下步骤： 曰21.校正曝光指数，获取图像斜率Slop和偏移量Offset;a22.读入原始图像； a23.分割出乳腺组织区域R0I，计算ROI的像素均值Mean和像素点个数Count;a24.对乳腺组织区域进行二次分割，分割出组织区域中的致密区域secondROI，计算 secondROI的像素均值SecondMean和像素点个数secondCount; 曰25.利用公式计算乳腺组织区域的感兴趣值Value,所述公式为： Value= (secondMean*secondCount+Mean*(Count-secondCount))/Count;a26.利用公式EI= 10*(Slop*Value+0ffset)计算曝光指数EI，并将EI写入Dicom 文件。4. 根据权利要求3所述的一种数字乳腺X射线成像系统曝光剂量的评价方法，其特征 在于，步骤a21具体包括如下步骤： A. 将剂量仪放置于平板探测器的胸壁侧，确定RQM8福射质量； B. 分别采集6不同mAs的图像，并记录每幅图像对应的剂量Dose; C. 计算每幅图像的均值Mean; 化计算Mean与Dose间的线性关系，通过线性拟合公式Dose=Slop*Mean+Offset确 定斜率Slop和偏移量Offset,将Slop和Offset写入配置文件或数据库，供计算曝光指数 EI时调用。5. 根据权利要求1所述的一种数字乳腺X射线成像系统曝光剂量的评价方法，其特征 在于，所述步骤b为：使用公式计算偏差指数DI，所述公式为将生成的DI与入Dicom文件。
【专利摘要】本发明公开一种数字乳腺X射线成像系统曝光剂量的评价方法，包括如下步骤:制定目标曝光指数；计算偏差指数DI。使用本发明的技术方案，根据曝光指数与目标曝光指数的偏差来评价乳腺图像曝光剂量高低的方法，从而方便医生后续调整曝光参数，获得曝光剂量适当的图像，算法简单，使用方便。
【IPC分类】A61B6/00
【公开号】CN105193436
【申请号】CN201410280940
【发明人】郭朋, 秦启兴, 鄢照龙, 李勇, 赵明
【申请人】深圳市蓝韵实业有限公司
【公开日】2015年12月30日
【申请日】2014年6月20日