CN103971387B - Ct图像重建方法 - Google Patents
Ct图像重建方法 Download PDFInfo
- Publication number
- CN103971387B CN103971387B CN201310032989.XA CN201310032989A CN103971387B CN 103971387 B CN103971387 B CN 103971387B CN 201310032989 A CN201310032989 A CN 201310032989A CN 103971387 B CN103971387 B CN 103971387B
- Authority
- CN
- China
- Prior art keywords
- data
- image
- projection
- uproar
- carried out
- 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
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
本发明提供一种CT图像重建的方法，包括如下步骤：根据CT扫描获得原始数据；对原始数据进行校正获得校正数据；对校正数据进行除噪处理获得除噪数据；对除噪数据进行重排获得重排数据；对重排数据进行卷积获得卷积数据；对卷积数据进行反投影处理获得全扫描视野图像与临床选择视野图像；根据全扫描视野图像获得误差图像；根据临床选择视野图像与误差图像获得CT图像。本发明的方法通过全扫描视野图像获得误差图像，并根据临床选择视野图像与误差图像获得CT图像，可以保证重建的CT图像的分辨率，并且可以大大减少重建过程中的计算量，提高图像重建的速度。
Description
【技术领域】
本发明涉及医学影像处理的技术领域，尤其涉及一种CT图像重建的方法。
【技术背景】
在对被扫描对像进行CT扫描的时候会有因被扫描对象携带有金属物品而导致根据CT扫描获得的数据进行重建的图像在与含有金属物品的扫描部位对应的区域产生条状、带状的伪影(即金属伪影)，进而影响医生对被扫描对象的疾病的诊断，甚至出现误诊。
现有的一种CT图像重建中去除金属伪影的步骤如下：根据CT扫描获得原始生数据；根据原始生数据进行全扫描视野图像重建获得全扫描视野图像；对全扫描视野图像进行正投影获得原始正投影数据；对全扫描视野图像进行分割获得含有金属区域的图像；对含有金属区域的图像进行正投影获得金属投影数据并根据金属投影数据确定金属投影区域；对金属投影数据进行线性插值处理获得修正的金属投影数据；根据修正的金属投影数据进行图像重建获得修正图像；对修正图像进行图像建模获得模型图像；对模型图像进行正投影处理获得模型投影数据；在金属投影区域内用模型投影数据对原始正投影数据进行修正获得修正后的投影数据；根据修正后的投影数据进行图像重建获得CT图像。
该CT图像重建中去除金属伪影的方法是一种纯图像后处理的方法，只在图像域进行去除金属伪影处理，虽然比较容易实现去除金属伪影的目的，但由于该方法重建的图像是在金属投影区域内用模型投影数据对原始正投影进行修正后再进行图像重建获得的，因对原始正投影数据进行滤波反投影处理获得的全扫描视野的图像的分辨率远不如根据原始生数据重建的全扫描视野的图像的分辨率，造成了重建CT图像分辨率的损失，进而降低了最终重建的CT图像的图像质量。
目前还有一种前处理的CT图像重建去除金属伪影的方法，该方法的步骤与上述CT图像重建中去除金属伪影的步骤大致相同，区别之处在于：该前处理的方法是在金属投影区域内用与原始生数据对应的模型投影数据对原始生数据进行修正后再进行图像重建的，而不是在金属投影区域内采用模型投影数据对原始正投影数据进行修正后再进行图像重建的。该前处理的CT图像重建去除金属伪影的方法因在金属投影区域内采用与全扫描视野图像相同几何的模型图像进行正投影获得的模型投影数据对原始生数据进行修正的，虽然可以保证不损失重建的CT图像的分辨率，保证重建的CT图像的图像质量，但要获得与全扫描视野图像相同几何的含有金属区域的图像、与全扫描视野图像相同几何的模型图像的精度要求比较高，实现起来比较困难，并且该方法需要在投影域和数据域之间来回切换，也会大大降低CT图像重建的速度。另外，该方法如果用在CT螺旋扫描时也会因为边缘图像很难获得准确的正投影数据；如果用在3D正投影时也需要获得较多的图像才能进行正投影。
因此，确有必要提供一种CT图像重建方法，用于克服现有技术存在的缺陷。
【发明内容】
本发明的目的在于提供一种CT图像重建的方法，可以在保证重建图像的图像质量的基础上，提高图像重建的速度，并可以比较容易地去除图像中的金属伪影。
为达到上述目的，本发明是通过如下技术方案实现的：一种CT图像重建的方法，包括如下步骤：根据CT扫描获得原始数据；对原始数据进行校正获得校正数据；对校正数据进行除噪处理获得除噪数据；对除噪数据进行重排获得重排数据；对重排数据进行卷积获得卷积数据；对卷积数据进行反投影处理获得全扫描视野图像与临床选择视野图像；根据全扫描视野图像获得误差图像；根据临床选择视野图像与误差图像获得CT图像。
获得所述误差图像的过程包括如下步骤：对全扫描视野图像进行正投影获得原始投影数据；对全扫描视野图像进行分割获得金属图像；对金属图像进行正投影获得金属投影数据，并根据金属投影数据确定金属投影区域；根据原始投影数据与金属投影数据获得非金属投影数据；对非金属投影数据中与金属投影区域对应的位置进行修正获得首次修正的投影数据；使用全扫描视野对首次修正的投影数据进行重建获得修正图像；根据修正图像进行图像建模获得模型图像；对模型图像进行正投影获得模型投影数据；使用模型投影数据对原始投影数据进行修正获得二次修正的投影数据；根据非金属投影数据与二次修正的投影数据获得误差投影数据；使用临床选择视野和临床重建中心对误差投影数据进行重建获得误差图像。
优选地，所述CT图像是通过从临床选择视野图像中去除误差图像的方式获得的。
优选地，所述非金属投影数据是通过从所述原始投影数据中去除金属投影数据的方法获得的。
优选地，所述误差投影数据是通过从所述二次修正的投影数据中去除所述非金属投影数据的方法获得的。
优选地，还包括步骤：在获得金属图像时对金属图像之外的区域进行置零处理。
优选地，所述除噪数据是通过对所述校正数据进行非线性滤波除噪处理获得的。
优选地，所述除噪数据是通过对校正数据进行基于全变分的图像域去噪处理获得的。
优选地，所述除噪数据n(γ,ξ)是满足公式
时所对应的除噪数据n(γ,ξ)，其中，minimize{}是最小值函数，即minimize{}括号中的部分取得最小值，γ是探测器的通道数，ξ是探测器的排数，是全变分，n0(γ,ξ)代表校正数据，(δn(γ,ξ))x代表对校正数据n0(γ,ξ)的噪声估计，β是结合噪声估计用来控制去噪力量的参数，x为2或2.5或3。
优选地，所述除噪数据是通过对所述校正数据进行线性滤波除噪处理获得的。
本发明的CT图像重建的方法通过全扫描视野图像获得误差图像，并从临床选择视野图像中去除误差图像获得CT图像，可以保证重建的CT图像的分辨率，并且能大大减少计算量，提高图像重建的速度。
【附图说明】
图1为本发明CT图像重建的方法的流程示意图。
图2为本发明CT图像重建方法中获得误差图像的流程示意图。
图3示出采用本发明的方法去除金属伪影获得的CT图像与未进行金属伪影去除处理重建的CT图像的对比图。
【具体实施方式】
以下结合附图和具体实施例对本发明的CT图像重建方法作进一步详细说明。根据下面说明和权利要求书，本发明的优点和特征将更清楚。需说明的是，附图均采用非常简化的形式且均使用非精准的比率，仅用于方便、明晰地辅助说明本发明实施例的目的。
本发明提供一种CT图像重建方法通过全扫描视野图像获得误差图像，并从临床选择视野图像中去除误差图像获得CT图像，可以保证重建的CT图像的分辨率，并且能大大减少计算量，提高图像重建的速度。
图1为本发明CT图像重建方法的流程示意图。该CT图像重建方法包括以下步骤：
S11、根据CT扫描获得原始数据；
S12、对原始数据进行校正获得校正数据，所述校正数据是由原始数据依次进行空气校正、焦点校正、非线性校正及CT值校正获得的；
S13、对校正数据进行除噪处理获得除噪数据，所述除噪数据是通过对校正数据进行非线性滤波除噪处理获得的，具体地，除噪数据是通过对校正数据进行基于全变分的图像域去噪处理获得的；
S14、对除噪数据进行重排获得重排数据，所述重排数据是通过对除噪数据依次进行角向重排、径向重排获得的；
S15、对重排数据进行卷积获得卷积数据；
S16、对卷积数据进行反投影处理获得全扫描视野图像与临床选择视野图像；
S17、根据全扫描视野图像获得误差图像；
S18、根据临床选择视野图像与误差图像获得CT图像，所述CT图像是通过从临床选择视野图像中去除误差图像的方式获得的。
具体地，步骤S13中的除噪数据n(γ,ξ)是满足公式 时所对应的除噪数据n(γ,ξ)，其中，γ是探测器的通道数，ξ是探测器的排数，是全变分，n0(γ,ξ)为校正数据，(δn(γ,ξ))x为对校正数据n0(γ,ξ)的噪声估计，β是结合噪声估计用来控制去噪力量的参数，x优选参数为2或2.5或3。
所述除噪数据还可以通过对校正数据进行线性滤波除噪处理获得。
图2为图1的步骤S17获得误差图像的流程示意图。该获得误差图像的过程包括如下步骤：
S20、对全扫描视野图像进行正投影获得原始投影数据；
S21、对全扫描视野图像进行分割获得金属图像并对金属图像之外的区域进行置零处理；
S22、对金属图像进行正投影获得金属投影数据，并根据金属投影数据确定金属投影区域；
S23、根据原始投影数据与金属投影数据获得非金属投影数据，所述非金属投影数据是通过从所述原始投影数据中去除金属投影数据的方法获得的；
S24、对非金属投影数据中与金属投影区域对应的位置进行修正获得首次修正的投影数据；
S25、使用全扫描视野对首次修正的投影数据进行重建获得修正图像；
S26、根据修正图像进行图像建模获得模型图像；
S27、对模型图像进行正投影获得模型投影数据；
S28、使用模型投影数据对原始投影数据进行修正获得二次修正的投影数据；
S29、根据非金属投影数据与二次修正的投影数据获得误差投影数据，所述误差投影数据是通过从所述二次修正的投影数据中去除所述非金属投影数据的方法获得的；
S30、使用临床选择视野和临床重建中心对误差投影数据进行重建获得误差图像。
图3示出采用本发明的方法重建的CT图像与未采用本发明的方法重建的CT图像的对比图。本实施例中，最大扫描视野即全扫描视野为500mm，临床选择的视野为350mm，探测器的排数为16，探测器的通道数为816。图3a为未进行金属伪影去除处理重建的CT图像。图3b为采用本发明的方法重建的CT图像。从图3a与图3b图中可以很明显的看出，图3b中采用本发明的方法重建的CT图像具有较高的图像分辨率，并且不具有条状、带状的伪影(即金属伪影)。
综上，以上仅为本发明的较佳实施例而已，不应以此限制本发明的范围，即凡是依本发明的权利要求书及本发明说明书内容所作的简单的等效变化与修饰，均应仍属本发明专利涵盖的范围内。
Claims (9)
1.一种CT图像重建的方法，其特征在于，包括如下步骤：
根据CT扫描获得原始数据；
对原始数据进行校正获得校正数据；
对校正数据进行除噪处理获得除噪数据；
对除噪数据进行重排获得重排数据；
对重排数据进行卷积获得卷积数据；
对卷积数据进行反投影处理获得全扫描视野图像与临床选择视野图像；
对全扫描视野图像进行正投影获得原始投影数据；
对全扫描视野图像进行分割获得金属图像；
对金属图像进行正投影获得金属投影数据，并根据金属投影数据确定金属投影区域；
根据原始投影数据与金属投影数据获得非金属投影数据；
对非金属投影数据中与金属投影区域对应的位置进行修正获得首次修正的投影数据；
使用全扫描视野对首次修正的投影数据进行重建获得修正图像；
根据修正图像进行图像建模获得模型图像；
对模型图像进行正投影获得模型投影数据；
使用模型投影数据对原始投影数据进行修正获得二次修正的投影数据；
根据非金属投影数据与二次修正的投影数据获得误差投影数据；
使用临床选择视野和临床重建中心对误差投影数据进行重建获得误差图像；
根据临床选择视野图像与误差图像获得CT图像。
2.如权利要求1所述的CT图像重建的方法，其特征在于：所述CT图像是通过从临床选择视野图像中去除误差图像的方式获得的。
3.如权利要求2所述的CT图像重建的方法，其特征在于：所述非金属投影数据是通过从所述原始投影数据中去除金属投影数据的方法获得的。
4.如权利要求3所述的CT图像重建的方法，其特征在于：所述误差投影数据是通过从所述二次修正的投影数据中去除所述非金属投影数据的方法获得的。
5.如权利要求2至4中任意一项所述的CT图像重建的方法，还包括步骤：在获得金属图像时对金属图像之外的区域进行置零处理。
6.如权利要求5所述的CT图像重建的方法，其特征在于：所述除噪数据是通过对所述校正数据进行非线性滤波除噪处理获得的。
7.如权利要求6所述的CT图像重建的方法，其特征在于：所述除噪数据是通过对校正数据进行基于全变分的图像域去噪处理获得的。
8.如权利要求7所述的CT图像重建的方法，其特征在于：所述除噪数据n(γ,ξ)是满足公式时对应的除噪数据n(γ,ξ)，其中，γ是探测器的通道数，ξ是探测器的排数，是全变分，n0(γ,ξ)代表校正数据，(δn(γ,ξ))x代表对校正数据n0(γ,ξ)的噪声估计，β是结合噪声估计用来控制去噪力量的参数，x为2或2.5或3。
9.如权利要求5所述的CT图像重建的方法，其特征在于：所述除噪数据是通过对所述校正数据进行线性滤波除噪处理获得的。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310032989.XA CN103971387B (zh) | 2013-01-29 | 2013-01-29 | Ct图像重建方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310032989.XA CN103971387B (zh) | 2013-01-29 | 2013-01-29 | Ct图像重建方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103971387A CN103971387A (zh) | 2014-08-06 |
| CN103971387B true CN103971387B (zh) | 2017-10-27 |
Family
ID=51240836
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310032989.XA Active CN103971387B (zh) | 2013-01-29 | 2013-01-29 | Ct图像重建方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103971387B (zh) |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104574292B (zh) * | 2014-11-26 | 2018-06-26 | 沈阳东软医疗系统有限公司 | 一种ct图像的校正方法和装置 |
| CN106548453B (zh) * | 2015-09-21 | 2021-03-16 | 上海联影医疗科技股份有限公司 | Pet图像重建方法及系统 |
| US10049449B2 (en) | 2015-09-21 | 2018-08-14 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for image reconstruction |
| CN105748105A (zh) * | 2016-01-18 | 2016-07-13 | 刘焕涛 | 一种乳腺检查系统 |
| CN105852969A (zh) * | 2016-03-29 | 2016-08-17 | 鞠克丰 | 一种新型骨科手术导航系统 |
| CN106175809A (zh) * | 2016-06-30 | 2016-12-07 | 王礼泉 | 一种乳腺检查系统 |
| CN107066690A (zh) * | 2017-03-02 | 2017-08-18 | 平顶山天安煤业股份有限公司 | 一种基于计算机的图像计算钻孔深度控制系统 |
| CN107424140A (zh) * | 2017-03-02 | 2017-12-01 | 平顶山天安煤业股份有限公司 | 一种基于全景遥视成像和钻孔轨迹测量控制系统 |
| CN107978003B (zh) * | 2017-10-30 | 2021-03-30 | 东软医疗系统股份有限公司 | 一种ct图像金属伪影的处理方法及装置 |
| WO2019090533A1 (en) | 2017-11-08 | 2019-05-16 | Shenzhen United Imaging Healthcare Co., Ltd. | System and method for diagnositic and treatment |
| CN109448070B (zh) * | 2018-11-06 | 2020-10-09 | 明峰医疗系统股份有限公司 | 一种ct宽展视野的重建方法 |
| CN110084866B (zh) * | 2019-04-24 | 2023-02-17 | 上海联影医疗科技股份有限公司 | 一种计算机断层成像方法和装置 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101675884A (zh) * | 2008-09-19 | 2010-03-24 | Ge医疗系统环球技术有限公司 | 设置影像重建中心参数的方法及其ct扫描设备 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE10036142B4 (de) * | 2000-07-25 | 2004-04-29 | Siemens Ag | Röntgen-Computertomographieeinrichtung |
-
2013
- 2013-01-29 CN CN201310032989.XA patent/CN103971387B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101675884A (zh) * | 2008-09-19 | 2010-03-24 | Ge医疗系统环球技术有限公司 | 设置影像重建中心参数的方法及其ct扫描设备 |
Non-Patent Citations (2)
| Title |
|---|
| A self-adapting fast reconstruction method for high resolution images in cone beam CT;Mingjun Li等;《The Fifth International Conference on Image and Graphics》;20100830;第23-27页 * |
| 扇束CT几何伪影的校正方法;王亮等;《电子学报》;20110515;第39卷(第5期);第1143--1149页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103971387A (zh) | 2014-08-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103971387B (zh) | Ct图像重建方法 | |
| CN106683144B (zh) | 一种图像迭代重建方法及装置 | |
| JP7330703B2 (ja) | 医用画像処理装置及びｘ線ｃｔシステム | |
| US9349198B2 (en) | Robust artifact reduction in image reconstruction | |
| CN106600614B (zh) | 基于凹凸性的sd-oct视网膜图像cnv分割方法 | |
| CN105225208B (zh) | 一种计算机断层成像金属伪影校正方法及装置 | |
| EP2783344B1 (en) | Image domain de-noising | |
| CN102737395B (zh) | 一种医用x射线系统中的图像处理方法及装置 | |
| La Riviere et al. | Reduction of noise-induced streak artifacts in X-ray computed tomography through spline-based penalized-likelihood sinogram smoothing | |
| CN103136773B (zh) | 一种稀疏角度x射线ct成像方法 | |
| CN105046675B (zh) | 基于mvct图像消除kvct图像中金属伪影的方法 | |
| CN105528766B (zh) | Ct金属伪影处理方法和装置 | |
| CN106023200A (zh) | 一种基于泊松模型的x光胸片图像肋骨抑制方法 | |
| CN111223156B (zh) | 一种用于牙科锥束ct系统的金属伪影消除方法 | |
| JP2014517713A5 (zh) | ||
| CN103116879A (zh) | 一种基于邻域加窗的非局部均值ct成像去噪方法 | |
| CN110310346A (zh) | 一种ct和cbct图像中金属伪影校正方法 | |
| CN106530236B (zh) | 一种医学图像处理方法及系统 | |
| CN105528800A (zh) | 一种计算机断层成像伪影校正方法及装置 | |
| CN106296615A (zh) | Ct图像中校正金属伪影的方法和系统 | |
| CN105761226B (zh) | 一种超视野ct扫描图像的补偿重建方法 | |
| EP2096604A2 (en) | Computed tomography reconstruction from truncated scans | |
| CN112330565A (zh) | 基于改善的U-net的低剂量CT投影域中图像去噪方法 | |
| KR20170073640A (ko) | 다중 수집을 통해 콘트라스트를 정규화하는 방법 및 시스템 | |
| Li et al. | Metal artifact reduction in CT based on adaptive steering filter and nonlocal sinogram inpainting |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 201815 No. 1180 Xingxian Road, Jiading Industrial Zone, Jiading District, Shanghai Patentee after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201815 No. 1180 Xingxian Road, Jiading Industrial Zone, Jiading District, Shanghai Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
| CP02 | Change in the address of a patent holder | ||
| CP02 | Change in the address of a patent holder |
Address after: 201807 2258 Chengbei Road, Jiading District, Shanghai Patentee after: Shanghai Lianying Medical Technology Co.,Ltd. Address before: 201815 No. 1180 Xingxian Road, Jiading Industrial Zone, Jiading District, Shanghai Patentee before: Shanghai Lianying Medical Technology Co.,Ltd. |