CN103020928A - 锥束ct***的金属伪影校正方法 - Google Patents

锥束ct***的金属伪影校正方法 Download PDF

Info

Publication number
CN103020928A
CN103020928A CN2012104754831A CN201210475483A CN103020928A CN 103020928 A CN103020928 A CN 103020928A CN 2012104754831 A CN2012104754831 A CN 2012104754831A CN 201210475483 A CN201210475483 A CN 201210475483A CN 103020928 A CN103020928 A CN 103020928A
Authority
CN
China
Prior art keywords
image
metal
metal part
cone
res
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2012104754831A
Other languages
English (en)
Other versions
CN103020928B (zh
Inventor
胡战利
郑海荣
夏丹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Institute of Advanced Technology of CAS
Original Assignee
Shenzhen Institute of Advanced Technology of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenzhen Institute of Advanced Technology of CAS filed Critical Shenzhen Institute of Advanced Technology of CAS
Priority to CN201210475483.1A priority Critical patent/CN103020928B/zh
Publication of CN103020928A publication Critical patent/CN103020928A/zh
Application granted granted Critical
Publication of CN103020928B publication Critical patent/CN103020928B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明涉及一种锥束CT***的金属伪影校正方法，包括下述步骤：从原始正投影图像f(x，y)中分离出金属投影图像M(x，y)；对金属投影图像M(x，y)进行重建，得到金属部分的CT重建图像XMetal；将原始正投影图像f(x，y)减去金属投影图像M(x，y)，即得到不含金属部分的投影图像fres(x，y)；对不含金属部分的投影图像fres(x，y)进行重建，得到不含金属部分的CT图像Xres；及将金属部分CT图像XMetal与不含金属部分的CT图像Xres相加，即得到金属伪影校正后的最终CT重建图像Xcorrection。上述校正方法由于不进行坐标变换，所以不会损失图像空间分辨率，提高了图像质量。

Description

锥束CT***的金属伪影校正方法
技术领域
本发明涉及医学图像处理领域，尤其涉及一种锥束CT***的金属伪影校正方法。
背景技术
锥束CT(Computed Tomography)以扫描速度快、射线利用率高、剂量低、空间分辨率各向同性等优势在口腔疾病诊断、人体局部器官成像、工业检测等领域发展迅速。但是，CT成像时高衰减系数物质的存在(如人体内的金属整形器件、病人必须携带的活检剂、芯片中金属等)会导致严重的金属伪影，表现为条纹状亮纹、暗带或阴影，严重破坏了CT图像的真实性，极大的限制了其在医学与工业领域的高精度检测。因此去除或最大程度的减少金属伪影是十分必要的。
中国专利申请号201010616842.1提出的“一种去除CT图像中金属伪影的图像后处理方法”。在该专利中，将CT的原始图像由直角坐标图像转换为极坐标图像，在极坐标图像中确定金属投影区域，在极坐标图像中建立模型，采用上述的模型进行模型修正，修正上述模型修正中引入的正反投影误差，将极坐标图像转换成直角坐标图像。由于要进行坐标变换，所以会大大损失图像空间分辨率，导致图像质量下降。
中国专利申请号201110123647.X提出的“用于减少CT图像数据中的图像伪影、尤其是金属伪影的方法”。在该专利中，在不同的平均X射线能量时产生两个CT图像数据组。通过将这两个CT图像数据组加权地组合来计算出新的CT图像数据组。在此在加权的组合中所采用的加权因子被选择为，使得在该新的CT图像数据组中的图像伪影与在两个原始CT图像数据组中的图像伪影相比明显减少。由于要在不同能量时产生两个CT图像数据组，需要重复扫描两次，增加了病人所承受的辐射剂量，而且由于病人扫描时间加长，可能会带来潜在的运动伪影。
发明内容
基于此，有必要针对上述金属伪影校正方法存在的缺陷，提供一种图像质量高的锥束CT***的金属伪影校正方法。
一种锥束CT***的金属伪影校正方法，包括下述步骤：
从原始正投影图像f(x，y)中分离出金属投影图像M(x，y)；
对所述金属投影图像M(x，y)进行重建，得到金属部分的CT图像XMetal；
将所述原始正投影图像f(x，y)减去所述金属投影图像M(x，y)，即得到不含金属部分的投影图像fres(x，y)；
对所述不含金属部分的投影图像fres(x，y)进行重建，得到不含金属部分的CT图像Xres；及
将所述金属部分CT图像XMetal与所述不含金属部分的CT图像Xres相加，即得到金属伪影校正后的最终CT重建图像Xcorrection。
在本发明提供的实施例中，其中，从原始正投影图像f(x，y)中分离出金属投影图像M(x，y)，包括下述步骤：
步骤S11：定义初始阀值I0，所述初始阀值为I0＝(Imax+Imin)/2，其中，Imax为所述原始正投影图像f(x，y)中像素最大灰度值，Imin为所述原始正投影图像f(x，y)中最小灰度值；
步骤S12：基于所述初始阀值I0，将所述原始正投影图像f(x，y)分为两组图像T1(x，y)、T2(x，y)，其中，T1(x，y)的所有像素灰度值大于I0，T2(x，y)的所有像素灰度值均小于或等于I0；
步骤S13：计算所述图像T1(x，y)内所有像素的均值及图像T2(x，y)内所有像素的均值，并分别记为E1和E2；
步骤S14：定义新阀值I，所述新阀值I＝(E1+E2)/2；
步骤S15：判断|I-I0|是否大于ε，若“是”则进行下一步，若“否”则将所述I赋予I0，并返回步骤S12，其中，所述ε为0.5；
步骤S16：将所述新阀值I赋予所述初始阀值I0；及
步骤S17：所述金属投影图像M(x，y)为原始正投影图像f(x，y)中的像素灰度值大于I0的像素。
在本发明提供的实施例中，其中，对所述金属投影图像M(x，y)进行重建，得到金属部分的CT重建图像XMetal，是采用下述构造公式进行重建：
Figure BSA00000809111400031
其中，AXMetal＝M(x，y)，XMetal≥0
上述式中，A表示***矩阵，XMetal表示金属部分CT图像，M(x，y)表示分离出的金属部分投影图像，
Figure BSA00000809111400032
表示XMetal的l1范数，
Figure BSA00000809111400033
表示需满足XMetal取最小值。
在本发明提供的实施例中，其中，对所述不含金属部分的投影图像fres(x，y)进行重建，得到不含金属部分的CT图像Xres，是采用FDK重建方法进行重建。
在本发明提供的实施例中，其中，对所述不含金属部分的投影图像fres(x，y)进行重建，得到不含金属部分的CT图像Xres，是采用BPF重建方法进行重建。
在本发明提供的实施例中，其中，对所述不含金属部分的投影图像fres(x，y)进行重建，得到不含金属部分的CT图像Xres，是采用FBP重建方法进行重建。
本发明提供的锥束CT***的金属伪影校正方法，首先从原始正投影图像f(x，y)中分离出金属投影图像M(x，y)；对金属投影图像M(x，y)进行重建，得到金属部分的CT图像XMetal；并将原始正投影图像f(x，y)减去所述金属投影图像M(x，y)，即得到不含金属部分的投影图像fres(x，y)；对不含金属部分的投影图像fres(x，y)进行重建，得到不含金属部分的CT图像Xres；最后将金属部分CT图像XMetal与不含金属部分的CT图像Xres相加，即得到金属伪影校正后的最终CT重建图像Xcorrection。上述校正方法由于不进行坐标变换，所以不会损失图像空间分辨率，提高了图像质量。
另外，由于上述校正方法提供的从原始正投影图像f(x，y)中分离出金属投影图像M(x，y)的方法，采用自适应阀值分割方法，所以能够适用于金属结构复杂的情况。
同时，上述校正方法提供的对金属投影图像M(x，y)的重建方法，运算量小，速度快，易于实用化。
附图说明
图1为本发明实施例提供的锥束CT***的金属伪影校正方法的步骤流程图。
图2为本发明实施例提供的从原始正投影图像f(x，y)中分离出金属投影图像M(x，y)的步骤流程图。
具体实施方式
请参阅图1，为本发明实施例提供的锥束CT***的金属伪影校正方法的步骤流程图，包括下述步骤：
步骤S10：从原始正投影图像f(x，y)中分离出金属投影图像M(x，y)。
请参阅图2，为本发明实施例提供的从原始正投影图像f(x，y)中分离出金属投影图像M(x，y)的步骤流程图，步骤S10具体为：
步骤S11：定义初始阀值I0。其中，初始阀值为I0＝(Imax+Imin)/2，上述式中，Imax为原始正投影图像f(x，y)中像素最大灰度值，Imin为原始正投影图像f(x，y)中最小灰度值。
步骤S12：基于初始阀值I0，将原始正投影图像f(x，y)分为两组图像T1(x，y)及T2(x，y)。其中，图像T1(x，y)的所有像素灰度值大于I0，图像T2(x，y)的所有像素灰度值均小于或等于I0。
步骤S13：计算图像T1(x，y)内所有像素的均值及图像T2(x，y)内所有像素的均值，并分别记为E1和E2。即E1＝(1/n)∑T1，E2＝(1/n)∑T2，上述式中，T1、T2分别为图像T1(x，y)及图像T2(x，y)的像素，n为像素的数量。
步骤S14：定义新阀值I。其中，新阀值I＝(E1+E2)/2。
步骤S15：判断|I-I0|是否大于ε，若“是”则进行下一步，若“否”则将新阀值I赋予初始阀值I0，并返回步骤S12。其中，ε为预先设定的参数，可根据不同金属物体的属性(如人体内的金属整形器件、病人必须携带的活检剂、芯片中金属等)，依照经验选取。在本发明提供的实施例中，ε优选为0.5。可以理解，当|I-I0|大于ε时，执行下一步；当|I-I0|小于或等于ε时，返回步骤S12，并将此时I值赋予I0值，执行后续步骤。
步骤S16：将新阀值I赋予初始阀值I0。可以理解，基于步骤S15，当|I-I0|大于ε时，将新阀值I赋予初始阀值I0。
步骤S17：金属投影图像M(x，y)为原始正投影图像f(x，y)中的像素灰度值大于I0的像素。基于步骤S16中的I0，金属投影图像M(x，y)为原始正投影图像f(x，y)中的像素灰度值大于I0的像素。
上述实施例提供的从原始正投影图像f(x，y)中分离出金属投影图像M(x，y)的方法，由于采用自适应阈值分割方法，尤其适用于金属结构复杂的情况。
步骤S20：对金属投影图像M(x，y)进行重建，得到金属部分的CT图像XMetal。在本发明提供的实施例中，对于分离出的金属部分投影图像M(x，y)，采用下述构造公式进行重建，得到金属投影图像M(x，y)的CT重建图像XMetal
Figure BSA00000809111400051
其中，AXMctal＝M(x，y)，XMetal≥0
上述式中，A表示***矩阵，XMetal表示金属部分CT图像，M(x，y)表示分离出的金属部分投影图像，
Figure BSA00000809111400061
表示XMetal的l1范数，表示需满足XMetal取最小值。
上述对金属投影图像M(x，y)重建的方法，采用稀疏数据重建方法，运算量小，速度快，易于实用化。
步骤S30：将原始正投影图像f(x，y)减去步骤S10得到的金属投影图像M(x，y)，即得到不含金属部分的投影图像fres(x，y)，即不含金属部分的投影图像fres(x，y)＝f(x，y)-M(x，y)。
步骤S40：对不含金属部分的投影图像fres(x，y)进行重建，得到不含金属部分的CT图像Xres。
在本发明提供的实施例中，对不含金属部分的投影图像fres(x，y)进行重建，优选采用FDK重建方法进行重建。具体为，对投影图像fres(x，y)序列集进行反log操作，即通过下述公式：
-log{(物体图像-平均暗场图像)/(平均亮场图像-平均暗场图像)}
得到处理后的投影图像。对投影图像中的数据进行加权，用以修正像素到射线源的距离及角度引起的误差，然后对不同投影角度的投影数据进行水平方向上的一维滤波，最后进行三维反投影得到重建的图像。其中，重建的体素值是通过该体素的所有投影角度的射线的贡献之和。
在本发明提供的实施例中，对不含金属部分的投影图像fres(x，y)进行重建，还可以采用BPF重建方法进行重建。具体为，对投影图像fres(x，y)序列集进行反log操作，即通过下述公式：
-log{(物体图像-平均暗场图像)/(平均亮场图像-平均暗场图像)}
得到处理后的投影图像。在投影图像中，对每个角度下的锥束投影数据进行求导，并由求导后的投影数据对螺旋轨道内每条PI线上的投影图像进行加权反投影。对于加权反投影得到的PI线上的数据，沿着该条PI线求Hilbert变换(希尔伯特变换)，得该条PI线上的重建结果。最后对螺旋轨道内所有PI线上的重建结果进行重采样，得到重建的图像。
在本发明提供的实施例中，对不含金属部分的投影图像fres(x，y)进行重建，还可以采用FBP重建方法进行重建。具体为，对投影图像fres(x，y)序列集进行反log操作，即通过下述公式：
-log{(物体图像-平均暗场图像)/(平均亮场图像-平均暗场图像)}
得到处理后的投影图像。再对投影图像进行加权、卷积滤波、加权反投影后获得重建的图像。
步骤S40：将金属部分CT图像XMetal与不含金属部分的CT图像Xres相加，即得到金属伪影校正后的最终CT重建图像Xcorrection。
上述实施例提供的锥束CT***的金属伪影校正方法，首先从原始正投影图像f(x，y)中分离出金属投影图像M(x，y)；对金属投影图像M(x，y)进行重建，得到金属部分的CT图像XMetal；并将原始正投影图像f(x，y)减去所述金属投影图像M(x，y)，即得到不含金属部分的投影图像fres(x，y)；对不含金属部分的投影图像fres(x，y)进行重建，得到不含金属部分的CT图像Xres；最后将金属部分CT图像XMetal与不含金属部分的CT图像Xres相加，即得到金属伪影校正后的最终CT重建图像Xcorrection。上述校正方法由于不进行坐标变换，所以不会损失图像空间分辨率，提高了图像质量。
可以理解，本发明上述实施例提供的金属伪影校正方法，除应用于锥束CT***外，做适当变后，也可应用于PET、SPECT***。
以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，虽然本发明已以较佳实施例揭露如上，然而并非用以限定本发明，任何熟悉本专业的技术人员，在不脱离本发明技术方案范围内，当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例，但凡是未脱离本发明技术方案内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属于本发明技术方案的范围内。

Claims (6)

1.一种锥束CT***的金属伪影校正方法，其特征在于，包括下述步骤：
从原始正投影图像f(x，y)中分离出金属投影图像M(x，y)；
对所述金属投影图像M(x，y)进行重建，得到金属部分的CT图像XMetal；
将所述原始正投影图像f(x，y)减去所述金属投影图像M(x，y)，即得到不含金属部分的投影图像fres(x，y)；
对所述不含金属部分的投影图像fres(x，y)进行重建，得到不含金属部分的CT图像Xres；及
将所述金属部分CT图像XMetal与所述不含金属部分的CT图像Xres相加，即得到金属伪影校正后的最终CT重建图像Xcorrection。
2.根据权利要求1所述的锥束CT***的金属伪影校正方法，其特征在于，其中，从原始正投影图像f(x，y)中分离出金属投影图像M(x，y)，包括下述步骤：
步骤S11：定义初始阀值I0，所述初始阀值为I0＝(Imax+Imin)/2，其中，Imax为所述原始正投影图像f(x，y)中像素最大灰度值，Imin为所述原始正投影图像f(x，y)中最小灰度值；
步骤S12：基于所述初始阀值I0，将所述原始正投影图像f(x，y)分为两组图像T1(x，y)、T2(x，y)，其中，T1(x，y)的所有像素灰度值大于I0，T2(x，y)的所有像素灰度值均小于或等于I0；
步骤S13：计算所述图像T1(x，y)内所有像素的均值及图像T2(x，y)内所有像素的均值，并分别记为E1和E2；
步骤S14：定义新阀值I，所述新阀值I＝(E1+E2)/2；
步骤S15：判断|I-I0|是否大于ε，若“是”则进行下一步，若“否”则将所述I值赋予I0值，并返回步骤S12，其中，所述ε为0.5；
步骤S16：将所述新阀值I赋予所述初始阀值I0
步骤S17：所述金属投影图像M(x，y)为原始正投影图像f(x，y)中的像素灰度值大于I0的像素。
3.根据权利要求1所述的锥束CT***的金属伪影校正方法，其特征在于，其中，对所述金属投影图像M(x，y)进行重建，得到金属部分的CT重建图像XMetal，采用下述构造公式进行重建：
Figure FSA00000809111300021
其中，AXMetal＝M(x，y)，XMetal≥0
上述式中，A表示***矩阵，XMetal表示金属部分CT图像，M(x，y)表示分离出的金属部分投影图像，
Figure FSA00000809111300022
表示XMetal的l1范数，
Figure FSA00000809111300023
表示需满足XMetal取最小值。
4.根据权利要求1所述的锥束CT***的金属伪影校正方法，其特征在于，其中，对所述不含金属部分的投影图像fres(x，y)进行重建，得到不含金属部分的CT图像Xres，是采用FDK重建方法进行重建。
5.根据权利要求1所述的锥束CT***的金属伪影校正方法，其特征在于，其中，对所述不含金属部分的投影图像fres(x，y)进行重建，得到不含金属部分的CT图像Xres，是采用BPF重建方法进行重建。
6.根据权利要求1所述的锥束CT***的金属伪影校正方法，其特征在于，其中，对所述不含金属部分的投影图像fres(x，y)进行重建，得到不含金属部分的CT图像Xres，是采用FBP重建方法进行重建。
CN201210475483.1A 2012-11-21 2012-11-21 锥束ct***的金属伪影校正方法 Active CN103020928B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210475483.1A CN103020928B (zh) 2012-11-21 2012-11-21 锥束ct***的金属伪影校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210475483.1A CN103020928B (zh) 2012-11-21 2012-11-21 锥束ct***的金属伪影校正方法

Publications (2)

Publication Number Publication Date
CN103020928A true CN103020928A (zh) 2013-04-03
CN103020928B CN103020928B (zh) 2016-06-29

Family

ID=47969495

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210475483.1A Active CN103020928B (zh) 2012-11-21 2012-11-21 锥束ct***的金属伪影校正方法

Country Status (1)

Country Link
CN (1) CN103020928B (zh)

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103700124A (zh) * 2013-12-27 2014-04-02 公安部第一研究所 一种x射线双能ct成像抑制金属效应方法
CN104778667A (zh) * 2015-04-14 2015-07-15 南京邮电大学 一种基于水平集的锥束ct中杯状伪影的校正方法
CN105046675A (zh) * 2015-08-24 2015-11-11 重庆邮电大学 基于mvct图像消除kvct图像中金属伪影的方法
CN105118030A (zh) * 2015-08-11 2015-12-02 上海联影医疗科技有限公司 医学图像金属伪影的校正方法及装置
CN105528766A (zh) * 2014-12-04 2016-04-27 沈阳东软医疗***有限公司 Ct金属伪影处理方法和装置
CN105528771A (zh) * 2016-01-19 2016-04-27 南京邮电大学 一种使用能量函数方法的锥束ct中杯状伪影的校正方法
CN105631908A (zh) * 2015-12-16 2016-06-01 沈阳东软医疗***有限公司 一种pet图像重建方法和装置
CN105701778A (zh) * 2016-01-11 2016-06-22 赛诺威盛科技（北京）有限公司 一种从ct图像中去除金属伪影的方法
CN105719245A (zh) * 2016-01-12 2016-06-29 南方医科大学 一种利用投影数据去除ct探测元故障造成的环状伪影的方法
CN105787973A (zh) * 2014-12-19 2016-07-20 合肥美亚光电技术股份有限公司 Ct***中投影图像的重建方法及装置
WO2017063569A1 (en) * 2015-10-14 2017-04-20 Shanghai United Imaging Healthcare Co., Ltd. System and method for image correction
CN106909927A (zh) * 2015-12-23 2017-06-30 通用电气公司 Ct投影数据中可信数据的提取装置及方法、ct成像***
CN106960429A (zh) * 2017-02-16 2017-07-18 中国科学院苏州生物医学工程技术研究所 一种ct图像金属伪影校正方法及装置
CN107016655A (zh) * 2017-03-30 2017-08-04 中国人民解放军信息工程大学 锥束cl几何全参数迭代校正方法
CN107714069A (zh) * 2016-08-10 2018-02-23 南京普爱医疗设备股份有限公司 一种基于cbct影像的金属伪影的测试方法
WO2018126434A1 (zh) * 2017-01-06 2018-07-12 深圳先进技术研究院 Ct图像阴影校正方法、装置及电子设备
CN108596993A (zh) * 2018-02-26 2018-09-28 上海奕瑞光电子科技股份有限公司 校正图像不饱和伪影的***及校正方法
CN109920020A (zh) * 2019-02-27 2019-06-21 西北工业大学 一种锥束ct病态投影重建伪影抑制方法
CN111542853A (zh) * 2017-10-31 2020-08-14 皇家飞利浦有限公司 数据采集期间的运动伪影预测
CN113729747A (zh) * 2021-08-09 2021-12-03 中科超精(南京)科技有限公司 一种球形金属标记的锥束ct金属伪影去除***及去除方法
CN117830456A (zh) * 2024-03-04 2024-04-05 中国科学技术大学 用于校正图像金属伪影的方法、装置及电子设备

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20010028696A1 (en) * 2000-04-07 2001-10-11 Yosihiro Yamada Image processing method of X-ray CT, X-ray CT and X-ray CT image-taking recording medium
CN101777177A (zh) * 2009-12-29 2010-07-14 上海维宏电子科技有限公司 基于衰减滤波的ct图像去金属伪影混合重建法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20010028696A1 (en) * 2000-04-07 2001-10-11 Yosihiro Yamada Image processing method of X-ray CT, X-ray CT and X-ray CT image-taking recording medium
CN101777177A (zh) * 2009-12-29 2010-07-14 上海维宏电子科技有限公司 基于衰减滤波的ct图像去金属伪影混合重建法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈豫: "CT图像金属伪影校正算法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》, no. 07, 15 July 2010 (2010-07-15), pages 27 - 39 *

Cited By (41)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103700124A (zh) * 2013-12-27 2014-04-02 公安部第一研究所 一种x射线双能ct成像抑制金属效应方法
CN103700124B (zh) * 2013-12-27 2017-11-14 公安部第一研究所 一种x射线双能ct成像抑制金属效应方法
CN105528766B (zh) * 2014-12-04 2018-08-03 沈阳东软医疗***有限公司 Ct金属伪影处理方法和装置
CN105528766A (zh) * 2014-12-04 2016-04-27 沈阳东软医疗***有限公司 Ct金属伪影处理方法和装置
CN105787973A (zh) * 2014-12-19 2016-07-20 合肥美亚光电技术股份有限公司 Ct***中投影图像的重建方法及装置
CN104778667A (zh) * 2015-04-14 2015-07-15 南京邮电大学 一种基于水平集的锥束ct中杯状伪影的校正方法
CN104778667B (zh) * 2015-04-14 2019-05-03 南京邮电大学 一种基于水平集的锥束ct中杯状伪影的校正方法
CN105118030B (zh) * 2015-08-11 2018-08-03 上海联影医疗科技有限公司 医学图像金属伪影的校正方法及装置
CN105118030A (zh) * 2015-08-11 2015-12-02 上海联影医疗科技有限公司 医学图像金属伪影的校正方法及装置
US10964072B2 (en) 2015-08-24 2021-03-30 Chongqing University Of Posts And Telecommunications Methods, systems, and media for noise reduction in computed tomography images
CN105046675B (zh) * 2015-08-24 2017-09-29 重庆邮电大学 基于mvct图像消除kvct图像中金属伪影的方法
CN105046675A (zh) * 2015-08-24 2015-11-11 重庆邮电大学 基于mvct图像消除kvct图像中金属伪影的方法
US10521886B2 (en) 2015-10-14 2019-12-31 Shanghai United Imaging Healthcare Co., Ltd. System and method for image correction
WO2017063569A1 (en) * 2015-10-14 2017-04-20 Shanghai United Imaging Healthcare Co., Ltd. System and method for image correction
US11756164B2 (en) 2015-10-14 2023-09-12 Shanghai United Imaging Healthcare Co., Ltd. System and method for image correction
US11232543B2 (en) 2015-10-14 2022-01-25 Shanghai United Imaging Healthcare Co., Ltd. System and method for image correction
GB2548303A (en) * 2015-10-14 2017-09-13 Shanghai United Imaging Healthcare Co Ltd System and method for image correction
GB2548303B (en) * 2015-10-14 2018-02-21 Shanghai United Imaging Healthcare Co Ltd System and method for image correction
CN105631908A (zh) * 2015-12-16 2016-06-01 沈阳东软医疗***有限公司 一种pet图像重建方法和装置
CN105631908B (zh) * 2015-12-16 2018-11-20 沈阳东软医疗***有限公司 一种pet图像重建方法和装置
CN106909927A (zh) * 2015-12-23 2017-06-30 通用电气公司 Ct投影数据中可信数据的提取装置及方法、ct成像***
CN105701778A (zh) * 2016-01-11 2016-06-22 赛诺威盛科技（北京）有限公司 一种从ct图像中去除金属伪影的方法
CN105719245B (zh) * 2016-01-12 2018-06-19 南方医科大学 一种利用投影数据去除ct探测元故障造成的环状伪影的方法
CN105719245A (zh) * 2016-01-12 2016-06-29 南方医科大学 一种利用投影数据去除ct探测元故障造成的环状伪影的方法
CN105528771B (zh) * 2016-01-19 2018-09-18 南京邮电大学 一种使用能量函数方法的锥束ct中杯状伪影的校正方法
CN105528771A (zh) * 2016-01-19 2016-04-27 南京邮电大学 一种使用能量函数方法的锥束ct中杯状伪影的校正方法
CN107714069A (zh) * 2016-08-10 2018-02-23 南京普爱医疗设备股份有限公司 一种基于cbct影像的金属伪影的测试方法
WO2018126434A1 (zh) * 2017-01-06 2018-07-12 深圳先进技术研究院 Ct图像阴影校正方法、装置及电子设备
CN106960429A (zh) * 2017-02-16 2017-07-18 中国科学院苏州生物医学工程技术研究所 一种ct图像金属伪影校正方法及装置
CN106960429B (zh) * 2017-02-16 2019-08-27 中国科学院苏州生物医学工程技术研究所 一种ct图像金属伪影校正方法及装置
CN107016655A (zh) * 2017-03-30 2017-08-04 中国人民解放军信息工程大学 锥束cl几何全参数迭代校正方法
CN111542853B (zh) * 2017-10-31 2024-05-14 皇家飞利浦有限公司 数据采集期间的运动伪影预测
CN111542853A (zh) * 2017-10-31 2020-08-14 皇家飞利浦有限公司 数据采集期间的运动伪影预测
CN108596993B (zh) * 2018-02-26 2022-07-12 上海奕瑞光电子科技股份有限公司 校正图像不饱和伪影的***及校正方法
CN108596993A (zh) * 2018-02-26 2018-09-28 上海奕瑞光电子科技股份有限公司 校正图像不饱和伪影的***及校正方法
CN109920020B (zh) * 2019-02-27 2022-10-18 西北工业大学 一种锥束ct病态投影重建伪影抑制方法
CN109920020A (zh) * 2019-02-27 2019-06-21 西北工业大学 一种锥束ct病态投影重建伪影抑制方法
CN113729747A (zh) * 2021-08-09 2021-12-03 中科超精(南京)科技有限公司 一种球形金属标记的锥束ct金属伪影去除***及去除方法
CN113729747B (zh) * 2021-08-09 2024-04-26 中科超精(南京)科技有限公司 一种球形金属标记的锥束ct金属伪影去除***及去除方法
CN117830456A (zh) * 2024-03-04 2024-04-05 中国科学技术大学 用于校正图像金属伪影的方法、装置及电子设备
CN117830456B (zh) * 2024-03-04 2024-05-28 中国科学技术大学 用于校正图像金属伪影的方法、装置及电子设备

Also Published As

Publication number Publication date
CN103020928B (zh) 2016-06-29

Similar Documents

Publication Publication Date Title
CN103020928B (zh) 锥束ct***的金属伪影校正方法
WO2021159948A1 (zh) 一种基于深度学习的低剂量pet三维重建方法
CN109146994B (zh) 一种面向多能谱x射线ct成像的金属伪影校正方法
JP2020036877A (ja) 反復的画像再構成フレームワーク
CN103065279B (zh) 圆轨道锥束计算机断层摄影（ct）中用于大幅度地减轻伪影的方法以及***
US10867375B2 (en) Forecasting images for image processing
CN107481297A (zh) 一种基于卷积神经网络的ct图像重建方法
CN105118039B (zh) 实现锥束ct图像重建的方法及***
CN101917906A (zh) 通过将物体的围绕物用作动态约束在断层扫描中减少剂量并增强图像
CN101980302A (zh) 投影数据恢复导引的非局部平均低剂量ct重建方法
CN104574416A (zh) 一种低剂量能谱ct图像去噪方法
CN102024251A (zh) 用于双源ct的基于多图像虚拟非对比图像增强的***和方法
US20120275673A1 (en) Truncation correction imaging enhancement method and system
EP3847623A1 (en) A method of generating an enhanced tomographic image of an object
CN113870138A (zh) 基于三维U-net的低剂量CT图像去噪方法及***
CN103793890A (zh) 一种能谱ct图像的恢复处理方法
Jin et al. Bone-induced streak artifact suppression in sparse-view CT image reconstruction
Michielsen et al. Patchwork reconstruction with resolution modeling for digital breast tomosynthesis
Kim et al. Convolutional neural network–based metal and streak artifacts reduction in dental CT images with sparse‐view sampling scheme
CN107316282A (zh) 一种图像降噪的***和方法
CN114511497A (zh) 应用于锥束ct稀疏采样的成像方法及装置
Zhu et al. Iterative CT reconstruction via minimizing adaptively reweighted total variation
CN110176045A (zh) 一种由单能ct图像生成双能ct图像的方法
US20220414832A1 (en) X-ray imaging restoration using deep learning algorithms
EP4330862A1 (en) System and method for medical imaging

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant