一种去除X射线图像条纹噪声的方法与流程

本发明涉及图像噪声处理研究领域，特别涉及一种去除x射线图像条纹噪声的方法。

背景技术：

x射线图像已经广泛应用于工业探伤、货物检测等领域。但是在使用x射线线阵探测器时，由于线阵增益校准不当或机械故障等原因，在获得的x射线图像上会出现条纹噪声，具体表现为在一幅应该是灰度均匀的图像中某些列的灰度明显高于或低于平均像素灰度的竖条纹。这些不均匀的竖条纹会作为背景出现在实际采集的x射线图像中，破坏图像质量，并且如果基于这种有条纹噪声的图像进行目标识别检测或测量将会出现严重的偏差。

目前，对于条纹噪声的去除方法主要分为以下两类，即直接插值法或基于傅里叶变换的方法。直接插值法即利用双三次插值等方法，在检测出条纹的位置后，根据其两侧的像素值对条纹上的像素重新赋值。这种方法对较细的条纹效果较好，但是对于较粗的条纹，该方法效果不明显。而基于傅里叶变换的方法不能去除间距无规律的条纹，并且计算的复杂度较高，对正常部分也会带来一定的画质损伤。

为此，针对位置不规律、条纹粗细不一的噪声进行有效去除，具有重要的研究意义和实用价值。

技术实现要素：

为了克服现有技术的上述缺点与不足，本发明的目的在于提供一种去除x射线图像条纹噪声的方法，该方法利用了线阵采集图像每一列具有统一的增益系数的特性，可将位置不规律、粗细不一的条纹噪声进行有效去除，使图像画质得以大大提高，另外对其他常见噪声也有一定的抑制作用。

本发明的目的通过以下的技术方案实现：一种去除x射线图像条纹噪声的方法，包括步骤：

首先，对原始噪声图像建立高斯金字塔和拉普拉斯金字塔，对得到的拉普拉斯金字塔细节函数进行高频滤波，利用滤波后的细节函数重建图像，得到平滑后的图像；

其次，将原始噪声图像和平滑后的图像作点除运算，取结果中每一列的中值得到行向量，将行向量作为对应于原始噪声图像每一列的附加增益值；

最后，将原始噪声图像的每一列乘以对应的附加增益值，即完成去噪。

本发明根据线阵采集图像每一列具有统一的增益系数的特性，利用图像金字塔，计算每一列的噪声图像的附加增益，通过将原始噪声图像乘以对应的附加增益值调整图像各列的灰度，有针对性的消除条纹噪声，在基本不损伤原图的画质的情况下大大提高去噪效果。

优选的，对噪声图像建立高斯金字塔的方法是：

按行对原始噪声图像i1进行一维中值滤波得到滤波后的图像i2；

对i2建立高斯金字塔g，g的层数为12，原始噪声图像作为第1层输入，g中第2层到第6层是其下一层对列进行一维下采样后得到的，而第7层到第12层尺寸不变。

优选的，对噪声图像建立拉普拉斯金字塔的方法是：

根据建立的高斯金字塔g，由g建立对应于i2的拉普拉斯金字塔l，显然l中仅有关于列的细节函数l1～l11；建立l的具体步骤如下：

设g中第i层矩阵为gi，将g12利用双三次插值法恢复到g11的尺寸，得到g′11；

计算l11＝g11-g′11；

对g中除g1以外的所有元素重复上述两个步骤，得到l11～l1。

优选的，对拉普拉斯金字塔中的细节函数l1～l11进行高频滤波，方法是：

将l11的全部元素取绝对值，得到|l11|；

取|l11|的中值元素m0，令m1＝4*m0；

对|l11|中的每个元素值p，计算v＝p/m1-1；若v＞1，则将v置1，若v＜0，将v置0，由此得到矩阵v；

计算新的细节函数l′11＝l11.*v，其中“.*”符号表示矩阵对应位置元素相乘；

对l中剩下的每一层均进行上述相同的操作，得到新的拉普拉斯金字塔l′。

优选的，利用滤波后的细节函数重建图像，方法是：借助高频滤波后的拉普拉斯金字塔l′，由高斯金字塔g的最顶层重建原始尺寸的图像，得到平滑后的图像i′。

优选的，将原始噪声图像i1的每一列乘以对应的附加增益值y，得到初始去噪图像i3(i)，对于计算的第i列的附加增益值y(i)，若|y(i)-1|＞p，则i3(i)进行中值滤波，滤波后图像为最终的结果图像。

本发明与现有技术相比，具有如下优点和有益效果：

1、针对噪声位置不规律，形成的噪声条纹粗细不一的问题，本发明利用图像金字塔，计算噪声图像的附加增益，从而去除条纹噪声，去噪的效果好。

2、本发明方法具有去噪效果好且去噪效果稳定的优点，能有效处理不同粗细，不同分布的条纹噪声。算法几乎不会对原图画质造成损伤。

附图说明

图1为本发明实施例去除x射线线阵图像条纹噪声的算法流程图。

图2为本发明实施例带有条纹噪声的x射线原始噪声图像。

图3为图2所示原始噪声图像通过一维中值滤波后的图像。

图4为本实施例经过去噪处理后由金字塔重建的平滑图像。

图5为本实施例附加增益的函数图像。

图6为本实施例最终结果。

具体实施方式

下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。

实施例

在使用x射线线阵探测器时，由于线阵增益校准不当或机械故障等原因，可能产生图2所示的带有条纹噪声的x射线图像，本实施例特别针对该类图像进行去噪方法的研究，当然，该方法处理对象不仅仅限于条纹噪声，对于其他常见噪声也有一定的抑制作用，这里不再展开说明。

以图2所示图像为例，结合图1的流程图，对本实施例一种去除x射线图像条纹噪声的方法的具体步骤说明如下。

s1、按行对原始噪声图像i1进行一维中值滤波，滤波器的阶数取7，得到滤波后的图像i2，图3给出了本实施例中i2的图像，可见仅用滤波或插值的方法不能有效的去除条纹噪声，也不能保证画质。

s2、对i2建立高斯金字塔g，g的层数为12，原图像作为第1层输入。需要注意的是，g中第2层到第6层是其下一层对列进行一维下采样后得到的，采样的步长为2；而第7层到第12层尺寸不变，采样的步长为1。采样用到的核f大小为1*5，f的元素值大小可以根据情况进行调整，这里取f＝[0.05,0.25,0.4,0.25,0.05]。

在本实施例中，i2的分辨率为2048*1536，而g中第12层的分辨率为2048*48。g只对列进行了下采样。

s3、由g建立对应于i2的拉普拉斯金字塔l，显然l中仅有关于列的细节函数l1～l11。建立l的具体过程如下：

(3-1)设g中第i层矩阵为gi，i＝1,2,……,n；n＝12。将gn利用双三次插值法恢复到gn-1的尺寸；得到g′n-1；

(3-2)ln-1＝gn-1-g′n-1；

(3-3)对g中除g1以外的所有元素进行步骤(3-1)至(3-2)中的操作，得到ln-1～l1。

s4、对l中的每一层进行去噪，包括以下几个步骤：

(4-1)将ln-1的全部元素取绝对值，得到|ln-1|，

(4-2)取|ln-1|的中值元素m0，令m1＝4*m0，

(4-3)对|ln-1|中的每个元素值p，计算v＝p/m1-1；若v＞1，则将v置1，若v＜0，将v置0。由此得到矩阵v；

(4-4)计算新的细节函数l′n-1＝ln-1.*v，其中“.*”符号表示矩阵对应位置元素相乘；

(4-5)对l中剩下的层进行步骤(4-1)到(4-4)中相同的操作，得到新的拉普拉斯金字塔l′。

s5、借助l′，由g的最顶层重建原始尺寸的图像，得到平滑后的图像i′。图4给出了实施例对应的i′图像，可以看到图像上已经完全滤除了明显的高频信号，即条纹噪声，但因为对细节函数的修改并不精细，此时图像i′的有效部分相比原图有效部分存在复原不准确和过度平滑的问题，因此i′不适合直接作为最终结果用于后续测量等操作。

在本发明所提出的算法中，i′被当作i1在增益设置正确的情况下，所应该呈现的结果，即期望图像。基于i′的平滑效果和原图每一列存在的统一的增益系数，将在后续步骤中对i1按列进行图像采集以后的第二次校准。

s6、计算矩阵y＝i′./i1，其中“./”符号表示矩阵对应位置元素相除。取y中每一列的中值构成行向量y。此步骤即选取附加的增益值，对应于本实施例的行向量y如图5所示。

s7、按列遍历i1，对i1中第i列的元素i1(i)，计算i1(i)*y(i)，由此得到新的矩阵i3。

s8、遍历y(i)，若|y(i)-1|＞0.02，0.02为根据经验设定的阈值，对i3的第i列的元素i3(i)进行中值滤波。滤波器的阶数为7。由此得到矩阵i4。

s9、i4即为本发明实施例方法的最终结果，对应于实施例的i4如图6所示，从图6可看到，相较于原始噪声图像，图像中的条纹噪声得到有效去除，噪声之外的信息得到保留，基本没有损伤，去噪效果很好，能够满足工程应用的需要。

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。