实验报告Ⅰ-CT..docxVIP

计算机断层成像实验摘要： 计算机断层成像（Computed Tomography,简称CT）是计算机技术、电子技术、数字化图像重建技术和核技术等相结合的产物。它是在无损状态下获得被检测物体断层的灰度图像，以其灰度来分辨被检测断面内部的几何结构、材质情况、缺陷种类等，且清晰、准确的展示出来。本文就重点就CT的数学、物理基础以及其图像重建原理等进行阐述与说明。这其中就包括一些目前常用的图像重建算法。关键词：断层成像 CT图像 图像重建引言：1895年，Wilhelm C. Roentgen 发现了 X 射线，并为夫人拍下了世界上第一张 X 片 —— 戴戒指的手掌照片。1967年，Godfrey N. Hounsfield 发明了第一台 CT 设备，能够从多个角度摄片，采集被摄物体的三维信息，在不破坏物体的情况下观察其内部结构。1970年代，医院开始使用CT诊断疾病。数十年来，这一伟大技术已经广泛应用于各种领域，例如医学（组织器官、生理代谢过程成像）、药学（药效检测、新药开发）、材料学（新材料的开发）、工业（各种器件的质检和探伤）、农业（木材和种子的质检和分析）、工程（建筑材料内部孔隙度、连通度和渗透性分析）、珠宝（真伪识别和最佳切割方案设计）、考古（化石的结构和成分分析）等领域。实验目的了解工业CT机的基本结构以及组成部分理解计算机断层成像的基本原理掌握CT教学实验仪的使用方法实验原理CT的物理原理和数学基础CT的物理原理是基于射线与物体的相互作用的原理。一束射线穿过物质并与物质相互作用后，射线会因吸收或散射而减弱。如图1。减弱的程度与物质的厚度和组成成分（或吸收系数）有关。 图1 射线穿过衰减系数为U长度为t的物体当入射强度为I0的单色射线穿过厚度为l，线性衰减系数为μ的物体时，探测器检测到的射线强度I服从郎伯-比尔定理（Lambert-Beers）确定，即：。对于一般的物体，设f（）为物体被测物体某断面的线性衰减系数的分布函数，当强度为的单色射线穿过物体时，探测器检测到的射线强度，两边取对数，并令，称为射线穿过物体的投影值。可见射线穿过物体的投影值是被检物体某断面的线性衰减系数分布函数沿射线路径的线积分值。所谓图像重建就是如何由各个方向不同位置的投影值来得到线性衰减系数分布函数。CT图像重建原理(1) 重建的概念断层图像重建是以扫描检测所得的衰减系数线积分数据集为基础，经过必要数据校正，按一定的图像重建算法，得到衰减系数具体的二维分布，再将其以灰度形式显示，从而获得可在计算机屏幕上观察的断层图像。图像重建的关键是重建算法，既要考虑图像质量，又要注意运算速度。重建算法多种多样，各有特色，主要有两类：变换法及级数展开法。 (2) 图像重建算法原理A、反投影法 反投影法是第一个用于医学CT的成像算法，该方法是将射线穿过被检测工件断面所得到的投影值，沿原方向“反投影”到重建区域内沿射线路径上的每一点，并将经过各点的射线值加起来得到某像素点的值。这就是反投影重建法。该算法实现简单，但从它的思想原理可以看出起潜在的局限就是没有消除噪声的能力，同时伴有严重的点扩散效应。 将射线穿过断层所检测到的数据称为投影，儿把射线路径对应于图像上所有像素点赋以相同的投影值则称为反投影，反投影以灰度表示将形成一个图案或图形。对断层各个方向上投影完成反投影并形成相应的图形或图案，将多有的反投影图形或图案叠加，则得到反投影法重建的断层图像。B、卷积反投射原理CT成像原理在不少刊物上有详细介绍。所谓CT图像重建是按照采集后的数据，求解图像矩阵中像素，然后重新构造图像的过程；而图像矩阵的求解由计算机完成，在CT图像重建的解析法中，当前最常用的是采用卷积运算的滤波反投影法，也称卷积反射影法。卷积反投影重建图像时，先把由检测器上获得的原始数据与一个滤波函数进行了卷积运算，得到各方向卷积的投影函数；然后再把它们从各方向进行反投影，即按其原路径平均分配到每一矩阵元上，进行叠加后得到每一矩阵元的CT值；再经过适当处理后就可以得到被扫描物体的断层图像，卷积反投影可消除单纯的反投影产生的边缘失锐效应，补偿投影中的高频成分和降低投影中心密度，并保证重建图像边缘清晰和内部分布均匀。卷积反投影法的数学表达式可由雷当投影定理推导出；下面直接给出表达式。如图1所示，取物体某断层面为坐标系，射线管发出的单能窄束射线束穿过衰减系数为的物体断层，沿路径线束的投影可写为：　式中是线管输出强度，是检测器测得的线束强度。　图1　被测物体断面　 由上式可知，投影值就是衰减系数在断层平面的二维积分，若测出投影值，可求出该断层平面上衰减系数的二维分布。在X线束的平移-旋转扫描中线的投影值总是与线束路径S有关的，在图1中，引入一个新坐标系R-θ来描述线束路径S在x-y坐标系中的位置；取R为线束路径S到坐