基于扇束拼接的工业CT重建方法研究.pdfVIP

中国体视学学会CT理论与应用分会、中国工程物理研究院 2008年7月 应用电子学研究所、《cT理论与应用研究》编辑部论文集 基于扇束拼接的工业CT重建方法研究 张立凯，魏东波，王倩妮，王世洋 (北京航空航天大学机械工程及自动化学院，北京100083) 摘要：针对工业CT(ComputerizedTomography)大尺寸构件的快速重建问题，提出了一种利 用2次偏置扫描扩大扫描视场的方法，并推导了扇束拼接、投影数据预处理和滤波反投影重建 算法，实现了将探测器在2个位置获得的投影数据首先重排成一组等价的扇束投影数据，然后 将该投影数据重排成一组完整的平行束投影数据，从而使滤波反投影重建得以进行。计算机仿 真结果表明，该方法在重构大尺寸构件时，获得了满意的质量。分析表明，其有效扫描视场达 到标准三代工业cT的4．2倍。该方法突破了扇束角大小的限制，且与现有三代工业CT兼容， 工程实现容易。 关键词：工业cT：滤波反投影；扇束；大视场 工业CT是以x射线扫描、探测器采集的数字投影序列重建扫描区域内被检试件横截面 射线衰减系数分布的映射图像。据此图像，可对被检试件的结构组成、材质分布、密度变 化、试件有无缺陷以及缺陷大小、性质和位置等内部信息进行判读和计量，是目前国际无 损检测界常用无损检测技术之一。 滤波反投影算法速度快，是工业CT中常用的算法。但受探测器和扇束角大小的限制， 二维CT面临的一个突出问题是大尺寸构件的快速重建。目前，检测大尺寸构件的方式主要 有：二代扫描方式，可以实现对大尺寸构件的精确重建，但费时效率低；基于三代扫描方 式的有投影拼接法u’51、伪数据拟合法H1和冗余数据挖掘法¨1。投影拼接法通过扇束拼接可 以实现较大的扫描视场，效率较高，但投影拼接部分需要进行插值。伪数据拟合法利用对 扫描构件的先验知识，通过对缺失投影进行估计和填充实现对构件的重建。冗余数据挖掘 法利用三代扫描方式存在的投影数据冗余特征，通过改进扫描方式扩大扫描视场，但扫描 视场相对依然较小。难以满足大尺寸构件的层析需求。 为此，本文提出了一种通过转台2次偏置扫描，充分利用投影拼接法和冗余数据挖掘法 各自的优点来扩大扫描视场的二维ICT成像方法，根据投影的几何特征，推导了扇束拼接 的函数映射关系，基于三代ICT投影数据存在的冗余特征，推导了投影数据重排预处理方 法和滤波反投影重建算法，并通过计算机仿真验证了其正确性。分析表明，该方法的有效 扫描视场达到标准三代ICT方法的4．2倍左右，扫描效率较二代扫描方式大为提高。 1 扫描原理 大尺寸构件扇束扫描结构如图1所示。扫描时，转台首先在位置0做一次3600回转扫 描，然后在位置1再做一次回转扫描，即可完成扫描层析。可检构件回转直径为R，而三 代cT可检构件回转直径最大为足，基于三代CT的偏置扫描的冗余数据挖掘法可检构件回 全国射线数字成像与CT新技术研讨会 转直径最大为是。图l中，射线源到探测器的距离为D，探测器的大小为Ⅳo 图1大尺寸构件扇束cT扫描几何结构 2重建算法 2．1重建算法思路 本文提出的扫描方式的重建算法思路是：首先将探测器在2个位置获得的扇束投影拼 接成一幅等价的扇束投影，然后根据扇束投影数据的冗余特征，将其重排成平行束投影， 最后运用平行束滤波反投影重建公式进行重建。平行束几何结构如图2所示，重建公式如 下： ‘ 口(，，≯)2 Oj办(7)·尸(f，口)l=””伊一”d口 (1)… ‘ 2Jg(f，≯)I，l，。os(州)dO 0 影重建位置，口为投影角度。 射线源 探潺蕃 图Z二维平行束