1、第九章第九章 图象重建图象重建9.1 概述9.2 基本原理9.3 滤波逆投影法9.4 代数投影变换法9.5 图象重建技术的应用第九章第九章 图象重建图象重建9.1 9.1 概述概述 由物体截面投影来重建该截面图象是近年来发展起来并获得广泛应用的图象处理技术。图象重建最典型的应用是医学上的计算断层摄影技术(CT)。它用于人体头部、腹部等内部器官的无损伤诊断,其基本方法就是根据人体截面投影,经过计算机处理来重建截面图象。计算机发展后才出现的一个分支,不同于传统处理与图形 X 光透视图像 CT断层扫描图像 投影几何对规则形状已有系统办法,圆柱(认为内部是均匀的)视图增加,即无限多视图可以解决任意三维
2、物体原图。Compute Aided Tomography(CAT)正视图 侧视图问题:能否从投影中恢复原图问题:能否从投影中恢复原图?答复是肯定的。yxstcossinsincos 一条射线沿S方向穿透物体,投影轴与X轴夹角为,建立s、t坐标系,(t,s)与(x,y)关系如下式:xtysP(t)X射线 dsststfyxfdsyxftpcossin,sincos),(,角位置之平行投影在是射线沿射线积分组成投影:物理上X射线到人体有个衰减过程:u(x,y)为x,y点的衰减 Nin :入射X射线(光子)强度Nd :X射线穿透物体后被检测到的射线强度u(x,y):反映了人体各部组织的性质,在空间
3、上的分布就形成了人体 各部组织的图象,所以u(x,y)实质上反映了图象灰度分布f(x,y)射dsyxuNNind,exp射得:dsyxuNNdin,ln 上述为一个切面,一系列切面构成三维物体,两种方法:直接、间接。dxdyvyuxjyxfvuF2exp,9.2 9.2 基本原理:基本原理:对图像函数f(x,y)付氏变换:uwysxtdtwtjtpwStP,02exp时当的付氏变换 dxuxjxpuS2exp00当 v=0 时,dxuxjdyyxfdxdyuxjyxfuF2exp,2exp,0,wSuSdxuxjxpuF0002exp0,即故有:此结果推广到一般情况下:此结果推广到一般情况下:
4、每一个下P(t)付氏变换后对应于F(u,v)在对应下的剖面值。只要有足够多的P(t)n个对应的S(u),就是足够多个F(u,v)的剖面,近而可以逼近F(u,v),反变换即可求得f(x,y)。xyuv由付氏变换旋转不变性:得:得:S S (w)=F(w,(w)=F(w,)=F(u,v)=F(u,v)(一般的S(w)=F(u,v)的证明)证:证:f(t,s)是f(x,y)在t,s坐标上为函数 dsstftpyxst,cossinsincosxtysuv yxstyxsyxtdtedsstfdtetpwSwtjwtj,cossinsincos,22由代入这里 ,)sin vcos u (,2sinc
5、os2wFvuFwSdxdyeyxfwwdxdyeyxfwSvyuxjyxwj 则:在约束条件下:实现流程:极坐标 直角坐标 yxfvuFwFtp,9.3 9.3 滤波滤波逆投影法逆投影法 极坐标F的付氏反变换:dwwdyxwjwFwdwdwtjwFyxtwdwdwtjwFyxf sincos2exp,2exp,sincos2exp,0000200这里:F(,)vu,wFwF对极坐标:vu ddwewwSdwdwtjwwFddwwwFwdwwFdwdwwFwdwwFwtj 0200000002exp,)(exp,exp,exp,exp,dyxQyxfsincos,0即:dwewwStQwtj2
6、)(令下面的关键问题是如何求得:下面的关键问题是如何求得:)(tQ之贡献。对下在某一离散形式:yxQyxQkyxfkiii,sincos,1 wwSFdwwtjwwStQ12exp下面讨论:tPFdtetPwStgtfsGsFFwtj21注意:dthPtQwFth111则:若令 这里:P(t)已知(可实验测出),只要求得h1(t)即可求Q(t)和f(x,y)了。|w|w|看作另一函数的付氏变换,其反变换是什么?看作另一函数的付氏变换,其反变换是什么?有如下几种近似办法:有如下几种近似办法:1、若取|w|1,即相当加入了|1/w|滤波,(故|1/w|w|=1)|1/w|加入低通滤波,图象变模糊了
7、。2、找一个函数其结果近似|w|:用G(w)=|w|R(w)R(w)=1|w|N例:对128128图象,150个投影,每个投影150条射线。aij 有15021282=3.68108个实际上方程无法解,一般用迭代法逐次近似求解。,修正通过加大,修正通过减少非正确解时,当给定初值fgfafgfaffiTiiTi000 iTiiiTiaaagfaff001 反复迭代即可求解。9.5 9.5 图象重建技术的应用图象重建技术的应用 主要在放射医疗、工业检测设备中,显示人体(对象)各个部位断层图象。历史:历史:理论源于1917年,奥地利数学家(Radon)所发表的论文。证明了二维或三维物体能够通过其无限
8、多个投影来确定,但限于当时技术条件没能实现。60年代初,计算机技术发展,重建受到重视,不少学者做了卓有成效的研究。其中:英国EMI公司中央研究所工程师Housfield,经四年努力,在1972年研制成诊断头颅用的第一台计算机X线断层摄影装置,这一新设备在1974年5月蒙特利尔(Montreal)召开第一次国际CT会议上,正式命名为“电子计算机断层摄影技术”简称CT。1979年EMI公司又研制出全身CT。1979年这项技术获诺贝尔奖。几代产品(按扫描方式的改进)几代产品(按扫描方式的改进):第一代:第一代:单束扫描方式,X射线管与探侧器同步水平直线运动,时间长,一周需4分钟左右。第二代:第二代:窄角扇束扫描,张角1020度,2030个探测器相配合。扫描时间20秒左右。第三代:第三代:广角扇束扫描方式,张角为30度左右,探测器增加到250350个,射线源和探测器同步旋转扫描,扫描时间可缩短2.5秒。第四代第四代:三代基础上发展,检测器1500个左右,布满整个360。固定不动,X线源旋转扫描。约2秒。第五代第五代:28射线管,0.01秒,1秒内重复60次。目前拓展、超声CT、放射性同位素正电子CT、质子CT。CT其它领域:电子天文学、电子显微镜。