[发明专利]基于切片轮廓重投影的锥束CT射束硬化校正方法

说明书

技术领域

本发明涉及一种锥束CT射束硬化校正方法，属于CT技术领域。

背景技术

锥束CT(Cone-Beam Computed Tomography，CBCT)利用锥形束射线源和面阵探测器采集被测物体的投影数据，是近年发展起来的一种三维CT技术。与传统的二维CT相比，锥束CT具有很高的扫描速度，重建出的切片图像记录了物体内部各点的材质、密度等物性参数分布，具有切片序列连续、切片内和切片间的空间分辨率相同、精度高等特点，在逆向工程与工业内视等领域已显示出广阔的应用和发展前景。

在锥束CT系统中，射线源发出的X射线具有一定范围的能量分布(此种射线称为多色射线)，当射线与物质相互作用时，由于低能量光子的衰减量大于高能量光子，造成穿过一定厚度的物体后射线的平均能量增大，此时射线平均能量对应的线性衰减系数也不再是常数而是逐渐减小，导致投影值与穿越长度也不再呈线性关系。而CT重建算法是基于X射线是单能谱假设的，重建时直接用多色投影代替单色投影，导致在重建的切片图像上呈现杯状伪影(Cupping Artifacts)，严重时图像会产生变形，使得该切片内的结构、尺寸、密度、成分等物理化学性质无法准确地判读和计量，这种现象就称为射束硬化(Beam Hardening)。

射束硬化是锥束CT实际应用中必须解决的一个重大问题。目前，CT的射束硬化校正方法主要分为单能法和双能法两大类。由于操作上的复杂性，双能法在工程实践中很少被采用。单能法易于实现，实际应用效果也比较好，因此被广泛研究。杨民、路宏年、路远等人在《光学技术》(2003，29(2)：177-182)的文章“CT重构中射线硬化的校正研究”中提出的方法就是一种典型的单能校正法，其校正思路是：利用楔形模体来获取射线贯穿物体长度与多色投影值之间的关系曲线，再对该曲线进行拟合，然后从坐标原点对该曲线做切线，以该切线建立多色数据与单色数据的函数关系，从而达到硬化校正的目的。该方法实现起来简单，但要求有与被检试件相同材质的模体，这就影响了该方法的应用灵活性。

另外，目前的射束硬化校正方法多数是针对二维CT的，尽管有的方法通过改造可以推广到锥束CT，但实际开展的研究仍然较少。

发明内容

为了克服现有技术应用不够灵活、不能有效适用于锥束CT的不足，本发明提供了一种基于切片轮廓重投影的锥束CT射束硬化校正新方法，以对锥束CT应用于实际无损检测中出现的射束硬化问题进行有效校正。

本发明针对工业锥束CT无损检测中的射束硬化问题提出一种切片轮廓重投影的锥束CT射束硬化校正新方法，其特征在于包括下列步骤：

(1)对被检测零件进行圆轨迹锥束CT扫描，从平板探测器获取一组连续的多色投影图像数据，这些图像在获取过程中已经过必要的暗场校正、坏像素校正和增益校正，校正手段可采用平板探测器厂方配套程序进行，也可以自行根据公知技术开发相应程序进行；

(2)对多色投影进行锥束CT重建得到一个中心切片图像；若中心切片没有包含射线穿越物体的最大长度，即所有投影图像的最大灰度值像素不位于中心切片位置，则在所有投影图像的最大灰度值像素位置附近重建若干连续的切片图像，并利用这些切片图像根据剖切体数据的公知方法生成一个过最大灰度位置的斜截面图像；

(3)对中心切片图像采用数字图像处理中的相关算法进行轮廓提取，获得若干封闭的中心切片零件截面轮廓；若有斜截面图像，则进行相同处理，获取若干封闭的零件斜截面轮廓；

(4)对中心切片零件截面轮廓按实际锥束CT扫描成像几何关系进行重投影，得到射线穿越长度-多色投影灰度间对应的数据；若有零件斜截面轮廓，则进行相同处理，并将得到射线穿越长度-多色投影灰度间对应的数据与由中心切片轮廓得到的数据合并成一组数据；

(5)采用指数函数L＝A·exp(B·G^C)-A来拟合射线穿越长度-多色投影灰度间对应的数据，其中L是射线穿越物体长度，G是多色投影灰度值，A、B、C是拟合系数；

(6)将射线穿越长度-多色投影灰度间对应的数据靠近原点的一部分进行线性拟合，得到一个过原点的校正直线函数；

(7)根据指数函数和校正直线函数进行锥束CT射束硬化校正计算：对各幅多色投影图像的每个像素，先将其灰度代入指数函数计算得到对应长度，然后将该长度代入校正直线函数得到校正后的近似单色投影灰度值。

上述的基于配准模型仿真的锥束CT射束硬化校正方法，不仅可以应用于锥束CT，还可以应用于扇束CT。