[发明专利]一种基于扇束X光CT滤波反投影重建的混合滤波方法

权利要求书

1.一种基于扇束X光CT滤波反投影重建的混合滤波方法，其特征在于，包括以下步骤：

（1）获取扇束X光CT投影数据；

（2）对投影数据进行混合滤波处理；

（3）沿着X射线方向进行反投影重建得到断层图像。

2.根据权利要求1所述一种基于扇束X光CT滤波反投影重建的混合滤波方法，其特征在于，步骤（1）中，获取扇束X光CT投影数据具体为：沿着圆形扫描轨迹每隔一个角度，采集一个投影数据，一共获取360个采样数据。

3.根据权利要求1所述一种基于扇束X光CT滤波反投影重建的混合滤波方法，其特征在于，步骤（2）中，对投影数据进行联合滤波处理的具体方法为：

（2-1）对获取的X光CT投影数据进行离散化，设定离散化后的投影数据为p(s,β)；

其中，采样点s=nd，采样间隔d=1，采样角度β=1～2π；

（2-2）对离散化后的投影数据p(s,β）进行混合滤波处理，具体为：

设定h(s)为混合滤波器，根据p(s,β)=p(s,β)**h(s)对投影数据p(s,β)进行卷积混合滤波，得到滤波后的数据

其中，**号表示卷积。

4.根据权利要求3所述一种基于扇束X光CT滤波反投影重建的混合滤波方法，其特征在于，步骤（2）中对投影数据进行混合滤波处理，其混合滤波器为：

h(nd)=λ₁h_SL(nd)+λ₂h_L(nd)；

加权系数：λ₁≥0，λ₂≥0，λ₁+λ₂=1，用户根据实际需要确定；

SL滤波器：n=0，±1，±2，±3……；

Lewitt滤波器：

5.根据权利要求1所述一种基于扇束X光CT滤波反投影重建的混合滤波方法，其特征在于，步骤（3）具体为：

（3-1）设定待重建断层图像内任一点的位置为(r,θ)，反投影地址为s₁，X光射线源到重建中心的距离为D，根据以下公式计算反投影地址：

s1=Drcos(β-θ)D+rsin(β-θ)]]>

（3-2）设定待重建断层图像为f(r,θ)，滤波反投影数据为根据以下公式进行累加反投影重建：

f(r,θ)=∫12πp~(s1,β)dβ.]]>