专利名称：高能双能ct系统的图象重建方法
技术领域：
本发明涉及辐射成像，具体涉及一种高能X射线双能CT图像重 建方法。
背景技术：
计算机断层成像技术(CT)在医疗诊断和工业无损检测领域应用 广泛。随着社会的发展，它在公共安全和社会保障的需求也逐渐增加， 其中包括了大量以X射线成像为基础的CT成像系统，简称XCT。这 些XCT系统按使用的成像技术可以分为单能和双能CT。虽然多能X 射线成像技术已经有了一定的发展，但是在实际的应用中还是以单、 双能XCT为主，而且技术相对成熟。单能XCT通过重建物体断层的 衰减系数图像，得到物质内部的结构和物理信息，但是无法准确分辨 是何种物质。双能XCT不仅能够得到物质内部衰减系数的信息，也可 以通过重建得到物质组成的信息， 一种典型的结果就是物质的有效原 子序数和等效特征密度。这样就可以进行准确的物质识别，为公共安 全领域提供一种高效的检查手段。
双能XCT技术己经在医学成像和小型物件的安全检查领域得到 了大量的应用，其技术已经相对比较成熟。但是我们看到，在大部分 的双能系统中，都采用的低能X射线进行成像， 一般在200KeV以下。 选择这样的能量，其原因有以下几条低能X射线能较为容易通过X 光管产生，其辐射防护也相对较为简便，而且物质在这个能段的衰减 系数差异较大，所以物质的图像区分性也比较好；另外由于扫描的物 体往往体积较小，物体对X射线的衰减较少，所以低能X射线双能技 术在这些领域广泛应用。
但是在大型货物的安全检查领域，这样能量范围的X射线的穿透力是远远不够，无法得到清晰有效的投影数据用于图像重建。对于大
型货物的X射线透视成像， 一般需要MeV量级的X射线，能量范围 一般是1 10MeV。在这个能量范围内，用于传统低能双能重建方法不 再适用，其根本原因在于这些方法都是基于一个事实，即X射线与物 质的作用只有两种效应光电效应和康普顿散射，而不包含电子对效 应。但是在高能双能领域，由于X射线的能量往往大于电子对效应发 生的最低能量1.02MeV，所以在两种效应基础上提出的原有方法都无 法继续使用，而必须开发新的技术来满足这样的要求。
发明内容
本发明的目的是提出一种基于高能X射线(大于1MeV)的双能 CT重建方法，以便解决大型货物双能XCT检査的问题。根据本发明 实施例的方法，能够准确有效地得到原子序数和特征密度断层图像从 而进行材料识别，为大型货物的安全检查提供了一种高效的手段。
本发明的CT双能系统使用可以获取双能信息的射线源和探测 器，可以采用标准的CT圆轨道或者其他数据釆集方式得到投影数据， 并利用这些数据重建断层图像。
在本发明的一个方面，提供了一种在高能双能CT系统中的图像 重建方法，包括步骤利用高能双能射线扫描物体，以获得高能双能 投影值；根据预先创建的査找表或者通过解析求解方程组的方法，计 算双能投影值对应的基材料系数投影值；基于基材料的投影值，得到 基材料系数的分布图像。
优选地，所述高能双能射线的能量大于1MeV。
优选地，所述查找表是通过以下方法创建的选定两种基材料； 计算双能射线通过这两种材料不同厚度组合下的投影值，按照高低能 投影值和不同厚度组合的关系，得到查找表。
优选地，解析方程组的方法包括利用实际得到的高能双能投影 值，通过求解基材料分解下的高低能投影方程组，得到相应的厚度组 合。
优选地，对于物体包含混合物或化合物的情况，所述原子序数和特征密度是等效物理量。
优选地，所述的图像重建方法还包括步骤基于基材料系数分布 图像计算被检物体的原子序数图像。
优选地，所述的图像重建方法还包括步骤基于基材料系数分布 图像计算被检物体的特征密度图像。
优选地，所述的图像重建方法还包括步骤基于基材料系数分布
图像计算被检物体的衰减系数图像。
根据上述方法，解决了高能双能情况下重建的问题，提供了更 为有效的物质识别和违禁品检査手段，大大提高了安全检査的精度和 效率。
从下面结合附图的详细描述中，本发明的上述特征和优点将更明 显，其中
图l可应用本发明的CT扫描系统结构示意图，图中系统采用扇 束圆轨道扫描方式。
图2本发明重建方法的计算流程图。其中连接线上的数据表示从 上一步得到的结果，其作为下一步流程的输入数据。
3A和3B示出了模拟重建得到的原子序数图和特征密度图，采用 单一的石墨模型；以及
图3C和3D示出了重建的原子序数图和特征密度图与理论值的比较。
具体实施例方式
下面，参考附图详细说明本发明的优选实施方式。在附图中，虽 然示于不同的附图中，但相同的附图标记用于表示相同的或相似的组 件。为了清楚和简明，包含在这里的已知的功能和结构的详细描述将 被省略，以防止它们使本发明的主题不清楚。
* CT数学原理将二维分布W(X,力沿着某个方向0求线积分，便得到一维的函数
Pe(O，该函数称为"(x，力在《角度的投影。如果能够得到各个方向的投 影内(O，那么可以根据Radon变换精确计算得到二维分布"(;c,力。从 投影得到二维分布的过程称为重建。这便是CT的数学原理。
实际应用中，X光机和探测器围绕物体旋转一圈，便测量得到物 体的某个切片的线衰减系数分布在各个方向的投影，从而可以根据CT 原理重建得到物体切片的衰减系数二维分布。
基材料分解模型
物质对X射线的线衰减系数，可以通过三种基本效应的线衰减
系数的和来表示-
=^ + A + A =《(五)+ "2/ v (￡) + A/e (￡) 其中^表示物质对X射线的线性衰减系数，//p，A，A分别表示光电、
康普顿、电子对效应所对应的线性衰减系数，而每一项又可以近似的 表示成为两项的乘积，其中系数a与物质的原子序数和密度相关，而
/(^)则与X射线的能量相关。
由物质衰减系数的理论分析和实验结果可以得到，在高能情况 下，光电效应的贡献 相比其他两项非常小，在计算时是可以忽略的。
由此得到-
A = A + A = (五)+ "3,e (￡)
由于每种材料的线衰减系数都可以被两个系数al和a2唯一确 定，因此可以选取两种基材料，比如碳和铝，用基材料的线衰减系数 的线性组合表示其它所有材料的线性衰减系数，也就是
其中^为任意一种材料的线衰减系数，A，仏为基材料的线性衰减系 数，而&，62是基材料系数。这也就是基材料分解方法的核心表达式。 按照下式定义a2，A:
。z ， z2
<32 = 2p — ，a3 = — 其中p为物质的密度，Z为物质的原子序数，^为物质的原子质量数。
6则可以得到通过基材料分解方法计算有效原子序数和特征密度的公 式<formula>formula see original document page 7</formula>(2)
*基材料投影模型
从x光管或者加速器中产生的x射线都是连续谱的，在x射线 投影过程中，可以将射线源能谱和探测器能谱合并为d(e)，便于计算， 其满足归一化条件
fm举)必=1
那么由此，对于投影方程，也可以改写成连续谱的形式<formula>formula see original document page 7</formula>
将基材料分解模型带入上式，则双能的投影也可以表示成
<formula>formula see original document page 7</formula>
其中/^,P2为高能和低能两个能量下的投影，^，z)2为高能和低能下的 x射线系统能谱，而a，&为基材料厚度，定义为
利用査找表的方式，或者直接求解方程组(3)的方法，得到基材料投影 系数&，采用通用的滤波反投影重建算法得到卜<formula>formula see original document page 7</formula>M2(,)=化(病W
0 —w ;r +00
其中(P,。表示为投影值的径向和角向坐标。
利用公式(2)得到原子序数Zeff和特征密度A，利用公式(l)可以得 到任意能量下的线衰减系数图像。
图1是根据本发明实施例的双能CT系统的示意图。如图1所示， 射线源100在控制器的控制下按照预定的时序产生能量连续分布的双
能X射线。被检物体200放在承载机构300上。承载机构300在控制 器500的控制下可以匀速旋转，也可以上升和下降。探测器阵列400 设置在与射线源IOO相对的位置，在控制器的控制下接收透射过被检 物体200的透射射线，得到针对第一能量的探测信号和针对第二能量 的探测信号。由探测器阵列400检测的信号被转换成数字信号，存储 在计算机中，用于后续的重建操作。
如图l所示，根据
具体实施例方式
(1) X射线源IOO采用高能双能加速器射线源，其可以快速交替产生 两种高压下的X射线。由于大型货物的体积较大，使用加速器 射线源以产生高功率的射线，以获得更清晰重建图像。
(2) 承载机构例如是转动平稳的载物平台。
(3) 线阵探测器400阵列垂直于X射线源100和承载机构300中心
的轴线，水平放置。
(4) 整个CT系统的机电控制、数据传输、图像重建都由计算机工作 站完成。工作站进行断层图像的重建，在显示器上以二维或三维 的方式显示。
(5) 为了能够精确的重建图像，CT系统精确地测量或标定以下系统 参数X射线源点到探测器的距离D， X射线源到转台旋转轴的 距离仏X射线源在探测器上映射位置c，探测器的像素尺寸A
8以及所有探测器单元的准确几何位置Xi，转台旋转角度S。
(6) 系统的扫描方式，采用标准的扇束圆轨道扫描方式，即射线源和 探测器固定在某一高度，物体随载物台旋转一周，得到双能的
CT投影数据。
(7) 图像重建过程采用上述的方法，从双能的CT投影数据中，通过 计算机实现，得到扫描断层物质的原子序数和特征密度图像，从 而可以为后续的物质识别和判定提供依据。
下面结合附图2描述本发明实施例的图像重建方法的详细过程。 在步骤Sll，利用双能射线扫描物体，以获得双能投影值； 在步骤S12,根据预先创建的查找表或者通过解析求解方程组的 方法，计算双能投影值对应的基材料系数投影值。创建査找表的方法 是，选定两种基材料，计算双能射线通过这两种材料不同厚度组合下 的投影值，按照高低能投影值和不同厚度组合的关系，得到查找表。 而解析求解方法组的方法，则利用实际得到的高低能投影值，通过求 解基材料分解下的高低能投影方程组，得到相应的厚度组合，这种方 法精度高，但是计算速度慢，在实际应用中优先使用査找表的方法。 在步骤S13由基材料的投影值，就可以得到基材料系数的分布图像。
在步骤S14，从基材料的系数分布图得到被检物体的原子序数、 特征密度图像，和任意能量下被检物体衰减系数图像。这一步骤在硬 件系统得以实施。本发明建立了高能X射线下的基材料分解双能CT 重建方法，可以较为准确的从双能投影中，得到断层的原子序数和特 征密度图像。对于物质中包含化合物和混合物的情况，得到的原子序
数和特征密度都是等效物理量。
图3A到3D示出了部分实验结果，它是利用单一物质石墨模型进 行的模拟实验结果，其中图3C和3D分别表示原子序数图像和特征密 度图像的水平方向的剖线图，虚线表示重建值，实线表示理论值。
本发明作为一种通用的高能双能图像重建方法，在一般的以高 能X射线为射线源的双能CT系统中，都有其适用性，具有广泛的应
用前景。上面的描述仅用于实现本发明的实施方式，本领域的技术人员应 该理解，在不脱离本发明的范围的任何修改或局部替换，均应该属于 本发明的权利要求来限定的范围，因此，本发明的保护范围应该以权 利要求书的保护范围为准。
权利要求
1、一种在高能双能CT系统中的图像重建方法，包括步骤利用高能双能射线扫描物体，以获得高能双能投影值；根据预先创建的查找表或者通过解析求解方程组的方法，计算双能投影值对应的基材料系数投影值；基于基材料的投影值，得到基材料系数的分布图像。
2、 如权利要求1所述的图像重建方法，其中所述高能双能射线的 能量大于1MeV。
3、 如权利要求1所述的图像重建方法，其中所述查找表是通过以 下方法创建的 选定两种基材料；计算双能射线通过这两种材料不同厚度组合下的投影值，按照高 低能投影值和不同厚度组合的关系，得到查找表。
4、 如权利要求1所述的图像重建方法，其中求解方程组的方法包括利用实际得到的高能双能投影值，通过求解基材料分解下的高低 能投影方程组，得到相应的厚度组合。
5、 如权利要求1所述的图像重建方法，其中对于物体包含混合物或化合物的情况，所述原子序数和特征密度是等效物理量。
6、 如权利要求l所述的图像重建方法，还包括步骤 基于基材料系数分布图像计算被检物体的原子序数图像。
7、 如权利要求l所述的图像重建方法，还包括步骤 基于基材料系数分布图像计算被检物体的特征密度图像。
8、 如权利要求l所述的图像重建方法，还包括步骤-基于基材料系数分布图像计算被检物体的线衰减系数图像。
全文摘要
公开了一种在高能双能CT系统中的图像重建方法，包括步骤利用高能双能射线扫描物体，以获得高能双能投影值；根据预先创建的查找表或者通过解析求解方程组的方法，计算双能投影值对应的基材料系数投影值；基于基材料的投影值，得到基材料系数的分布图像。解决了高能下双能重建的问题，提供了更为有效的物质识别和违禁品检查手段，大大提高了安全检查的精度和效率。
文档编号A61B6/03GK101647706SQ20081011830
公开日2010年2月17日 申请日期2008年8月13日 优先权日2008年8月13日
发明者刘以农, 康克军, 丽 张, 李元景, 段新辉, 胡海峰, 邢宇翔, 陈志强, 黄清萍 申请人:清华大学;同方威视技术股份有限公司