CN103065342A - 一种ct图像的重建方法 - Google Patents
一种ct图像的重建方法 Download PDFInfo
- Publication number
- CN103065342A CN103065342A CN2012105607929A CN201210560792A CN103065342A CN 103065342 A CN103065342 A CN 103065342A CN 2012105607929 A CN2012105607929 A CN 2012105607929A CN 201210560792 A CN201210560792 A CN 201210560792A CN 103065342 A CN103065342 A CN 103065342A
- Authority
- CN
- China
- Prior art keywords
- image
- program
- centerdot
- formula
- subcycle
- 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
Links
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Description
技术领域
本发明属于医学和工业领域计算机断层扫描系统领域,尤其涉及一种CT图像的重建方法。
背景技术
计算机断层成像(CT)是通过无损方式获取物体内部结构信息的一种重要成像手段,它拥有高分辨率、高灵敏度以及多层次等众多优点,是我国装机量最大的医疗影像诊断设备之一,被广泛应用于各个医疗临床检查领域。然而,由于CT扫描过程中需要使用X射线,随着人们对辐射潜在危害的逐步了解,CT辐射剂量问题越来越受到人们的重视。合理使用低剂量(As Low AsReasonably Achievable,ALARA)原则要求在满足临床诊断的前提下,尽量降低对患者的辐射剂量。因此,研究和开发新的低剂量CT成像方法,既能保证CT成像质量又减少有害的辐射剂量,对于医疗诊断领域具有重要的科学意义和应用前景。
现有技术提出了一种稀疏角度CT图像的重建方法。该方法是:获取CT机的系统参数和一周期内的稀疏角度投影数据,并将所获得的稀疏角度投影数据按采集的时间顺序等分成若干组;应用重建模型采用辅助函数法求解得到要求的迭代运算方法进行循环迭代运算,当循环次数达到预设的次数时即停止迭代运算,并将当前所得到的迭代运算结果作为重建图像。
然而,上述现有技术方案由于要进行全采样,导致CT扫描时间长;由于全采样所采集的数据量较大,导致图像重建速度慢;由于扫描时间长,导致出现病人运动所引起的伪影;病人所承受的CT辐射剂量高。
发明内容
本发明的目的在于提供一种CT图像的重建方法,旨在解决现有技术的稀疏角度CT图像的重建方法由于要进行全采样,导致CT扫描时间长;由于全采样所采集的数据量较大,导致图像重建速度慢;由于扫描时间长,导致出现病人运动所引起的伪影;病人所承受的CT辐射剂量高的问题。
本发明是这样实现的,一种CT图像的重建方法,所述方法包括以下步骤;
步骤S101,获取CT扫描所采集的投影数据;
步骤S102,应用成像模型方程式G=MX,采用 公式计算出重建图像;其中,G为CT扫描所采集的投影数据,M为系统矩阵,X为重建的CT图像,Xi表示第i次迭代后得到的重建图像,λ表示收敛系数。
在本发明中,本发明在CT扫描时,进行高欠采样数据采集可以有效缩短CT扫描时间、所采集的数据量较小,从而提高了图像重建速度、由于扫描时间短,从而减少病人运动所引起的伪影,能在保证CT图像质量的前提下大幅降低病人所承受的CT辐射剂量。
附图说明
图1是本发明实施例提供的CT图像的重建方法的实现流程示意图。
具体实施方式
为了使本发明的目的、技术方案及有益效果更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
请参阅图1,为本发明实施例提供的CT图像的重建方法的实现流程,其包括以下步骤;
在步骤S101中,获取CT扫描所采集的投影数据;
具体的CT图像的重建方法如下:
1、定义成像模型方程如公式(1),其中G为CT扫描所采集的投影数据,M为系统矩阵,X为重建的CT图像。定义参数a=0.2;b=0.9,b一般取0-1之间的数值。特殊的,b也可直接取值为0,以用来简化程序;c=0.96,λ=1,d=0.993。
G=MX (1)
2、将重建图像X的初始值赋为0,即X=0;
3、进入主循环程序,将程序执行i=1,L,istepnum次,一般取istepnum为20到100之间的整数;
4、将X的值赋给X0,即X0=X;
5、进入子循环程序I,按公式(2)计算重建图像X。其中Xi表示第i次迭代后得到的重建图像;λ表示收敛系数,一般取小于1的正数,按实际经验一般在0.9-1.0范围内取值。
子循环程序I结束;
6、进入子循环程序II,对重建图像X内的所有点依次进行如下操作
子循环程序II结束;
7、由公式(4),计算此时X与步骤4中X0的差异参数η
8、进行如下判断操作,即如果主循环程序执行的是第一次循环,则给ω赋初值ω=a×η;
9、将X的值赋给X0,即X0=X;
10、进入子循环III,
子循环程序III结束;
11、再次给ω和λ赋值
ω=ω×c (11)
λ=λ×d (12)
12、返回主循环程序。
综上所述,本发明实施例在CT扫描时,进行高欠采样数据采集可以有效缩短CT扫描时间、所采集的数据量较小,从而提高了图像重建速度、由于扫描时间短,从而减少病人运动所引起的伪影,能在保证CT图像质量的前提下大幅降低病人所承受的CT辐射剂量。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,所述的程序可以存储于一计算机可读取存储介质中,所述的存储介质,如ROM/RAM、磁盘、光盘等。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (7)
2.如权利要求1所述的方法,其特征在于,步骤S102具体包括以下步骤:步骤1、定义成像模型方程如公式(1),其中G为CT扫描所采集的投影数据,M为系统矩阵,X为重建的CT图像;
G=MX (1)
步骤2,将重建图像X的初始值赋为0,即X=0;
步骤3、进入主循环程序,将程序执行i=1,L,istepnum次;
步骤4、将X的值赋给X0,即X0=X;
步骤5、进入子循环程序I,按公式(2)计算重建图像X;其中Xi表示第i次迭代后得到的重建图像;λ表示收敛系数;
子循环程序I结束;
步骤6、进入子循环程序II,对重建图像X内的所有点依次进行如下操作
子循环程序II结束;
步骤7、由公式(4),计算此时X与步骤4中X0的差异参数η;
步骤8、进行如下判断操作,即如果主循环程序执行的是第一次循环,则给ω赋初值ω=a×η;
步骤9、将X的值赋给X0,即X0=X;
步骤10、进入子循环III,
子循环程序III结束;
步骤11、再次给ω和λ赋值
ω=ω×c (11)
λ=λ×d (12)
步骤12、返回主循环程序。
3.如权利要求2所述的方法,其特征在于,定义参数a=0.2;b取0-1之间的数值;c=0.96,λ=1,d=0.993。
4.如权利要求3所述的方法,其特征在于,所述b=0.9。
5.如权利要求2所述的方法,其特征在于,所述取istepnum为20到100之间的整数。
6.如权利要求2所述的方法,其特征在于,所述λ取小于1的正数。
7.如权利要求6所述的方法,其特征在于,所述λ在0.9-1.0范围内取值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210560792.9A CN103065342B (zh) | 2012-12-21 | 2012-12-21 | 一种ct图像的重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210560792.9A CN103065342B (zh) | 2012-12-21 | 2012-12-21 | 一种ct图像的重建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103065342A true CN103065342A (zh) | 2013-04-24 |
CN103065342B CN103065342B (zh) | 2016-10-19 |
Family
ID=48107958
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210560792.9A Active CN103065342B (zh) | 2012-12-21 | 2012-12-21 | 一种ct图像的重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103065342B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104183286A (zh) * | 2013-05-23 | 2014-12-03 | 环境保护部核与辐射安全中心 | 用于堆芯熔融物状态监测的图像重建方法、装置和系统 |
CN104933744A (zh) * | 2015-06-19 | 2015-09-23 | 深圳先进技术研究院 | Ct图像重建方法和系统 |
CN104992457A (zh) * | 2015-06-19 | 2015-10-21 | 深圳先进技术研究院 | Ct图像重建方法和系统 |
WO2018133003A1 (zh) * | 2017-01-19 | 2018-07-26 | 深圳先进技术研究院 | Ct三维重建方法及系统 |
CN109717886A (zh) * | 2017-10-30 | 2019-05-07 | 上海交通大学 | 一种低辐射剂量的ct扫描方法 |
CN111260647A (zh) * | 2020-03-12 | 2020-06-09 | 南京安科医疗科技有限公司 | 一种基于图像检测的ct扫描辅助方法、计算机可读存储介质、ct扫描装置 |
CN111724452A (zh) * | 2020-06-17 | 2020-09-29 | 深圳高性能医疗器械国家研究院有限公司 | 一种低剂量ct图像重建方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060050838A1 (en) * | 2004-09-01 | 2006-03-09 | Ge Medical Systems Global Technology Company, Llc | CT image reconstruction method, CT apparatus, and program |
CN102063728A (zh) * | 2010-11-23 | 2011-05-18 | 南方医科大学 | 基于标准剂量图像冗余信息的低剂量ct图像重建方法 |
CN102663790A (zh) * | 2012-05-08 | 2012-09-12 | 南方医科大学 | 一种稀疏角度ct图像的重建方法 |
-
2012
- 2012-12-21 CN CN201210560792.9A patent/CN103065342B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060050838A1 (en) * | 2004-09-01 | 2006-03-09 | Ge Medical Systems Global Technology Company, Llc | CT image reconstruction method, CT apparatus, and program |
CN102063728A (zh) * | 2010-11-23 | 2011-05-18 | 南方医科大学 | 基于标准剂量图像冗余信息的低剂量ct图像重建方法 |
CN102663790A (zh) * | 2012-05-08 | 2012-09-12 | 南方医科大学 | 一种稀疏角度ct图像的重建方法 |
Non-Patent Citations (2)
Title |
---|
LI JIAN ET AL.: "Rotating polar-coordinate ART applied in industrial CT image reconstruction", 《NDT&E INTERNATIONAL》 * |
李毅等: "一种基于数据外插改进的ART迭代算法", 《CT理论与应用研究》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104183286A (zh) * | 2013-05-23 | 2014-12-03 | 环境保护部核与辐射安全中心 | 用于堆芯熔融物状态监测的图像重建方法、装置和系统 |
CN104933744A (zh) * | 2015-06-19 | 2015-09-23 | 深圳先进技术研究院 | Ct图像重建方法和系统 |
CN104992457A (zh) * | 2015-06-19 | 2015-10-21 | 深圳先进技术研究院 | Ct图像重建方法和系统 |
CN104992457B (zh) * | 2015-06-19 | 2018-07-27 | 深圳先进技术研究院 | Ct图像重建方法和系统 |
CN104933744B (zh) * | 2015-06-19 | 2018-11-20 | 深圳先进技术研究院 | Ct图像重建方法和系统 |
WO2018133003A1 (zh) * | 2017-01-19 | 2018-07-26 | 深圳先进技术研究院 | Ct三维重建方法及系统 |
CN109717886A (zh) * | 2017-10-30 | 2019-05-07 | 上海交通大学 | 一种低辐射剂量的ct扫描方法 |
CN111260647A (zh) * | 2020-03-12 | 2020-06-09 | 南京安科医疗科技有限公司 | 一种基于图像检测的ct扫描辅助方法、计算机可读存储介质、ct扫描装置 |
CN111724452A (zh) * | 2020-06-17 | 2020-09-29 | 深圳高性能医疗器械国家研究院有限公司 | 一种低剂量ct图像重建方法 |
CN111724452B (zh) * | 2020-06-17 | 2024-05-28 | 深圳高性能医疗器械国家研究院有限公司 | 一种低剂量ct图像重建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103065342B (zh) | 2016-10-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103065342A (zh) | 一种ct图像的重建方法 | |
US20230210477A1 (en) | Positron emission tomography imaging system and method | |
US9155514B2 (en) | Reconstruction with partially known attenuation information in time of flight positron emission tomography | |
CN101499173B (zh) | 一种pet成像中卡尔曼滤波图像重建方法 | |
CN107095691B (zh) | 一种pet成像方法及系统 | |
CN106232009B (zh) | 断层扫描设备和由断层扫描设备重构断层扫描图像的方法 | |
Cheng et al. | Applications of artificial intelligence in nuclear medicine image generation | |
Wu et al. | Spatial-temporal total variation regularization (STTVR) for 4D-CT reconstruction | |
CN105147312A (zh) | Pet图像获取方法及系统 | |
CN103282941A (zh) | 利用门控信号检测和校正列表模式pet数据中的运动的方法和装置 | |
US8913710B2 (en) | Truncation correction imaging enhancement method and system | |
CN105608717A (zh) | 一种ct系统和ct图像重建方法 | |
CN103106676B (zh) | 一种基于低剂量投影数据滤波的x射线ct图像重建方法 | |
CN101917906A (zh) | 通过将物体的围绕物用作动态约束在断层扫描中减少剂量并增强图像 | |
US20240177377A1 (en) | Spect imaging prediction model creation method and apparatus, and device and storage medium | |
CN102663790B (zh) | 一种稀疏角度ct图像的重建方法 | |
CN112568923A (zh) | X射线相位衬度图像提取方法、装置、终端及存储介质 | |
CN103810734A (zh) | 一种低剂量x射线ct投影数据恢复方法 | |
Ma et al. | An encoder-decoder network for direct image reconstruction on sinograms of a long axial field of view PET | |
JP2015066445A (ja) | コンピュータ断層撮影(ct)画像再構成におけるアーチファクトを低減するための方法および装置 | |
Xue et al. | A 3D attention residual encoder–decoder least-square GAN for low-count PET denoising | |
CN103793890A (zh) | 一种能谱ct图像的恢复处理方法 | |
Gavriilidis et al. | Positron range corrections and denoising techniques for gallium-68 PET imaging: a literature review | |
CN101268950B (zh) | 基于cell宽频引擎的螺旋ct精确重建系统 | |
Santarelli et al. | CZT detectors-based SPECT imaging: how detector and collimator arrangement can determine the overall performance of the tomograph |
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 |