利索能及
我要发布
收藏
专利号: 2011103766108
申请人: 中国科学院深圳先进技术研究院
专利类型:发明专利
专利状态:已下证
更新日期:2024-12-09
缴费截止日期: 暂无
联系人

摘要:

权利要求书:

1.一种CT图像的环形伪影校正方法,所述CT图像来自以m个扫描投影角度θi扫描所得的m幅分辨率为n×n的图像,其中,0°<θi<360°,1≤i≤m且为正整数,所述CT图像的环形伪影校正方法包括以下步骤:将所有的CT图像变换到以探测器像元通道n为横坐标、投影角度θ为纵坐标的直角坐标系中,获得n×m个分辨率为n×1的第一中间图像R1(n,θi);

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

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

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

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

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

4.根据权利要求1所述的CT图像的环形伪影校正方法,其特征在于,所述将所有的CT图像变换到以探测器像元通道n为横坐标、投影角度θ为纵坐标的直角坐标系中,获得n×m个分辨率为n×1的第一中间图像R1(n,θi)的步骤包括:将二维面阵探测器探测的数据转换成一系列的一维线阵探测器探测的数据;

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

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

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

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

6.一种CT图像的弧形伪影校正方法,所述CT图像来自以m个扫描投影角度θi扫描所得的m幅分辨率为n×n的图像,其中,0°<θi<360°,1≤i≤m且为正整数,所述CT图像的弧形伪影校正方法包括以下步骤:将所有的CT图像变换到以探测器像元通道n为横坐标、投影角度θ为纵坐标的直角坐标系中,获得n×m个分辨率为n×1的第一中间图像R1(n,θ);

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

将所述第一中间图像R1(n,θ)n列的像素分别叠加求和后各自减去所述数据E1(n),即可得到校正部分弧形伪影后的第二中间图像R2(n,θ);

将第二中间图像R2(n,θ)分为m个子集 对每个子集 依次进行n列的像素各自叠加求和、高通滤波提取高频成分、n列的像素叠加求和后取平均值等数据处理得到数据将所述m个子集 各行依次减去对应的平均值 即可得到完全校正弧形伪影后第三中间图像R3(n,θ);

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

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

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

9.根据权利要求6所述的弧形伪影校正方法,其特征在于,所述将所有的CT图像变换到以探测器像元通道n为横坐标、投影角度θ为纵坐标的直角坐标系中,获得n×m个分辨率为n×1的第一中间图像R1(n,θ)的步骤包括:将二维面阵探测器探测的数据转换成一系列的一维线阵探测器探测的数据;

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

10.根据权利要求6所述的CT图像的弧形伪影校正方法,其特征在于,所述对第一中间图像R1(n,θ)依次进行n列的像素各自叠加求和、n列的像素叠加求和后取平均值等数据处理后得到数据E1(n)的步骤包括:对第一中间图像R1(n,θ)的n列的像素各自进行叠加求和,得到数据S(n);

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

对所述数据S(n)n列叠加求和后取平均值,得到数据E1(n)。

11.根据权利要求6所述的CT图像的弧形伪影校正方法,其特征在于,所述将第二中间图像R2(n,θ)分为m个子集 对每个子集 依次进行所有行的像素各自叠加求和、高通滤波提取高频成分、所有行的像素叠加求和后取平均值等数据处理得到数据的步骤包括:将第二中间图像R2(n,θ)按总的投影角度θ等分为m个子集对 每个 内 进 行所 有 行 的 像 素 各自 叠 加 求 和,得到 m个 数 据对所述数据Ai(n)进行高通滤波提取高频成分,得到数据对所述数据 所有行的像素叠加求和后取平均值,得到数据