X-射线CT卷积后投影再建算法中滤波的选取.pdfVIP

R‰‘)、 x-射线GT卷积后投影再建算法中滤波核的选取 段；／9 李海芳 李梅 太艨理J^擎， ：太原理工^学) ■薹露营．X-qT太多散秉用巷识后投影再建算法，投髟尊菠植的选取对再建质量弗响艰大．折囊地袁粗棱的高攫段 盛较好地挣村警声．提高再整质量。阜天就此饿，操1}． 关键词毒根舌投雾． 宙囊重建算法．藏最拄 计算机断层成像的目的是获得人体内一个断面中各个位置上物体性质的信。息．图像再琏 算法是cr的关键．而再建算法中巷积后投影苹法运算谴度决，再建质量好，在卷积后投影算 法中滤波核和内插法的选取是很重要的，因此本文主要对巷积后投影算法中滤波核的选取 进行7讨论，为改进cT成像质量提高成像速度做j一些尝试 X-CT得到的是人体一个横截面由组织对x射线衰减系数的分布，组织对x秆线的吸收 取决干它的密度，原子组成和x射线的能量，我们用函数i’(X、y)来描述所研窀截面钓衰减事 发分币(衰减／单位长度j．点源发出的射线穿过人体可得到 ，)(目．『1=f厂lTj)c，v (1) ．H 上式中n(o，t)称为函数“x，y)沼射线的线积分、又称为投黟，它可以通过测量得到，所谓臣 傈童建尉是由一系列f0，￡)的口(0，￡j解得式中雌一r(J．y) 在x—cT中常用的图像重建算法是变换法，式I1)可看作函数f的一种至换，昧为雪警鼍 按，明算于R表示，有： p(e．，)=月(厂j=f，(r，J-nLs(二) 腑 变换法就是寻找一种新的算子．隶投影的雷当逆雯抉满足：R。R(I、)=I从而得到再建像 r(x，n 卷积后投影算法就是一种变换法，它是宙当逆变抉【内近似完成，它是让投影P(o，【)先通 过一个特定的滤波，然后对变量t作一维内插，再求后投影． 付立叶切片定理是卷积后投影算法的根据，付立叶切片定理指出：在一个角度0所得到 与u轴成角0的一个切片，如图所示，p(0，t)对t的一维付立叶变换给出了F(u，v)在直线BB 上的值，从许多方向照射物体即0从0到2 n雯化，得到f(x，y)的全部投影p(o，1)，频域宅 问(u，V)就可以完全被填满，通过直接求F(u，vj的付氏逆变换就可以再建像，设s(e．∞)是 p(e，t)对t的一维付立叶变换，Ff∞，0)是fky)二维付立叶变换的极坐标表示，司得到： ‘ F(0，【|))=S{0．∞) 卷积后投影算法是在一系列e角(ot，o一…oJ上用一组平行射束作投影而取得 ·S87· 0龆蒜： V B p(0 稳◇ U ／啦 } I 添7 片定理示意盈 p¨，。，1)．．．．．．．，Dl u。，t)，对这些投影作爵敢付立叶变换获得Sf0．，“j)．．． 【{j。．∞)，另一方面，我们知道二维襄减密庋挣有函数I’(x，yj的二维付氏变换满足 ，l』，j，1=I ／(∥．Ⅲ)Ⅲc‘…。，他，口 =i…}^f矽、m)Ⅲ】c，“c／cⅨ／O 戈(3)可以分成二步来完成，首先今