首页 [重点]线性代数方法建模6 CT图像重建--数学建模案例分析

[重点]线性代数方法建模6 CT图像重建--数学建模案例分析

[重点]线性代数方法建模6 CT图像重建--数学建模案例分析[重点]线性代数方法建模6 CT图像重建--数学建模案例分析 ?6 CT图像重建 CT是Computed Tomography的简称，即计算机断层成像技术，也称为计算机辅助断层扫描CAT—Scanner。为什么通过CT扫描能够比较清楚地了解被扫描物体断层的组织结构呢,它与数学又有什么样的联系, 拍X光片是将三维对象(立体)显示在二维的胶片或荧光屏上，待检测物体与胶片平行，X射线垂直投射到胶片上，这样，在深度方向的信息重叠在一起，混淆不清。另外，由于胶片的密度分辩力低，不能区分软组织的细节，只能区分密度差别大的内...

[重点]线性代数方法建模6 CT图像重建--数学建模案例分析 ?6 CT图像重建 CT是Computed Tomography的简称，即计算机断层成像技术，也称为计算机辅助断层扫描CAT—Scanner。为什么通过CT扫描能够比较清楚地了解被扫描物体断层的组织结构呢,它与数学又有什么样的联系, 拍X光片是将三维对象(立体)显示在二维的胶片或荧光屏上，待检测物体与胶片平行，X射线垂直投射到胶片上，这样，在深度方向的信息重叠在一起，混淆不清。另外，由于胶片的密度分辩力低，不能区分软组织的细节，只能区分密度差别大的内脏器官，影响了诊断的效力。CT的创立，解决了这个问题。它不同于传统的X射线，它的X射线束则位于待检测物体的横截面内，X射线源发射出极细的笔束X射线，在其对面放置一检测器，测量出X射线源发出的射线的强度，以及经过物体衰减后达到检测器的X射线强度，然后，将X射线源与检测器在观测平面II0 内不断同步改变位置(平移或旋转)，得到关于X射线强度I的若干组数据(可以是几万组甚至几十万组)。如果物体是均匀的，物体对X射线的衰减系数为常数。设强度为的射线在物体中I0行进距离I后衰减至，由Beer定理: x I,,,,x0IIex (1)(或,ln),,,,0I,, 但若物体在待检测的平面内是不均匀的，则。此时X射线在某一方向沿某一路径xy,,,(x,y) L的总衰减可以用线积分表示: I,,0dl (2),ln,,,,LI,, 称(2)为射线投影。若未指明路径，只指明方向，即，称为投影，投影是一组射线投影,dl, I的集合。(1)与(2)中的I和可由实验测得。 0 II 我们的任务是:根据测得的一系列和来求。求得的仅是一个离散的,,,(x,y),(x,y)0 二元函数，通过转换成CT数(这是相对于水的线性衰减系数，它比易于分辩)，以数模,(x,y)转换器转换成图像信号，由电视屏以不同灰度等级或色彩显示出来或拍摄下来。这就是投影重建图像，所建立的图像为扫描平面的物体断层图像。 CT是G.N.Hounsfield在1971年提出并在英国EMI Ltd.生产了第一台CT，1979年他与A.M.Cormack一起获得了Nobel医学奖。CT一出现，就得到了广泛的应用，不仅在医学方面，在工业无损检测、农林业、生态环境检测、地球物理等方面都得到了广泛的应用，现已由第一代发展到第五代。我们这里讨论的是如何求,(x,y)的问题。投影重建图像有许多方法，如反投影重建算法、滤波重建算法和迭代重建算法等。这里仅介绍迭代重建算法。 , 将待检测物体的截面分成许多边长为的小正方形，每个小正方形称为一个像元。设一束宽 ,度为的射线，平行于像元的边，整个穿边像元。组成X射线的光子以一定的速率被像元内的组织吸收，其速率正比于组织的X射线密度(线性衰减系数)。记第个像元的X射线密度为，xjj定义进入第j个像元的光子数目x,lnj离开第j个像元的光子数目 ,,ln(穿过第j个像元而未被吸收的光子的百分比) kii如果第束X射线穿过由个像元组成的一行像元，那么第束X射线经过这行像元未被吸收的 i百分比=第束X射线经过待检测物体的截面而未被吸收的百分比。记第i束射线未经过待检测截面进入检测器的光子数目b,lni 第i束射线经过待检测截面进入检测器的光子数目 ,,ln(第i束射线经过待检测截面后未被吸收的光子的百分比) i称为第束X射线的密度，从而 (3)x，x，？，x,b12ki 式中可由测量得到，是未知的。 bx,x,？x12ki 2NiN,n 将待检测的截面分成个像元，它们的标号由1到，第束X射线穿过个像元组n成的一行像元，这些像元的标号为，则(3)可写成 j,j,？j12n x，x，？，x,b (4)jjji12n 1,若j,j,j,？j,12n,a令 ,ij0,其它, 则(4)可写成线性方程组 ax，ax，？，ax,b,i,1,2,？,M (5)i11i22iNNi iM其中为射线的数目，称(5)为第束X射线的方程，矩阵称为投影矩阵。，，A,aijM，N a 然而，一束射线不一定沿平行于像元的方向穿过像元，因而需要将(5)中的加以修正，ij常用的方法如下: (1)像元中心法: 1,若第i束射线经过第j个像元的中心,a, (6),ij0,其它, (2)中心线法: 第束射线的中心线位于第个像元内的长度ijl (7),,aij,第个像元内的长度j (3)面积法: 第i束射线位于第j个像元内的面积 (8)a,ij第i束射线若平行于像元边穿过第j个像元时位于第j个像元内的面积 N 若总共有束射线，个像元，利用上述三种方法中的任何一种选择，则可得到一个含MaijN有个未知数，M个方程的线性方程组: N (9)ax,b,i,1,2,？,M,ijji,1j NNM一般地，可等于，也可以大于或小于。在实际的CT设备中，等，N,80，80;320，320 等。方程组(9)可能有解，也可能无解;有解时，可能有无穷多组解。客M,28800;184300 观上，由实际问题的意义，(9)应该有解。一般地，应该有无穷多组解。由于模型是由实际问题经过若干近似抽象而来的，数据又是测量得到的，不可避免地有误差，因此也会出现无解的情形。所以，我们将方程组(9)写成: Ax，e,b (10) A在式中，CT中用迭代法求解。它的思想是:先建立一个判据，即接受的标准，选择一个解的初始值，为误差向量。另外，由于矩阵的元素很多，又有许多元素为零，用消元法来解很费时，若它符合判据，则接受为解，若达不到要求，则按一定的方法对原值加以修正，通过不断的迭代，直到可以接受为解为止。e

                    本文档为【[重点]线性代数方法建模6 CT图像重建--数学建模案例分析】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
                    图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。 
 该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。

                    [版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。

                    本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。

                    网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
                

下载需要：免费已有0 人下载

立即下载

[重点]线性代数方法建模6 CT图像重建--数学建模案例分析

你可能还喜欢