拉冬变换.ppt.Convertor.docVIP

拉东（Radon）变换 Radon变换是数学上的一种变换。 本章先介绍拉东（Radon）变换（RT）及其反变换（IRT）; 介绍以Radon变换为理论基础的CT（Computer Tomography）技术，以及它们在地球物理勘探中的一些应用。 Radon变换原理 一个二维函数的Radon变换（RT）是由沿多条直线的积分所组成的，这个运算同医学、结晶学、地球物理学等诸多领域中的大量物理实验相对应，一般将它们称之为投影试验。 在一些地震学的问题中也会涉及到RT或它在曲线积分上的推广，这里我们将讨论层析成像在地震学中的各种应用。从层析的意义看，沿着射线路径传播的信号（至少在高频极限上可以这样近似地讲）累加起来构成了模型的某些性质——如慢度或慢度异常、衰减等等，当多道射线路径从许多方向上穿经了该模型时，就可以提供出足以重建出该模型的信息。 Radon变换原理 若 是二维位置矢量，那么函数 的Radon变换定义为: 式中直线 所对应的方程表达式为: 其中 分别为直角坐标和极坐标的表达。 Radon变换原理 直射线投影中的几何分布与变量。直角坐标 和极坐标 定义了模型中的位置矢量 投影角 决定了 坐标系。投影试验是在S 坐标的方向上沿着线L 对模型做积分，它给出了 坐标系中的投影数据。 CT技术---扇形束投影示意图 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域中直接运算时间是非常大的，为了降低运算时间，可以将t-x域中求逆转换到F-X域中。 在F-X域拉当变换对为： （2.5） 即我们把原来的t-x域数据通过拉当变换转化到τ-ρ域去噪的方法，改为把F-X域地震数据通过拉当变换转化到f-ρ域去噪。 Radon变换原理 2.3 参数的确定 参数包含斜率参数ρ和斜率参数的采样间隔Δρ。 1、采样间隔Δρ的选择 设 是原始数据中的最高有效频率， 和 分别是最大炮检距和最小炮检距。 Radon变换原理 F-X域拉当变换中斜率参数ρ采样率应该满足的条件为： 2.12 Radon变换原理 由此我们可以知道F-X域拉当变换中斜率参数ρ的临界采样率Δρc为： 2.13