生物医学工程博士论文答辩.pptVIP

浙江大学博士学位论文答辩 答辩人： 徐茂林 导师：汪元美 教授 浙江大学生物医学工程与仪器科学学院 2004.3.11 下午2：00 论文题目 透射及散射断层成像中若干反演算法的研究 A New Approach to Reconstruction Algorithms for Transmission and Scattering Tomography 说明： 本课题承蒙国家自然科学基金（No.资助。 文中所用的实测数据由中国计量科学研究院提供。 下图是中国计量科学研究院研制的工业CT机。 论文主要内容框架图： CT的数学基础： 记被测物体对X射线的线性衰减系数为： 反演算法 变换法 用给定的一系列解析式（算子或函数变换）反演公式(A)来重建图像 。常用Randon求逆公式的“离散化”形式。 级数展开法 直接对公式(A) “离散化”，得到一大型稀疏线性方程组： Radon变换示意图： 一.滤波反投影算法的改进 及其它应用 新滤波函数及其时频特性分析 各滤波函数主瓣、旁瓣的幅值和宽度 讨论： 滤波函数对重建误差的影响主要有以下三个因素：主瓣，临近旁瓣及远处旁瓣。主瓣越高越窄，空间分辨率越好；旁瓣越小，数值精度越高，通常密度分辨率就越好。从时域特性图中及上表中可以看出，新滤波函数的主瓣高而窄，尤其是远处旁瓣几乎为零，能够在保证空间分辨率的同时，相对地提高密度分辨率。 3.数字实验结果： 模拟数据 实测数据： 局部重建 1.新滤波函数用于局部重建的特性 时域曲线局部区域放大比较图： 2.局部重建原理 选择适当的滤波函数，使得滤波函数的旁瓣迅速衰减，譬如在某一有限区间外很小，以至于几乎为零，则可以通过如下的公式作重建： 3.局部重建实验结果 误差估计 结论： CBP算法是CT中标准的方法，关于它的理论和算法已有很多研究。同时，人们认为CBP没有局部性，似乎不可能直接用来作局部重建。本文证明，无论从理论或数据实验来看，直接用CBP实现图像局部重建是可能的。所谓直接用CBP是指，就用CBP的计算公式、算法，以至于计算程序。但是，要恰当选取其中的滤波函数。本文选取一种新的滤波函数，它的旁瓣衰减很快，经证明在可能选取的滤波函数中，它的衰减是最快的。 滤波反投影算法的其它应用 二.级数展开法的若干研究 经典ART和SIRT方法的迭代格式 对称块迭代算法 重建误差及时间比较： 4.实测数据 三.背散射成像的多准则神经网络优化重建 Compton散射成像原理 四.总结与展望 由投影重建图像是CT成像技术的理论基础，本文立足于二维由投影重建图像理论，对透射及散射断层成像问题，从成像系统模型的角度，抽象出其数学的共同点，作了统一的研究和分类。主要是： 1．提出了一种新的滤波函数，其主要部分是泛函分析中的磨光算子，它有许多优异的性质，特别是利用旁瓣的迅速衰减性，用滤波反投影算法简便快速地进行了图像局部重建，重建图像空间分辨率高，且能较好地克服Gibbs效应。 对称块迭代 局部重建 由Radon变换似可进一步揭示一般区域上的Green函数的形式和性质 Compton背散射成像 三维CT技术： 螺旋CT、核磁共振（MRT）、三维正电子成像（3D PET）、三维单电子反射成像（SPECT）等。??? 本文所提出的算法思想推广至三维成像，是否有相应的理论和结果，同时是否会产生新的理论或实验结果，这些都有待于展开进一步的研究工作。 Thanks 是第i扫描行第k个能级的光子计数率， 是探测器第k个能级的能量。 表示第(i，j)个象素的电子密度。N是被测截面图象在一维方向的象素个数(设分成N×N的网格)； 是探测器在第i扫描行所能分辩的能级数。G是一个常数，它是由系统的几何形状和探测器的效率决定的。 是Kronecker符号，表示第(i，j)个象素对第k个能级有贡献时为1，否则为0。 是第(i，j)个象素的有效散射体积。 是探测器对第(i，j)个象素所张的立体角；若 以表示探测器表面在散射方向的投影值， 表示探测器到第(i，j)个象素的距离，则 ＝ / 。 现将第i扫描行的电子密度及能谱分别记为向量 及 ，则可写为如下矩阵形式： 若记 ， ，即 是由 ( )构成