生物医学工程基础 生医图象2-林江莉.ppt

1、2020/8/13,四川大学材料学院生物医学工程系,1,生物医学图象(Biomedical Image),生物医学工程基础（四） X射线断层成象 林江莉 （2007.11）,Imaging technologies are changing the way science is done (R.P. Crease, Science, Vol. 261, July 1993),2020/8/13,四川大学材料学院生物医学工程系,2,第一部分：医学成像,1.1 医学成像过程的组成部分及成像系统模型 1.2 基本概念 1.3 主要的成像系统及其应用 1.4 评价图像系统与图像质量 1.5 X射线断层

2、成象 1.6 核磁共振成象 1.7 超声医学成象,2020/8/13,四川大学材料学院生物医学工程系,3,1.5 X射线断层成象,2020/8/13,四川大学材料学院生物医学工程系,4,The History,Johan Radon (1917) showed how a reconstruction from projections was possible. Cormack (1963,1964) introduced Fourier transforms into the reconstruction algorithms. Hounsfield (1972) invented the

3、X-ray Computer scanner for medical work (for which Cormack and Hounsfield shared a Nobel prize). EMI Ltd (1971) announced development of the EMI scanner which combined X-ray measurements and sophisticated algorithms solved by digital computers.,2020/8/13,四川大学材料学院生物医学工程系,5,X-CT的发展历程,1917 奥地利数学家Radon：

4、从投影重建图像的原理 1920-60 美国物理学家A.M.Cormack：可行性证明、试验研究、仿真 1972 G.N.Housfield：装置和断层成像结果,2020/8/13,四川大学材料学院生物医学工程系,6,X-CT成像的演变,1901年(投影） 1975年（断层） 2000年（断层）,2020/8/13,四川大学材料学院生物医学工程系,7,第一幅脑CT图象,2020/8/13,四川大学材料学院生物医学工程系,8,一.投影X线成像技术,1. X线的产生 X线管 加热的灯丝发射出电子 30200千伏高压作用 灯丝射出电子 吸引到阳极靶子上 电子与靶内的原子相互作用产生X射线(少于1%的入

5、射电子能量转换成了X射线，旋转阳极) X射线穿过管壁发射出来,2020/8/13,四川大学材料学院生物医学工程系,9,2.X线的性质,(a)穿透作用: 波长短，能量大，能穿透一般光线不能穿透的物质，0.50.001A，25Kev1Mev。 (b)荧光作用: 当X线照射某些物质时可产生荧光，利用这一性质，可以在荧光屏上直接观察X线图像。 (c)电离作用:具有足够能量的X线光子不仅能击脱物质原子轨道上的电子，使该物质产生一次电离，而且脱离原子的电子又与其它原子相碰，还会产生二次电离。可根据气体分子电离电荷的多少来测定X线的剂量。 (d)生物效应:X线是一种电离辐射。生物细胞经一定量的X线照射后会受

6、到损害甚至坏死。放射治疗破坏肿瘤组织；也会导致一定正常组织的损伤,2020/8/13,四川大学材料学院生物医学工程系,10,3. X线与人体组织的作用,当X射线穿过物体时，能量 物质的吸收，转换成其它形式的能量 散射，改变了原来的传播方向 透射，穿过被探查物体沿原方向继续向前传播，称为透射分量。,2020/8/13,四川大学材料学院生物医学工程系,11,最早的X线成像方法是靠投影成像 投影成像分为： 透视： 荧光透视成像系统 摄影： 胶片摄影系统,二、投影X线成像系统,2020/8/13,四川大学材料学院生物医学工程系,12,第一代的荧光透视接收器,平板荧光屏 优点：观察组织形态、位置外，还可

7、以观察脏器的运动 缺点：屏的亮度比较低，观察起来比较吃力。放射科医生一般要在黑暗环境中待15分钟左右才能使自己的眼睛适应，可观察到的信息也较少）,2020/8/13,四川大学材料学院生物医学工程系,13,采用影像增强管 影像增强管的引入是透视X射线成像系统的一项重大改进。,解决方案,使亮度的总增益达到5000倍左右,2020/8/13,四川大学材料学院生物医学工程系,14,现代的投影X射线成像设备,采用影像增强管-电视系统 包含影像增强管、光学图像分配系统和一个包括摄像机、监视器的闭路视频系统及辅助电子设备 视频系统的引入是荧光透视成像系统的又一项重大改进,2020/8/13,四川大学材料学院

8、生物医学工程系,15,用摄影胶片代替透视的荧光屏 X射线在胶片上形成潜影，经过显影、定影，将影像固定在胶片上 临床中使用屏-胶片系统：它是由涂上感光乳胶的胶片和与胶片紧密接触的一个或两个荧光增强屏组成的。 X射线的能量由增强屏吸收，并将大约520%转变为光线，使胶片曝光。这样胶片曝光所需的实际X射线辐射剂量大幅度地降低。但产生一定的模糊 数字X射线摄影,胶片摄影系统,2020/8/13,四川大学材料学院生物医学工程系,16,数字电视X射线摄影系统,2020/8/13,四川大学材料学院生物医学工程系,17,重叠：三维物体结构投影在一个二维平面上 办法：设法获取某一指定断面的图像，或者是人为突出人

9、体中某种结构（如血管、骨骼）的形象。,投影X线成像的问题？,2020/8/13,四川大学材料学院生物医学工程系,18,三、 数字减影技术DSA,最早应用于血管系统的研究 在病人的血管中注入造影剂（造影剂对X线的衰减系数大于人体的骨骼和软组织） 摄下同一部位注入造影剂前后的两帧图像。原则上讲，只要将这两帧图像相减，就能消除血管外其余的组织形象，而只保留造影血管的图像,2020/8/13,四川大学材料学院生物医学工程系,19,四、运动断层摄影,为了获取某一深度的图像 只要X线源与检测器按指定规律运动，就可能使聚焦平面中A、B点的影像在整个移动过程中始终落在胶片的固定位置上，而其它深度处的结构的影像

10、则从胶片的一端移动到另一端。最后图像将突出聚焦平面所在深度上的断面结构。,2020/8/13,四川大学材料学院生物医学工程系,20,五. X线计算机断层摄影（X-CT）,X射线发现后的七八十年，始终没有解决两个根本性的问题 一是：器官重叠，读片困难； 二是：固有的分辨率差 投影X射线成像系统对软组织的分辨能力比较差。,XCT 从根本上克服了上述困难，是80多年来X射线诊断学上的一次重大突破。,2020/8/13,四川大学材料学院生物医学工程系,21,X射线被准直后成为一条很窄的射线束。当X射线管沿一个方向平移时，与之相对应的检测器也跟着作平移运动。这样，射线束就对整个感兴趣的平面进行了一次扫描

11、，检测器接收到了与脏器衰减系数直接相关的投影数据。,1.基本原理,检测器作同步平移运动，并旋转进行扫描来获得投影数据,2020/8/13,四川大学材料学院生物医学工程系,22,平移旋转扫描,2020/8/13,四川大学材料学院生物医学工程系,23,假设用一种单一能量的X线照射厚度为d的物体，其入射射线强度为Io，透射后的强度为I，则有： 称为线性衰减系数，表示特定能量的X射线照射某种特定的物质时单位距离上的衰减分数。,衰减系数 VS 投影数据,2020/8/13,四川大学材料学院生物医学工程系,24,式中d为像素的宽度。,2020/8/13,四川大学材料学院生物医学工程系,25,1)能鉴别出较

12、小的衰减系数差，从而提高了对软组织的诊断能力。从CT图像上能识别出与周围组织的衰减系数只差0.5%的病灶。 2)能精确地测定出组织的X射线衰减系数值，从而对组织性质做出判断。 3)可做出人体任意部位的断面图像，并利用计算机图像处理技术构造出人体结构的三维图像。,2. XCT成像的突出优点,2020/8/13,四川大学材料学院生物医学工程系,26,根据投影数据计算出原始图像 迭代法和反投影法 (1)迭代法:用逐次近似法来求解联立方程组 目的：寻找二维分布度函数f(x,y) 具体做法为：先假设一个最初的密度分布(如假设所有各点的值为0)，根据这个假设得出相应的投影数据，然后与实测到的数据进行比较。

13、,3.图像重建技术,2020/8/13,四川大学材料学院生物医学工程系,27,如果不符，就根据所使用的迭代程序进行修正，得出一个修正后的分布。这就是第一次迭代过程。以后，就可以把前一次迭代的结果作为初始值，进行下一次迭代。在进行了一定次数的迭代后，如果认为所得结果已足够准确，则图像重建过程就到此结束。,2020/8/13,四川大学材料学院生物医学工程系,28,右图给出了一幅由四个象素组成的图像，若四个象素的值分别为5、7、6、2。则可以分别获得6个投影数据，包括两个水平方向，两个垂直方向和两个对角线方向，分别是7、11、9、13、12和8。这就是所能得到的所有的已知数据。,例：代数重建技术,2

14、020/8/13,四川大学材料学院生物医学工程系,29,迭代开始，先令所有的重建单元的值为0，第一步计算出垂直方向的投影值，分别都是0，如右上图。把这个计算值与实测值11和9相比较后，将其差值除以2以后分别加到相应的单元上去，就可得到垂直方向的迭代结果，如右下图。,2020/8/13,四川大学材料学院生物医学工程系,30,在此基础上可以再做水平方向的迭代，此时有计算值10、10，实测值为12、8，将它们比较后求出差值，分别加到有关的象素上去，结果是：,2020/8/13,四川大学材料学院生物医学工程系,31,最后再作对角线方向的迭代，就得到了所要求的真实数据。实际上，要求重建的矩阵很大，因此选

15、代算法是非常耗费时间的。,2020/8/13,四川大学材料学院生物医学工程系,32,基本做法： 把每次测得的投影数据原路反投影到射线经过的各个象素上。也就是说，指定投影线上所有各象素点的值等于所测得的投影值。,(2)直接反投影法,2020/8/13,四川大学材料学院生物医学工程系,33,2020/8/13,四川大学材料学院生物医学工程系,34,将测得的投影数据与一个“核函数”作卷积运算，然后再用所得的结果作反投影。 选择“核函数” 序列-1/3，1，-1/3，用这个序列与测量到的投影数据0，1，0做卷积运算后可得到一个新的序列-1/3，1，-1/3，用这个序列数据作反投影,(3)卷积反投影法,

16、2020/8/13,四川大学材料学院生物医学工程系,35,卷积：设序列x(n),h(n),它们的卷积和y(n)定义为：,n: -1， 0， 1 h(n):-1/3，1，-1/3 x(n) 0， 1， 0 y(n): -1/3，1，-1/3,2020/8/13,四川大学材料学院生物医学工程系,36,2020/8/13,四川大学材料学院生物医学工程系,37,(4)傅立叶滤波反投影法,Radon 变换 中心面片理论 滤波反投影法,2020/8/13,四川大学材料学院生物医学工程系,38,Fourier Transform,傅里叶1768年生于法国 Fourier的两个最主要的贡献 1807年提出 “

17、周期信号都可表示为成谐波关系的正弦信号的加权和”傅里叶的第一个主要论点 “非周期信号都可用正弦信号的加权积分表示”傅里叶的第二个主要论点,2020/8/13,四川大学材料学院生物医学工程系,39,周期函数的Fourier级数,直流分量 n =0,基波分量 n =1,谐波分量 n1,2020/8/13,四川大学材料学院生物医学工程系,40,直流系数,余弦分量 系数,正弦分量 系数,周期函数的Fourier级数的系数,2020/8/13,四川大学材料学院生物医学工程系,41,1D Fourier Transform,定义：若函数f（x）满足狄里赫赖条件： 1）具有有限个间断点； 2）具有有限个极值

18、点； 3）绝对可积 则把变换称为：,Fourier正变换：,Fourier反变换：,2020/8/13,四川大学材料学院生物医学工程系,42,1D Discrete Fourier Transform（1D DFT）,Fourier正变换：,Fourier反变换：,求反变换：,2020/8/13,四川大学材料学院生物医学工程系,43,2D Fourier Transform,2020/8/13,四川大学材料学院生物医学工程系,44,2D DFT,2020/8/13,四川大学材料学院生物医学工程系,45,FT性质,对称性和叠加性 奇偶虚实性 尺度变换特性 时移特性和频移特性 微分和积分特性 卷积

19、定理 Paseval定理,2020/8/13,四川大学材料学院生物医学工程系,46,2020/8/13,四川大学材料学院生物医学工程系,47,Fourier变换,2020/8/13,四川大学材料学院生物医学工程系,48,谱的表示,2020/8/13,四川大学材料学院生物医学工程系,49,The function is known as the Radon transform of the function f(x,y).,Radon TransformLine Integrals and Projections,2020/8/13,四川大学材料学院生物医学工程系,50,Fourier Slic

20、e Theorem,2020/8/13,四川大学材料学院生物医学工程系,51,Fourier Reconstruction,x,y,g(s,),(x,y),wx,wy,s,1D-Fourier transform Fg(s),2020/8/13,四川大学材料学院生物医学工程系,52,Fourier Reconstruction,1D Fourier transform Fg(s,) of one projection,2D Inverse Fourier transform,2020/8/13,四川大学材料学院生物医学工程系,53,Fourier Reconstruction,x,y,g(s,

21、1,2),(x,y),wx,wy,s,1D-Fourier transform Fg(s),2020/8/13,四川大学材料学院生物医学工程系,54,Fourier Reconstruction,1D Fourier transforms Fg(s,1,2) of two projections,2D Inverse Fourier transform,2020/8/13,四川大学材料学院生物医学工程系,55,Fourier Reconstruction,x,y,wx,wy,s,g(s,1,2,3,4),1D-Fourier transform Fg(s),(x,y),2020/8/13,四川

22、大学材料学院生物医学工程系,56,Fourier Reconstruction,1D Fourier transforms Fg(s,1,2,3,4) of four projections,2D Inverse Fourier transform,2020/8/13,四川大学材料学院生物医学工程系,57,Fourier Reconstruction,x,y,wx,wy,s,g(s,1,2,8),1D-Fourier transform Fg(s),(x,y),2020/8/13,四川大学材料学院生物医学工程系,58,Fourier Reconstruction,1D Fourier tran

23、sforms Fg(s,1,2,8) of eight projections,2D Inverse Fourier transform,2020/8/13,四川大学材料学院生物医学工程系,59,Influence of # of Projections,64 projections,16 projections,8 projections,1 projection,2 projections,4 projections,2020/8/13,四川大学材料学院生物医学工程系,60,例：反投影,One Projection,2020/8/13,四川大学材料学院生物医学工程系,61,例：反投影（完全

24、：未滤波）,2020/8/13,四川大学材料学院生物医学工程系,62,滤波反投影,Filtered Sinogram,Sinogram,2020/8/13,四川大学材料学院生物医学工程系,63,|w| - Filters,Shepp-Logan,Hamming,Lowpass Cosine,RAM - LAK,2020/8/13,四川大学材料学院生物医学工程系,64,例：滤波反投影,Filtered Sinogram,Reconstructed Image,2020/8/13,四川大学材料学院生物医学工程系,65,discrete implementation:,Filtered Backpr

25、ojection,2020/8/13,四川大学材料学院生物医学工程系,66,图像重建是X-CT中的核心技术，也是众多科技人员研究的热门课题。它是确保获得高质量图像的关键。,2020/8/13,四川大学材料学院生物医学工程系,67,五代CT,2020/8/13,四川大学材料学院生物医学工程系,68,First Generation,EMI Mark I (Hounsfield), “pencil beam” or parallel-beam scanner (highly collimated source) excellent scatter rejection, now outdated 1

26、80 - 240 rotation angle in steps of 1 Used for the head 5-min scan time, 20-min reconstruction Original resolution: 80 80 pixels (ea. 3 3 mm2), 13-mm slice,2020/8/13,四川大学材料学院生物医学工程系,69,Second Generation,Hybrid system: Fan beam, linear detector array (30 detectors) Translation and rotation, however R

27、educed number of view angles scan time 30 s Slightly more complicated reconstruction algorithms because of fan-beam projection,2020/8/13,四川大学材料学院生物医学工程系,70,Third Generation,Wide fan beam covers entire object 500-700 detectors (ionization chamber or scintillation detector) No translation required sca

28、n time seconds (reduced dose, motion artifacts) Reconstruction time seconds Pulsed source (reduces heat load, radiation dose),2020/8/13,四川大学材料学院生物医学工程系,71,Fourth Generation,Stationary detector ring (600 4800 scintillation detectors) Rotating x-ray tube (inside or outside detector ring) Scan time, re

29、construction time seconds Source either inside detector ring or outside (rocking, nutating detectors),2020/8/13,四川大学材料学院生物医学工程系,72,Comparison of 3rd and 4th Generation,Both designs currently employed, neither can be considered superior 3rd Generation (GE, Siemens): Fewer detectors (better match, cheaper) Good scatter rejectio