[发明专利]CT图像的环形伪影及弧形伪影校正方法

权利要求书

1.一种CT图像的环形伪影校正方法，所述CT图像来自以m个扫描投影角度θ_i扫描所得的m幅分辨率为n×n的图像，其中，0°＜θ_i＜360°，1≤i≤m且为正整数，所述CT图像的环形伪影校正方法包括以下步骤：

将所有的CT图像变换到以探测器像元通道n为横坐标、投影角度θ为纵坐标的直角坐标系中，获得n×m个分辨率为n×1的第一中间图像R₁(n，θ_i)；

对第一中间图像R₁(n，θ)依次进行n列的像素各自叠加求和、n列的像素叠加求和后取平均值等数据处理后得到数据E₁(n)；

将所述第一中间图像R₁(n，θ)n列的像素分别叠加求和后各自减去所述数据E₁(n)，即可得到校正环形伪影后n×m个分辨率为n×1的第二中间图像R₂(n，θ)；

对所述第二中间图像R₂(n，θ)进行锥形束CT图像重建，即可得到校正环形伪影后的CT图像。

2.根据权利要求1所述的CT图像的环形伪影校正方法，其特征在于，所述扫描投影角度θ_i通过等分圆周得到。

3.根据权利要求2所述的CT图像的环形伪影校正方法，其特征在于，所述扫描投影角度θ_i通过将圆周等分360份得到。

4.根据权利要求1所述的CT图像的环形伪影校正方法，其特征在于，所述将所有的CT图像变换到以探测器像元通道n为横坐标、投影角度θ为纵坐标的直角坐标系中，获得n×m个分辨率为n×1的第一中间图像R₁(n，θ_i)的步骤包括：

将二维面阵探测器探测的数据转换成一系列的一维线阵探测器探测的数据；

将每个一维线阵探测器在扫描过程中获得的所有CT图像变换到直角坐标系中。

5.根据权利要求1所述的CT图像的环形伪影校正方法，其特征在于，所述对第一中间图像R₁(n，θ)依次进行n列的像素各自叠加求和、n列的像素叠加求和后取平均值等数据处理后得到数据E₁(n)的步骤包括：

对第一中间图像R₁(n，θ)n列的像素各自进行叠加求和，得到数据S(n)；

对所述数据S(n)进行高通滤波；

对所述数据S(n)进行n列的像素叠加求和后取平均值，得到数据E₁(n)。

6.一种CT图像的弧形伪影校正方法，所述CT图像来自以m个扫描投影角度θ_i扫描所得的m幅分辨率为n×n的图像，其中，0°＜θ_i＜360°，1≤i≤m且为正整数，所述CT图像的弧形伪影校正方法包括以下步骤：

将所有的CT图像变换到以探测器像元通道n为横坐标、投影角度θ为纵坐标的直角坐标系中，获得n×m个分辨率为n×1的第一中间图像R₁(n，θ)；

对第一中间图像R₁(n，θ)依次进行n列的像素各自叠加求和、n列的像素叠加求和后取平均值等数据处理后得到数据E₁(n)；

将所述第一中间图像R₁(n，θ)n列的像素分别叠加求和后各自减去所述数据E₁(n)，即可得到校正部分弧形伪影后的第二中间图像R₂(n，θ)；

将第二中间图像R₂(n，θ)分为m个子集对每个子集依次进行n列的像素各自叠加求和、高通滤波提取高频成分、n列的像素叠加求和后取平均值等数据处理得到数据

将所述m个子集各行依次减去对应的平均值即可得到完全校正弧形伪影后第三中间图像R₃(n，θ)；

对所述第三中间图像R₃(n，θ)进行锥形束CT图像重建，即可得到校正弧形伪影后的CT图像。

7.根据权利要求6所述的CT图像的弧形伪影校正方法，其特征在于，所述扫描投影角度θ_i通过等分圆周得到。