CN109919868B - 一种锥束ct射束硬化曲线侦测及投影加权校正方法 - Google Patents

Abstract

本发明公开了一种锥束CT射束硬化曲线侦测及投影加权校正方法，利用射线采样矩阵信息先验知识和预先重建的切片图像，以射线与二值化图像进行求取交线长度，侦测得到射束硬化曲线。通过对拟合函数进行投影加权修正，实现多能投影的硬化校正，改善因多能投影重建引起的杯状伪影，提高图像质量。本发明提供的锥束CT射束硬化曲线侦测及投影加权校正方法，适用于任意复杂度的被测物体的硬化曲线侦测及杯状伪影校正，方法的可靠性、稳定性、抗噪性好，可在很大程度上减少锥束CT硬化投影杯状伪影对图像的干扰与影响，明显改善锥束CT图像质量。

Description

一种锥束CT射束硬化曲线侦测及投影加权校正方法

技术领域

本发明属于锥束CT应用相关的医学成像和工业无损检测领域，涉及一种锥束CT射束硬化曲线侦测及投影加权校正方法。

背景技术

锥束CT(Cone Beam Computed Tomography,CBCT)作为一种先进的医学成像和工业无损检测技术，在不破坏物体情况下，以二维或三维断层图像的形式清晰、准确、直观地展现被检测物体的内部结构，定量地提供物体内部缺陷位置和尺寸。

实际工业射线检测中使用的射线不是单色射线，而是宽束、连续谱的多色射线。且工业 CT扫描能量较高，由于不同能量射线衰减程度不同，射线穿越物体后能谱发生变化。射束硬化伪影是影响重建图像质量的重要因素。在多色X射线谱中，低能光子更容易被吸收，随着射线穿透深度的增加，高能光子所占比例增加，射线平均能量升高，重建时依然以单能等效能量进行计算，导致重建图像出现杯状伪影，各点灰度值同实际物体的吸收系数之间出现偏差，降低了图像质量。

目前，X射线的射束硬化校正方法主要分为单能法和双能法。多项式拟合法是一种单能射束硬化校正方法，虽然实现起来简单，但灵活性差。Zou X等人在Journal ofMedical Imaging &Health Informatics(2016,6(7):1701-1707)的文章“TV-BasedCorrection for Beam Hardening in Computed Tomography”中提到的后处理法是基于重建-滤波-投影的迭代重建校正方法。一般先用前投影法校正，然后对重建出的切片进行门限值滤波，然后再进行投影、重建。该方法对各部分的材料吸收系数差别较大的物体有一定的效果，但实施起来复杂、费时。Abbema J等人在Physica Medica(2012,28(1):25-32)的文章“Feasibility and accuracy of tissue characterization with dual sourcecomputed tomography”中提出的双能量法是理论上消除射束硬化的最佳方法，它是对引起射束硬化现象的原因进行建模，实现物***置函数和能量函数的分离，但该方法需要知道X射线的能谱分布。

上述不同方法虽然取得了一定的射束硬化伪影抑制效果，但对于锥束CT医学成像要求和工业无损检测需求，已有方法的通用性较差，实际应用中往往受到诸多限制。

发明内容

针对多能X射线衰减系数随着射线穿越路径的变长而减小，出现投影灰度随射线路径长度的非线性变化，重建时产生的杯状伪影问题，本发明提供一种锥束CT射束硬化曲线侦测及投影加权校正方法。利用射线采样矩阵信息先验知识和预先重建的切片图像，以射线与二值化图像进行求取交线长度，侦测得到射束硬化曲线。通过对拟合函数进行投影加权修正，实现多能投影的硬化校正，改善因多能投影重建引起的杯状伪影，提高图像质量。

本发明解决其技术问题所采用的技术方案包括以下步骤：

(1)设置投影采样几何参数；

(2)对多色投影P_p进行预重建，获取中心切片图像S对应的二值化图像B；

(3)获取射线与二值化图像的交线长度，拟合得到侦测硬化曲线；

(4)对侦测硬化曲线进行投影加权及修正；

(5)以修正后的函数模型对多能射线进行校正，三维重建获取切片。

在上述步骤(1)中，投影采样几何参数包括：射线源到探测器中心距离、射线源到旋转中心距离、重建分辨率及体素大小，采样幅数，采样的起始位置。

在上述步骤(1)中，获取射线采样的矩阵信息的具体步骤包括：

1)将待重建图像离散化，以M×N表示离散网格图像，称之为D，M和N分别表示沿 X方向和Y方向的网格数目；

2)构建平面坐标系，设置射线源起始坐标及探测器探元坐标，建立射线源与探测器像元之间的射线方程；

3)计算射线与二维离散网格图像D中每个像素的交线长度，以此交线长度作为对应像素位置ID的权值系数，不相交的像素位置，记为0；

4)遍历每条射线，以每条射线与二维离散网格图像D中每个像素的交线长度信息作为先验知识，存入采样矩阵。

在上述步骤(2)中，对多色投影预重建之前，需要先对投影进行散射校正，散射校正完成后需要对校正投影进行对数变换然后重建，重建时只需重建中心层切片即可，重建方法可以选择滤波反投影算法，也可以选择迭代类算法，然后对中心切片图像S根据OTSU大律法得到二值化图像B。

在上述步骤(3)中，获取射线与二值化图像的交线长度的具体步骤包括：

1)构建平面坐标系，设置射线源起始坐标及探测器探元坐标，建立射线源与探测器像元之间的射线方程，根据平面坐标系，将步骤2中得到的二值化图像B离散化，以M×N表示离散网格图像，称之为B1，M和N分别表示沿X方向和Y方向的网格数目；

2)选择不同角度下的射线，射线索引，记为R，求射线方程与二值化图像的离散网格图像B1的交点，将单个网格内的线段长度作为射线与该网格的交线长度，并记录该网格索引为

以此类推，获取同一射线下获得的全部网格索引

并对

对应的交线长度求和，作为该条射线穿越路径长度。

在上述步骤(3)中，拟合得到侦测硬化曲线的具体步骤包括：

1)首先根据射线索引R，提取射线索引R在探测器上投影的灰度值

获取射线索引R 对应的射线穿越路径长度Ls；

2)建立射线穿越路径长度Ls和投影灰度

的对应关系，即硬化离散数据点

3)选择函数L＝At^B+Ct，以最小二乘法拟合离散数据点

得到射线穿越长度Ls 关于投影灰度

变化的射束硬化侦测曲线，以A、B、C作为拟合系数，L表示射线穿越物体长度，t为多色投影灰度值。

在上述步骤(4)中，对侦测硬化曲线进行修正及投影加权的具体步骤包括：

1)以拟合系数C对函数模型就行修正，得到更新模型

2)对更新模型增加投影加权补偿项λt，得到最终模型

其中λ表示加权系数，计算方法为：

在上述方法中，如果检测对象的扫描方式及扫描参数与上次检测相同或批量检测时，则步骤(3)中侦测的硬化曲线只需要获取一次，后面的不同检测对象可以直接应用，无需再次获取。

本发明的有益效果是：本发明提供的锥束CT射束硬化曲线侦测及投影加权校正方法，适用于任意复杂度的被测物体的硬化曲线侦测及杯状伪影校正，方法的可靠性、稳定性、抗噪性好，可在很大程度上减少锥束CT硬化投影杯状伪影对图像的干扰与影响，明显改善锥束CT图像质量。

下面结合附图和具体实施方式对本发明进行详细说明。

附图说明

图1为本发明算法流程图。

图2为重建图像杯状伪影校正前后相同位置的线性灰度比较。

具体实施方式

通过现有的工业锥束CT设备(X射线源为Comet的MXR-451HP/11，平板探测器为PerkinElmer的XRD 1621AN15 ES，并具备扫描机构、***控制及计算用计算机)，对两个零件进行投影采样，应用本发明方法对锥束CT射束硬化曲线侦测及投影加权校正方法，执行以下步骤：

(1)通过多能谱射线源工业锥束CT设备，选择射线源电压420kV和电流0.75mA，扫描几何参数为：射线源到探测器距离1241.832886mm，射线源到旋转中心距离965.323242mm；重建分辨率为512×512，圆周扫描得到检测对象锥束CT投影360幅。

1)图像网格离散化成512×512，图像网格与投影采样中的分辨率设置相同；

2)构建平面坐标系，设置射线源起始坐标及探测器探元坐标，建立射线源与探测器像元之间的射线方程；

3)计算射线与二维离散网格图像D中每个像素的交线长度，以此交线长度作为对应像素位置ID的权值系数，不相交的像素位置，记为0；

4)遍历每条射线，以每条射线与二维离散网格图像D中每个像素的交线长度信息作为先验知识，构建采样矩阵M,大小为360×262144(512×512＝262144)。

(2)对多色投影预重建之前，需要先对投影进行散射校正，散射校正完成后需要对校正投影进行对数变换然后重建，重建时只需重建中心层切片即可，重建方法可以选择滤波反投影算法，也可以选择迭代类算法，然后对中心切片图像S根据OTSU大律法得到二值化图像B。

(3)获取射线与二值化图像的交线长度，拟合得到侦测硬化曲线，具体步骤包括：

1)构建平面坐标系，设置射线源起始坐标及探测器探元坐标，建立射线源与探测器像元之间的射线方程，根据平面坐标系，将步骤2中得到的二值化图像B离散化成512×512的 B1图像；

2)选择不同角度下的射线，射线索引，记为R，求射线方程与二值化图像的离散网格图像B1的交点，将单个网格内的线段长度作为射线与该网格的交线长度，并记录该网格索引为

3)以此类推，获取同一射线下获得的全部网格索引

并对

对应的交线长度求和，作为该条射线穿越路径长度；

4)根据射线索引R，提取射线索引R在探测器上投影的灰度值

获取射线索引R对应的射线穿越路径长度Ls，建立射线穿越长度和对应投影灰度关系，得到硬化离散数据点；

5)选择函数L＝At^B+Ct，以最小二乘法拟合得到射线穿越长度关于多色投影灰度变化的射束硬化侦测曲线，以A、B、C作为拟合系数，L表示射线穿越物体长度，t为多色投影灰度值。

(4)对侦测硬化曲线进行投影加权及修正，具体步骤包括：

1)以拟合系数C对函数模型就行修正，得到更新模型

2)对更新模型增加投影加权补偿项λt，得到最终模型

其中λ表示加权系数，计算方法为：

(5)以修正后的函数模型对多能射线进行校正，三维重建获取切片。

本实施例中，基于多能谱射线硬化投影的锥束CT射束硬化曲线侦测及投影加权校正方法特点在于：

(1)根据检测对象的投影信息结合先验知识侦测得到射束硬化曲线，描述并表示了具体对象的射束硬化行为；

(2)通过对硬化曲线模型拟合完成了模型修正，经过投影加权补偿项，使得硬化校正模型更精准的建立多能投影到单能投影之间的映射；

(3)以最终模型完成多能谱硬化投影校正，抑制杯状伪影，获得高质量图像。

图2为重建图像杯状伪影校正前后相同位置的线性灰度比较，可见本发明方法可使硬化投影重建图像的杯状伪影进行抑制，使得图像轮廓对比度和清晰度得到显著提高。

Claims (3)

1.一种锥束CT射束硬化曲线侦测及投影加权校正方法，其特征在于包括下述步骤：

步骤1：设置投影采样几何参数；

步骤2：对多色投影P_p进行预重建，获取中心切片图像S对应的二值化图像B；

步骤3：获取射线与二值化图像的交线长度，拟合得到侦测硬化曲线；

步骤4：对侦测硬化曲线进行投影加权及修正；

步骤5：以修正后的函数模型对多能射线进行校正，三维重建获取切片；

在所述步骤1中，投影采样几何参数包括：射线源到探测器中心距离、射线源到旋转中心距离、重建分辨率及体素大小，采样幅数，采样的起始位置；

在所述步骤2中，对多色投影预重建之前，需要先对投影进行散射校正，散射校正完成后需要对校正投影进行对数变换然后重建，重建时只需重建中心层切片即可，重建方法可以选择滤波反投影算法，也可以选择迭代类算法，然后对中心切片图像S根据OTSU大律法得到二值化图像B；

在所述步骤3中，获取射线与二值化图像的交线长度的具体步骤包括：

1)构建平面坐标系，设置射线源起始坐标及探测器探元坐标，建立射线源与探测器像元之间的射线方程；

2)根据平面坐标系，将步骤2中得到的二值化图像B离散化，以M×N表示离散网格图像，称之为B1，M和N分别表示沿X方向和Y方向的网格数目；

3)选择不同角度下的射线，射线索引，记为R，求射线方程与二值化图像的离散网格图像B1的交点，将单个网格内的线段长度作为射线与该网格的交线长度，并记录该网格索引为

4)以此类推，获取同一射线下获得的全部网格索引

并对

对应的交线长度求和，作为该条射线穿越路径长度；

在所述步骤3中，拟合得到侦测硬化曲线的具体步骤包括：

1)首先根据射线索引R，提取射线索引R在探测器上投影的灰度值

2)获取射线索引R对应的射线穿越路径长度Ls；

3)建立射线穿越路径长度Ls和投影灰度

的对应关系，即硬化离散数据点

4)选择函数L＝At^B+Ct，以最小二乘法拟合离散数据点

得到射线穿越长度Ls关于投影灰度

变化的射束硬化侦测曲线，以A、B、C作为拟合系数，L表示射线穿越物体长度，t为多色投影灰度值；

在所述步骤4中，对侦测硬化曲线进行修正及投影加权的具体步骤包括：

1)以拟合系数C对函数模型就行修正，得到更新模型

2)对更新模型增加投影加权补偿项λt，得到最终模型

其中λ表示加权系数，计算方法为：

2.根据权利要求1所述的一种锥束CT射束硬化曲线侦测及投影加权校正方法，其特征在于：在所述方法中，如果检测对象的扫描方式及扫描参数与上次检测相同或批量检测时，则步骤3中侦测的硬化曲线只需要获取一次，后面的不同检测对象可以直接应用，无需再次获取。

3.根据权利要求1所述的一种锥束CT射束硬化曲线侦测及投影加权校正方法，其特征在于：基于多能谱射线硬化投影的锥束CT射束硬化曲线侦测及投影加权校正方法中：

(1)根据检测对象的投影信息结合先验知识侦测得到射束硬化曲线，描述并表示了具体对象的射束硬化行为；

(2)通过对硬化曲线模型拟合完成了模型修正，经过投影加权补偿项，使得硬化校正模型更精准的建立多能投影到单能投影之间的映射；

(3)以最终模型完成多能谱硬化投影校正，抑制杯状伪影，获得高质量图像。

Images