CT形状过滤器缺陷的双能矫正方法与流程

本发明涉及一种ct形状过滤器缺陷的双能矫正方法，主要可用于人体计算机断层扫描的医用ct机中的形状过滤器缺陷的矫正。

背景技术：

形状过滤器(例如，bowtiefilter)是医用ct机中的核心部件之一，用于x射线射束硬化校正以获得良好、稳定的图像质量，为了获得形状过滤器几何尺寸的精确数据，或更准确地讲，为了获得x射线在从源达到各探测器单元之前在形状过滤器中所穿行的路径长度或总路径长度，一种直截了当的方法是使用高品质材料和精加工工艺使得实际的形状过滤器与设计的形状过滤器高度一致，但这样做无疑会提高产品成本，并且能够达到的精度还受到加工技术的限制；另一种方式是通过后续复杂的矫正技术并配合额外的模体在投影数据中进行补偿，这同样在很大程度上导致最终用户成本的增加以及后期维护的难度。

技术实现要素：

为了克服现有技术的上述缺陷，本发明提供了一种ct形状过滤器缺陷的双能矫正方法，利用这种方法可以方便地计算获得形状过滤器的实际几何尺寸或相关数据，以便在数据处理中对形状过滤器缺陷进行矫正，进而减小用户成本和维护难度。

本发明所采用的技术方案：一种ct形状过滤器缺陷的双能矫正方法，其以同一ct机，在安装好形状过滤器状态下，在x射线源球管电压分别取两个不同值的两种状态下分别进行空气扫描，依据探测器单元在所述两种状态下获得的测量数据，推算出达到该探测器单元的x射线在形状过滤器中所穿行的路径长度。

可以采用下列公式计算到达探测器单元n的x射线在形状过滤器中所穿行的路径长度lb(n)：

其中，

in(kv1)为在x射线源球管电压为kv1的状态下扫描，探测器单元n被x射线照射后产生的光电流测量值；

in(kv2)为在x射线源球管电压为kv2的状态下扫描，探测器单元n被x射线照射后产生的光电流测量值；

itube，1为在x射线源球管电压为kv1的状态下扫描所用的x射线源球管电流；

itube，2为在x射线源球管电压为kv2的状态下扫描所用的x射线源球管电流；

tview，1为在x射线源球管电压为kv1的状态下扫描一个视角所需要的时间；

tview，2为在x射线源球管电压为kv2的状态下扫描一个视角所需要的时间；

e为x射线源球管产生的光子的能量；

为到达探测器单元n的光子能量为e的x射线在传播过程中因形状过滤器之外的其他所有材料导致的衰减率；

nkv(e)为x射线源球管电压为kv时x射线源球管的光子数能谱，nkv(e)de代表每个打在x射线源球管阳极上的电子能产生的能量在e到e+de范围内的光子数，kv＝kv1，kv2；

μb(e)为形状过滤器的材料的线性衰减系数；

ηd(e)为探测器或探测器单元的能量探测效率；

n为探测器单元编号，编号为n的探测器单元称为探测器单元n。

可以采用下列公式计算到达探测器单元n的x射线在形状过滤器中所穿行的路径长度lb(n)：

其中，

in(kv1)为在x射线源球管电压为kv1的状态下扫描，探测器单元n被x射线照射后产生的光电流测量值；

in(kv2)为在x射线源球管电压为kv2的状态下扫描，探测器单元n被x射线照射后产生的光电流测量值；

e为x射线源球管产生的光子的能量；

为到达探测器单元n的光子能量为e的x射线在传播过程中因形状过滤器之外的其他所有材料导致的衰减率；

是x射线源球管电压为kv时归一化的光子数能谱，nkv(e)为x射线源球管电压为kv时x射线源球管的光子数能谱，nkv(e)de代表每个打在x射线源球管阳极上的电子能产生的能量在e到e+de范围内的光子数，kv＝kv1，kv2；

μb(e)为形状过滤器的材料的线性衰减系数；

ηd(e)为探测器或探测器单元的能量探测效率；

iref(kv1)为在x射线源球管电压为kv1的状态下扫描，参考探测器单元被x射线照射后产生的光电流测量值；

iref(kv2)为在x射线源球管电压为kv2的状态下扫描，参考探测器单元被x射线照射后产生的光电流测量值；

为到达参考探测器单元的光子能量为e的x射线在传播过程中因形状过滤器之外的其他所有材料导致的衰减率；

ηref(e)为参考探测器单元的能量探测效率；

lb(nref)为已知的到达参考探测器单元的x射线在形状过滤器中所穿行的路径长度；

n为探测器单元编号，编号为n的探测器单元称为探测器单元n。

本发明的有益效果：能够方便地实现对形状过滤器的矫正，在不改变现有工艺、材料品质以及设备附件、也不需要拆卸形状过滤器的前提下，只需要在不同x射线源球管电压的两个状态下进行空气扫描，由于两种状态下测量数据的差异部分源于形状过滤器的设置，由此通过两种状态下的测量数据之间的比对，就能够揭示出形状过滤器对测量数据的影响，计算出x射线在形状过滤器中所穿行的路径长度，进而获得有效的形状过滤器的几何尺寸数据或将相应的路径长度数据用于其他相关数据的运算，实现对形状过滤器缺陷的双能矫正，相应路径长度数据的计算可以依据x射线的传播规律以及实际中的各影响因素，采用任意适宜的现有技术或其他可能的技术实现，由此无需要求将形状过滤器加工得与设计完全一致，省略了现有技术下矫正补偿模体的相关设计、加工及使用和维护，由此大幅度降低了形状过滤器的加工难度和制备成本，方便了使用，降低了维护费用，并且计算方式和计算过程也较为简便。另外，能够实现的补偿精度与ct机的测量精度相关，不受形状过滤器加工技术和加工精度的限制。

附图说明

图1给出了在已知光子数能谱的情形下，计算得到的一种形状过滤器的计算路径长度曲线；为了比较，图中也给出了形状过滤器在整个探测器范围的实际路径长度曲线和理论设计路径长度曲线；在该图的显示分辨率下，计算路径长度曲线和实际路径长度曲线几乎完全重叠；

图2是与图1对应的形状过滤器在探测器中心区域范围的计算路径长度曲线、实际路径长度曲线和理论设计路径长度曲线，在该图的显示分辨率下，计算路径长度曲线和实际路径长度曲线基本重叠；

图3给出了在已知沿边缘探测器单元(参考探测器单元)光路上穿过的各种材料的厚度的情形下，计算得到的一种形状过滤器的计算路径长度曲线；为了比较，图中也给出了形状过滤器在整个探测器范围的实际路径长度曲线和理论设计路径长度曲线；在该图的显示分辨率下，计算路径长度曲线和实际路径长度曲线几乎完全重叠；

图4是与图3对应的形状过滤器在探测器中心区域范围的计算路径长度曲线、实际路径长度曲线和理论设计路径长度曲线，在该图的显示分辨率下，计算路径长度曲线和实际路径长度曲线基本重叠；

图5给出了归一化的x射线源球管电压为80kvp和140kvp的实验光谱和模型光谱。

具体实施方式

本发明用于计算x射线在形状过滤器中所穿行的路径长度lb(n)的公式可以通过下列方式推导出：

根据探测器单元n被x射线照射后光电流的产生机理，可知：

其中αn是公式推导过程中引入的只与探测器性质有关的比例系数，单位可为纳安·千电子伏^-1；

ad为相应探测器单元在垂直于x射线入射方向的投影面积；

dsd为x射线源球管到相应探测器单元的距离。

在已知光子数能谱和的前提下，依据公式(1)和(2)可得

如果已知到达某个特定探测器单元的x射线在形状过滤器中所穿行的路径长度，或者已知沿某个特定探测器单元光路上穿过的各种材料的厚度进而可以计算获得到达该特定探测器单元的x射线在形状过滤器中所穿行的路径长度，可以以该特定探测器单元作为参考探测器单元，根据参考探测器单元被x射线照射后光电流的产生机理，可知：

其中αref是公式推导过程中引入的只与参考探测器单元性质有关的比例系数，单位可为纳安·千电子伏^-1。

aref为参考探测器单元在垂直于x射线入射方向的投影面积；

dref为x射线源球管到参考探测器单元的距离。

依据公式(1)-(2)和(4)-(5)可得

公式(3)和(6)中积分的范围为实际中有效的光子能量e的范围，可以依据具体的计算精度要求通过试验或理论分析确定。利用公式(3)或(6)求解获得路径长度lb(n)，即可用于涉及形状过滤器各种计算，实现对其缺陷的矫正。

进行空气扫描时，优选以等视角间距进行360度扫描。

在此情形下，可以各视角下的光电流测量值的平均值用作计算路径长度lb(n)。

具体地，在已知光子数能谱的情形下计算路径长度lb(n)时，所述空气扫描通常可优选为等视角间隔的360度扫描，在此扫描下，分别获得各视角下探测器单元n被x射线照射后产生的光电流测量值，并以各状态下探测器单元n在各视角的光电流测量值的平均值作为相应状态下的光电流测量值in(kv1)或in(kv2)。

具体可以采用下列步骤：

1)在x射线源球管电压为kv1的状态下进行一次空气扫描，获得检测数据in(kv1)；

2)用同一ct机在x射线源球管电压为kv2的状态下进行一次空气扫描，获得检测in(kv2)；

3)利用x射线源球管的已知相关参数，计算和

4)利用ct机的已知相关参数，分别计算和ηd(e)；

将上述in(kv1)、in(kv2)、ηd(e)、itube，1、itube，2、tview，1、tview，2以及形状过滤器材料的线性衰减系数μb(e)代入相应方程，求解出达到探测器单元n的x射线在形状过滤器中穿行的路径长度lb(n)。

在已知x光在达到参考探测器单元之前在形状过滤器中所穿行的路径长度为lb(nref)的情形下计算路径长度lb(n)时，所述空气扫描通常可优选为等视角间隔的360度扫描，在此扫描下，分别获得各视角下探测器单元n被x射线照射后产生的光电流测量值，并以各状态下探测器单元n在各视角的光电流测量值的平均值作为相应状态下的光电流测量值in(kv1)和in(kv2)；分别获得各视角下参考探测器单元被x射线照射后产生的光电流测量值，并以各状态下参考探测器单元在各视角的光电流测量值的平均值作为相应状态下的光电流测量值iref(kv1)和iref(kv2)。

具体可以采用下列步骤：

1)在x射线源球管电压为kv1的状态下进行一次空气扫描，获得检测数据in(kv1)和iref(kv1)；

2)用同一ct机在x射线源球管电压为kv2的状态下进行一次空气扫描，获得检测in(kv2)和iref(kv2)；

3)利用x射线源球管的已知相关参数，计算x射线源球管的归一化光子数能谱和

4)利用ct机的已知相关参数，分别计算和以及ηd(e)和ηref(e)；

5)将上述in(kv1)、iref(kv1)、in(kv2)、iref(kv2)、ηd(e)、ηref(e)、lb(nref)和形状过滤器材料的线性衰减系数μb(e)代入相应方程，求解出达到探测器单元n的x射线在形状过滤器中穿行的路径长度lb(n)。

图1、图2、图3和图4显示了在相应情形下(所用x射线源球管电压为80kvp和140kvp，x射线源球管电流均为50ma)，依据本发明计算出的计算路径长度曲线与实际路径长度曲线和理论设计路径长度曲线之间的对比。计算时假设了，如图5所示，模型光谱与x射线源球管电压为80kvp和140kvp的实验光谱只有2％的相对误差。结果表明，用本发明的方法，依据方程(3)计算得出的x射线在形状过滤器中的路径长度与真实路径长度的平均误差为0.0125mm，最大误差为0.0214mm；依据方程(6)计算得出的x射线在形状过滤器中的路径长度与真实路径长度的平均误差为0.0095mm，最大误差为0.0188mm。由此可以看出，本发明能够准确地计算出达到各探测器单元的x射线在形状过滤器中的路径长度。

关于本发明计算公式涉及的部分参数的进一步说明：

lb(n)为x光在从x射线源球管达到探测器单元n之前在形状过滤器中所穿行的路径长度，为本发明所要计算的路径长度，在ct检测的相关数据处理中，这些数据作为到达相应的各检测器单元的x光经过形状过滤器的实际路径长度或者对设计路径长度进行补偿，能够有效地消除因形状过滤器实际路径长度与设计路径长度不一致导致的误差；

in(kv1)和in(kv2)分别为在x射线源球管电压为kv1和kv2的状态下进行空气扫描，探测器单元n被x射线照射后产生的光电流测量值。当采用多个不同视角进行空气扫描时，可以以各视角下的光电流测量值的平均值(通常可以为算术平均)作为光电流测量值。探测器单元n被x射线照射后产生的光电流可依据探测器单元n上的读数，单位可采用纳安(10^-9安培)；

iref(kv1)和iref(kv2)分别为在x射线源球管电压为kv1和kv2的状态下进行空气扫描，参考探测器单元被x射线照射后产生的光电流测量值。当采用多个不同视角进行空气扫描时，可以以各视角下的光电流测量值的平均值(通常可以为算术平均)作为光电流测量值。参考探测器单元被x射线照射后产生的光电流可依据参考探测器单元上的读数，单位可采用纳安(10^-9安培)；

itube，1和itube，2分别为在x射线源球管电压为kv1和kv2的状态下扫描所用的x射线源球管电流，单位可采用毫安；

tview，1和tview，2分别为在x射线源球管电压为kv1和kv2的状态下扫描一个视角所需要的时间；

e代表x射线源球管产生的光子的能量，单位可为千电子伏；

和分别为代表到达探测器单元n和参考探测器单元的x射线在传播过程中在除形状过滤器外的其他材料中经历的衰减率(以e为自变量的函数)，这些函数是ct机或ct系统的已知函数，由ct机除形状过滤器外的系统特性确定，依赖于变量e，这些函数也可以依据理论分析和/或实验数据确定；

nkv(e)为光子数能谱，nkv(e)de代表每个打在x射线源球管阳极上的电子能产生的能量在e到e+de范围内的光子数；

为归一化的光子数能谱，依据前面所述的公式计算；

μb(e)为形状过滤器的材料的线性衰减系数，可通过理论计算和/或实验获得；

ηd(e)和ηref(e)分别为探测器单元n和参考探测器单元的能量探测效率，可通过理论计算和/或实验获得；

lb(nref)为x光在从x射线源球管达到参考探测器单元之前在形状过滤器中所穿行的路径长度；

kv1和kv2的是进行空气扫描时使用的两个不同的x射线源球管电压；

n为探测器单元编号，编号为n的探测器单元称为探测器单元n，可以顺序编号，也可以采用其他编号方式。

本发明涉及x射线的路径或路径长度均指相应x射线的实际路径(光路)或路径长度。

本发明所称x射线源球管泛指ct机所用的用于产生x射线的任意形式的x射线源，可以采用常见的球管形式，也可以采用其他形式。

本发明公开的各优选和可选的技术手段，除特别说明外及一个优选或可选技术手段为另一技术手段的进一步限定外，均可以任意组合，形成若干不同的技术方案。