图象恢复重建for 地质工程物.pptVIP

图象恢复与重建简介 图象恢复：也称图象复原。是图象处理的一大类技术。 和图象增强的相同之处：为了得到（某种意义上）改进的图象，或者说是改进图象的视觉质量。不同之处：图象恢复处理退化或恶化（品质下降）的图象。需要根据相应的退化模型和知识重建或恢复原始的图象。 图象恢复技术是要将图象退化的过程模型化，并据此采取相反的过程以得到原始的图象。它要根据一定的图象退化模型来进行。 例：(老书上P103页 　　　　　 新书上P185页，) 退化模型示例。 图象重建 投影重建：指从一个物体的多个（轴向）投影图重建目标图象的过程。 从投影重建图象是一类特殊的图象处理方法。 通过投影重建可以直接看到原来被投影物体某些特性的空间分布。 X-CT中的图象重建方法很多，常见的有以下几种： 1.联立方程法 2.逐次近似法（掌握） 3.反投影法 4.卷积法（卷积反投影法） 5. 傅立叶变换法 代数重建法 ＡＲＴ ＳＩＲＴ 逐次近似法　　计算实例： 设Ｘ线从６个方位射入，测得结果如图所示，求各单元值？ 逐次近似法 原则：　（原值－新值）／２＝修正值 　　　　　　前次值＋修正值　＝近似值 ①先设各单元初始值均为零 ②选择方向，填入平均值（选Ｃ，Ｄ方向，则C/2=5.5,D/2=4.5) ③选取另一方向，计算修正值.[ 选Ａ，Ｂ，计算修正值(12-10)/2=1, (8-10)/2=-1 ] ④再选取另外方向进行修正（选Ｅ，Ｆ，修正值 　(7-10)/2=-1.5 , (13-10)/2=1.5 计算图例过程： 傅立叶变换投影（定理） （要求会叙述,证明傅立叶变换投影定理：） 可表述为：图像　f (x,y)　在与Ｘ轴成θ角的直线上投影的傅立叶变换是　f (x,y)　的傅立叶变换在朝向角θ上的一个截面。 　或：f (x,y)是以θ角进行投影的傅立叶等于　f (x,y)的傅立叶变换在傅立叶空间（Ｒ， θ）处的值。 投影重建示意图 上图中： 发射源到接收器直线：（由①它与原点的距离s ②与Ｙ轴的夹角θ 确定）。 沿直线（s, θ)对f(x,y)的积分为： 定理的数学描述为： 设Ｇ（R, θ)是g(s, θ)对应第一个变量s的１－Ｄ傅立叶变换，即： 证明过程如下： 证： 方法２ 证明：不失一般性的，可取θ＝０ 　　设图像f(x,y)在Ｘ轴上的投影为gy（x),则有： 众所周知，医学CT技术（Computerized Tomography）是获得巨大成功的一种反演技术。 其成功的原因： ● 采用了不同的探测方法 成像清晰 ● 全方位的数据采集 以 X-CT采集为例 CT技术在医学上发展很快，应用成熟，主要是由于有好的探测手段（如X线等），全方位的数据采集方式，而且人体组织成分相对单一，因此现有的成像技术已经能取得令人满意的图像。 地震CT的与医学CT相比，两个难点： ⑴ 数据不全问题（有限角观测） ⑵ 非线性问题（即积分与路径有关） 地震层析成像是在曲线簇上的积分（它是广义Radon变换），因此无论从理论上来说，还是从实际应用来看，CT技术研究者有这样一种说法：“医学CT是昨天的CT，地震CT是明天CT”。它正在逐步完善，走向成熟。 地震层析成像的广泛应用 90年以来，渗入到应用地球物理界各个部门 ● 在工程地质和灾害地质中 三峡工程，两侧地质构造、工程危体 ● 在地球科学研究中 了解地幔结构、研究板块运动 ● 在石油的开发中 开发地震学（层析技术、四维地震） 进行油气储集层描述和动态监测 层析成像的基本理论 ●层析成像是根据物体外部的测量数据，依照一定的物理和数学关系反演物体内部物理量的分布，最后又以图象的形式表现。 这就是说，人们可以在不切开物体的情况下，对其内部结构成像。对物体进行逐层剖析成像，相当于把物体切成片，然后照相。它是在非破坏性的条件下，通过测定不同入射角的射线(X射线、γ射线、地震波、电磁波)等的走时、衰减或阻抗值——它们分别是待测物体中慢度、吸收系数或者密度对射线路径的线积分(称为投影函数)。 层析成像是成像技术发展到目前的最高阶段。-----由投影重建图像。 总之，层析成像不仅利用波作为探测手段，而且要同时利用各种不同方向投射的波场信息。 ● 拉当变换----层析成像的理论基础： 1917年，数学家拉当证明：已知所有入射角?的