1.一种锥束CT系统的金属伪影校正方法,其特征在于,包括下述步骤:从原始正投影图像f(x,y)中分离出金属投影图像M(x,y);
对所述金属投影图像M(x,y)进行重建,得到金属部分的CT图像XMetal;
将所述原始正投影图像f(x,y)减去所述金属投影图像M(x,y),即得到不含金属部分的投影图像fres(x,y);
对所述不含金属部分的投影图像fres(x,y)进行重建,得到不含金属部分的CT图像Xres;
及
将所述金属部分CT图像XMetal与所述不含金属部分的CT图像Xres相加,即得到金属伪影校正后的最终CT重建图像Xcorrection,其中,从原始正投影图像f(x,y)中分离出金属投影图像M(x,y),包括下述步骤:步骤S11:定义初始阈值I0,所述初始阈值为I0=(Imax+Imin)/2,其中,Imax为所述原始正投影图像f(x,y)中像素最大灰度值,Imin为所述原始正投影图像f(x,y)中最小灰度值;
步骤S12:基于所述初始阈值I0,将所述原始正投影图像f(x,y)分为两组图像T1(x,y)、T2(x,y),其中,T1(x,y)的所有像素灰度值大于I0,T2(x,y)的所有像素灰度值均小于或等于I0;
步骤S13:计算所述图像T1(x,y)内所有像素的均值及图像T2(x,y)内所有像素的均值,并分别记为E1和E2;
步骤S14:定义新阈值I,所述新阈值I=(E1+E2)/2;
步骤S15:判断|I-I0|是否大于ε,若“是”则进行下一步,若“否”则将所述I值赋予I0值,并返回步骤S12,其中,所述ε为0.5;
步骤S16:将所述新阈值I赋予所述初始阈值I0;
步骤S17:所述金属投影图像M(x,y)为原始正投影图像f(x,y)中的像素灰度值大于I0的像素。