1、CT的图像重建的的图像重建的数学模型*MRI图像处理的数学模型*咖玛刀图像处理的数学模型3、1989：单方向连续旋转型CT机螺旋式扫描CT第五代CT机1、CT：Computer Tomography计算机体层摄影，计算机断层成像2、1972：Allan M.Cormack (USA)Godfrey N.Hounsfield (UK)1979:Nobel 奖一、简介：二、基本原理和概念?X射线在穿过物质时其强度按指数形式衰减，因此X射线穿过均匀物质后的强度 Iout与入射强度Iin的关系为loutinIIe?其中 为X射线在均匀物质中的传播距离，为物质对X射线的衰减系数?lexp(d)outin
2、LIIl?因此对于不均匀物质，其出射强度和入射强度的关系为lllinoutneeeII?21X和入射强度的关系为的模块后，其出射强度，而衰减系数不同射线穿过一组长度相同当但如果用入射强度相同的射线穿过两个不同的模块时，其出射强度也可能相同，尽管两个物质的内部结构不同。水水该物质值某物质的?1000CT人们关心的是各组织密度间的差异。此时就需要把整个模块分成若干小的 基本单元。这些小单元称为体素(Voxel)，我们要求的就是这些小单元的CT值，它反映了物质的密度。CTCT值反映物质的密度，物质的值越高，物质的密度越高10000CT?水水水水值010001000CT?水空气水值三、三、CTCT的结
3、构与原理的结构与原理平面平面CTCT的成像原理的成像原理体素、矩阵和象素体素、矩阵和象素体素：将选定层面分成若干个体积相同的立方体数字矩阵：每个体素的 X线衰减系数排列成矩阵像素：一幅CT的图象由许多按矩阵排列的小单元组成，这些组成图象的基本单元称为像素在CT中，用专门的计算机将收集到的原始数据经过复查的运算而得到一个显示数据的矩阵，由原始数据一显示数据过程称为重建，CT的本质就是重建图象，CT的重建图象克服了常规 X射线设备线积分测量的局限性1917:拉东（Radon）变换:()(,)dLP fF x y sLs?L其中 为平面中的一条直线，表示沿L的弧长Q0QdF1F()fP f?L（q）
4、（Q）q其中Q为平面中任意一点，（q）表示沿所有与点Q相距q（q0)的直线L的的平均值.将待测物体的截面分成若干个体素，每个体素为边长为的正方形，设有一束射线的宽度为四、模型的建立?X射线的强度以一定的速率被体素内的组织吸收而衰减，衰减速率正比于该体素的衰减系数jjXx记第 个体素的射线密度为线性方程个立束射线，我们就可以建如果有，则射线离开时的强度束第射线进入时的强度束第个体素，记束射线穿过连续如果第射线的强度个体素的离开第射线的强度个体素的进入第射线密度为个体素的定义第MMbxxxkixxXjikijj?21XiXilnbXjXjln1121212121 1221,9.91,0,(9.9)
5、(9.10)NiikiikjikjjjikijiiiNNiNNikjjjxikjjjxxxbjjjjaa xa xa xbM?LLLLL将待检测的截面分成个体素，并将它从 到 编号（如图）设第 束穿过 个体素，这些体素的编号为则设第 束穿过 个体素，这些体素的编号为则（）当令其他所以式可写为若总共有束射线，个体素1,1,2,.Nijjija xbiM?L，最后得到如下方程组X1X2X3X4CT图像重建1234这样CT图像重建的数学模型为：1,1,2,.,NijjijNa xbiMAx bxRMN?L也就是线性方程组：这里因此这是有个方程，个变量的线性方程组X射线不一定沿平行于体素边穿过，因而对
6、aij有如下三种修正。(1)981,0(2)99(3)9 10ijijiijijaijijlajijAaijA?中心法（图）当第 束射线经过第 个体素中心时，当第 束射线不经过第 个体素中心时中心线法（图）第 束射线的中心线经过第 个体素内的长度第 个体素的边长面积法（图）第 束射线位于第 个体素内的面积第 束射线若平行于体素边长时位于第 个体素内的面积五、模型的求解矩阵有许多的元素是零，是大型稀疏矩阵，我们采用迭代法求解，有雅可比(Jacobi)迭代法，高斯-赛德而(Guass-Seidel)迭代法,超松驰迭代法等。1、Jacobi 迭代法迭代法:NMijaA?)(11()(1,2,0,1,
7、2,)kkiiijjj iiixba xin ka?LL111()kkkAx bADL UxDL U xD bBxf?2、Guass-Seidel 迭代法111111()(1,2,0,1,2,)inkkkiiijjijjjj iiixba xa xain k?LL11111()()kkkkkkAx bADL UDxLxUxbxDLxDLbGxf?3、超松驰迭代法：11111()(1,2,0,1,2,)inkkkkiiiijjijjjj ixxba xa xin k?LL1()kkkkkkAx bCIAxCxbxbAxxr?1()(1,2,0,1,2,)kkkiiiijjj ixxba xin
8、k?LL11111()(1,2,0,1,2,)inkkkkiiiijjijjjj ixxba xa xin k?LL1(1,2,)nijijjiiiiabxinaa?L考虑方程组：11111()(1,2,0,1,2,)inijijkkkkiiijjjj iiiiiiiaabxxxxaaain k?LL11111()(1,2,0,1,2,)021:inijijkkkkiiijjjj iiiiiiiaabxxxxaaain kGauss Seidel?LL迭代公式六、含有测量误差的处理六、含有测量误差的处理1、无解的情形：12(1,2,)minnnijjijx Ra xbinAx bAx bAx b?L即无解此时处理方式为：2、无穷解的情形：1(1,2,)min.nijjijTa xbinAx bx xAx bstAx b?L无穷多解此时处理方式为：实践与思考题：P114，习题1
