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

摘要:

权利要求书:

1.一种CT图像的重建方法,其特征在于,所述方法包括以下步骤;

步骤S101,获取CT扫描所采集的投影数据;

步骤S102,应用成像模型方程式G=MX,采用公式计算出重建图像;其中,G为CT扫描所采集的投影数据,M为系统矩阵,X为重建的CT图像,Xi表示第i次迭代后得到的重建图像,λ表示收敛系数。

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结束;

0

步骤7、由公式(4),计算此时X与步骤4中X 的差异参数η;

步骤8、进行如下判断操作,即如果主循环程序执行的是第一次循环,则给ω赋初值ω=a×η;

0 0

步骤9、将X的值赋给X,即X =X;

步骤10、进入子循环III,

子循环程序III结束;

其中: 为先验图像,一般为全扫描重建得到的CT图像;而||Xs,t,v||可由公式(9)得到, 可由公式(10)得到, 可由公式(10)得到,步骤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范围内取值。