TJΘ(??u Ρ)=G Θ 3(??u Ρ????u Ρ)=G Θ 3(??u Ρ??u Ρ) 该构造张量反映的局部方位是沿第一本征向量的方向, 即相干性最强的方向, 一般决定于连续同相轴控 000100-10 3103 ------------------------------------------------------------------------------------------------ —————————————————————————————————————— 30-(11) 32 30-3 的一阶导数算子, 该模板的近似旋转不变性要比 Sobel 滤波中的通用模板好得多。理想的一阶一维ζζ导数滤波变换函数由 iΠk 给出(k 为归一化的波数)。 6 54 石油地球物理勘探 2004 年制的信息部分。如果数据中存在边缘, 由于边缘特征在构造张量中往往不沿连续同相轴控制的方位, 因此在扩散过程中很容易被摧毁。为了保持尖峰边缘, 就应当设法在构造张量中增大边缘信息的贡献量。注意到尖峰处二阶导数值往往是局部极大值, 而在其他较平滑的位置, 二阶导数值往往较小, 因此可以在构造张量中加入二阶导数的贡献, 这样更有利于保持边缘。构造张量改写为 TTJΘ(??u Ρ)=G Θ 3??u Ρ??u Ρ+(u Ρ)xx(u Ρ)xx+ (12)+(u Ρ)yy(u Ρ)Ty 由式(12) 可知, 在构造张量中加入二阶导数的贡献后, 对于原先平滑的数据位置影响不大, 而对于边缘位置则增加了贡献量。因此新构造张量确定的方位对边缘有很好的保护, 使其不至于在扩散中被平滑掉。构造张量经过改造之后, 算法的计算流程基本不变, 的计算。导算子。干扰或人为等因素( 采集、处理手段), 常常会发生时断时续的现象, 而非地层特征引起。滤波平滑的目的就是要使时断时续的同相轴连续
全文预览
