拉冬变换数学基础PPT.ppt

模型测试与分析 对这个地震数据进行F-X域Radon变换，变换后的信号在f-ρ域如图六所示。 图六 原始信号F-X域Radon变换 在图六中我们可以看到图中存在一个脉冲，由于在x-t域共炮点道集是有限的，做Radon变换会引起畸变端点效应，即能够看到端点发散效应，变换到f-ρ域是一个能量团，这与理论是一致的。 模型测试与分析 我们把f-ρ域的信号返回到F-X域，信号谱图如图七所示，经过反傅立叶变换我们将Radon变换后的F-X信号变换到t-x域，信号图如图八所示。 图七 f-ρ域映射到F-X域的频谱图 图八 F-X域Radon变换后返回到t-x域 的信号图 模型测试与分析 为了更加明确的看到Radon变换前后信号是否相近，我们取Radon变换前后t-x域中的第5道数据进行比较，如图九所示。 图九 F-X域Radon前后t-x域信号比较 拉东（Radon）变换 Radon变换是数学上的一种变换。 本章先介绍拉东（Radon）变换（RT）及其反变换（IRT）; 介绍以Radon变换为理论基础的CT（Computer Tomography）技术，以及它们在地球物理勘探中的一些应用。 Radon变换原理 一个二维函数的Radon变换（RT）是由沿多条直线的积分所组成的，这个运算同医学、结晶学、地球物理学等诸多领域中的大量物理实验相对应，一般将它们称之为投影试验。 在一些地震学的问题中也会涉及到RT或它在曲线积分上的推广，这里我们将讨论层析成像在地震学中的各种应用。从层析的意义看，沿着射线路径传播的信号（至少在高频极限上可以这样近似地讲）累加起来构成了模型的某些性质——如慢度或慢度异常、衰减等等，当多道射线路径从许多方向上穿经了该模型时，就可以提供出足以重建出该模型的信息。 Radon变换原理 若 是二维位置矢量，那么函数 的Radon变换定义为: 式中直线 所对应的方程表达式为: 其中 分别为直角坐标和极坐标的表达。 Radon变换原理 直射线投影中的几何分布与变量。直角坐标 和极坐标 定义了模型中的位置矢量 投影角 决定了 坐标系。投影试验是在S 坐标的方向上沿着线L 对模型做积分，它给出了 坐标系中的投影数据。 Radon变换原理 2.1 拉当变换的数学原理 自1917年Radon先生提出这个变换以后,拉当变换在医学、物理学、天文学等许多领域都已得到了广泛的应用。 设函数y=g(x)连续可导,而且其反函数是单值的，d(x,t)满足可积,则定义： (2.1) 为拉当正变换的连续公式, (2.2) 为拉当反变换的连续公式。 Radon变换原理 其中是d(x,t)拉当反变换的结果， ， 称为Hilbert算符。根据g(x)的不同,可以把Radon变换分为线性拉东变换和非线性拉当变换：①如果g(x)=x，则我们定义的拉当变换就是线性Radon变换，既τ-ρ变换，该变换把t-x域中的一条直线映射成τ-ρ域中的一个点；②如果g(x)是其它的非线性函数，则我们定义的拉当变换就是非线性Radon变换,或称为广义Radon变换。 Radon变换原理 拉当变换有明显的物理意义，它是将时间、空间域t-x的一条直线t=τ+ρx映射到τ-ρ域上的一个点，如图1所示。 （a）t-x域一条直线 （b）由t-x域映射到τ-ρ域中的一个点 图1 t-x域一条直线与τ-ρ域中一个点的关系 Radon变换原理 在二维连续空间-时间域的Radon正反变换对： 2.3 Radon变换原理 在计算机实现中，由于在时间域和空间域的离散采样，不能应用连续函数方程，因此用离散的累加来代替连续域的积分运算；为了消除离散采样的有限孔径的影响，利用最小平方法计算离τ-ρ变换。二维离散时间、空间域的拉当正反变换对： 2.4 Radon变换原理 2.2 F_X域拉当变换的数学原理 由于在t-x域中直接运算时间是非常大的