l基于几何迭代ART的断层成像问题.docxVIP

基于几何迭代的断层成像问题 李金臣 邵春雨 2013.7前言ART算法可以追溯到柏拉图从投到墙上的影子重构问题，一直到柯马克的计算X射线断层成像问题。所有的问题都是要求从某个中重构物体。这一个投影可以看作是物体的切片的平均。在本文中我们介绍ART算法并且提出断层成像相关的正问题和反问题,同时讨论所提问题的性质。摘要本文用ART算法研究断层成像问题,讨论断层成像的解收敛性，并提出矛盾方程和不定方程收敛值的特点，同时讨论迭代到什么程度认为是所提问题的解。讨论不同情况下，上机的实验结果。关键词：ART算法，收敛性，超定解，不定解正文ART算法简介ART算法使用连续正交投影来近似方程的解。下面考虑k个未知数的单个方程 当k=2时，这一方程的解集为一直线，当k=3时，它是一个平面，而一般情况下解集叫做超平面。这一超平面由向量v的系数和标量u来决定，它可以表示为H={x}，其中（. , .）是通常的内积。向量v是超平面的法线，因为对于所有的x，z，（v，x-z）。对于给定x，y关于H的正交投影因此为唯一的向量Py标量t必须满足但是Py=y-tv，所以Py=y+因此，超平面的投影是容易实现的，特别是如果向量v有相对少的非零元时。由m个线性方程k个未知数构成的方程组的解属于m个超平面的交集，其中每个超平面由一个方程决定。而ART算法是通过连续投影到这m个超平面上来近似求解的。把这个方法应用于y，一个周期后得到新的 其中重复这个周期，给出进一步的近似式ART算法收敛性的证明及迭代步数的选取(1)射影定理：记H={x}，其中（. , .）是通常的内积，向量v是超平面的法线，对于给定x，y关于H的正交投影因此为唯一的向量Py，且满足其中 是欧几里得范数。（k）（2）ART算法收敛性的证明ART算法是通过连续投影到这m个超平面上来近似求解的，并设起精确解为x 。把这个方法应用于y，一个周期后得到新的逼近值 其中重复这个周期，给出进一步的近似式，并且 a为某个常数证明：由射影定理得 …………………….由以上m个式子即 ……………………………………………………………… {11}有界由数列单调有界定理得的极限存在，并且 a为某个常数如果当给定方程为不定方程，a=0；当给定方程为矛盾方程时，a非零。（3）.给定阈值由收敛性得当 我们取，ART算法在断层成像的应用及其正问题和反问题我们将ART算法应用到一个很简单的射线摄影术问题中来。考虑物体为n排列的“像素”。我们假定物体在每个像素内是一样的，但像素与像素之间可能存在差异，我们的目标是得到关于像素之间变化情况的图片。假设物体被通过像素中心的放射线所照射。我们从第一行开始从左到右排列这些像素，光线的观察向量为 V=[0，1，0，1，0，0，0，0，0]。假设材料的每个像素吸收射入它的放射线的一部分。这个部分（0,1）区间上的一个数，接近于0的值说明像素是透明的，而接近1的值表明像素是基本不透明的。我们称第j个像素吸收部分为该像素的吸收系数。假设一束光线的观察向量为v=[]，穿过物体的吸收系数为 []。如果光线穿过第j个像素，那么放射线的部分被吸收，而1-部分传递出去了。从物体中穿越出去放射线p为 P=，其中v为光线的观察向量，或等价地， 注意到每个ln（1-），lnp为负数，所以我们可以把这一关系写成 其中 为正数。因此这一X射线断层摄影术问题等价于寻找由m个观察向量和相关的右端项确定的m个超平面的交集中的一个点。此外，对于某一特定观察，其右端项是射线穿过的像素对应的给定m个观察向量和相应的值我们把这些观察值作为观察矩阵的行向量，其第ij个元素并且把与这些观察值相应的测量值作为m维行向量mu。给定维行向量表示初始近似值。正问题：给定观察向量v和物体n*n像素的内部结构及吸收程度,求相应的u值；反问题：给定观察向量v和相应的u，求物体n*n像素的内部结构及吸收程度；4.上机实验 （1）一条虫子在一个划分为5网格的苹果上打洞。穿过行（从上到下）的水平探测测量值分别为0,2,1,2,0，而沿列（从左到右）垂直探测测量值分别为0,1,3,1,0.主对角线（西北到东南）之和为1，而另一条对角线（西南到东北）之和为3.使用零向量为初始近似，调用程序“art1”和“displa”计算虫子的图像。迭代n=5 000.117500．0560 0.00700．06190.90450.96360．007000.01310.96070.013100.00700.96360.90450.06190.00700.056000.117500迭代n=10000.041100．0011 0.00600．00720.96101.00210．00600 00.9936000.00601.00210.96100.00720.00600