一种基于fbp的能谱ct迭代成像方法和成像系统的制作方法
【技术领域】
[0001] 本发明涉及计算机断层成像技术领域，尤其是一种基于FBP的能谱CT迭代成像方 法和成像系统。
【背景技术】
[0002] 传统的X射线CT成像系统，使用X光机在一个电压下扫描被成像物体，由探测器获 取穿过物体的X射线流强数据，并由该数据对物体进行图像重建。由于不同材料可能会有相 同的CT值，因此传统CT图像对某些物质无法区分。与传统CT成像系统不同，多能谱CT成像系 统能够获得被成像物体在多个不同能谱下的扫描数据，并采用特别的方法进行图像重建。 由重建的CT图像可以获得被成像物体更多物理特性参数的信息，比如：电子密度、等效原子 序数等。与传统CT相比，多能谱CT具有更好的物质区分能力，在医学诊断、无损检测以及安 全检查等应用中都重要应用前景。
[0003] 常见的多能谱CT成像系统是双能谱CT成像系统。获取双能谱CT数据的方法有多 种，常用的有双源双探测器扫描方法、双层探测器扫描方法和电压快速切换扫描方法等。还 有其他获得多能谱数据的方法，如电压慢变切换扫描方法，不同电压下多次扫描方法等。另 外，近年来，光子计数探测器技术迅速发展，利用具有能量分辨的光子计数探测器通过单次 CT扫描也可以获得多组不同能谱下的扫描数据。
[0004] 重建方法方面，目前研究较多的是双能谱CT重建。相关研究始于上个世纪七十年 代。已有的重建方法可分为：前处理重建方法、后处理重建方法和优化迭代重建方法。前处 理重建方法由实际获得的高低两组多能投影数据估计出基材料密度图像的线性投影值，然 后用传统方法进行重建，并获得基材料的密度图像，该类方法的缺点是要求获得的两组多 能投影数据是几何相容的，即沿每条射线路径采集高、低两个能谱的扫描数据。后处理重建 方法是先对实际获得的两组多能投影或由此得到的加权投影分别进行重建，然后对重建图 像进行线性组合以获得基材料的密度图像。通常，这类方法无法完全消除重建图像中由能 谱引起的硬化伪影。优化迭代方法对双能谱CT的数据获取过程进行优化建模，并构造相应 的迭代格式求解。该类方法的缺点是计算量大，速度慢，难以满足某些实时性要求较高的应 用需求。
[0005] 通过检索，发现一篇与本发明专利申请相关的专利公开文献：
[0006] 一种多能谱CT成像方法及成像系统(CN103900931A)，公开了一种多能谱CT成像方 法及成像系统，所述多能谱成像方法包括对被测物体扫描、重建基材料的密度图像最终获 取被测物体的电子密度、等效原子序数分布或被测物体对单能X射线的线性衰减系数分布。 其中多能谱CT重建方法包括:为待重建基材料的密度图像赋初值;对估计图像进行正投影 以获得投影估计值;估计重建图像的误差;利用重建图像误差的估计更新图像;重复上述步 骤，直到重建的密度图像收敛。相应的，所述多能谱CT成像系统包括扫描模块、重建模块和 统计模块，其中所述重建模块又包括初始化模块、投影模块、校正模块和更新模块。
[0007] 通过对比，本发明专利申请与上述专利公开文献存在本质的不同。
【发明内容】
[0008] 针对现有技术中多能谱CT成像方法的不足之处，提供一种基于FBP的能谱CT迭代 成像方法和成像系统，该方法能够由多能谱CT扫描数据重建出高质量、无硬化伪影的CT图 像，同时该方法对多能谱扫描数据无几何相容的要求，另外，该方法计算量小，重建图像速 度快，该系统能够扫描被成像物体，获得多能谱CT扫描数据并得到高质量、无硬化伪影的CT 图像。
[0009] 为了实现上述目的，本发明所采用的技术方案如下：
[0010] -种基于FBP的能谱CT迭代成像方法，该方法对多能谱CT成像系统采集到的多个 能谱CT扫描数据进行重建，得到多能谱CT图像。
[0011] 而且，每次由FBP算法重建残差图像，进而更新估计图像，步骤如下：
[0012] 以Pi，i = 1，2，…，N表示不同能谱下的多能投景多数据，N表示能谱的个数，f i，i = 1， 2，…，N表示待重建的基材料图像，Ai，i = 1，2，…，N表示基材料图像fi的线性投影；
[0013] 以口土二口土以!，A2，…，An)，i = 1，2，…，N和Ai = Ai(pi，p2,…，pn)，i = 1，2，…，N表示基 材料的线性投影与多能投影间的映射关系，该映射关系通过扫描典型模体直接获得，或者 利用X光机的能谱以及基材料的线性衰减系数获得解析的表达式。
[00M] 而且，具体步骤如下：
[0015] ⑴使用多能谱CT成像系统采集多个能谱下的多能谱CT扫描数据，对多能谱CT扫描 数据进行预处理得到不同能谱下的多能投影数据 Pl，i = 1，2，…，N;
[0016] ⑵为待重建基材料的密度图像赋初值，作为初始估计图像，i = 1，2，…，N; [0017]⑶对估计图像#11^〇4 = 1，2，~』进行线性正投影，获得估计的线性投影
[0018] ⑷根据线性投影与多能投影间的映射关系，计算多能投影数据
[0019] (5)计算残差投影Δ/f)=凡U = /V .
[0020] (6)对残差投影分别进行加权得到加权残差投影：
[0021]
〇4
[0022] 其中+表示Aj对Pi的一阶偏导数；
[0023] ⑵对上述加权残差投影利用传统的和成像系统图像重建算法进行重建得到加权 分量残差图像，S
>其中FBP表示滤波反投 影算子；
[0024] (8)将上述加权分量残差图像组合得到残差图像：
[0025]
[0026] (9)将Af')加到估计图像负⑴上，得到新的估计图像负(1{+1)，即 (k),i = l,2,---,N；
[0027] (10)重复步骤⑶-(9)，直到重建图像满足收敛条件，即得多能谱CT成像。
[0028] 而且，所述步骤⑵中选择基于投影或者基于图像的重建方法重建的图像作为初始 估计图像，或者以〇作为初始图像，即图像中每个像素值均为〇。
[0029] 而且，所述步骤⑶中正投影得到的估计线性投影.= …对应的路径应包 括所有实际多能投影数据Pl，i = 1，2，…，N对应的路径；
[0030] 扫描数据Pi，i = 1，2，…，N几何不相容时，同样有效。
[0031 ] 而且，具体步骤如下：
[0032] ⑴使用多能谱CT成像采集多个能谱下的多能谱CT扫描数据，对多能谱CT扫描数据 进行预处理得到不同能谱下的多能投影数据 ?14 = 1，2，'"，1投影数据?14 = 1，2，一少几 何相容；
[0033] ⑵为待重建基材料的密度图像赋初值，作为初始估计图像fi(())，i = 1，2，…，N;
[0034] ⑶对估计图像5(1°汴2〇4 = 1，2，~川进行线性正投影，获得估计的线性投影 4k\i=ia,--,N',
[0035] ⑷根据线性投影与多能投影间的映射关系，计算多能投影数据
[0036] (5)计算残差投影 Apf4 =只-尸,(A 1, i = 1，2, · ", i￥.
[0037] (6)对残差投影进行加权求和，得到组合残差投影：
[0038]
(；A
[0039] 其中f表示Aj对pi的一阶偏导数； cp;
[0040] (7)对上述组合残差投影Δ4('，利用传统的FBP图像重建算法进行重建得到残差图 ??
[0041] (8)将加到估计图像上，得到新的估计图像负~+1)，即 (k),i = l,2,---,N；
[0042] (9)重复步骤⑵_(8)，直到重建图像满足收敛条件，即得多能谱CT成像。
[0043] 而且，所述步骤(7)中加权分量残差图像或残差图像再经如下处理：
[0044] 对加权分量残差图像做降噪处理。
[0045] 一种利用如上所述的基于FBP的能谱CT迭代成像方法的CT成像系统，其特征在于： 所述系统包括载物台、X射线放生装置、探测器和控制和数据处理装置，所述载物台能够承 载被成像物体，所述X射线放生装置能够发出X射线照射被测物体，所述探测器能够接收X射 线放生装置发出的X射线，获得扫描数据;所述控制和数据处理装置能够对所述的载物台、