基于重建前后图像的环形伪影校正方法、CT控制系统技术方案

本发明专利技术属于影像处理技术领域，公开了一种基于重建前后图像的环形伪影校正方法、CT控制系统；通过对投影数据进行预处理；计算沿正弦图中所有视角的每个探测器元素的平均值得到一维校正矢量与正弦图在所有视角上逐行相乘；计算沿正弦图中每个角度的平均值得到一维校正向量与正弦图逐列相乘；对重建后原始CT图像与获得的校正CT图像之间的差异图像应用高斯滤波；滤波后的差异图像被添加到校正的CT图像实现对比度补偿。本发明专利技术根据所使用的基于光子计数探测器的能谱CT系统所存在的信号不稳定的缺陷，引入多步算法，更适合能谱CT图像数据的图像特点并使图像环形伪影得到了很大的改善。

Ring artifact correction method and CT control system based on reconstructed images

The invention belongs to the field of image processing technology, and discloses a ring artifact correction method and a CT control system based on reconstructed images; preprocessing projection data; calculating the average value of each detector element along the sinusoidal graph to get one-dimensional correction vector multiplied by the sinusoidal graph row by row in all views; calculating the average value of each angle along the sinusoidal graph. It is worth multiplying one-dimensional correction vectors and sinusoids one by one. Gauss filtering is applied to the difference image between the reconstructed original CT image and the corrected CT image. The filtered difference image is added to the corrected CT image to realize contrast compensation. According to the defect of unstable signal in the energy spectrum CT system based on photon counting detector, the multi-step algorithm is introduced, which is more suitable for the image characteristics of the energy spectrum CT image data and greatly improves the annular artifact of the image.

【技术实现步骤摘要】
基于重建前后图像的环形伪影校正方法、CT控制系统
本专利技术属于影像处理
，尤其涉及一种基于重建前后图像的环形伪影校正方法、CT控制系统。
技术介绍
目前，业内常用的现有技术是这样的：环形伪影通常出现在计算机断层扫描(CT)系统采集的图像中并以同心圆的形式存在，严重降低了获得的图像质量。环形伪影可能由多种原因导致，其中包括存在缺陷的探测器像素、像素的非线性响应、探测器增益的不一致性、探测器和X射线源的结构模式、随时间和温度条件的系统性变化等因素。关于环形伪影校正方法已有大量的研究，可以在正弦域和CT图像域上进行环状伪影校正。在正弦域中，像素之间响应的不一致性表现为竖条纹，各个角度间整个探测器的响应的不一致性表现为横条纹，这些条纹的存在会直接造成重建后的环形伪影，影响图像的质量。因此大多数正弦图像阈方法估计利用正弦图中的条纹的不均匀性，此方法对环状伪影有不错的抑制效果，但是校正不彻底，并且容易带来新的伪影。在重建后的图像域中用到了坐标变换，无论是从极坐标系转换到直角坐标系，还是其逆变换，其中都要用到插值运算，插值算法本身就是通过周围数据估算的不准确数据，这样不仅在坐标变换上就已经导致了两次丢失原始图像信息的情况，图像分辨率必然会下降。如果在极坐标中平滑滤波时采用滤波模板，也会对所要校正的环形伪影个数有一定的局限性。关于锥束CT二维投影图像上的环形伪影校正方法很少。近期也有研究者利用重建前后图像的环形伪影校正方法，采用两步的方法进行环形伪影的校正，因为光子计数型探测器由于工艺的不成熟性，会导致像素间相应不一致性更加明显，每个角度间整个探测器的相应特性不平稳，在正弦域中存在大量的竖条纹的基础上又存在大量的横条纹，严重影响了重建后图片的质量。此种方法对于基于光子计数探测器的能谱CT效果欠佳。综上所述，现有技术存在的问题是：(1)在正弦域中，不均匀性表现为条纹，因此大多数正弦图像阈方法估计利用正弦图中的条纹的不均匀性。(2)利用重建前后图像的环形伪影校正方法，采用两步的方法进行环形伪影的校正，对于基于光子计数探测器的能谱CT效果欠佳。解决上述技术问题的难度和意义：如果在正弦域中进行校正，由于光子计数型探测器各个角度间整个探测器的响应的不一致性所表现出来的大量横条纹，不能针对光子计数型探测器的能谱CT图像得到很好的改善，并且此种方法校正不彻底，并且容易带来新的伪影。如果在重建后的图像域中用到了坐标变换，无论是从极坐标系转换到直角坐标系，还是其逆变换，其中都要用到插值运算，插值算法本身就是通过周围数据估算的不准确数据，这样不仅在坐标变换上就已经导致了两次丢失原始图像信息的情况，图像分辨率必然会下降。如果在极坐标中平滑滤波时采用滤波模板，也会对所要校正的环形伪影个数有一定的局限性。因此，如果能够结合重建前的正弦域图像和重建后的图像域图像的特点，再结合光子计数型探测器所存在的特征，设计出一种能够有效的基于光子计数型的能谱CT的环形伪影校正方法，也能推动基于光子计数型探测器的能谱CT的发展和应用。
技术实现思路
针对现有技术存在的问题，本专利技术提供了一种基于重建前后图像的环形伪影校正方法、CT控制系统。本专利技术是这样实现的，一种基于重建前后图像的环形伪影校正方法，所述基于重建前后图像的环形伪影校正方法包括：步骤一，对投影数据进行预处理，纠正投影域中值接近零或饱和的缺陷像素，平稳投影域中由于光子计数探测器的每一个小探测器间的相应不一致性所造成的信号浮动，去除光子计数探测器的每一个小探测器间间存在的高亮缝隙；步骤二，对于对应于一个切片的单行检测器，计算沿正弦图中视角的每个探测器元素的平均值，以获得一维平均向量，通过取平均向量的逆来计算一维校正向量，并将校正矢量与正弦图在所有视角上逐行相乘；根据探测器相应的不稳定性，计算沿正弦图中每个角度的所有行探测器元素平均值，获得一个一维平均向量，然后通过取平均向量的逆也计算一个一维校正向量，并将校正矢量与正弦图逐列相乘；步骤三，对原始CT图像与上一步获得的校正CT图像之间的差异图像应用高斯滤波，使用沿行方向的中值滤波器对差异图进行模糊；滤波后的差异图像被添加到校正的CT图像以补偿由于正弦图中的对比度变化出现的可能的对比度异常。进一步，对投影图像进行特殊的预处理包括：(1)对采集到的投影图像进行有效区域裁剪，仿体实验中采取663×256裁剪大小；(2)对裁剪后的投影数据进行暗场校正和空场校正；(3)对暗场校正和空场校正后的投影数据进行坏点坏线校正；(4)对高亮缝隙进行补偿，利用图像高亮缝隙的右侧一列像素的非正常特性，采取右2补偿右1，高亮像素由左1和右1取平均来补偿；图像数据为I0，处理后的图像数据为Ie，则对于缝隙t处的校正方法如下公式：I'e(t+1,:)＝I0(t+2,:)；Ie(t,:)＝(Ie(t-1,:)+Ie(t+1,:))*0.5；(5)对每个小块间的相应不一致性所导致的块与块之间的偏差采取整流的方法进行补偿，利用空场校正的思想对每一个小模块取平均值，每个小模块中的所有像素除以自身小模块的平均值，再乘整张图片的平均值以达到整流的效果。具体实现方法如下公式，其中:第t小块，m:第t小块的所有像素的平均值，M：整张人图片中所有像素的平均值：Iz＝I(t)/m*M；(6)对各个图片间存在的相应不一致性，将校正方法融入到环形伪影校正中，通过在正弦域的图像数据中，计算沿正弦图中每个角度的所有行探测器元素平均值，也获得一个一维平均向量，然后通过取平均向量的逆也计算一个一维校正向量，并将校正矢量与正弦图逐列相乘消除角度间图片的相应不一致造成的偏差，再利用后期的对比度补偿来补偿由于平均后所造成的对比度损失。进一步，利用分步的方法分别在正弦域和重建后的图像域分别进行处理，在正弦域进行半校正处理，在重建后的图像域进行对比度补偿。进一步，在正弦域中对于对应于一个切片的单行检测器，计算沿正弦图中视角的每个探测器元素的平均值，以获得一维平均向量，然后通过取平均向量的逆来计算一维校正向量，并将校正矢量与正弦图在所有视角上逐行相乘；根据探测器相应的不稳定性，计算沿正弦图中每个角度的所有行探测器元素平均值，也获得一个一维平均向量，然后通过取平均向量的逆也计算一个一维校正向量，并将校正矢量与正弦图逐列相乘。进一步，在正弦图域中，X射线探测器逐行2D进行校正，在探测器的第j行处，通过在不同视角处从每个投影图像采取相同的行来制作正弦图Sj(i,k)＝S’(i,j,k)k＝1,2,3...其中S’是坏点坏线校正后的投影图像；然后，在所有视角上计算第i个像素处的平均值以获得1D平均矢量Mj；在所有视角上计算第i个像素处的平均值以获得1D平均矢量：校正向量Mj',Mk'如下：经过校正投影域数据：在所有列上计算第k个角度处的平均值以获得1D平均矢量，并计算校正向量对进行校正；进一步，对重建后的图像进行差分图的计算，选择高斯滤波器以平滑差分图像，高斯滤波器大小为15×15，标准偏差为10个像素，具体实现方法如下:D＝D1-D2；以平滑差高斯滤波差分图像和半校正图像的总和得到最终图像：本专利技术的另一目的在于提供一种包含所述基于重建前后图像的环形伪影校正方法的锥束CT成像平台。本专利技术的另一目的在于提供一种包含所述基于重本文档来自技高网...

【技术保护点】
1.一种基于重建前后图像的环形伪影校正方法，其特征在于，所述基于重建前后图像的环形伪影校正方法包括：步骤一，对投影数据进行预处理，纠正投影域中值接近零或饱和的缺陷像素，平稳投影域中由于光子计数探测器的每一个小探测器间的相应不一致性所造成的信号浮动，去除光子计数探测器的每一个小探测器间间存在的高亮缝隙；步骤二，对于对应于一个切片的单行检测器，计算沿正弦图中视角的每个探测器元素的平均值，以获得一维平均向量，通过取平均向量的逆来计算一维校正向量，并将校正矢量与正弦图在所有视角上逐行相乘；根据探测器相应的不稳定性，计算沿正弦图中每个角度的所有行探测器元素平均值，获得一个一维平均向量，然后通过取平均向量的逆也计算一个一维校正向量，并将校正矢量与正弦图逐列相乘；步骤三，对原始CT图像与上一步获得的校正CT图像之间的差异图像应用高斯滤波，使用沿行方向的中值滤波器对差异图进行模糊；滤波后的差异图像被添加到校正的CT图像以补偿由于正弦图中的对比度变化出现的可能的对比度异常。

【技术特征摘要】
1.一种基于重建前后图像的环形伪影校正方法，其特征在于，所述基于重建前后图像的环形伪影校正方法包括：步骤一，对投影数据进行预处理，纠正投影域中值接近零或饱和的缺陷像素，平稳投影域中由于光子计数探测器的每一个小探测器间的相应不一致性所造成的信号浮动，去除光子计数探测器的每一个小探测器间间存在的高亮缝隙；步骤二，对于对应于一个切片的单行检测器，计算沿正弦图中视角的每个探测器元素的平均值，以获得一维平均向量，通过取平均向量的逆来计算一维校正向量，并将校正矢量与正弦图在所有视角上逐行相乘；根据探测器相应的不稳定性，计算沿正弦图中每个角度的所有行探测器元素平均值，获得一个一维平均向量，然后通过取平均向量的逆也计算一个一维校正向量，并将校正矢量与正弦图逐列相乘；步骤三，对原始CT图像与上一步获得的校正CT图像之间的差异图像应用高斯滤波，使用沿行方向的中值滤波器对差异图进行模糊；滤波后的差异图像被添加到校正的CT图像以补偿由于正弦图中的对比度变化出现的可能的对比度异常。2.如权利要求1所述的基于重建前后图像的环形伪影校正方法，其特征在于，对投影图像进行特殊的预处理包括：(1)对采集到的投影图像进行有效区域裁剪，仿体实验中采取663×256裁剪大小；(2)对裁剪后的投影数据进行暗场校正和空场校正；(3)对暗场校正和空场校正后的投影数据进行坏点坏线校正；(4)对高亮缝隙进行补偿，利用图像高亮缝隙的右侧一列像素的非正常特性，采取右2补偿右1，高亮像素由左1和右1取平均来补偿；图像数据为I0，处理后的图像数据为Ie，则对于缝隙t处的校正方法如下公式：I'e(t+1,:)＝I0(t+2,:)；Ie(t,:)＝(Ie(t-1,:)+Ie(t+1,:))*0.5；(5)对每个小块间的相应不一致性所导致的块与块之间的偏差采取整流的方法进行补偿，利用空场校正的思想对每一个小模块取平均值，每个小模块中的所有像素除以自身小模块的平均值，再乘整张图片的平均值以达到整流的效果；具体实现方法如下公式，其中:第t小块，m:第t小块的所有像素的平均值，M：整张人图片中所有像素的平均值：Iz＝I(t)/m*M；(6)对各个图片间存在的相应不一致性，将校正方法融入到环形伪影校正中，通过在正弦域的图像数据中，计算沿正弦图中每个角度的所有行探测器元素平均值，也获得一个一维平均向...

【专利技术属性】
技术研发人员：朱守平，闫秋娟，时煜，孟凡珍，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西,61