结合图像平滑处理的CT图像环形伪影抑制方法与流程

本发明属于医学图像处理技术领域，涉及结合图像平滑处理的ct图像环形伪影抑制方法。

背景技术：

图像伪影是ct图像在重建过程中受到的一些干扰信号产生的或者是原来物体中所不存在的影像。ct图像的伪影常常表现为是环状伪影。图像质量的下降正是由于这些环状伪影的出现会导致，从而造成一些误判的发生。虽然环形伪影的出现只是图像中的一个干扰，但是它如果出现在一些研究方向的重要位置上就可能会产生比较严重的问题。尤其在医学上，环状伪影的出现会影响ct图像中的一些病灶处的影像，甚至还会掩盖一些有效区域，从而导致医生对疾病的诊断，造成一些重大的医疗事故。而且ct设备是一个很复杂的系统，由于使用不当、保养不及时或未及时利用标准模型校准等均可引起一些伪影。由这些因素所导致的环形伪影通常很容易被看出，但如果所造成的伪影不明显或是由于一些机器伤影而产生的时，很容易把它当成图像本身的特点而造成误判,如当环形伪影出现在脊髓中央时，会导致诊断错误等。此外，在环形伪影中有一个比较特殊表现形式，即它在图像中表现为环的中央模糊，这种环形伪影一般表现在小部分范围内的。又由于它的表现形态与病理结构很相似，从而可能导致误判。

技术实现要素：

本发明针对现有的技术存在的上述问题，提供结合图像平滑处理的ct图像环形伪影抑制方法，本发明所要解决的技术问题是：如何提供一种可以结合图像平滑处理ct图像以对ct图像中的环形伪影进行抑制的ct图像环形伪影抑制方法。

本发明的目的可通过下列技术方案来实现：

结合图像平滑处理的ct图像环形伪影抑制方法，包括如下步骤：

s1：将直角坐标下的第一环形图像转化为极坐标下的第一条形图像，通过图像纹理平滑处理方法将第一条形图像处理为第二条形图像；

s2：将极坐标下的第二条形图像转化为直角坐标下的第二环形图像，通过图像纹理平滑处理方法将第二环形图像处理为环形伪影抑制后的第三环形图像。

优选的，步骤s1中第一环形图像转化为第一条形图像，第一环形图像中的第一环形伪影转化为第一条形伪影，通过对第一条形图像进行图像纹理平滑处理以将第一条形图像中的第一条形伪影消除。

优选的，步骤s2中第二条形图像转化第二环形图像，第二条形图像中的第二条形伪影转化为第二环形伪影，通过对第二环形图像进行图像纹理平滑处理以将第二环形图像中的第二环形伪影消除。

优选的，图像纹理平滑处理过程具体包括：将图像窗口的固有变差和图像窗口总变差结合起来，形成结构-纹理分解正则化器，结构-纹理分解正则化器的目标函数为，其中使得输入和输出的结果图像不发生变化，正则项为，正则项为图像窗口的相对总变差，为权重值，像素级窗口总变差分别为、，窗口固有变化分别为、，为根据空间相关关系定义的权重函数，，为图像窗口的固有变化，大于0。

优选的，步骤s1中直角坐标下的第一环形图像的边长为，极坐标下的第一条形图像的极径为，第一条形图像的角度为。

优选的，步骤s1中直角坐标下的第一环形图像转化为极坐标下的第一条形图像时对于超出第一环形图像的采样点通过最邻近插值法确定其超出预设范围的值。

优选的，步骤s2中极坐标下的第二条形图像转化为直角坐标下的第二环形图像时对于超出第二条形图像的位置通过最邻近插值法进行填补。

优选的，步骤s1第一环形图像为同心圆，第一环形图像转化为第一条形图像时，第一条形图像的中心设置于极坐标原点，直角坐标下的同心圆经过变换为极坐标下的一条平行于轴的直线。

优选的，参数的大小为0~0.05之间。

优选的，参数的值大于0.09。

本发明中首先将直角坐标下的第一环形图像转化为极坐标下的第一条形图像，通过图像纹理平滑处理方法将第一条形图像处理为第二条形图像，第一环形图像中的部分环形伪影的程度得到减弱，接着将极坐标下的第二条形图像转化为直角坐标下的第二环形图像，通过图像纹理平滑处理方法将第二环形图像处理为环形伪影抑制后的第三环形图像，未能被完全抑制的环形伪影能够有效得到抑制，能够在很好的保护图像细节和图像清晰度的情况下，能够对图像中的环形伪影有比较好的抑制作用。

附图说明

图1是本发明的流程示意图。

具体实施方式

以下是本发明的具体实施例并结合附图，对本发明的技术方案作进一步的描述，但本发明并不限于这些实施例。

请参阅图1，本实施例中的结合图像平滑处理的ct图像环形伪影抑制方法，包括如下步骤：

s1：将直角坐标下的第一环形图像转化为极坐标下的第一条形图像，通过图像纹理平滑处理方法将第一条形图像处理为第二条形图像；

s2：将极坐标下的第二条形图像转化为直角坐标下的第二环形图像，通过图像纹理平滑处理方法将第二环形图像处理为环形伪影抑制后的第三环形图像。

此处，首先将直角坐标下的第一环形图像转化为极坐标下的第一条形图像，通过图像纹理平滑处理方法将第一条形图像处理为第二条形图像，第一环形图像中的部分环形伪影的程度得到减弱，接着将极坐标下的第二条形图像转化为直角坐标下的第二环形图像，通过图像纹理平滑处理方法将第二环形图像处理为环形伪影抑制后的第三环形图像，未能被完全抑制的环形伪影能够有效得到抑制，能够在很好的保护图像细节和图像清晰度的情况下，能够对图像中的环形伪影有比较好的抑制作用。

步骤s1中第一环形图像转化为第一条形图像，第一环形图像中的第一环形伪影转化为第一条形伪影，通过对第一条形图像进行图像纹理平滑处理以将第一条形图像中的第一条形伪影消除。

步骤s2中第二条形图像转化第二环形图像，第二条形图像中的第二条形伪影转化为第二环形伪影，通过对第二环形图像进行图像纹理平滑处理以将第二环形图像中的第二环形伪影消除。

图像纹理平滑处理过程可以具体包括：将图像窗口的固有变差和图像窗口总变差结合起来，形成结构-纹理分解正则化器，结构-纹理分解正则化器的目标函数为，其中使得输入和输出的结果图像不发生变化，正则项为，正则项为图像窗口的相对总变差，为权重值，像素级窗口总变差分别为、，窗口固有变化分别为、，为根据空间相关关系定义的权重函数，，为图像窗口的固有变化，大于0，是一个正数，为了避免被除数出现为0的情况。只包含纹理的窗口的结果的通常比包含结构边缘的窗口的结果的小。

步骤s1中直角坐标下的第一环形图像的边长可以为，极坐标下的第一条形图像的极径可以为，第一条形图像的角度为。

步骤s1中直角坐标下的第一环形图像转化为极坐标下的第一条形图像时对于超出第一环形图像的采样点可以通过最邻近插值法确定其超出预设范围的值。

步骤s2中极坐标下的第二条形图像转化为直角坐标下的第二环形图像时对于超出第二条形图像的位置可以通过最邻近插值法进行填补。

步骤s1第一环形图像可以为同心圆，第一环形图像转化为第一条形图像时，第一条形图像的中心设置于极坐标原点，直角坐标下的同心圆经过变换为极坐标下的一条平行于轴的直线。

参数的大小可以为0~0.05之间，这样可以保护图像中想细节部分。

参数的值可以大于0.09，这样可以保护图像中的细节部分。在对数据图像处理时，参数在大于0.025和参数在小于0.09时，在对极坐标下处理后的图像再在直角坐标下处理中，图像会变得很模糊，甚至图像中的细节部分也消失了，反之，图像中的环形伪影处理效果不是特别的好。对数据图像处理时，参数在大于0.05和参数在小于0.03时，在对极坐标下处理后的图像再在直角坐标下处理中，图像会变得很模糊，甚至图像中的细节部分也消失了，反之，图像中的环形伪影处理效果不是特别的好。

本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。