数字X射线尘肺图像的处理方法及系统技术方案

本发明专利技术实施例公开了一种数字X射线尘肺图像的处理方法及系统，所述方法包括：灰度拉伸步骤：统计采集的受检者的尘肺图像A中灰度值的最小值和最大值，对图像A进行灰度拉伸，调整灰度的动态范围，得到图像B；对比度调整步骤：统计图像B中灰度值的最大值，调整图像B的对比度，得到图像C；窗宽窗位计算步骤：统计图像C中灰度值的最小值和最大值，计算最优窗宽和最优窗位；输出显示步骤：根据窗宽窗位显示图像C。本发明专利技术得到的图像动态范围广、对比度好、视觉效果佳，粉尘小阴影及其细微结构显示清晰，有效提高尘肺病的检出率，对尘肺病的早期诊断具有重要的临床意义。

【技术实现步骤摘要】
数字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射线的吸收能力随着管电压的升高而变小，从而能获得较大的曝光宽容度，提高图像的对比度，改善图像质量；本专利技术实施例的曝光时间要小于10本文档来自技高网...

【技术保护点】
1.一种数字X射线尘肺图像的处理方法，其特征在于，包括：/n灰度拉伸步骤：统计采集的受检者的数字X射线尘肺图像A中灰度值的最小值MinGrayA和最大值MaxGrayA，对图像A进行灰度拉伸，调整灰度的动态范围，得到图像B；/n对比度调整步骤：统计图像B中灰度值的最大值MaxGrayB，调整图像B的对比度，得到最终的处理后图像C；/n窗宽窗位计算步骤：统计图像C中灰度值的最小值MinGrayC和最大值MaxGrayC，计算最优窗宽WW和最优窗位WL，改善图像显示时的视觉效果；/n输出显示步骤：根据得到的WW和WL输出显示图像C。/n

【技术特征摘要】
1.一种数字X射线尘肺图像的处理方法，其特征在于，包括：
灰度拉伸步骤：统计采集的受检者的数字X射线尘肺图像A中灰度值的最小值MinGrayA和最大值MaxGrayA，对图像A进行灰度拉伸，调整灰度的动态范围，得到图像B；
对比度调整步骤：统计图像B中灰度值的最大值MaxGrayB，调整图像B的对比度，得到最终的处理后图像C；
窗宽窗位计算步骤：统计图像C中灰度值的最小值MinGrayC和最大值MaxGrayC，计算最优窗宽WW和最优窗位WL，改善图像显示时的视觉效果；
输出显示步骤：根据得到的WW和WL输出显示图像C。


2.如权利要求1所述的数字X射线尘肺图像的处理方法，其特征在于，灰度拉伸步骤中，采用下式对图像A进行灰度拉伸：



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


3.如权利要求1所述的数字X射线尘肺图像的处理方法，其特征在于，对比度调整步骤中，采用下式调整图像B的对比度：


















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


4.如权利要求1所述的数字X射线尘肺图像的处理方法，其特征在于，窗宽窗位计算步骤中，采用下式计算最优窗宽WW和最优窗位WL：
WW＝WWGain*(MaxGrayC-MinGrayC+1)；



其中，WWGain表示窗宽系数，取值范围为0.5～1.5；WLGain表示窗位系数，取值范围为0.5～1.5。


5.如权利要求1所述的数字X射线尘肺图像的处理方法，其特征在于，所述灰度拉伸步骤前还包括：
图像采集步骤：受检者保持胸部后前位，源像距为180cm，X射线对准受检者的第六胸椎；根据受检者的体型...

【专利技术属性】
技术研发人员：郭朋，李勇，鄢照龙，郑杰，
申请(专利权)人：深圳蓝韵医学影像有限公司，
类型：发明
国别省市：广东;44