CN113017662B - 一种ct图像的混叠伪影去除方法及系统、ct扫描仪 - Google Patents
一种ct图像的混叠伪影去除方法及系统、ct扫描仪 Download PDFInfo
- Publication number
- CN113017662B CN113017662B CN202110115941.XA CN202110115941A CN113017662B CN 113017662 B CN113017662 B CN 113017662B CN 202110115941 A CN202110115941 A CN 202110115941A CN 113017662 B CN113017662 B CN 113017662B
- Authority
- CN
- China
- Prior art keywords
- ray
- rays
- weight
- current
- detector
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 239000000523 sample Substances 0.000 claims description 6
- 230000000694 effects Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Veterinary Medicine (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及一种CT图像的混叠伪影去除方法,包括:S1、计算每条射线进行反投影的权重;S2、采集满足条件的所有目标射线,并将所有目标射线与当前射线按照在XY平面内的投影方向分为两类,与当前射线同向的射线为正向射线,与当前射线反向的射线为共轭射线;S3、分别计算正向射线的权重和、共轭射线的权重和以及两者之和得到的总权重,根据正向射线的权重和、总权重得到正向射线的权重比,根据当前射线的权重、总权重得到当前射线的权重比;S4、根据正向射线的权重比以及当前射线的权重比,重新计算当前射线对体素进行反投影的反投影值;S5、根据重新计算得到的每条射线的反投影值对CT图像进行重建。本发明的混叠伪影去除效果佳。
Description
技术领域
本发明属于医学影像技术领域,具体涉及一种CT图像的混叠伪影去除方法及系统、CT扫描仪。
背景技术
混叠伪影是CT图像重建中常见的一种伪影,该伪影产生的原因一般为帧采样不足或者探测器采样不足,其在图像上的表现为图像边缘或者高密度物质的边缘出现放射状条纹。
为了解决帧采样不足导致的混叠伪影,常用方法为提高投影方向的采样率。
为了解决探测器采样不足导致的混叠伪影,在系统硬件上常用的方法有以下两种:
一种是使用1/4探测器偏移,通过使用对向角度采集的数据来提高分辨率,减少混叠伪影,但是使用该技术时只能保证探测器视野中心的位置有较好的效果;
另一种是飞焦点技术,在采集投影时焦点在两个固定位置上来回摆动,该技术也可以有效地提高分辨率,较少混叠,但是对球管的要求较高。
除了上述在系统硬件上使用1/4探测器偏移、飞焦点等技术外,还可以通过各种图像处理技术来消除混叠伪影。常见的图像处理方法包括以下两类:
一类是在投影域对数据进行处理,一般是通过使用平滑的卷积核来实现的,使用平滑的卷积核可以明显地压制投影数据中的高频成分,减少混叠伪影,但该方法容易导致图像细节丢失,变得更模糊;
另一类是在图像域进行处理,一般会对图像进行傅里叶变换,结合空间域和频域的相关信息对图像中混叠伪影的特征进行提取,但需要能够准确地识别混叠的特征。
发明内容
基于现有技术中存在的上述缺点和不足,本发明的目的之一是至少解决现有技术中存在的上述问题之一或多个,换言之,本发明的目的之一是提供满足前述需求之一或多个的一种CT图像的混叠伪影去除方法及系统、CT扫描仪。
为了达到上述发明目的,本发明采用以下技术方案:
一种CT图像的混叠伪影去除方法,包括以下步骤:
S1、计算每条射线进行反投影的权重;
S2、采集满足以下条件的所有目标射线:射线在XY平面内的投影与当前射线的投影平行、射线穿过体素、射线能被探测器接收,并将所有目标射线与当前射线按照在XY平面内的投影方向分为两类,与当前射线同向的射线为正向射线,与当前射线反向的射线为共轭射线;
S3、分别计算正向射线的权重和、共轭射线的权重和以及两者之和得到的总权重,根据正向射线的权重和、总权重得到正向射线的权重比,根据当前射线的权重、总权重得到当前射线的权重比;
S4、根据正向射线的权重比以及当前射线的权重比,重新计算当前射线对体素进行反投影的反投影值;
S5、根据重新计算得到的每条射线的反投影值对CT图像进行重建。
作为优选方案,所述步骤S1中,计算每条射线进行反投影的权重为:
作为优选方案,所述步骤S3,包括:
所述正向射线的权重和、共轭射线的权重和分别记为wdirc、wconj;
所述总权重为:wdirc+wconj;
作为优选方案,所述步骤S4,包括:
判断rdirc是否大于0.5;
若是,则当前射线对体素进行反投影的反投影值val为:
若否,则当前射线对体素进行反投影的反投影值val为:
val=rcurr×proj(t,v);
其中,proj(t,v)为当前射线入射至探测器位置对应的像素坐标(t,v)对应的投影值,t为沿着探测器通道对应的X方向,v为沿着探测器排对应的Z方向;proj(t-1,v),proj(t+1,v)为当前射线入射至探测器位置对应的相邻两像素坐标对应的投影值。
本发明还提供一种CT图像的混叠伪影去除系统,包括:
采集模块,用于采集满足以下条件的所有目标射线:射线在XY平面内的投影与当前射线的投影平行、射线穿过体素、射线能被探测器接收;
分类模块,用于将所有目标射线与当前射线按照在XY平面内的投影方向分为两类,与当前射线同向的射线为正向射线,与当前射线反向的射线为共轭射线;
计算模块,用于计算每条射线进行反投影的权重,还用于计算正向射线的权重和、共轭射线的权重和以及两者之和得到的总权重,还用于根据正向射线的权重和、总权重计算得到正向射线的权重比,还用于根据当前射线的权重、总权重计算得到当前射线的权重比,还用于根据正向射线的权重比以及当前射线的权重比重新计算当前射线对体素进行反投影的反投影值;
重建模块,用于根据重新计算得到的每条射线的反投影值对CT图像进行重建。
作为优选方案,计算每条射线进行反投影的权重为:
作为优选方案,所述正向射线的权重和、共轭射线的权重和分别记为wdirc、wconj;
所述总权重为:wdirc+wconj;
作为优选方案,混叠伪影去除系统,还包括:
判断模块,用于判断判断rdirc是否大于0.5;
若是,则计算模块计算当前射线对体素进行反投影的反投影值val为:
若否,则计算模块计算当前射线对体素进行反投影的反投影值val为:
val=rcurr×proj(t,v);
其中,proj(t,v)为当前射线入射至探测器位置对应的像素坐标(t,v)对应的投影值,t为沿着探测器通道对应的X方向,v为沿着探测器排对应的Z方向;proj(t-1,v),proj(t+1,v)为当前射线入射至探测器位置对应的相邻两像素坐标对应的投影值。
本发明还提供一种CT扫描仪,应用如上任一方案所述的混叠伪影去除方法。
本发明还提供一种CT扫描仪,搭载如上任一方案所述的混叠伪影去除系统。
本发明与现有技术相比,有益效果是:
本发明的CT图像的混叠伪影去除方法及系统,通过在反投影时对经过每个体素的所有射线按照特定方式计算权重,以此计算相应的反投影值进行图像重建,从而消除CT图像的混叠伪影。
附图说明
图1是本发明实施例的CT图像的混叠伪影去除方法的流程图;
图2是本发明实施例的光源穿过体素后入射至探测器上的位置及相应权重分布的示意图;
图3是本发明实施例的CT图像的混叠伪影去除方法去除前后的CT图像对比图;
图4是本发明实施例的CT图像的混叠伪影去除系统的构架图。
具体实施方式
以下通过具体实施例对本发明的技术方案作进一步解释说明。
本发明通过在反投影时对经过每个体素的所有射线按照特定方式计算权重的方法来消除混叠伪影。具体原理为:在将某一条射线反投影到一个体素时,先找出所有在XY平面内的投影与该射线的投影相平行的射线,之后按照射线在XY平面上的方向分为正向射线和反向(即共轭)射线,与当前射线方向相同的射线为正向射线,不同的为共轭射线;然后根据正向射线权重之和占总投影权重的比例,重新计算每条射线的权重,最后使用新的权重对每条射线进行反投影,根据相应的反投影值对CT图像进行重建。
具体地,如图1所示,本发明实施例的CT图像的混叠伪影去除方法,包括以下步骤:
S1、计算每条射线进行反投影的权重;
如图2所示,光源穿过体素后打在探测器上的位置及相应权重分布的示意图,其中,每条射线进行反投影的权重的计算公式为:
其中,Q为权重调节因子,取值为0~1;dz为探测器在Z方向的总长度;
S2、在对体素voxelj进行反投影时,采集满足以下条件的所有目标射线:射线在XY平面内的投影与当前射线的投影平行、射线穿过体素voxelj、射线能被探测器接收,并将所有目标射线与当前射线按照在XY平面内的投影方向分为两类,与当前射线同向的射线为正向射线,与当前射线反向的射线为共轭射线;
S3、分别计算正向射线的权重和、共轭射线的权重和以及两者之和得到的总权重,根据正向射线的权重和、总权重得到正向射线的权重比,根据当前射线的权重、总权重得到当前射线的权重比;
具体地,正向射线的权重和、共轭射线的权重和分别记为wdirc、wconj;
总权重为:wdirc+wconj;
正向射线的权重比rdirc为:
当前射线的权重比rcurr为:
其中,wcurr为当前射线的权重;
S4、根据正向射线的权重比以及当前射线的权重比,重新计算当前射线对体素进行反投影的反投影值;
具体地,根据正向射线的权重比rdirc使用不同的方式进行反投影:
判断正向射线的权重比rdirc是否大于0.5;
若是,则当前射线对体素进行反投影的反投影值val为:
若否,则当前射线对体素进行反投影的反投影值val为:
val=rcurr×proj(t,v);
其中,proj(t,v)为当前射线入射至探测器位置对应的像素坐标(t,v)对应的投影值,t为沿着探测器通道(channel)对应的X方向,v为沿着探测器排(slice)对应的Z方向;proj(t-1,v),proj(t+1,v)为当前射线入射至探测器位置对应的相邻两像素坐标对应的投影值。
S5、根据重新计算得到的每条射线的反投影值对CT图像进行重建。
如图3所示,通过本发明实施例的上述混叠伪影去除方法前后的图像对比可知,采用本发明实施例的混叠伪影去除方法之后,CT图像的混叠伪影明显减少。
如图4所示,本发明实施例还提供CT图像的混叠伪影去除系统,包括:
采集模块,用于采集满足以下条件的所有目标射线:射线在XY平面内的投影与当前射线的投影平行、射线穿过体素、射线能被探测器接收;
分类模块,用于将所有目标射线与当前射线按照在XY平面内的投影方向分为两类,与当前射线同向的射线为正向射线,与当前射线反向的射线为共轭射线;
计算模块,用于计算每条射线进行反投影的权重,还用于计算正向射线的权重和、共轭射线的权重和以及两者之和得到的总权重,还用于根据正向射线的权重和、总权重计算得到正向射线的权重比,还用于根据当前射线的权重、总权重计算得到当前射线的权重比,还用于根据正向射线的权重比以及当前射线的权重比重新计算当前射线对体素进行反投影的反投影值;
重建模块,用于根据重新计算得到的每条射线的反投影值对CT图像进行重建;
具体地,每条射线进行反投影的权重计算公式为:
其中,Q为权重调节因子,取值为0~1;dz为探测器在Z方向的总长度;
另外,正向射线的权重和、共轭射线的权重和分别记为wdirc、wconj;
总权重为:wdirc+wconj;
正向射线的权重比rdirc为:
当前射线的权重比rcurr为:
其中,wcurr为当前射线的权重;
另外,本发明实施例的混叠伪影去除系统还包括:
判断模块,用于判断正向射线的权重比rdirc是否大于0.5;
若是,则计算模块计算当前射线对体素进行反投影的反投影值val为:
若否,则计算模块计算当前射线对体素进行反投影的反投影值val为:
val=rcurr×proj(t,v);
其中,proj(t,v)为当前射线入射至探测器位置对应的像素坐标(t,v)对应的投影值,t为沿着探测器通道对应的X方向,v为沿着探测器排对应的Z方向;proj(t-1,v),proj(t+1,v)为当前射线入射至探测器位置对应的相邻两像素坐标对应的投影值。
最后,根据重新计算得到的每条射线的反投影值对CT图像进行重建,得到混叠伪影去除后的CT图像,如图3所示。
本发明实施例还提供一种CT扫描仪,应用上述混叠伪影去除方法以及搭载上述混叠伪影去除系统。
以上所述仅是对本发明的优选实施例及原理进行了详细说明,对本领域的普通技术人员而言,依据本发明提供的思想,在具体实施方式上会有改变之处,而这些改变也应视为本发明的保护范围。
Claims (4)
1.一种CT图像的混叠伪影去除方法,其特征在于,包括以下步骤:
S1、计算每条射线进行反投影的权重;
S2、采集满足以下条件的所有目标射线:射线在XY平面内的投影与当前射线的投影平行、射线穿过体素、射线能被探测器接收,并将所有目标射线与当前射线按照在XY平面内的投影方向分为两类,与当前射线同向的射线为正向射线,与当前射线反向的射线为共轭射线;
S3、分别计算正向射线的权重和、共轭射线的权重和以及两者之和得到的总权重,根据正向射线的权重和、总权重得到正向射线的权重比,根据当前射线的权重、总权重得到当前射线的权重比;
S4、根据正向射线的权重比以及当前射线的权重比,重新计算当前射线对体素进行反投影的反投影值;
S5、根据重新计算得到的每条射线的反投影值对CT图像进行重建;
所述步骤S1中,计算每条射线进行反投影的权重为:
所述步骤S3,包括:
所述正向射线的权重和、共轭射线的权重和分别记为wdirc、wconj;
所述总权重为:wdirc+wconj;
所述步骤S4,包括:
判断rdirc是否大于0.5;
若是,则当前射线对体素进行反投影的反投影值val为:
若否,则当前射线对体素进行反投影的反投影值val为:
val=rcurr×proj(t,v);
2.一种CT图像的混叠伪影去除系统,其特征在于,包括:
采集模块,用于采集满足以下条件的所有目标射线:射线在XY平面内的投影与当前射线的投影平行、射线穿过体素、射线能被探测器接收;
分类模块,用于将所有目标射线与当前射线按照在XY平面内的投影方向分为两类,与当前射线同向的射线为正向射线,与当前射线反向的射线为共轭射线;
计算模块,用于计算每条射线进行反投影的权重,还用于计算正向射线的权重和、共轭射线的权重和以及两者之和得到的总权重,还用于根据正向射线的权重和、总权重计算得到正向射线的权重比,还用于根据当前射线的权重、总权重计算得到当前射线的权重比,还用于根据正向射线的权重比以及当前射线的权重比重新计算当前射线对体素进行反投影的反投影值;
重建模块,用于根据重新计算得到的每条射线的反投影值对CT图像进行重建;
计算每条射线进行反投影的权重为:
所述正向射线的权重和、共轭射线的权重和分别记为wdirc、wconj;
所述总权重为:wdirc+wconj;
CT图像的混叠伪影去除系统,还包括:
判断模块,用于判断rdirc是否大于0.5;
若是,则计算模块计算当前射线对体素进行反投影的反投影值val为:
若否,则计算模块计算当前射线对体素进行反投影的反投影值val为:
val=rcurr×proj(t,v);
3.一种CT扫描仪,其特征在于,应用如权利要求1所述的混叠伪影去除方法。
4.一种CT扫描仪,其特征在于,搭载如权利要求2所述的混叠伪影去除系统。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110115941.XA CN113017662B (zh) | 2021-01-28 | 2021-01-28 | 一种ct图像的混叠伪影去除方法及系统、ct扫描仪 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110115941.XA CN113017662B (zh) | 2021-01-28 | 2021-01-28 | 一种ct图像的混叠伪影去除方法及系统、ct扫描仪 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113017662A CN113017662A (zh) | 2021-06-25 |
CN113017662B true CN113017662B (zh) | 2022-06-14 |
Family
ID=76459395
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110115941.XA Active CN113017662B (zh) | 2021-01-28 | 2021-01-28 | 一种ct图像的混叠伪影去除方法及系统、ct扫描仪 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113017662B (zh) |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4703424A (en) * | 1985-01-03 | 1987-10-27 | General Electric Company | Method and means of correcting for a shift in the center of rotation of a rotating fan beam CT system |
US5708691A (en) * | 1996-07-05 | 1998-01-13 | Kabushiki Kaisha Toshiba | X-ray computed tomographic imaging device and x-ray computed tomographic method |
US6490334B1 (en) * | 2001-06-29 | 2002-12-03 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for high pitch helical computed tomography image reconstruction |
US6873676B2 (en) * | 2003-03-05 | 2005-03-29 | Ge Medical Systems Global Technology Company, Llc | Convolution reconstruction algorithm for multi-slice CT |
JP2005087592A (ja) * | 2003-09-19 | 2005-04-07 | Hitachi Ltd | X線計測装置 |
US8824760B2 (en) * | 2009-10-20 | 2014-09-02 | University Of Utah Research Foundation | Modification and elimination of back projection weight during the CT image reconstruction |
JP5600271B2 (ja) * | 2010-06-30 | 2014-10-01 | 富士フイルム株式会社 | 放射線撮影装置および方法並びにプログラム |
AU2012207195B2 (en) * | 2011-01-21 | 2016-09-15 | Mercury Systems, Inc. | Optimized implementation of back projection for computed tomography (CT) |
CN102521853A (zh) * | 2011-12-08 | 2012-06-27 | 苏州生物医学工程技术研究所 | 一种反投影权重锥束ct重建算法 |
CN103472074B (zh) * | 2013-06-19 | 2016-01-20 | 清华大学 | Ct成像系统和方法 |
WO2016111016A1 (ja) * | 2015-01-09 | 2016-07-14 | 株式会社島津製作所 | 放射線画像生成方法および画像処理装置 |
CN109949411B (zh) * | 2019-03-22 | 2022-12-27 | 电子科技大学 | 一种基于三维加权滤波反投影和统计迭代的图像重建方法 |
CN111476860B (zh) * | 2020-04-22 | 2023-10-24 | 上海联影医疗科技股份有限公司 | 图像重建方法、装置、计算机设备和存储介质 |
-
2021
- 2021-01-28 CN CN202110115941.XA patent/CN113017662B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113017662A (zh) | 2021-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | Denoising 3D MR images by the enhanced non-local means filter for Rician noise | |
US8761540B2 (en) | Method and system for estimating noise level | |
US11727609B2 (en) | Limited-angle CT reconstruction method based on anisotropic total variation | |
US20070053477A1 (en) | Method and apparatus of global de-noising for cone beam and fan beam CT imaging | |
Abu Anas et al. | Comparison of ring artifact removal methods using flat panel detector based CT images | |
US9401009B2 (en) | Method and apparatus for enhancing quality of 3D image | |
CN110717956A (zh) | 一种有限角投影超像素引导的l0范数最优化重建方法 | |
CN101980302A (zh) | 投影数据恢复导引的非局部平均低剂量ct重建方法 | |
CN109920020A (zh) | 一种锥束ct病态投影重建伪影抑制方法 | |
CN109523458B (zh) | 一种结合稀疏诱导动态引导滤波的高精度稀疏角度ct重建方法 | |
Jeevakala | Sharpening enhancement technique for MR images to enhance the segmentation | |
CN111161182A (zh) | Mr结构信息约束的非局部均值引导的pet图像部分容积校正方法 | |
CN107464219A (zh) | 连续图像的运动检测及降噪方法 | |
CN114708352A (zh) | 一种基于无监督学习的低剂量ct图像恢复方法及系统 | |
CN113017662B (zh) | 一种ct图像的混叠伪影去除方法及系统、ct扫描仪 | |
CN105844678A (zh) | 基于全广义变分正则化的低剂量x射线ct图像重建方法 | |
Arata et al. | Correction of organ motion in SPECT using reprojection data | |
Vani et al. | Multi focus and multi modal image fusion using wavelet transform | |
CN114359126A (zh) | 三维图像的几何校正方法、装置、电子设备及存储介质 | |
CN108734674B (zh) | 一种改进nas-rif的oct图像盲复原方法 | |
Xu et al. | Hessian-based vessel enhancement combined with directional filter banks and vessel similarity | |
CN116029934A (zh) | 一种低剂量dr图像和ct图像去噪方法 | |
Saoji et al. | Speckle and rician noise removal from medical images and Ultrasound images | |
CN112785520A (zh) | Ct图像伪影去除方法 | |
Daniel et al. | Retinal image enhancement using wavelet domain edge filtering and scaling |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |