• / 20
  • 下载费用:10 金币  

Surfer中网格化方法的选取探究.ppt

关 键 词:
Surfer中网格化方法的选取探究.ppt
资源描述:
Surfer中网格化方法的选取探究曹创华2010年 11月• 绘制地球物理数据等值线图时,需要根据客观环境特征和数据本身的特点,选择合适的网格化方法。• Surfer中的插值方法有 12种,常用的有加权反距离插值法、克里格法、最小曲率法、最近邻点法、多项式回归法、径向基函数法、带线性插值的三角剖分法等,本文将结合实例将常用网格化方法的选取方法、适用范围及参数设置等使用技巧做了简单的探讨。前言• 网格化方法的特征及应用条件– 网格化 概念 —是指通过一定的插值方法,将稀疏的、不规则分布的数据插值加密为规则分布的数据,以适合绘图的需要。– 网格化方法 —加权反距离法 ( Inverse distance to a power)、克里格法( Kriging)、最小曲率法( Minimum curvature)、改进谢别德法( Modified shepard’s method)、自然邻点法( Natural neighbor)、最近邻点法(Nearest neighbor)、多项式回归法( Polynomial regression)、径向基函数法( radial basis function)、带线性插值的三角剖分法( triangulation/liner interpolation )、移动平均法( moving average )、数据度量法( data metrics )和局部多项式方法( local polynomial)。• 克里格法( Kriging)– 最初是由南非金矿地质学家克里格根据南非金矿的具体情况提出的计算矿产储量的方法:按照样品与待估块段的相对空间位置和相关程度来计算块段品位及储量,并使估计误差为最小。后来,法国学者马特隆对克里格法进行了详细的研究,使之公式化和合理化。– 克里格方法的基本原理是根据相邻变量的值(如若干样品元素含量值),利用变差函数所揭示的区域化变量的内在联系来估计空间变量数值。– 该方法总是尽可能地去描述原数据所隐含的趋势特征,以区域化变量理论为基础,以变差函数为主要工具,在保证研究对象的估计值满足无偏性条件和最小方差条件的前提下求得估计值。– 对于高值数据点会使之沿某一 “脊 ”分布,而不围绕该点孤立插值,不形成 “公牛眼 ”等值线。– 克里格法极为灵活,广泛地应用于各个科学领域,适于各种类型的离散数据,网格化精度高,是极佳的网格化方法。• 最小曲率法( Minimum curvature)– 采用迭代的方法逐次求取网格节点数据,其插值面类似一个薄的、线性 —弹性形变板,该 “板 ”经过所有的数据点,且每个数据点具有最小曲率。– 由于最小曲率法采用全区的数据进行网格化,因而比较适合于数据分布不均匀的情况。在尽可能体现原数据的同时,最小曲率法产生尽可能的光滑曲面,绘制的图件比较美观。– 使用最小曲率法需要用最大偏差参数和最大循环次数参数来控制最小曲率的收敛标准,且要求至少有 4个点。– 该方法速度快,适合于大量( 1000个以上)数据的网格化,由于其主要考虑曲面的光滑性,不能达到精确的插值结果,容易超出最大值和最小值的范畴。• 径向基函数法( radial basis function)– 又称距离基函数,是多个数据插值方法组合的一种多形式网格化方法。– 其基函数是由单个变量的函数构成的,通过选择不同的基本函数来定义不同的加权方法,进行不同方式的网格化。– 所有径向基函数插值法都是准确的插值器,它们都能尽量适应数据。若要生成一个更圆滑的曲面,对所有这些方法都可以引进一个圆滑系数。– 径向基函数插值方法具有很强的拟合数据点、产生光滑曲面的能力,其适应范围也类似克里格法。• 加权反距离法 ( Inverse distance to a power)– 首先是由气象学家和地质工作者提出的,后来由于 D. Shepard的工作被称为谢别德法方法。– 基本原理是设平面上分布一系列离散点,已知其位置坐标和属性值, P( x, y) 为任一网格点,根据周围离散点的属性值,通过距离加权插值求 P点属性值。– 实质是待插值点邻域内已知散乱点属性值的加权平均,权的大小与待插值点的邻域内散乱点之间的距离有关,是距离n次方的倒数。– 加权反距离插值法认为任何一个观测值都对邻近的区域有影响,且影响的大小随距离的增大而减小。– 该方法的优点是可以通过权重调整空间插值等值线的结构,但是其计算值容易受到数据点集群的影响,计算结果中常出现孤立点数据明显高于周围数据点的现象。• 最近邻点法( Nearest neighbor)– 是荷兰气象学家 A. H. Thiessen提出的一种分析方法。– 最初用于从离散分布气象站的降雨量数据中计算平均降雨量, GIS和地理分析中多采用其进行快速赋值。– 最近邻点插值的一个隐含的假设条件是任意网格点的属性值都是用距离它最近的位置点的属性值,用每一个网格节点的最邻点值作为待求的节点值。– 该方法适合对规则分布的数据进行网格化;或者大多数数据点位于网格节点上;或者在一个完整的数据文件中,只有少数点无值,可以采用该方法来填充无值的数据点。– 总之,最近邻点插值法是均质无变化的,更适合于均匀间隔的数据插值,可以有效填充无值数据区域。• 三角剖分法( triangulation/liner interpolation )– 是一种严格的插值方法,使用最佳的 Delaunay三角形,通过直线连接各数据点形成一系列三角形,并且所有的三角形互不相交,每个三角形内的网格节点值由该三角平面决定。– 由于采用所有的数据点去构造三角形,因而原数据能得到很好的体现,给定三角形内的全部节点都要受到该三角形的表面限制。– 该方法速度快,适合中等数量、均匀分布的数据的网格化,地图上稀疏区域将会形成截然不同的三角面。– 当数据量足够时,该方法对断线的保留具有其他方法不可比拟的优势。• 多项式回归法( Polynomial regression)– 严格地说并不是一种真正的插值方法,它仅仅通过定义趋势面类型来表明原数据的大状态趋势,并不增加未知的网格节点值,实际上是一种趋势面分析作图程序,可用来确定数据的大规模趋势和图案。– 使用该方法需要考虑两方面问题:一是趋势面数学表达式的确定;二是拟合精度的确定。– 该方法具有速度快的特点,但也去掉了原数据中的局部细节,不利于资料的详细分析。类型 特征 应用条件克里格法 根据相邻变量的值,利用变差函数所揭示的区域变量的内在联系来估计空间变量数值,网格化精度高数量小于 250 个点数据的网格化,对于 2 50~1000 个数据点,效果也不错最小曲率法 采用迭代的方法逐次求取网格节点数据 方法速度快,适合于大量( 1000 个以上)数据的网格化径向基函数法 多个数据插值方法组合的、多形式的方法 适应范围类似克里格法加权反距离插值法认为任何一个观测值都对邻近的区域有影响,且影响的大小随距离的增大而减小可以通过权重调整空间插值等值线的结构,计算值容易受到数据点集群的影响,常出现孤立点数据明显高于周围数据点最近邻点法 采用距离网格节点最近的数据点的值来表明网格节点的值适合规则分布、或者大多数数据点位于网格节点上的数据,更适合于均匀间隔的数据插值,可以有效填充无值数据区域三角形剖分法 通过直线连接各数据点形成一系列互不相交的三角形,每个三角形内的网格节点值由该三角平面决定方法速度快,适合中等数量、均匀分布的数据网格化多项式回归法 仅仅通过定义趋势面类型来表明原数据的大状态趋势,并不增加未知的网格节点值实际上是一种趋势面分析作图程序,可用来确定数据的大规模趋势和图案。被广泛应用于地质科学。该方法具有速度快特点、然而其去掉了原数据中的局部细节,不利于资料的详细分析表一:地球物理中常用的 8种网格化方法及其特征1、高密度电法温纳装置实例• 实例A—三角剖分B—克里格C—反距离加权D—最小曲率E—径向函数法F—最近邻点法图 1 高密度不同网格化方法视电阻率等值线图G—多项式回归此例子的高密度数据点 553个。高密度电阻率法实测数据,取相同的参数,采用不同的网格化方法绘制视电阻率等值线断面(图 1A~F),可见效果不同。采用克里格法网格生成的图 1B,整个断面呈矩形,网格化的结果是扩大了实测数据边界,没有数据的区域插值产生,呈现 2个低阻异常区域,由一些渐变的异常点组成,有利于异常区的圈定和解释。采用三角剖分法网格生成的图 1A,整个断面呈倒梯形,网格化的结果是严格控制了实测数据边界,可清晰分辨出局部高阻异常,有利于局部异常区的圈定和解释。采用加权反距离法网格生成的图 1C,整个断面呈矩形,网格化的结果是扩大了实测数据边界,没有数据的区域插值产生,低阻和高阻异常的分界面很清晰,有的地方呈现串珠状高阻异常,形成一些孤立的异常点,不利于异常区的圈定和解释。采用最小曲率法网格生成的图 1D,与图 2B克里格插值方法的效果基本相同,有 2个低阻异常区域,由一些渐变的异常点组成,有利于异常区的圈定和解释。采用径向基函数法网格生成的图 1E,整个断面呈矩形,网格化的结果是扩大了实测数据边界,没有数据的区域插值产生,可分辨出局部高阻异常,因此有利于局部异常区的圈定和解释。但图形左侧和右侧等值线杂乱,表示插值效果不好。采用最近邻点法网格生成的图 1F,整个断面呈矩形,网格化的结果是扩大了实测数据边界,没有数据的区域插值产生,呈现2个孤立的低阻异常区域,异常区等值线稀疏,边界呈矩形。采用多项式回归法对本案例进行网格化生成的图 1G,效果最差,基本找不到异常点,网格化的结果是扩大了实测数据边界,没有数据的区域插值产生,仅显示了上下 2个渐变的异常区,分辨不出局部高阻异常,因此不利于局部异常区的圈定和解释。2、 EH4电导率成像实例加权反距离 克里格图 2 EH4电导率成像不同网格化方法视电阻率等值线图上图为今年九月份我们在青海多龙恰柔 EH4测深 L4线测深地电断面图,可以看出虽然在本测区用加权反距离、克里格法网格化绘制此等值线图,可清晰显示异常,满足绘图目的,但在实际工作过程中,还是应根据研究内容的环境特征和离散数据原有的特点,选择合适的网格化方法。在此等值线下部,克里格法绘制等值线具有良好的闭合性,对于高值数据点会使之沿某一 “脊 ”分布,而不围绕该点孤立插值,不形成 “公牛眼 ”等值线;而加权反距离则局部有 “公牛眼 ”等值线的存在。这与此矿区的实际地质情况应该比较吻合。• 总结通过对以上几种网格化方法的分析比较,可以看出,用克里格法、径向基函数法、最小曲率法网格化绘图效果较好,同时结合野外实际数据的特点,所以在此选择克里格方法、径向基函数法或最小曲率法进行网格化是切实可行的,这样既保证了网格化的精度,产生的等值线也较为光滑,符合野外实际情况。克里格插值法解决了空间连续变化的属性非常不规则时数据的插值,用反距离加权法可以有效提高所预测数据的精度。总之,在针对实际数据插值网格化时,需要根据客观环境特征和数据本身的特点,选择合适的网格化方法,进行相应的数据分析,科学、合理地选择正确有效的插值方法,才能绘制出准确、有意义的等值线图。
展开阅读全文
  微传网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

暂无评论,赶快抢占沙发吧。

关于本文
本文标题:Surfer中网格化方法的选取探究.ppt
链接地址:https://www.weizhuannet.com/p-11246923.html
微传网是一个办公文档、学习资料下载的在线文档分享平台!

网站资源均来自网络,如有侵权,请联系客服删除!

 网站客服QQ:80879498  会员QQ群:727456886

copyright@ 2018-2028 微传网络工作室版权所有

     经营许可证编号:冀ICP备18006529号-1 ,公安局备案号:13028102000124

收起
展开