基于投影匹配的X射线双能计算机层析成像投影分解算法_李保磊
基于投影匹配的X射线双能计算机层析成像投影分解算法李保磊 层析成像
基于投影匹配的X射线双能计算机层析成像投影分解算法李保磊层析成像第31卷第3期xx年3月光学学报V ol . 31, N o . 3March , xx基于投影匹配的X 射线双能计算机层析成像投影分解算法李保磊1, 2张耀军11(公安部第一研究所, 北京100048; 2北京大学数学科学学院, 北京100871)摘要 X 射线双能计算机层析成像(CT ) 技术是安全检查领域一种重要的材料探测与识别手段。
双能C T 投影分解是双能CT 预处理重建算法的核心内容和关键步骤。
针对现有投影分解算法的不足, 提出了一种基于投影匹配的双能CT 投影分解算法。
依据系统能谱和基材料线性衰减系数曲线, 通过求解投影积分方程组建立高低能投影查找表。
对于给定的高低能投影, 在查找表中寻找最佳匹配点, 进而获取基材料分解投影。
该算法避免了现有算法复杂的迭代优化过程, 简化了系统的标定过程, 分解精度取决于查找表的设定步长。
相对现有算法该算法有实现过程简单、易于并行计算的优点。
仿真实验结果验证了算法的有效性。
关键词 X 射线光学; 双能计算机层析成像; 投影匹配; 投影分解; 预处理重建 T P391 A do i :10. 3788/AOS xx31. 0311002Projection De position Algorithm of X -Ray Dual -Ene rgy Compute d Tomography Base d on Projection MatchingLi Baolei 1, 2 Zhang Yaojun 11The First R esearch Instit ut e of Minist ry of Public Security , Bei __g 100048, China2School of Mathematical Sciences , Peking University , Bei __g 100871, ChinaAbstract X -ray dual -energy puted tomography imaging technique is an important material detec tion and recognition method in the field of security inspection . The projection deposition is the nuclear content and key technique in the pre -rec onstruction algorithm of dual -energy puted tom ography . Ac cording to the disadvantages of the c urrent algorithms , a projection deposition algorithm based on projection matching is proposed . Firstly , based on energy spectrum of the system and the linear attenuation coefficient curve of the ba sic materials , the high -and low -energy projection lookup table can be got by solving the projection integral equations set . For a given dual -energy projection , find the best match point in the lookup table and then obtain the deposition projection of basic m ateria ls . The proposed algorithm avoids the processof plex iteration and optimization and simplifies the proc ess of system calipation . The deposition ac curacy depends on the step setting of lookup table . Compared to the c urrent algorithms , the proposed algorithm ′s realization is more simple and easy for parallel putation . The feasibility of the a lgorithm is va lidated by the results of sim ulation experiment . Key wo rds X -ray optic s ; dual -energy puted tomography ; projection matching ; projection deposition ; pre -reconstructionOCIS co des 110. 0110; 340. 0340; 020. 0020; 250. 02501 引言X 射线双能计算机断层成像(C T ) 作为一种重要的成像技术, 在安全检查领域起着越来越重要的作用[1～3]包括单能X 射线透视、双能X 射线透视、多视角、背散射和X 射线CT 技术等) , 双能C T 技术具有较高的探测能力, 它能够较准确的获取材料的电子密度、有效原子序数二维信息, 而利用此二维信息能够实。
层析成像的重建方法[发明专利]
![层析成像的重建方法[发明专利]](https://img.taocdn.com/s3/m/b9dcdaea03d8ce2f01662338.png)
专利名称：层析成像的重建方法专利类型：发明专利
发明人：孔明明,李冬东,赵辉
申请号：CN201410211463.2申请日：20140519
公开号：CN105096357A
公开日：
20151125
专利内容由知识产权出版社提供
摘要：本发明是关于层析成像的重建方法。
根据本发明的一实施例，一层析成像的重建方法包含：获取投影数据；图形处理单元使用开放计算语言对该投影数据进行滤波；及该图形处理单元使用该开放计算语言对所滤波的投影数据进行反投影，对经该反投影而重建的容积图像进行归一化处理。
本发明实施例所提供的层析成像的重建方法通过协调CPU与GPU之间的任务分担，巧妙利用开放计算语言的结构特点而减少数据传输占用的额外时间，可大幅度提高重建的速度。
另外，本发明实施例还通过加权及归一化等处理提高成像的质量。
申请人：锐珂(上海)医疗器材有限公司
地址：上海市浦东新区浦东新金桥路27号
国籍：CN
代理机构：北京律盟知识产权代理有限责任公司
代理人：章蕾
更多信息请下载全文后查看。
星载反射式成像系统计算光学设计方法综述
第44卷第6期航天返回与遥感2023年12月SPACECRAFT RECOVERY & REMOTE SENSING1星载反射式成像系统计算光学设计方法综述高金铭1,2,3郭劲英3孙越4卞殷旭1,2于钱5彭吉龙5贺融6罗茂7匡翠方1,2司徒国海3,*（1 浙江大学极端光学技术与仪器全国重点实验室，杭州310027）（2 浙江大学杭州国际科创中心，杭州311215）（3 中国科学院上海光学精密机械研究所，上海201800）（4 中国科学院南京天文光学技术研究所，南京210042）（5 北京环境卫星工程研究所，北京100094）（6 重庆嘉陵华光光电科技有限公司，重庆400700）（7 重庆嘉陵特种装备有限公司，重庆400032）摘要星载反射式遥感成像系统的结构集成度高、光学面型复杂，面临光学初始结构难以求解和像质优化难以收敛的设计难题。
综述着重探讨星载反射式成像系统中的计算光学设计方法。
星载反射式成像系统计算光学设计，包括复杂曲面反射镜、曲面反射型光栅等硬件设计，以及计算成像图像解码等算法设计。
面对星载遥感成像的不同应用，文章从深度学习光学结构设计、合成孔径计算成像、景深延拓计算成像、主动光学计算像差补偿和曲面光栅计算光谱成像等几个方面进行系统性地分类讨论。
本综述结论为：基于计算光学的设计方法，无论在光学结构与面型求解方面，还是图像信号非线性逆问题求解方面，都具有强大的设计能力。
在发展趋势上，计算光学在航天光学系统的设计潜力刚刚被挖掘。
在“人工智能时代”的软、硬件算力支持下，计算光学方法将大大提升星载反射式成像系统设计的便捷性。
关键词遥感成像光学设计人工智能计算光学中图分类号: O439文献标志码: A 文章编号: 1009-8518(2023)06-0001-11DOI: 10.3969/j.issn.1009-8518.2023.06.001Overview of Computational Optics Design for Space-Borne ReflectiveImaging SystemGAO Jinming1,2,3GUO Jinying 3SUN Yue 4BIAN Yinxu 1,2YU Qian 5PENG Jilong 5HE Rong 6LUO Mao 7KUANG Cuifang 1,2SITU Guohai3,*（1 State Key Laboratory of Extreme Photonics and Instrumentation, College of Optical Science and Engineering, ZhejiangUniversity, Hangzhou 310027, China）（2 ZJU-Hangzhou Global Scientific and Technological Innovation Center, Zhejiang University, Hangzhou 311215, China）（3 Shanghai Institute of Optics and Fine Mechanics, Chinese Academy of Sciences, Shanghai, 201800, China）（4 Nanjing Institute of Astronomical Optics & Technology, Chinese Academy of Sciences, Nanjing 210042, China）（5 Beijing Environmental Satellite Engineering Institute, Beijing 100094, China）（6 Chongqing Jialing Huaguang Optoelectronic Technology Co. Ltd., Chongqing 400700, China）（7 Chongqing Jialing Special Equipment Co. Ltd., Chongqing 400032, China）Abstract The space-borne reflective remote sensing imaging system are facing challenges to solve the收稿日期：2023-07-26基金项目：国家自然科学基金项目（61991452，62061136005，12104472，62005120）引用格式：高金铭, 郭劲英, 孙越, 等. 星载反射式成像系统计算光学设计方法综述[J]. 航天返回与遥感, 2023, 44(6): 1-11.GAO Jinming, GUO Jinying, SUN Yue, et al. Overview of Computational Optics Design for Space-Borne ReflectiveImaging System[J]. Spacecraft Recovery & Remote Sensing, 2023, 44(6): 1-11. (in Chinese)2航天返回与遥感2023年第44卷optical initial structure and optimize the image quality, which are with high structural integration and complex optical facets. This review focuses on computational optical design methods in space-borne reflective imaging systems. These include hardware design of complex curved mirrors and curved reflective gratings, as well as algorithm design of computational signal decoding. This review systematically discusses several aspects such as deep learning optical structure design, synthetic aperture computational imaging, depth-of-field extension computational imaging, active optics computational aberration compensation, and curved grating computational spectral imaging. In conclusion, the design methods based on computational optics have powerful design capabilities, both solving optical structure and optimizing image signal nonlinear inverse problem. As trend, the potential of computational optics in aerospace optical systems has just been explored. With the support of computing power in the "artificial intelligence era", the computational optics method will greatly improve the convenience of designing satellite-based reflective imaging systems.Keywords remote imaging; optical design; artificial intelligence; computational optics0 引言目前，大口径、宽波段、长焦距、轻量化以及极短波长探测等是空间遥感光学系统的发展趋势。
一种基于信息融合的图像分割方法
Abstract : Image segmentation is one of the basic steps in machine vision and the base of image analysis. Many studies have been done in this field in China and abroad. However , due to the influences of light source , tar2 get shape , image theory and so on , especially the former factor , the segmenting method has a poor compatibili2 ty. Comparing with the limitation of the previous image segmenting methods , this paper , taking hue and tex2 ture as well as the method of mathematical morphology into account , proposes a segmenting method based on information fusion. Experiments show that in the condition of no obvious change of light source , the segmenting method has a better compatibility with different homogeneous objects. Key words :information fusion ; hue ;texture ; image segmentation
结合全变差最小化的双能CT重建
ｃ ｎ ｔ ｃ ｉｎ ｉ ａ ａ ｌ ｆｇ ｔ ｎ ｅｔ ｒ ｕ ｌ ｙ ｉ ｇ ｓｏ ｔｒａｓ ｆ ｃｉ ｅａ ｏ ｃ ｎ ｍｂ ｒａ ｄ ｅｅ ｔ ｎ ｄ ｎ ｉ ，ａ ｄ ｉ ｒ ｏ ｓ ｔ ｓｃ ｐ ｂ ｅ ｏ ｅｔ ｇ ｂ ｔ ａｉ ｕ ｒ ｏ ｉ ｅｑ ｔ ｍａ ｅ ｆ ｍａ ｅ ｉｌ ’ｅｆ ｔ ｔｍｉ ｕ ｅ ｎ ｌｃｒ ｅ ｓ ｙ ｎ ｓｍｏ ｅ ｅ ｖ ｏ ｔ
展， 犯罪 团伙使用的犯罪技术 越来越 先进 ， 爆炸物 的种类越来
越 多样 化 ， 特别是 对于 那些 密度 相近但 是原 子序 数不 同 的材
料， 传统 的单能 Ｃ Ｔ就有些 力不 从心 了［ 。１７ ４ ９６年 Ａｖｒｚ等 ｊ ｌａｅ
人 提 出了双能 Ｃ Ｔ重建 的概念 ， 由于其能 检测 出物体 的有效 原子序数和电子密度 ， 以区别密度相近但原子序数不 同的物 可
ｃ ｎ ｔ ｃｉ ｎ ａ ｇ ｒ ｈ ，ａ ｄ ｔ ｅ ｏ ｌ ｈ ｆ ｃ ｉｅ ａｏ ｃ ｎ ｍｂ ｒａ ｄ ｅ ｅ ｔｏ ｅ ｓｔ ｍａ ｅ ． ｉ ａｌ ａ ｇ ｒ ｈ ｗｈｃ ｏ ｓｒ ｔ ｌｏ ｉ ｍ ｕ ｏ ｔ ｎ ｈ ｎ ｃ ｕｄ ｇ ｔｔ ｅ ｅｆ ｔ ｔ ｍｉ ｕ ｅ ｎ ｌｃ ｒｎ ｄ ｎ ｉ ｉ ｇ ｓ Ｆ ｎ ｌ ｅ ｖ ｙ ｙ， ｌｏ ｉ ｍ ｉｈ ｔ ｗａ ａ ｅ ｎ Ｔ ｎ ｍｉａｉｎ ｗ ｓ ｕ ｅ ｏ ｐ ｏ ｅ ｓｔ ｅ ａ ｏ ｅ ｄ ｔ ． ｈ ｘ ｅ ｍｅ ｔｌｒ ｓ ｌ ｈ ｗ ｔ ａ ｕ ｌ ｎ ｒ Ｔ ｃ ｎ ｏ — ｓｂ ｓ ｄ ｏ Ｖ ｍｉ ｉ ｚ ｔ ａ ｓ ｄ ｔ ｒ ｃ ｓ ｈ ｂ ｖ ａａ Ｔ ｅ ｅ ｐ ｒ ｎ ａ ｅ ｕ ｔ ｓ ｏ ｈ ｔ ａ ｅｇ Ｃ ａ ｂ ｏ ｉ ｓ ｄ ｅ ｙ
基于扇束双投影方向的改进迭代层析成像算法
基于扇束双投影方向的改进迭代层析成像算法
方静;程乐红;张玉萍;罗高鹏
【期刊名称】《光子学报》
【年(卷),期】2014(0)10
【摘要】提出一种适用于气体扩散分布重建的改进代数迭代算法.针对投影射线扇束构型及双投影方向,在算法迭代步骤中引入先验矩阵,对迭代结果间隔性地进行中值滤波.以高斯烟羽模型为测试对象,分析迭代次数、松弛因子及不同滤波器使用策略对重建结果的影响,结果表明:迭代次数为500次,松弛因子为2时,均方误差、峰值相对误差基本收敛,并将波峰控制在一定范围内.改进代数迭代算法的均方误差、峰值相对误差重建指标优于传统代数迭代重建算法;同时,改进后算法能够准确定位波峰,完成尾峰重建,具有较强的抗噪性.
【总页数】8页(P13-20)
【关键词】层析成像;迭代算法;中值滤波器;先验知识;遥感;双投影方向;扇束构型【作者】方静;程乐红;张玉萍;罗高鹏
【作者单位】合肥工业大学计算机与信息学院
【正文语种】中文
【中图分类】O439
【相关文献】
1.基于半扇束扫描的CT图像迭代重建算法 [J], 许娟;许乙凯
2.基于视觉灵敏度及改进方向距离彩色图像迭代滤波算法 [J], 李高西;曹军;张福元;
李华
3.有序表改进扇束滤波反投影重建算法求解速度的研究 [J], 王培;郑兆瑞;郝记秀
4.基于改进迭代反投影算法的遥感影像超分辨率重建 [J], 郭桐宇
5.基于改进梯度投影法的MIMO-NOMA系统迭代信号检测算法 [J], 李刚;胡哲;景小荣;周子钰
因版权原因，仅展示原文概要，查看原文内容请购买。
基于对称投影的2D-CT系统投影旋转中心自动确定方法
（８）
其中，Ｓ∈［Ｓａ，ｓｂ】代表物体投影区域，已，％代表正弦图中物体投影区域的最大左右边界。 值得注意，算法步骤中的投影Ｐ（ｆｌ，ｓ）的物理含义是物质衰减系数的线积分，需要对原 始信号取负对数运算，鉴于探测器探测到的射线场强信号与投影ｐ（ｆｌ，ｓ）一一对应关系，算 法也可以用射线场强信号，代替投影信号Ｐ（ｆｌ，ｓ），以免去求负对数的运算。
１．２算法特性分析和改进 经分析可知，本算法非但适用于图１） 示的情况，即射线源和旋转中心的连线垂 直于探测器，同样适用于射线源和旋转中 心的连线不严格垂直于探测器的情形（见
图２）。
在成像坐标系精确校准哺－１¨了以后， 射线源中心射束（垂直于探测器的射束） 坐标与旋转中心投影在探测器上的坐标应 该是重合的，然而在实际的２Ｄ—ｃＴ成像系 统中，很难保证上述两对坐标完全重合。 当两者不重合时等效于探测器有一定的偏 摆角，而实验表明，存在很小的偏角时 （±ｏ．２。）悔１对成像像质影响不大。 在图２所示的几何扫描布局中，仍有：
李保磊等：基于对称投影的２１卜ＣＴ系统投影旋转中心自动确定方法
３）该方案无须使用模体，亦无须知晓任何几何参数，只需要采集到的原始投影数据。 ４）该方案具有实时性，实现简便，且精度较高。 值得指出的是，该算法也有不足之处，在某些特殊情况下，通过ｍｅａｎ（ｓ）灰度曲线寻找 极值点时，极值点不唯一。如当被扫描物体是圆柱体且中心和转台旋转中心重合时，则 ｍｅａｎ（ｓ）灰度曲线理论上均为０值，故找不到极值点。实际上当圆柱体中心位于转台旋转中 心附近或者图像噪声过大时，均易产生极值不唯一的情况，导致算法定位误差增大。解决 此问题的方法是在标定系统时尽量不采用圆对称物体进行扫描或者将物体尽量不放在旋转 中心位置，通过这样的操作即可克服算法的不足之处。
超声CT层析成像方法
目 考虑空间所有 离散 点的走时和射线路径 的全局射 线 前，
追 踪方 法引起了研 究者 的广泛关注 。 这种方 法的理论基础 是 程 函方程、 ｅ ｍ ｔ Ｆ ｒ ａ原理和Ｈ ｙ ｅ 原理 ， ｕｇｎ 实际上各种全局算法 的
２波动层析成像
２ 世纪７ 年代末， Ｏ ０ 波动层析成像的研 究由于以下两个 因素
射线层析算法仅考虑声波 的折射问题 ， 在均匀介质或衍射 情况较弱时比较 有效。 当介质 内部粒子与入射波产生明显 的衍
信 息通 信 射 现 象 ， 该 方 法 确 定波 的传 播 路 径 并成 像 将 会 失效 。 时须 用 此
２ １年 ０１
第３ 期
总第 １３ 期 １
先验信息的情况下， 均能较好地收敛于问题 的真实解， 但是这 两种算法都要用到矩阵的奇异值分解技术， 所需的计算量和存 储量均较 大。 除了上述正则化方法外， 研究者们还根据不 同的 实际问题提 出了许 多种不同的正则化方法， 但这些方法 的正则
在误差 时， 可能会放 大噪声， 因此其数值稳定性不如Ｂ方法。 Ｉ 与 这类算法发展 的同时， 另一类算法则采用直接求解 非线性方程 的方法 ， Ｐｒ ｎ ｈ ｉ在１９年提 出的ＬＭ ￣Ｆ ａ ｃ ｏ ｓ ９ ７ — 算法和Ｊ ａ ｈ ｍ ｗ ｃ ｏ ｃ ｉ ｏ ｉｚ 在１ ９ 年提 出的Ｎ Ｋ ９８ — 算法 ， 它们 均是求解 非线性 问题Ｇ ｕ ｓ ａ ｓ— Ｎ ｗ ｏ Ｊ法 的变形 。 ｅ ｔ ｎ￣ － 近２ 年来， Ｏ 基于空间域的超声Ｃ 数迭代技术迅速 发展 。 玳 这 类迭代技术 困扰人们的主要问题是逆 散射方程 的不适定性 问题 ， 即解 的存在性 、 唯一性和稳定性 的满 足问题 。 于解的 对
方法 来反演料 内部结构， 提出了求解超声逆 散射成像 的 并
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第31卷 第3期光 学 学 报V ol .31,N o .32011年3月ACTA OPTICA SINICAMarch ,2011基于投影匹配的X 射线双能计算机层析成像投影分解算法李保磊1,2 张耀军1(1公安部第一研究所,北京100048;2北京大学数学科学学院,北京100871)摘要 X 射线双能计算机层析成像(CT )技术是安全检查领域一种重要的材料探测与识别手段。
双能C T 投影分解是双能CT 预处理重建算法的核心内容和关键步骤。
针对现有投影分解算法的不足,提出了一种基于投影匹配的双能CT 投影分解算法。
依据系统能谱和基材料线性衰减系数曲线,通过求解投影积分方程组建立高低能投影查找表。
对于给定的高低能投影,在查找表中寻找最佳匹配点,进而获取基材料分解投影。
该算法避免了现有算法复杂的迭代优化过程,简化了系统的标定过程,分解精度取决于查找表的设定步长。
相对现有算法该算法有实现过程简单、易于并行计算的优点。
仿真实验结果验证了算法的有效性。
关键词 X 射线光学;双能计算机层析成像;投影匹配;投影分解;预处理重建中图分类号 T P391 文献标识码 A do i :10.3788/AOS 201131.0311002Projection De composition Algorithm of X -Ray Dual -Ene rgy Compute d Tomography Base d on Projection MatchingLi Baolei 1,2 Zhang Yaojun 11The First R esearch Instit ut e of Minist ry of Public Security ,Beijing 100048,China2School of Mathematical Sciences ,Peking University ,Beijing 100871,ChinaAbstract X -ray dual -energy computed tomography imaging technique is an important material detec tion and recognition method in the field of security inspection .The projection decomposition is the nuclear content and key technique in the pre -rec onstruction algorithm of dual -energy c om puted tom ography .Ac cording to the disadvantages of the c urrent algorithms ,a projection decomposition algorithm based on projection matching is proposed .Firstly ,based on energy spectrum of the system and the linear attenuation coefficient curve of the ba sic materials ,the high -and low -energy projection lookup table can be got by solving the projection integral equations set .For a given dual -energy projection ,find the best match point in the lookup table and then obtain the decomposition projection of basic m ateria ls .The proposed algorithm avoids the process of complex iteration and optimization and simplifies the proc ess of system calibration .The dec omposition ac curacy depends on the step setting of lookup table .Compared to the c urrent algorithms ,the proposed algorithm ′s realization is more simple and easy for parallel computation .The feasibility of the a lgorithm is va lidated by the results of sim ulation experiment .Key wo rds X -ray optic s ;dual -energy computed tomography ;projection matching ;projection decomposition ;pre -reconstructionOCIS co des 110.0110;340.0340;020.0020;250.0250 收稿日期:2010-06-10;收到修改稿日期:2010-11-30基金项目:国家科技支撑计划(2009BA K 64B02)资助课题。
作者简介:李保磊(1980—),男,博士,主要从事X 射线计算机断层成像技术及其应用等方面的研究。
E -mail :bao lei li @yahoo .com .cn1 引 言X 射线双能计算机断层成像(C T )作为一种重要的成像技术,在安全检查领域起着越来越重要的作用[1～3]。
在常见的X 射线安全检查技术中(主要包括单能X 射线透视、双能X 射线透视、多视角、背散射和X 射线CT 技术等),双能C T 技术具有较高的探测能力,它能够较准确的获取材料的电子密度、有效原子序数二维信息,而利用此二维信息能够实光 学 学 报现对绝大部分材料的探测识别[4]。
同时CT型安检设备也是美国FAA唯一认定的EDS型安检设备。
可见双能C T技术在整个X射线安检技术中的地位。
双能CT的概念在20世纪70年代由Robert E. Alvarez等[5]首次提出并进行了研究。
最初该项技术主要用于医学射线图像的骨骼与软组织的分离和病灶诊断,之后其应用逐渐扩展到工业检测和安全检查领域,如利用双能技术对CT图像进行射束硬化校正等。
在安全检查领域,国外Reveal和Analogic公司均已推出自己的双能CT行包检查系统[6,7]。
双能CT安检系统的核心技术是双能CT重建算法。
双能CT重建算法大致可以分为三类。
1)迭代重建算法[8,9]。
这类方法通过迭代的过程,得到与能谱无关或者近似无关的CT图像,进而抑制由射束硬化效应带来的伪影,改善图像质量;或者通过双能迭代,得到材料分解图像。
不过此类方法计算速度慢,不易实际应用。
2)后处理重建算法[10～12]。
该方法首先利用高低能投影重建出高低能CT图像,然后再对高低能CT图像进行处理,获取物体断层的物理参数分布图像。
由于其自身固有的特点,此类方法不能从根本上消除射束硬化效应,因此精度难以保证。
3)基于投影分解的双能CT重建算法[13～16]。
对于这类方法,基于两种常见的物质线性衰减系数分解模型———双效应分解模型和基材料分解模型[13,17],对应着两种投影分解过程。
这类方法通过对投影数据的预处理分解,利用分解的参数完成双能C T重建进而对材料进行识别探测,因此又称为预处理重建算法。
这种算法理论上可以得到不受能谱影响的物质有效原子序数和电子密度图像,是目前双能CT 技术领域的主流重建算法。
在预处理重建算法中,投影分解是关键技术和核心内容,现有文献中已有多种投影分解方法,如多项式间接拟合法、多项式直接拟合法、分段多项式拟合法、等值线法和高次曲面函数拟合法等[16,18]。
这些方法在一定程度上促进了双能CT技术的发展,但是它们求解过程复杂,需要复杂的标定过程,计算速度慢,对噪声十分敏感,分解精度还有待提高。
2006年,Zhengro ng Ying等[19,20]提出了双能CT的优化投影分解方法,提高了分解算法的稳定度和收敛性能,但是仍然存在计算速度与初值的选取问题。
总之,现有投影分解方法还不能满足安检实际工作需求,需要研究简便实用的投影分解算法。
2 双能CT重建原理物质的线性衰减系数有两种常见的分解模型:μ(E)=a c f K N(E)+a p f p(E),(1)μ(E)=b1μ1(E)+b2μ2(E),(2)式中f p(E),f K N(E)为只与能量E有关而与材质无关的分解系数。
且有f p(E)=1E3,(3)f KN(α)=1+αα22(1+α)1+2α-1αln(1+2α)+12αln(1+2α)-1+3α(1+2α)2,(4)式中α=E/510.975keV,a p表示光电效应系数,a c 为康普顿散射效应系数,a p,a c是独立于能量只与材质有关的物理量,且有a p=l1ρZnA, a c=l2ρZA,n≈4～5,(5)式中l1,l2为两常数,ρ为物质密度,Z为原子序数,A 为原子量;μ1(E),μ2(E)分别为两种基材料的线性衰减系数。
b1,b2为对应两种基材料的分解系数,对于某一固定的物质,b1,b2是两个常数。
(1)式表示在一定的射线能量范围内,物质的衰减可由光电效应和康普顿散射两种作用共同组成。
(2)式表示任何一种物质的线性衰减系数都可由两种基材料的线性衰减系数线性叠加而成[13,17]。
依据这两种衰减系数分解模型,记A c=∫a c d l,A p=∫a p d l;B1=∫b1d l,B2=∫b2d l,根据宽能谱射线条件下的Beer定律可得P L=-ln{∫S L(E)ex p[-A c f K N(E)-A p f p(E)]d E}+ln∫S L(E)d EP H=-ln{∫S H(E)exp[-A c f KN(E)-A p f p(E)]d E}+ln∫S H(E)d E,(6) P L=-ln{∫S L(E)exp[-B1μ1(E)-B2μ2(E)]d E}+ln∫S L(E)d EP H=-ln{∫S H(E)ex p[-B1μ1(E)-B2μ2(E)]d E}+ln∫S H(E)d E,(7)李保磊等: 基于投影匹配的X射线双能计算机层析成像投影分解算法式中S L(E),S H(E),P L,P H分别为高低能谱和高低能投影。