[发明专利]一种利用投影数据去除CT探测元故障造成的环状伪影的方法

摘要

本发明涉及一种利用投影数据去除CT探测元故障造成的环状伪影的方法，该方法包括如下步骤：(1)导入投影数据得到矩阵A，检测矩阵A中故障区域的位置，用故障区域左右两列外的两列数据对故障区域所对应列进行线性插值得到矩阵B；(2)对矩阵B进行滤波反投影重建后再进行均值滤波、前向投影，得到投影正弦矩阵C；(3)用矩阵A减去矩阵C，再以差值矩阵中故障区域所对应的列外的左右列数据对故障区域所对应的列进行线性插值处理，得到矩阵D；(4)将矩阵C与矩阵D相加得到矩阵E；(5)返回到步骤(2)，不断循环步骤(2)～(4)，直至达到预设的迭代次数后，再对所得到的结果进行滤波反投影重建，得到去除伪影后的CT图像。本方法可消除CT探测元连续大面积故障所产生的环形伪影。 1

主权项

1.一种利用投影数据去除CT探测元故障造成的环状伪影的方法，该方法包括如下步骤：(1)导入圆轨道CT成像系统采集的投影数据，得到投影正弦矩阵A；以Sobel直线检测方法或Canny直线检测方法检测连续故障的CT探测元所对应的故障区域在矩阵A中的位置，并用其左右两列外的两列数据对故障区域所对应列进行线性插值，得到初步校正后的投影正弦矩阵B；(2)对所述的矩阵B进行滤波反投影重建得到矩阵image1，然后对矩阵image1进行均值滤波得到矩阵image2，再对矩阵image2进行前向投影，得到投影正弦矩阵C；其中，所述的均值滤波方法由下式(Ⅰ)所示：上式(Ⅰ)中，(i，j)表示矩阵image2中像素点位置；(i’,j’)表示满足以像素位置(i，j)为中心长宽分别为(2v+1)的矩形窗内的像素点位置，且，i’＝i‑v，…，i+v，j’＝j‑v，…，j+v，v为1～10的自然数，随伪影的宽度按正比例选取；w(i’,j’)表示image1(i’,j’)对image1(i,j)的权重；K表示以(i，j)为中心的邻域窗大小为(2v+1)*(2v+1)的范围内满足|image1(i',j')‑image1(i,j)|≤T像素点个数；T为预设阈值；(3)先用所述的矩阵A减去矩阵C，再以差值矩阵中故障区域所对应的列外的左右列数据对故障区域所对应的列进行线性插值处理，得到经差值运算和插值处理后的投影正弦矩阵D；(4)将所述的矩阵C与矩阵D相加，得到求和运算后的投影正弦矩阵E；(5)返回到步骤(2)，并用步骤(4)得到的矩阵E替代步骤(1)得到的矩阵B，不断循环步骤(2)～(4)，直至达到预设的迭代次数后，再对所得到的结果进行滤波反投影重建，得到去除伪影后的CT图像。