[发明专利]一种CT探测器像素响应不一致性的校正方法

权利要求书

1.一种CT探测器像素响应不一致性的校正方法，其特征在于，该方法包括：

将原始的投影正弦图在多个投影角内进行平均，获得初始平均响应曲线；

去除初始平均响应曲线中的轮廓线；

在去除轮廓线的平均响应曲线中探测响应不一致的像素点，并将响应不一致的像素点分为两类：单双点，以及连续三点及以上；

对于单双点的响应不一致像素点使用邻域中值滤波的方法估计偏差；

对于连续三点及以上的响应不一致像素点使用基于平滑性约束的填充方法估计偏差；

根据估计出的偏差补偿原始的投影正弦图，完成CT探测器像素响应不一致性的校正。

2.根据权利要求1所述的方法，其特征在于，所述将原始的投影正弦图在多个投影角内进行平均，获得初始平均响应曲线包括：

沿投影数n方向将原始的投影正弦图P(m,n)进行平均滤波，得到再将沿投影数n方向进行下采样得到初始平均响应曲线

其中，m表示探测器像素点坐标，n表示投影数。

3.根据权利要求2所述的方法，其特征在于，所述去除初始平均响应曲线中的轮廓线包括：

沿探测器像素点坐标m方向使用savitzky-golay滤波器进行滤波获取轮廓线Z(m,n)，表示为：

Z(m,n)=Σw=-LLP^(m-w,n)SG(w);]]>

其中，SG(w)表示使用savitzky-golay滤波器第w个系数进行滤波，平均滤波器的长度为M+1；

将获取轮廓线Z(m,n)从初始平均响应曲线中减去，再沿投影数n方向进行一次低频滤波，获得去除轮廓线的平均响应曲线R，其中的像素点的值记为R(m,n)，表示为：

R^(m,n)=P^(m,n)-Z(m,n);]]>

R(m,n)=Σi=n-M/2n+M/2R^(m,i).]]>