[发明专利]基于双能CT扫描数据重建投影图像检测方法

权利要求书

1.一种基于双能CT扫描数据重建投影图像检测方法，其特征是，包括下列步骤：

双能减影图像由下式表示：

IBone-canceledDE=ln(I0H)-wsln(I0L)-ρeWater·EPL{[fKN(EH)-wsfKN(EL)]+[1/(EH)3.2-ws/(EL)3.2]Cp(ZWater)3.8}---(1)]]>

公式（1）中的参数获得及说明如下:

减影图像中像素的灰度值

_I0^H，I₀^L：高能和低能入射X射线的强度，在这里这一项为常数，因此计算模拟图像时忽略这一项的值；

E^H,E^L:高能X线和低能X线的能量，单位为keV；

f_KN是Klein-Nishina方程：

fKN(E)=C0{1+EE2[2(1+E)1+2E-1Eln(1+2E)]+12Eln(1+2E)-1+3E(1+2E)2}---(2)]]>

其中C₀为常数，C₀=5.0×10^-25；

Z_Water:水的有效原子序数；

ρ_eWater：水的电子密度；

Cp:常数，Cp=9.8×10^-24；

w_s：消去骨性结构或软组织的权重因子,取经验值0.5；

EPL：等效路径长度，即除去被消除的结构，

EPL=∑_iΔd_iρ_ei                    （3）

i：X射线穿过组织中的第i种物质；

Δd_i：X射线穿过物质i路径的距离；

ρ_ei：物质i的相对电子密度。

2.如权利要求1所述的基于双能CT扫描数据重建投影图像检测方法，其特征是，EPL的具体计算方法如下：

1）使用CT对仿体或者病人扫描得到该仿体或病人的三维CT数据集。

2）使用射线追踪算法计算获得由射线源穿过仿体到达影像板中某像素的射线，这里射线追踪算法中需要确定的光源与影像板的位置，影像板的大小及分辨率等参数根据需要模拟的实际条件设定，通过射线追踪算法模拟通过三维CT数据集的射线，获得该射线经过的体素坐标以及该体素的CT数以及射线在该体素内穿过的距离，并且根据体素的CT数将射线穿过三维CT数据集中的骨性结构区分出来；

3）根据该CT的电子密度—CT数校正曲线，对射线经过的体素的CT数使用插值方法求得其对应的电子密度，由该体素的电子密度与射线在体素内经过的距离可以计算射线经过该体素的EPL，对整条射线中除骨性结构外的体素EPL求和即可得到由射线源到影像板内某一像素射线经过的EPL。