数字X射线尘肺图像的处理方法及系统与流程

本发明涉及图像处理技术领域，尤其涉及一种数字x射线尘肺图像的处理方法及系统。

背景技术：

尘肺病是由粉尘大量聚集于肺组织内引起的，造成肺组织的纤维化，形成进行性的损伤。尘肺病患者在早期并没有明显的症状，随着病变进展大多会出现咳嗽、咯血、胸痛、呼吸困难甚至更严重的并发症。尘肺病是我国最主要的职业病，是导致劳动能力降低、致残和影响寿命的疾病，也是国家和企业赔偿的主要职业病。因此，尘肺病的早期诊断和早期治疗对患者、企业和国家都有非常重要的意义。

尘肺病的诊断是以国家统一制作的尘肺标准片为对照，结合患者粉尘接触史，再由3名以上获得职业病尘肺诊断医师资格的诊断组进行诊断。随着计算机技术和图像处理技术的飞速发展，数字x射线成像技术得到迅速发展。在gbz70-2015尘肺病诊断标准中，对尘肺诊断标准片进行了调整和补充，并增加了附录f(规范性附录)“数字x射线胸片摄影的技术要求”来规范数字x射线摄影技术。与传统x线胸片相比，数字x线胸片包含的信息量更大，具有图像细腻、对比度好、辐射剂量小等优点。由于gbz70-2015尘肺病诊断标准规定不能对尘肺图像进行降噪和边缘增强等处理，因此常规的数字x射线图像处理方法不适用于数字x射线尘肺图像。

申请号为201710699009.x，发明名称为《用于尘肺病诊断的dr双能量减影高kv胸部摄片方法》的发明专利中，采用dr双能量减影联合高kv进行胸部后前位摄片，能清晰的显示粉尘小阴影及其细微结构阴影，提高一级片率和肺内小阴影的检出率，对尘肺病的早期诊断具有重要的临床意义。

目前公开的尘肺饼相关的资料基本都是关于尘肺图像的摄影方法和尘肺病诊断方法及标准，还没有明确的针对数字x射线尘肺图像的处理方法。

技术实现要素：

本发明实施例所要解决的技术问题在于，提供一种数字x射线尘肺图像的处理方法及系统，以更加清晰的显示粉尘小阴影及其细微结构，提高诊断的准确性和尘肺病的检出率。

为了解决上述技术问题，本发明实施例提出了一种数字x射线尘肺图像的处理方法，包括：

灰度拉伸步骤：统计采集的受检者的数字x射线尘肺图像a中灰度值的最小值mingraya和最大值maxgraya，对图像a进行灰度拉伸，调整灰度的动态范围，得到图像b；

对比度调整步骤：统计图像b中灰度值的最大值maxgrayb，调整图像b的对比度，得到最终的处理后图像c；

窗宽窗位计算步骤：统计图像c中灰度值的最小值mingrayc和最大值maxgrayc，计算最优窗宽ww和最优窗位wl，改善图像显示时的视觉效果；

输出显示步骤：根据得到的ww和wl输出显示图像c。

相应地，本发明实施例还提供了一种数字x射线尘肺图像的处理系统，包括：

灰度拉伸模块：统计采集的受检者的数字x射线尘肺图像a中灰度值的最小值mingraya和最大值maxgraya，对图像a进行灰度拉伸，调整灰度的动态范围，得到图像b；

对比度调整模块：统计图像b中灰度值的最大值maxgrayb，调整图像b的对比度，得到最终的处理后图像c；

窗宽窗位计算模块：统计图像c中灰度值的最小值mingrayc和最大值maxgrayc，计算最优窗宽ww和最优窗位wl，改善图像显示时的视觉效果；

输出显示模块：根据得到的ww和wl输出显示图像c。

本发明的有益效果为：本发明先用gamma曲线调整图像的动态范围；再用s曲线调整图像的对比度；最后计算图像的最优窗宽窗位，并用得到的最优窗宽窗位显示图像，用于诊断；本发明得到的图像动态范围广、对比度好、视觉效果佳，粉尘小阴影及其细微结构显示清晰，有效提高尘肺病的检出率，对尘肺病的早期诊断具有重要的临床意义。

附图说明

图1是本发明实施例的数字x射线尘肺图像的处理方法的流程示意图。

图2是本发明实施例的图像a的示意图。

图3是本发明实施例的图像c的示意图。

图4是本发明实施例的数字x射线尘肺图像的处理系统的结构示意图。

具体实施方式

需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互结合，下面结合附图和具体实施例对本发明作进一步详细说明。

本发明实施例中若有方向性指示(诸如上、下、左、右、前、后……)仅用于解释在某一特定姿态(如附图所示)下各部件之间的相对位置关系、运动情况等，如果该特定姿态发生改变时，则该方向性指示也相应地随之改变。

另外，在本发明中若涉及“第一”、“第二”等的描述仅用于描述目的，而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。

请参照图1，本发明实施例的数字x射线尘肺图像的处理方法包括灰度拉伸步骤、对比度调整步骤、窗宽窗位计算步骤、输出显示步骤。

灰度拉伸步骤：统计采集的受检者的数字x射线尘肺图像a中灰度值的最小值和最大值，分别记为mingraya和maxgraya，对图像a进行灰度拉伸，调整灰度的动态范围，得到图像b。

对比度调整步骤：统计图像b中灰度值的最大值，记为maxgrayb，调整图像b的对比度，得到最终的处理后图像c。

窗宽窗位计算步骤：统计图像c中灰度值的最小值和最大值，分别记为mingrayc和maxgrayc，计算最优窗宽ww和最优窗位wl，改善图像显示时的视觉效果。

输出显示步骤：根据得到的ww和wl输出显示图像c。

作为一种实施方式，灰度拉伸步骤中，采用下式对图像a进行灰度拉伸：

bitsstored＝1＜＜bits；

其中，i和j表示像素点的横纵坐标，bits表示图像的位数，m表示灰度拉伸力度，取值范围为0～1。优选地，m＝0.2。

作为一种实施方式，对比度调整步骤中，采用下式调整图像b的对比度：

其中，i和j表示像素点的横纵坐标，e为自然常数，e＝2.718281828，n表示对比度力度，取值范围为1～9。

作为一种实施方式，窗宽窗位计算步骤中，采用下式计算最优窗宽ww和最优窗位wl：

ww＝wwgain*(maxgrayc-mingrayc+1)；

其中，wwgain表示窗宽系数，取值范围为0.5～1.5；wlgain表示窗位系数，取值范围为0.5～1.5。优选地，wwgain＝0.9；wlgain＝0.8。

作为一种实施方式，所述灰度拉伸步骤前还包括：

图像采集步骤：受检者保持胸部后前位，源像距为180cm，x射线(中心线)对准受检者的第六胸椎；根据受检者的体型选择管电压和管电流，在受检者充分吸气后保持屏气状态时进行曝光，曝光10ms～30ms，得到图像a。

本发明实施例用高千伏(120kv—140kv)采集图像。人体组织对x射线的吸收能力随着管电压的升高而变小，从而能获得较大的曝光宽容度，提高图像的对比度，改善图像质量；本发明实施例的曝光时间要小于100ms，并且要求受检者在曝光前充分吸气后保持屏气状态，直到曝光结束。目的是避免曝光过程因受检者呼吸产生运动伪影而影响诊断。在曝光的过程中，病人的心跳和呼吸会在图像中产生运动伪影，影响图像质量，曝光时间越短，运动伪影的影响越小，因此要求曝光时间小于100ms。

示例：受检者保持胸部后前位，源像距(sid)为180cm，x射线对准受检者的第六胸椎；选择合适的管电压(120kv—140kv)，根据受检者的体型，适当调整管电压；选择小焦点(管电流≤160ma)，曝光时间控制在10ms—30ms之间。受检者充分吸气后保持屏气状态，曝光完成后，受检者可以正常呼吸，得到原始图像a，请参照图2；根据本发明实施例的数字x射线尘肺图像的处理方法对原始图像a进行处理，其中，m＝0.2，wwgain＝0.9；wlgain＝0.8，处理后输出显示的图像如图3所示。

请参照图4，本发明实施例的数字x射线尘肺图像的处理系统包括灰度拉伸模块、对比度调整模块、窗宽窗位计算模块及输出显示模块。

灰度拉伸模块：统计采集的受检者的数字x射线尘肺图像a中灰度值的最小值mingraya和最大值maxgraya，对图像a进行灰度拉伸，调整灰度的动态范围，得到图像b。

对比度调整模块：统计图像b中灰度值的最大值maxgrayb，调整图像b的对比度，得到最终的处理后图像c。

窗宽窗位计算模块：统计图像c中灰度值的最小值mingrayc和最大值maxgrayc，计算最优窗宽ww和最优窗位wl，改善图像显示时的视觉效果；

输出显示模块：根据得到的ww和wl输出显示图像c。

作为一种实施方式，灰度拉伸模块中，采用下式对图像a进行灰度拉伸：

bitsstored＝1＜＜bits；

其中，i和j表示像素点的横纵坐标，bits表示图像的位数，m表示灰度拉伸力度，取值范围为0～1。

作为一种实施方式，m＝0.2。

作为一种实施方式，对比度调整模块中，采用下式调整图像b的对比度：

其中，i和j表示像素点的横纵坐标，e为自然常数，e＝2.718281828，n表示对比度力度，取值范围为1～9。

作为一种实施方式，窗宽窗位计算模块中，采用下式计算最优窗宽ww和最优窗位wl：

ww＝wwgain*(maxgrayc-mingrayc+1)；

其中，wwgain表示窗宽系数，取值范围为0.5～1.5；wlgain表示窗位系数，取值范围为0.5～1.5。优选地，wwgain＝0.9；wlgain＝0.8。

尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同范围限定。