[发明专利]一种X射线图像中滤线栅伪影的抑制方法及装置

说明书

技术领域

本发明涉及一种医学图像处理技术，尤其是涉及一种X射线图像中滤线栅伪影的抑制方法，以及使用该抑制方法的装置。

背景技术

X射线成像系统中，不可避免的会存在散射线，散射线会使拍摄的图像整体对比度降低，影响图像的质量。滤线栅的作用是将散射线在到达图像前尽可能的去掉从而减少图像上的散射线量，改善图像的对比度，提高图像的影像质量。但当栅纹比率相对于所搭配的探测器分辨率来说较低时，往往会在图像上出现不同程度的肉眼可察觉的条状伪影，称之为滤线栅伪影。滤线栅伪影的存在，会影响正常的组织显示，进而影响医生对病症的诊断。目前能够搜索到的近似相关专利有：

中国专利申请CN201010154719.2提出一种滤线栅伪影检测及抑制方法及其装置，此方法利用傅里叶变换在频域中检测滤线栅伪影的频率并用插值的方法处理检测出的频率。此方法通过提取较小区域的背景图像虽然加快了傅里叶变换和检测滤线栅伪影的速度，但是对原图进行傅里叶变换时仍然是一项很费时的工作，并且检测滤线栅伪影频率和插值处理的精度对得到的抑制滤线栅伪影后图像的影响比较大，在抑制滤线栅伪影的同时容易造成新的伪影。

中国专利CN03115184.1提出一种像素式投影显示器的图像栅格去除方法及系统，此方法采用小尺度位移叠合去除像素间的栅格。此方法涉及到的栅格主要是有像素之间的间隔形成的，并且本方法对于系统器件有一定的要求，故本方法不适合应用于X射线成像系统。

目前常用的抑制滤线栅伪影的方法为：采用傅里叶变换到频域，在频域中检测出滤线栅伪影的频率，采用插值或直接滤除的方法对检测出的频率进行处理，然后傅里叶逆变换，得到抑制伪影后的图像。但是，傅里叶变换是非常费时的操作，并且频域处理对图像的整体影响很大，检测滤线栅伪影频率和对检测出的频率进行处理时出现微小的差错都会对整幅图像造成比较大的影响，甚至会出现新的伪影，因此，现有技术有待进一步的改进和完善。

发明内容

本发明提出一种X射线图像中滤线栅伪影的抑制方法及装置，可以有效抑制X射线图像中滤线栅伪影所形成的视觉干扰，从而方便医生诊断，提高诊断效率。

本发明采用如下技术方案实现：一种X射线图像中滤线栅伪影的抑制方法，其包括步骤：

A、在垂直于原始的X射线图像的滤线栅的栅纹方向上对原始的X射线图像进行一维高通滤波，得到第一栅纹图像；

B、去除第一栅纹图像中组织边缘的图像，得到第二栅纹图像；

C、根据连续栅纹的长度特征确定，在第二栅纹图像确定并保留滤线栅的栅纹，得到第三栅纹图像；

D、根据原始的X射线图像在像素点(x，y)处的灰度值确定系数矩阵T，将第三栅纹图像与系数矩阵T对应点相乘得到第四栅纹图像；

E、在原始的X射线图像中减去第四栅纹图像，得到抑制滤线栅伪影后的X射线图像。

其中，步骤B具体包括：

计算第一栅纹图像中像素点的灰度最大值Max，设定阈值TH1＝a*MAX，其中a为介于0～1之间的常数；

若第一栅纹图像中像素点(x，y)的灰度值B(x，y)＞TH1，则在第一栅纹图像中将其去除，得到第二栅纹图像。

其中，a＝0.2。

其中，步骤B具体包括：

根据提取出的栅纹长度特性，设定阈值TH2；

在第二栅纹图像中连续栅纹的长度记为L(x)，假如L(X)＜TH2，则在第二栅纹图像中将其去除，得到第三栅纹图像。

其中，TH2＝5。

其中，根据X射线图像A的像素点(x，y)处的灰度值，利用下式确定系数矩阵T：