一种基于滤波反投影‐迭代算法的CT系统图像重建方法技术方案

本发明专利技术提供一种基于滤波反投影‐迭代算法的CT系统图像重建方法，涉及CT系统图像重建领域，包括如下步骤：第一步：利用CT系统对待测物质测量，用于获取待测物体的投影值数据；第二步：利用滤波函数对获取到的数据进行滤波处理，并利用反投影算法求解出待测物体的衰减系数矩阵；第三步：将反投影算法求得的衰减系数矩阵作为迭代算法的初始迭代值，得到修正的衰减系数矩阵，完成待测物体的图像重建。其能快速成像，具有较好的抗噪声性能，并能保持较好的成像性能。

An image reconstruction method for CT system based on filtered backprojection iterative algorithm

The invention provides an image reconstruction method of CT system based on filter backprojection iterative algorithm, which involves the following steps: the first step: using the CT system to treat the measured material measurement to obtain the projection value data of the object to be measured; the second step is to use the filter function to carry out the obtained data. The attenuation coefficient matrix of the object to be measured is solved by the back projection algorithm. The third step: the attenuation coefficient matrix obtained by the back projection algorithm is used as the initial iteration value of the iterative algorithm, and the modified attenuation coefficient matrix is obtained to complete the image reconstruction of the object to be measured. It can fast imaging, has better anti noise performance, and can maintain good imaging performance.

【技术实现步骤摘要】
一种基于滤波反投影‐迭代算法的CT系统图像重建方法
本专利技术涉及CT系统图像重建领域，具体为一种基于滤波反投影‐迭代算法的CT系统图像重建方法。
技术介绍
目前CT系统图像重建算法主要包括反投影算法和迭代算法。反投影算法应用广泛，成像速度快，可以实时检测，但是成像性能较差，图像分辨率较低。在噪声较大的环境中，还需要加入滤波函数减小噪声对成像结果的影响。迭代算法的成像效果好，图像分辨率高，但其抗噪声性能差，对环境的变化极其敏感。而且不同迭代算法的收敛速度不同，其中有些迭代算法的收敛速度很慢，导致成像速度达不到实际中的应用要求。通过研究发现，迭代算法收敛速度较慢的原因主要有两种：1.环境噪声使迭代过程出现误差，影响算法收敛速度。2.迭代算法的初始值设置不当，导致算法的迭代次数增多。实际应用中迫切需要对迭代算法进行改良，既能保持较好的成像性能，同时具有抗噪声性能好、成像速度快的特点。
技术实现思路
本专利技术提供一种基于滤波反投影-迭代算法的CT系统图像重建方法，其能快速成像，具有较好的抗噪声性能，并能保持较好的成像性能。为实现上述技术目的，本专利技术采取的技术方案为，一种基于滤波反投影-迭代算法的CT系统图像重建方法，包括如下步骤：第一步：利用CT系统对待测物质测量，用于获取待测物体的投影值数据；第二步：利用滤波函数对获取到的数据进行滤波处理，并利用反投影算法求解出待测物体的衰减系数矩阵；第三步：将反投影算法求得的衰减系数矩阵作为迭代算法的初始迭代值，得到修正的衰减系数矩阵，完成待测物体的图像重建。作为本专利技术改进的技术方案，CT系统包括射线发射源与探测器，射线发射源与探测器分别设于待测物质两侧，并且射线发射源与探测器对应设置。作为本专利技术改进的技术方案，射线发射源有若干个，探测器有若干个，若干个射线发射源在待测物质一侧等间距排列，若干个探测器在待测物质另一侧等间距排列；并且每个射线发射源对应一个探测器。作为本专利技术改进的技术方案，利用CT系统对待测物质测量时，所有射线发射源始终保持在一条直线上，所有探测器始终保持在一条直线上，射线发射源所在直线与探测器所在直线平行。作为本专利技术改进的技术方案，射线发射源发射的射线垂直于探测器平面，每个探测器能看成一个接收点，用于接收穿过待测物体的射线。作为本专利技术改进的技术方案，利用CT系统对待测物质测量时，发射器和探测器的相对位置保持不变，发射器和探测器同时绕设定的固定的旋转中心逆时针旋转角度β的范围为0-180°，旋转步长α的范围为0°至2°。作为本专利技术改进的技术方案，滤波函数选用R-L滤波函数。有益效果因此，本专利技术引入滤波函数，将反投影算法和迭代算法结合起来，成为滤波反投影-迭代算法。先利用滤波函数滤除环境噪声，并将反投影算法的求解结果作为迭代算法的迭代初始值，可以有效提高算法收敛速度。本申请采用滤波反投影-迭代算法弥补了迭代算法抗噪声性能差，成像速度慢的缺点，能够在噪声较大的环境中快速稳定的获得分辨率较高的物体重建图像，在医疗、工程等领域可以得到广泛应用。具体实施方式为使本申请实施例的目的和技术方案更加清楚，下面将结合本申请实施例对本申请的技术方案进行清楚、完整地描述。显然，所描述的实施例是本申请的一部分实施例，而不是全部的实施例。基于所描述的本申请的实施例，本领域普通技术人员在无需创造性劳动的前提下所获得的所有其他实施例，都属于本申请保护的范围。本
技术人员可以理解，除非另外定义，这里使用的所有术语(包括技术术语和科学术语)具有与本申请所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是，诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义，并且除非像这里一样定义，不会用理想化或过于正式的含义来解释。一种基于滤波反投影-迭代算法的CT系统图像重建方法，包括如下步骤：第一步：利用CT系统对待测物质测量，用于获取待测物体的投影值数据；第二步：利用滤波函数对获取到的数据进行滤波处理，并利用反投影算法求解出待测物体的衰减系数矩阵；滤波函数选用R-L滤波函数；第三步：将反投影算法求得的衰减系数矩阵作为Richardson迭代算法的初始迭代值，得到修正的衰减系数矩阵，完成待测物体的图像重建。其中，CT系统，CT系统用于获取待测物体的投影值数据；滤波函数，滤波函数用于对投影值进行滤波，去除噪声；反投影算法，反投影算法用于计算衰减系数矩阵，为迭代算法提供初始值；迭代算法，迭代算法用于提高衰减系数矩阵的精确度，提高重建图像的分辨率。作为本专利技术改进的技术方案，CT系统包括射线发射源与探测器，射线发射源与探测器分别设于待测物质两侧，并且射线发射源与探测器对应设置。作为本专利技术改进的技术方案，射线发射源有若干个，探测器有若干个，若干个射线发射源在待测物质一侧等间距排列，若干个探测器在待测物质另一侧等间距排列；并且每个射线发射源对应一个探测器。作为本专利技术改进的技术方案，利用CT系统对待测物质测量时，所有射线发射源始终保持在一条直线上，所有探测器始终保持在一条直线上，射线发射源所在直线与探测器所在直线平行。作为本专利技术改进的技术方案，射线发射源发射的射线垂直于探测器平面，每个探测器能看成一个接收点，用于接收穿过待测物体的射线。作为本专利技术改进的技术方案，利用CT系统对待测物质测量时，发射器和探测器的相对位置保持不变，发射器和探测器同时绕设定的固定的旋转中心逆时针旋转角度β的范围为0-180°，旋转步长α的范围为0°至2°。理论基础：传统的图像重建代数模型，将待测图像分成若干个像素，不同像素可以对应不同的衰减系数，并合理的假定每个衰减系数为常数。对于m个像素(记j＝1,2,L,m)，n束射线(记i＝1,2,…,n)的情况，用I0表示射线的入射强度，μj表示像素j的衰减系数，Δlj表示射线在j中的穿行长度，可以得到射线衰减规律等式：其中，I为每个探测器的实际探测值，J(Li)表示射线束Li穿过的像素j的集合，λ为修正参数，由CT系统内部结构决定，可根据实际探测情况确定为一个常数。考虑到数据处理的难度，在实际应用中为了简单起见，常常将μj记为xj，将λln(I0/I)记为bi，aij代表射线束Li在像素j内的长度lij与像素边长之比。利用中心法可以将射线的宽度简化为0，当射线束Li穿过像素j时，记aij＝1，否则aij＝0。因此，可写成向量-矩阵的形式：AX＝b式二其中X＝(x1,x2,…,xm)T，A＝(aij)n×m，b＝(b1,b2,L,bn)T。实验数据获取：目前国内外CT系统探测方法主要采用等距平行束探测或者等角扇形束探测，二者的探测原理基本相同，本申请以平行束为例简要介绍探测数据的获取过程。在扫描过程中，射线发射源和探测器始终保持在一条直线上。平行入射的射线垂直于探测器平面，每个探测器可看成一个接收点，且等距排列。发射器和探测器的相对位置保持不变，整个CT系统绕某固定的旋转中心逆时针旋转β(0°≤β≤180°)。令旋转步长为α(0°<α<2°)，则可以获取βα组探测数据。滤波反投影算法过程：步骤：1.可设像素坐标为(x,y)，该像素的衰减系数为μ(x,y)。通过探测器可获得在旋转角度下测得的强度数据，将其作一维傅立叶变换后可得到本文档来自技高网...

【技术保护点】
1.一种基于滤波反投影‑迭代算法的CT系统图像重建方法，其特征在于，包括如下步骤：第一步：利用CT系统对待测物质测量，用于获取待测物体的投影值数据；第二步：利用滤波函数对获取到的数据进行滤波处理，并利用反投影算法求解出待测物体的衰减系数矩阵；第三步：将反投影算法求得的衰减系数矩阵作为迭代算法的初始迭代值，得到修正的衰减系数矩阵，完成待测物体的图像重建。

【技术特征摘要】
1.一种基于滤波反投影-迭代算法的CT系统图像重建方法，其特征在于，包括如下步骤：第一步：利用CT系统对待测物质测量，用于获取待测物体的投影值数据；第二步：利用滤波函数对获取到的数据进行滤波处理，并利用反投影算法求解出待测物体的衰减系数矩阵；第三步：将反投影算法求得的衰减系数矩阵作为迭代算法的初始迭代值，得到修正的衰减系数矩阵，完成待测物体的图像重建。2.根据权利要求1所述的一种基于滤波反投影-迭代算法的CT系统图像重建方法，其特征在于，CT系统包括射线发射源与探测器，射线发射源与探测器分别设于待测物质两侧，并且射线发射源与探测器对应设置。3.根据权利要求2所述的一种基于滤波反投影-迭代算法的CT系统图像重建方法，其特征在于，射线发射源有若干个，探测器有若干个，若干个射线发射源在待测物质一侧等间距排列，若干个探测器在待测物质另一侧等间距排列；并且每个射线发射源对应一个探测器。4.根...

【专利技术属性】
技术研发人员：林勇康，马意彭，赵江，毛诗涵，
申请(专利权)人：南京邮电大学，
类型：发明
国别省市：江苏,32