ct原理及发展历程

1、第一讲 概论目录第一讲 概论1一、CT简介21.1 什么是CT21.2 CT设备的构成31.3 CT分类及部分性能指标41.4 CT成像的特点51.5 CT应用领域61.6 CT发展简史91.7 检测成像工程中心简史13二、X射线成像物理学202.1 X射线的产生202.2 X射线与物质相互作用21三、CT成像的数学模型243.1 Beer定律243.2 CT成像的连续数学模型253.3 CT成像的离散数学模型27四、附录284.1 CT图书284.2 CT杂志294.3 CT 会议314.4 在线资源32一、CT简介1.1 什么是CTCT是Computerized Tomography或Co

2、mputed Tomography的缩写。“C”就是“计算机化的”或是“计算的”。我们看“T”或“Tomography”。Tomography=tomos(slice, section)+graphein(draw)所以，“Tomography”就是对断层进行绘图（成像）。Tomography是内视技术，即不打开物体就能看到物体内部的结构。具体的讲：l 目标：希望得到被测物体内部物理量B的分布。l 信息获取：用探测器探测一组物理量A。l 信息处理：根据物理量A与物理量B之间的物理关系A=F(B)，求解出物理量B的分布。CT就是计算机断层成像。广义上讲，以下技术均属于计算机断层成像技术：l X射

3、线（射线）断层成像技术（即狭义CT）l 激光CT成像（Optical Computer Tomography）l 核磁共振断层成像技术（Magnetic Resonance Imaging, MRI）l 超声成像技术（Ultrasonic Tomography, UT）l 电阻抗成像技术（Electrical Impedance Tomography）l 红外成像技术（Infrared Imaging）l 而本课程主要介绍X射线（射线）断层成像技术（即狭义CT，XCT）。即便是狭义的XCT，根据射线产生方式、获得信息内容、所关心物理过程等的不同，又分为l 透射式CT （Transmission

4、 CT）u 吸收衬度CT成像（Attenuation contrast CT, 即通常CT）u 相位衬度CT成像（Phase contrast CT）l 发射式CT（Emission CT）：PET，SPECTl 散射式CT（Scattering CT）本课程仅介绍吸收衬度CT。1.2 CT设备的构成前面提到，CT成像涉及“信息获取”和“信息处理”两个过程，根据需要，无论是工业CT或是医用CT，都包含以下组成部分：l 射线发生器l 射线探测器l 机械及控制系统l 计算机（控制/存储/处理）1.3 CT分类及部分性能指标可根据不同准则对CT进行分类：l 纳米焦点CT：对象10-30um，分辨率3

5、0-160nml 微焦点CT：对象1-30mm，分辨率0.5-5uml 医学CT：对象(人)300-800mm，分辨率0.3mm-3mml 工业CT：对象(构件)100-500mm，分辨率0.1-0.5mml 大型工业CT：构件300-1000mm，分辨率0.5-1.0mml 超大型工业CT：构件1000-2000mm，分辨率1.0mm左右l 工程CT：桥梁、隧道等，分辨率1米-几米l 地学CTl 天体CT著名的医学影像厂商：GPS（GE，Philips，Siemens），Toshiba著名的工业CT厂商：YXLON，GE，SkyScan，ZEISS(XRadia)，岛津1.4 CT成像的特点

6、l 非接触、无损成像l 无影像重叠l 密度和空间分辨率高l 数字图像，易于处理、存储、传输1.5 CT应用领域 l 医疗生物（医疗诊断、小动物成像）l 工业（国防：导弹、炮弹、飞机发动机叶片、火箭推进器、气缸.；民用：机械部件、电子器件、陶瓷、）l 材料（分子结构、岩心、管道、）l 安全检查（非金属武器、毒品、炸药 ）l 化石成像（古动植物三维结构成像）l CT原理也用于：天体成像、地学成像（地震、地质构造）、工程成像（桥梁、隧道、水坝、物探）编织材料泡沫材料芯片检测电子封装岩心分析火工品离体牙种子1.6 CT发展简史Wilhelm Conrad Rntgen (27 March 1845 1

7、0 February 1923) was a German physicist, who, on 8 November 1895, produced and detected electromagnetic radiation in a wavelength range known as X-rays or Rntgen rays, an achievement that earned him the first Nobel Prize in Physics in 1901. In honour of his accomplishments, in 2004 the International

8、 Union of Pure and Applied Chemistry (IUPAC) named element 111, roentgenium, a radioactive element with multiple unstable isotopes, after him.This 121-year-old picture of fingers is one of the first images ever made with x-rays. The hand belonged to Anna Bertha, wife of German physicist Wilhelm Rntg

9、en, the discover of x-rays. The black glob on the fourth finger is a ring made of gold, which absorbs x-rays.Johann Karl August Radon (16 December 1887 25 May 1956) was an Austrian mathematician.Sir Godfrey Hounsfield and the EMI header scannerThe first clinical CT scan on a patient took place on 1s

10、t October 1971 at Atkinson Morleys Hospital, in London, England. The patient, a lady with a suspected frontal lobe tumour, was scanned with a prototype scanner, developed by Godfrey Hounsfield and his team at EMI Central Research Laboratories in Hayes, west London. The scanner produced an image with

11、 an 80 x 80 matrix, taking about 5 minutes for each scan, with a similar time required to process the image data.国内情况1980s初期开始CT理论与应用研究较早开始CT研究的单位：上海交通大学（庄天戈）北京信息工程学院（邱佩璋）清华大学（安继刚、王经谨、张朝宗）东南大学，东北大学，重庆大学目前活跃的研究单位：清华大学（康克军、陈志强、李元景、刘以农、张丽、邢宇翔、赵自然等，威视），东北大学（郑全录、江根苗，东软），重庆大学（王珏、曾理，ICT中心），中北大学（韩炎、潘晋孝、桂志国等

12、），北京大学（姜明、杨建生、周铁），首都师范大学，中科院高能物理所（魏龙、姜晓明、朱佩平、单保慈、魏存峰等），上海交通大学（赵俊），西安交通大学（牟轩沁），北京交通大学（渠刚荣），北京信息科技大学（邱钧），大连理工大学（孙怡），北京航空航天大学（杨民、傅健），西北工业大学。还有一些搞PET、光CT的单位：中科院高能所、华中科技大学（骆清铭、谢庆国）等。1.7 检测成像工程中心简史l 2003年，张朋研究员调入首都师范大学，建设检测成像实验室l 2007年，完成实验室多功能锥束/扇束工业CT系统的搭建，国内率先实现GPU加速图像重建，对外做了大量的服务工作l 2009年，自主研发成功显微CT设备

13、，分辨率为10-20uml 2010年，获批检测成像北京高校工程中心l 2012年，与东营三英工程中心联合成功研制研发光耦合显微CT，分辨率到达亚微米l 2013年，成立天津三英精密仪器有限公司l 2013年，与中科天悦公司联合开发口腔CTl 2015年，自主研发成功国内首台双能谱显微CT设备，分辨率达到亚微米 二、X射线成像物理学2.1 X射线的产生X射线本质上是一种电磁波，波长介于0.001纳米到10纳米之间，X射线的能量与波长成反比。X射线通常由高速电子轰击高原子序数的金属靶（如钨、钼等）而产生。Error! Reference source not found.是由X光管产生X射线的示

14、意图。图 1 X光管产生X射线示意图高速运动的电子打到靶上（如钨、钼等），受靶中原子核库仑力的作用，速度骤减而发出X射线，称此种X 射线为轫致辐射。以这种方式产生的X 射线光子的能量分布从接近零的能量到电子的最高能量。若电子的能量足够大，还可以把靶面材料原子中的内层电子撞击出去，使原子处于不稳定的激发状态；为使原子恢复至稳定的低能态，邻近层的高能态电子立即自发地填补其空穴，该电子出发地能级扣除空穴处能级后多余的能量以辐射形式释放，即特征X射线。如上所述，加速电子轰击金属靶而产生的X射线由不同能量（或波长）的轫致辐射光子和特征光子组成。组成X射线的各个能量的光子个数的概率分布曲线称为X射线谱。E

15、rror! Reference source not found.是电压为100kV时电子以夹角13度轰击钨靶产生的X射线谱。图 2 X射线能谱示意图2.2 X射线与物质相互作用X射线穿过物质时，会与物质发生作用而衰减。在通常CT系统所用的X射线能量范围内，X射线与物质的相互作用主要有三种形式：光电效应（Photoelectric effect）、康普顿散射（Compton scatter）和电子对效应（Electron pair effect）。三种作用所占的比例与X射线的能量和被照射物体的原子序数有关，如Error! Reference source not found.所示。下面分别介绍

16、这三种作用的物理机制。图 3 X射线与物质作用中光电效应、康普顿散射和电子对效应起主导作用的范围1）光电效应：射线光子与成像介质的原子（作为整体）相互作用时，整个光子被原子吸收，其所有能量传递给原子中的一个电子（多发生于内层电子）。该电子获得能量后离开原子而被发射出来，称为光电子。光电子的能量等于入射X射线光子的能量减去电子的结合能。图 4 光电效应示意图2）康普顿散射：X射线光子与成像介质原子的外层电子（可视为自由电子）发生非弹性碰撞，光子只将部分能量传递给原子外层电子，使该电子脱离核的束缚从原子中射出，并按一定的概率分布产生一个与入射光子方向间夹角为的散射光子。散射角分布由0到180度。散

17、射光子还会以一定概率发生二次或多次散射。图 5 康普顿散射示意图3）电子对效应：能量大于1.02MeV的X射线光子从成像介质的原子核旁经过时，在原子核的库仑场作用下，光子转变成一个电子和一个正电子。所发射出的电子还能继续与介质产生激发、电离等作用；正电子与介质中的电子相结合湮没而产生两个能量为0.51MeV的X射线光子。图 6 电子对效应示意图X射线光子与物质的作用是随机的。通常我们用衰减系数（包括线性衰减系数和质量衰减系数）来刻画物质对X射线光子的阻挡本领。线性衰减系数表示光子穿过单位路径与物质发生作用的概率。量纲是长度的倒数（常用）。质量衰减系数是指每单位质量厚度光子和介质发生作用的概率。

18、常用单位是。质量衰减系数和线性衰减系数间的关系如下：其中为物质的密度。注意，通常物质的线性衰减系数和质量衰减系数还跟X射线光子的能量有关。三、CT成像的数学模型3.1 Beer定律按照线性衰减系数的定义，可以推导出著名的Beer定律。图 7 单色X射线穿过物体衰减示意图用能量为的单能（也叫单色）X射线照射材质均匀的物质。设物质对能量为的X射线的线性衰减系数为。X射线由一端射入物体，记初始强度为。设X射线穿过物质厚度后的强度为，如Error! Reference source not found.所示。则按照线性衰减系数的定义有经过简单的整理，有式中令，可得到如下微分方程求解该方程有由初始条件可

19、以计算出式中的。于是有式称为Beer定律。该定律是光吸收的基本定律，适用于所有的电磁辐射和所有的吸光物质。对于非均匀非单一物质，其在每一点的衰减系数与位置有关，记为。式可改写为若X射线不是单色的，其归一化能谱记为，同时衰减系数与X射线能量是有关的，记为。则上式进一步化为3.2 CT成像的连续数学模型在理想条件下，即假设X射线由单一能量的光子组成，忽略射线源焦点及探测器尺寸，考虑X射线穿过不均匀物体发生衰减。设物体的线性衰减系数分布为。以和表示X射线沿路径穿过物体前和穿过物体后的强度，则由Beer定律有上式两端同时除以，再经过对数变换和简单的整理有称为沿路径的投影。CT成像问题就是利用X射线探测

20、器探测沿不同路径的和，进而由这些值重建被测物体线性衰减系数分布的值或近似值。或者表述为：由一系列重建的值或近似值。图 8 Radon变换示意图考虑二维情况。平面内任意一条射线路径（直线），可以由原点到的有向距离和的法方向与轴的逆时针夹角表示，如Error! Reference source not found.所示。利用这两个变量，的方程可以表示为：其中。因此又可以写做式在数学上称为函数的二维Radon变换。下文中对的Radon变换也记为。固定角度，变化时，表示一系列平行射线，简称为平行束（parallel beam），称为平行束投影。于是，当射线用参数表示时，CT成像问题就是“由一组平行束投

21、影重建的问题”。3.3 CT成像的离散数学模型图 9 CT成像的离散模型以下以断层以下以断层CT成像为例介绍离散模型。如Error! Reference source not found.所示，设在圆域之外为零。将包含圆域的正方形区域等分成个小区域，每个小区域称为一个像素。设为的特征函数或示性函数，即在内为1，在外为0。令其中为的面积。称为的数字化图像。显然。设为与图像相交的射线，在沿射线的投影为式中即第条射线与第个像素的交线长。在离散CT问题中是未知量。习惯上以代替，记。于是得方程该式可以写成矩阵形式其中，是待求的的数字化图像的图像向量；是实测数据的投影向量；是投影矩阵。所谓离散CT问题就是

22、：已知投影矩阵和投影向量，求式的广义解。离散CT模型的求解方法可分为直接方法和迭代方法。直接方法有消元法、LU分解、奇异值分解方法等，其主要缺点是所需的存储量和计算量过大，例如图像像素数为1024x1024，投影角度数为720，每个角度的采样数为1024时，的行数为1024x1024，列数为720x1024。即使仅存储的非零元也需要约5.6GB。此外，直接方法还会遇到由于数据误差造成方程不相容的问题。因此，就现有计算机能力，直接方法很难求解离散CT模型。四、附录4.1 CT图书1) Herman G T. Fundamentals of computerized tomography: ima

23、ge reconstruction from projectionsM. Springer Science & Business Media, 2009.2) Kak A C, Slaney M. Principles of computerized tomographic imagingM. Siam, 1988.3) Hsieh J. Computed tomography: principles, design, artifacts, and recent advancesC. Bellingham, WA: SPIE, 2009.4) Natterer F. Mathematical

24、methods in image reconstructionM. Siam, 2001.5) Zeng G L. Medical Image ReconstructionM. Heidelberg: Springer, 2010.6) Buzug T M. Computed tomography: from photon statistics to modern cone-beam CTM. Springer Science & Business Media, 2008.7) Scherzer O. Handbook of Mathematical Methods in Imaging: V

25、ol. 1M. Springer Science & Business Media, 2011.8) 庄天戈. CT 原理与算法M. 上海交通大学出版社, 1992.9) 张朝宗, 郭志平, 张朋, 王贤刚. 工业 CT 技术和原理M. 科学出版社, 2009.10) 闫镔,李磊. CT图像重建算法M. 科学出版社, 2014.4.2 CT杂志1) IEEE Transactions on Medical Imaging2) International Journal of Biomedical Imaging3) Journal of X-Ray Science and Technology

26、4) Computerized Medical Imaging & Graphics5) Applied Radiology6) Medical Physics7) Academic Radiology8) Journal of Nuclear Medicine Technology9) Biomedical Engineering Online10) Journal of Biomedical Imaging11) IEEE Signal Processing Letters12) Medical imaging13) Journal of Computers and Mathematics

27、14) Journal of CT Theory and Applications15) Journal of Mathematical Analysis and Applications16) Physics in Medicine and Biology17) IEEE Transactions on Biomedical Engineering18) IEEE Transactions on Image Processing19) Biomedical Physics20) Journal of CT Theory and Applications21) Computed Tomography Theory and Application22) Electronic Imaging23) Applied Opt