[发明专利]一种基于图像金字塔的数字X射线影像校正方法

说明书

一种基于图像金字塔的数字X射线影像校正方法，包括以下步骤：步骤S1，输入图像A和触发行，其中所述触发行为N行；步骤S2，读取输入所述图像A，进行全局强度补偿，获得亮度均匀的图像B；步骤S3，对所述图像B中触发行附近的数据进行局部补偿得到图像C；步骤S4，基于图像金字塔对所述图像C中触发行附近数据进一步校正得图像D；步骤S5，对所述图像D中的坐标(x，y)包含于集合H的点，使用m_max插值，得到最终的图像E；步骤S6，输出校正后所述图像E。本发明可有效提升数字X射线系统AED模式得到图像的质量。

技术领域

本发明属于医学图像处理技术领域，具体地涉及一种基于图像金字塔的数字X射线影像校正方法。

背景技术

在便携DR或移动DR拍片使用自动曝光探测(Automatic Exposure Detection，AED)触发模式时，因探测射线所须的计算时间及平板暗电流存在，导致得到的图像有明显的触发行，图像模糊不均匀的问题。

发明内容

本发明的目的是提供一种基于图像金字塔的数字X射线影像校正方法以解决图像不均匀的问题。

为实现上述目的，本发明的技术方案如下：

一种基于图像金字塔的数字X射线影像校正方法，包括以下步骤：

步骤S1，输入图像A和触发行，所述触发行为N行；

步骤S2，读取输入所述图像A，并遍历每个像素点，得到过曝的坐标集合，记为集合H，进行全局强度补偿，获得亮度均匀的图像B，；

步骤S3，对所述图像B中触发行附近的数据进行局部补偿得到图像C；

步骤S4，基于图像金字塔对所述图像C中触发行附近数据进一步校正得图像D；

步骤S5，对所述图像D中的所有像素点，使用m_max替代坐标包含于所述步骤S2得到的坐标集合H的像素值，得到最终的图像E；

步骤S6，输出校正后所述图像E。

进一步地，所述步骤S2具体包括以下子步骤：

步骤S2.1，定义m_max为图像像素值的最大值；读取输入所述图像A，并遍历每个像素点，得到过曝的坐标集合，记为集合H；

步骤S2.2，计算全局校正所需参数：定义a为统计的行数；统计在图像A中N-a行到N+a行图像的属性；计算其表征图像属性值；所述表征图像属性值选用图像平均值，定义A_mean1为图像A的第N-a行到第N行图像的平均值，将N-a行到第N行图像的灰度值进行累加，后除N-a行到N行的像素个数可得图像均值A_mean1；A_mean2为图像A的第N行到第N+a行的平均值，将N行到第N+a行图像的灰度值进行累加，后除以N行到第N+a行的像素个数可得图像均值A_mean2。

步骤S2.3，根据步骤S2.2计算所得的所述参数，对所述图像A进行全局校正：定义全局校正因子θ，所述校正因子θ范围从0到1之间；则图像A中从0行到触发行N-1，校正的方法为A(x,y)＝A(x,y)-θ*(A_mean1–A_mean2)；从第N+1行到h行，校正方法为A(x,y)＝A(x,y)+θ*(A_mean1-A_mean2),校正后得到图像B

进一步地，将所述步骤S2.2中的所述表征图像属性值图像平均值替换为图像平均值、图像中位数、图像标准差、图像最大值和图像最小值中的一种。