基于加权平均的工业CT旋转中心校正方法.docxVIP

基于加权平均的工业CT旋转中心校正方法

摘要：旋转中心的精确定位是重建高质量CT图像的一个重要条件。本文提出一种基于加权平均的旋转中心校正方法，该方法首先将投影正弦图以180°为界限分成两组投影数据，计算出这两组投影数据的互相关系数及其对应的探测器探元坐标，然后使用最大与次最大的互相关系数对探测器探元坐标值进行加权平均得到实际旋转中心坐标。本文所提算法完全基于原始投影数据，无需确定CT测量的几何参数，不用提前对投影数据进行插值运算就可以将旋转中心精确到亚像素级，精度更高，实现简单。

关键词：CT图像；旋转中心校正；加权平均；正弦图；互相关

工业CT（industrialcomputerizedtomography，ICT）是工业计算机断层成像技术的简称，它能在对检测物体无损伤条件下，以二维断层图像或三维立体图像的形式，清晰、准确、直观地展示被检测物体的内部结构、组成、材质及缺损状况，被誉为当今最佳无损检测和无损评估技术[1-2]。获得清晰的重建图像是进行无损检测的前提，重建高质量图像需要获得准确的CT投影几何参数，旋转中心（centerofrotation，COR）是其中一项十分重要的参数[3-4]。重建过程中，COR的偏移会使CT重建图像产生伪影、图像精度变差，导致CT图像降质[5]，影响无损检测结果。COR能否准确定位对于CT图像成像质量好坏有举足轻重的影响。

现有的旋转中心的测量校正方法主要有正弦图中心法、相对角法、几何法、正弦图互相关法和点扩展函数法等。刘明进[6]根据系统旋转中心偏移量与点扩展函数极大值在一定范围内成高斯分布的关系，拟合点扩展函数的极值计算旋转中心，提出了点扩展函数法，该方法受工件形状影响较大。Liu等[7]提出了3种校正旋转中心的方法，分别为正弦图中心法、相对角法和几何法；正弦图中心法通过寻找正弦图的左右边界，求其平均值得到旋转中心，正弦图中心法要求射线源与旋转中心的连线垂直于探测器，但是在实际操作过程中，射线源与旋转中心的连线总会存在些许偏差，且该方法受噪声及边界像素定位精度影响较大；相对角法采用线模扫描，实现过程繁琐且精度受噪声影响较大；几何法虽然解决了射线源与旋转中心的连线必须垂直于探测器的情况，但是焦距的精确值和中心的精确位

置很难测得。Yang等[8]提出的基于图像互相关方法，需要预先将扇束投影转化为平行束投影，该方法计算效率偏低耗时较长。李保磊等[9]提出基于正弦图冗余信息的COR确定方法，有其不足之处，忽略了散射的影响。孟凡勇等[10]利用相差180°的投影会重复一次的性质，计算两组重复数据的相关程度获得旋转中心，提出正弦图互相关法（sinogramcross-correlation，SCC）结果较为准确，但是此方法计算的旋转中心的精度只能达到像素级。刘欢等[11]提出基于统计平均方法（statisticalaveragingmethod，SAM）的旋转中心校正，该方法先将相差180°的投影图像相减得到差值图像，然后对差值图像取绝对值作为样本，最后寻找样本中最小值的对应坐标作为旋转中心的坐标；该方法无需插值运算即可将旋转中心坐标精确到亚像素级，但是当探测器同时采集到载物台和被测物体的投影时，使用该方法来计算旋转中心会产生较大误差，导致重建的图像出现伪影，所以在实际工业CT成像系统中，使用该方法对不同工件进行旋转中心校正，需人为的对原始投影数据进行截取以获得精确的旋转中心。

由于实际CT系统与理想系统之间的差异使得旋转中心出现偏差，重建得到的CT图像不可避免地出现伪影[12]，严重降低了CT图像质量，对后续目标的识别和分析带来干扰。针对旋转中心偏移问题，本文提出一种加权平均校正方法（weightedaveragecorrection，WAC）。首先将采集到的投影正弦图以180°为界限分为两组数据；然后计算两组数据的互相关系数，查找并记录互相关系数的最大值和次最大值以及对应的探测器坐标；最后利用加权平均法对记录的探测器坐标计算得到实际旋转中心坐标。

经验证，本文所提算法不受工件尺寸大小的影响，而且无需提前对原始投影数据进行插值便可得到精确的旋转中心。通过本文算法进行旋转中心校正，校正精度可达亚像素级，同时所需的计算量小、实现简单。

1加权平均校正方法

1.1算法原理

锥束CT对工件进行全周扫描获取投影数据，旋转角度为阝和β+180。的投影数据具有几何对称关系，这两个角度的投影被称为对称投影[13]。不难发现，从射线源发射并经过旋

转中心的射线在相差180°的两个投影角度下经过的路径是相同