Surfer8_0在重力异常数据格网化中的应用
- 格式:pdf
- 大小:318.36 KB
- 文档页数:4
Surfer8.0线框图与曲⾯图绘图实例作业四⼀、实验内容:给定数据⽂件(作业⼆数据.dat),使⽤Surfer8.0软件绘制线框图与曲⾯图,并提交报告。
⼆、实现过程:1、数据⽹格化:执⾏Grid→Data,在打开的Open对话框中选择数据⼆⽂件,单击打开后,在Grid Data对话框中设置⽹格化参数,选择克⾥⾦插值法,在Grid Line Geometry中设置X、Y线的数量和间隔,如果线的数量设置的较⾼,则可以使线框图的⽹格更密,变得圆滑。
设置好各项参数之后确定,完成数据的⽹格化。
2、绘制线框图:(1)绘制线框图:执⾏Map→Wireframe命令,在Open Grid对话框中选择⽹格化之后确定,绘制出线框图。
(2)线框图参数设置:双击图区,在弹出的Map Wireframe Properties中设置参数。
General1)在Plot Lines of Constant 中勾选X、Y、Z,⽤于显⽰X、Y、Z三条线的任意组合,还可以设置每条线的颜⾊,线型,粗细程度。
X、Y等值线的数量在⽹格化数据时设置。
Z 等值线的数量在Z Levels中设置。
当设置好Z等值线的参数后,只显⽰Z等值线如右图所⽰:2)在Base中控制线框图的基准⾯、线框图与基准⾯间的垂线是否显⽰。
只有当把X和Y轴设为不可见后,基准⾯显⽰与否才能看清楚。
当把X、Y轴设置为不可见后基准⾯的显⽰如右图所⽰:Show Base 控制基准⾯是否显⽰Show Vertical Lines 控制是否显⽰基准⾯与线框图之间的垂线。
Elevation 设置基准⾯的⾼程。
Line设置基准线的线型、颜⾊、及粗细。
3)在Border中设置线框图表⾯的外边界是否显⽰。
有FrontOnly、All Sides和None三种显⽰效果。
当⽹格表⾯的X和Y线显⽰时,该外边线不可见。
当把X和Y线设置为不可见时,选择All Sides后线框图外边界的显⽰如右图所⽰:4)在Remove Hidden Lines中设置是否显⽰线框图上被其它线遮挡住的X和Y线,当线框图前⾯遮挡住后⾯的部分时,如果该选框被选中,则图形被遮挡住的部分将不画出来,否则,线框图将以透明的⽅式显⽰。
第3章Surfer8.0绘图软件的使用3.1 软件运行环境及特点Golden Software Surfer 8.0 (以下简称Surfer)是一款画三维图(等值线,image map,3d surface)的软件,是美国Golden Software公司的系列绘图软件之一。
该软件简单易学,可以在几分钟内学会主要内容,且其自带的英文帮助文件(help菜单)是相当完美且容易阅读的,对如何使用Surfer,解释的很详细,只要学过英语的人都可以很快上手。
Surfer的主要功能是绘制等值线图(contour map),是具有插值功能的绘图软件,因此,即使你的数据是不等间距的,依然可以用它作图。
此外它还可以绘制张贴图、分类张贴图、矢量图、影像图、线框图、3d surface map,等形式的图形,其功能是比较强大的。
Surfer的安装比较简单(目前,只有Windows操作系统下的版本,最为常用的是8.0版本),只要按其提示缺省安装即可。
其安装软件的大小不到30M,一般的计算机硬件基本能够顺利使用该软件。
安装好Surfer以后,其环境界面如图3-1所示。
命令菜单绘图命令目标管理窗口工作区状态栏图3-1 Surfer8.0软件界面3.2 软件界面及命令菜单Surfer软件的界面非常友好,继承了Windows操作系统软件的特点。
从图3-1中可以看到,其最上方为命令菜单,在命令菜单的下方是命令菜单中的快捷工具栏(共两行),左侧的空白区域为目标管理窗口,用来更加方便的管理绘制的各个图形要素,右侧的空白区域为工作区,用来绘制图形,最右侧的一个竖条工具栏是绘图命令的快捷方式。
下面详细介绍各个命令菜单的主要内容。
3.2.1文件菜单(F)“文件菜单”如图3-2所示,主要是对文件进行操作,如文件的建立、加载、打印设置等。
图3-2 文件菜单新建—用来新建一个工作窗口,点击后即出现图3-1界面。
打开—打开一个已经存在的Surfer可以识别的文件。
实习三Surfer8.0的初步使用(一)2005年9月30日实验完成课时:2学时(2005年10月10日星期一上午11:00 ~ 13:00)实验软件:Golden Software Surfer 8.0实验数据:.. \Program Files\Golden Software\Surfer8\Samples实验目标:Golden Software Surfer 8.0 (以下简称Surfer)是一款具有插值功能的二/三维图形(等高线、image map、3d surface)绘制软件。
Surfer的主要功能是绘制等高线图(contour map),此外它还可以绘制post map,classed post map,vector map,image map,wireframe map,3d surface map等形式的图形。
本实习实验目标包括:(一)等值线绘图软件Surfer软件的初步了解;(二)初步掌握Surfer的七种插值功能;(三)利用Surfer绘制等高线地形图;实验步骤:(一)由离散点的XYZ文件创建[.GRD] 文件为了产生等值线图或线网图, 或执行需要网格文件的操作如网格代数处理、计算体积和面积, 光滑或计算残差,必须预先对数据文件进行网格化处理。
由Grid→Data 打开所需的数据文件后,出现散点数据插值法对话框Scattered Data Interpolationd。
该对话框由4个Tab组成。
1 Date数据列组(Data Columns)设定X、Y和Z的值对应于数据文件内的哪一列。
相同数据点处理方式选择。
由下拉菜单Duplicates可以选择对有相同的X、Y坐标但有不同的Z值的两个或多个点Z值处理办法,包括取平均值(Average)、保留一个并删除其它值(Delete)、对各点Z值加和(Sum)等。
有些网格化方法还可以选择忽略这类数据点(Ignore)。
数据筛选(Filter)可以使用布尔表达式对原始数据点进筛选。
Golden Software Surfer 8.0 (以下简称Surfer)是一款画三维图(等高线,image map, 3d surface)的软件,该软件简单易学,可以在几分钟内学会主要内容,且其自带的英文帮助(help 菜单)对如何使用surfer解释的很详细,其中的tutorial 教程更是清晰的介绍了surfer的简单应用,应该说surfer 软件自带的帮助文件是相当完美且容易阅读的,只要学过英语的人都可以很快上手。
Surfer 是具有插值功能的绘图软件,因此,即使你的数据是不等间距的,依然可以用它作图。
但依据作者的经验,最好不使用Surfer 自带的插值功能,尤其是要精确确定等高线时。
由于surfer是美国的一个软件,它不提供对中文的支持,这可以算的上一个小的遗憾。
Surfer 的主要功能是绘制等高线图(contour map),此外它还可以绘制post map,classed post map,vector map,image map,wireframe map,3d surface map,等形式的图形。
其功能是比较强的,但没有各种投影变化是它的一大缺点。
尤其是在等高线领域,这不能不说是它的应用受到限制的地方。
由于surfer软件没有中文说明书,对一些初学者来说可能会存在上手较难的问题,鉴于此种需求,编写了这一初学者参考手册,希望对大家有所帮助。
二、等高线的绘制Surfer 的最主要的功能是绘制等高线图,但并不是我们具有了数据文件就可以直接绘制等高线,surfer要求绘制等高线的数据有特殊的格式要求,即首先要将数据文件转换成Surfer 认识的grd 文件格式,才能绘制等高线(当然,可以直接生成surfer 接受的ascii 码的grd 文件格式,这样就可以直接作图,此方法将在后面介绍,首先我们介绍常用的作图方法)。
假设你有三列数据分别为X,Y,Z,其中Z 为点(x,y)处的值,存在文件test.dat中(数据见附件),其中第X坐标,第Y坐标,第三x,y)上的值Z,则绘制等高线的步骤如下:步骤一:把数据文件转换成grd文件1. 打开菜单Grid | Data... ,在open 对话框中选择数据文件test.dat2. 这会打开“Grid Data”对话框。
一天学Surfer8.0作等值线图图心得从事煤炭资源勘查工作中,需要做很多等值线图,周围同事多是用传统的三角网插值方法,一个钻孔一个钻孔地量取之间的距离,然后在Excel表格中计算分点距离,再在CAD 中描点连线……虽然这是一种严密的插值法,但是感觉过程复杂,与传统手工作图无异,需耗费很大精力。
同学推荐用Surfer软件,我用了一天时间,基本掌握了这种作图方法。
这里,我要感谢水文07-1班的邓同学以及其他同学朋友给我提供的帮助。
Surfer做等值线图主要用到两种数据,一种是钻孔数据,包括钻孔坐标和需要描述的目标数据(如煤厚、灰分、硫分等),另一种是边界数据,主要是矿区边界,即拐点数据。
这两种数据分别对应着Surfer里的dat文件(也可以是其它可识别格式,如txt)和bln文件。
第一步:制作数据文件1、制作钻孔数据文件这一步可以在Excel、WPS、记事本中完成,也可以在Suerfer工作表中建立。
其格式如下图:数据录入完成后,保存为txt格式。
注意:我用的是WPS,另存为的时候须选择制表符分割的txt文件。
2、制作bln文件(边界数据)一般情况下,边界数据都是比较少的,可用Surfer工作表直接建立,保存为bln文件即可。
其格式如下图:第一行31代表边界拐点坐标的个数,后面的0值表示后面的等值线绘制在边界内部,如果写成1,表示等值线绘制在边界外部。
后面的数值为钻孔坐标,第一栏为横轴(Y值)。
保存为bln格式。
注意:我们工作区域只有30个拐点,写31的原因是最后一个钻孔坐标与第一个相同,形成一个封闭区域。
Ps:有时候拐点数较多,又没有拐点坐标文件,需要我们直接从图上量取,比较方便的方法是用工具展点。
比如,我用的CAD作图,作图比例为1:1000(保证了展点出来的坐标不变),矿区边界是多段线圈成,拐点较多,我就用到了Dxf2xyz2.0这个工具来展点。
把矿区边界带基点复制到空白文档中,并校正。
保存为dxf格式,据说dwg格式也行。
目录一.实验目的和要求 (2)二.硬件与软件环境 (2)三.实验内容 (3)3.1反距离加权插值法(Invme Distmeeto a Power) (3)3.2克里金插值法(Kriging) (11)3.3最小曲率法(Minimum Curvature) (14)3.4改进谢别德法(Modified Sheparfs Method) (16)3.5自然邻点插值法(Natural Neighbor) (18)3.6最近邻点插值法(Nearest Neighbor) (19)3.7多项式回归法(Polynomial Regression ) (21)3.8 径向基函数插值法(Radial Basis Function) (23)3.9 线性插值三角网法(Triangulation with Lin-ear Interpolation) (26)3.10 移动平均法(Moving Average) (28)3.11 数据度量法(Data Metries) (29)3.12 局部多项式法(Local Polynomia1) (31)四.结束语 (32)一.实验目的和要求1. 目的:✓掌握Surfer软件的基本操作技巧,如网格化、修改地图各种属性、导出图件等。
✓了解suffer自带的各种网格化方法的背景及原理。
✓使用不同的网格化方法,并对结果做对比分析。
2. 要求:✓熟悉Surfer软件的菜单,掌握Surfer软件的基本操作技巧。
✓根据数据文件(作业二.dat),对Surfer软件的各种网格化插值方法进行比较,并提交报告。
✓爱护实验室设备,独立完成实验报告。
二.硬件与软件环境1.硬件华硕笔记本K40IN2.软件平台操作系统:Windows XP软件:Suffer 8.0 中文版三.实验内容3.1反距离加权插值法(Invme Distmeeto a Power)3.1.1方法原理背景反距离加权插值法又被称为谢别德法(Shepard) 方法。
运用surfer8.0软件绘制物探平面剖面图胡江建(新疆维吾尔自治区有色地质勘查局地球物理探矿队 乌鲁木齐830011) 物探数据处理中的平面剖面图(简称平剖图),能够很好地反映数据的区域变化趋势及局部细节,对异常识别帮助很大,是物探数据处理、解释结果的一种重要的常规表现形式。
以往在进行物探数据处理时,往往在计算机上利用AutoCAD软件绘制出剖面观测值曲线及测线,然后手工填充上色。
由于某些物探测区大、测线多且密集,手工着色费时费力、效率低,不能适应快速高效、数字化、自动化的数据成图要求。
本人在研究了SU RFER8.0软件的数据的特点后,开发出利用SU RFER8.0软件的BL N文件进行绘制平剖图的成图技术。
实际应用表明,作为一种辅助成图手段,该方法可以充分利用SU RFER8.0软件强大的绘图能力,高效地绘制出美观、一致的平剖图,提高了物探数据处理的效率和质量。
1 运用SU RF ER8.0软件绘制平面剖面图的原理物探平剖图的绘制是将每条测线的测量数据沿测线按数据值的正、负大小,以一定的比例展布在测线两侧,形成色彩相对的平面图。
一般在绘制平剖图时,每条测线正观测值剖面的部分填充暖色,剖面位于垂直航迹线方向的上方;负观测值剖面的部分填充冷色,剖面位于垂直于航迹线方向的下方。
无论测线上测点数据值是大还是小,同为正(或负)观测值的剖面峰值的颜色,都为相同颜色。
我们可以利用SURFER8.0软件的BLN格式的文件将物探平剖图的观测值的边界表现出来,再通过SURFER8.0软件的张贴图将物探剖面的点、线在计算机中绘制出来,两图相叠加,即可绘制出物探平剖图。
SURFER8.0软件的BLN格式的文件称之为白化文件(Blanking File),用于表现边界的信息,也可以用于白化网格文件的某个多边形区域。
白化文件用逗号分隔,是个明码文件。
白化文件的格式是这样的:lengt h,flag“Pname1”x1,y1x2,y2...xn,ynlengt h,flag“Pname2”x1,y1x2,y2...xn,yn其中①lengt h是一个>0的整数,表示后面跟着的X,Y坐标对数目。
Surfer8.0学习手册一、简介这是一个非常不错的科学类绘图软件,是地质工作者必备的专业成图软件。
Golden Software Surfer 8.0 (以下简称Surfer)是目前最新的版本,可以轻松制作基面图、数据点位图、分类数据图、等值线图、线框图、地形地貌图、趋势图、矢量图以及三维表面图等;提供11种数据网格化方法,包含几乎所有流行的数据统计计算方法;提供各种流行图形图像文件格式的输入输出接口以及各大GIS软件文件格式的输入输出接口,大大方便了文件和数据的交流和交换;提供新版的脚本编辑引擎,自动化功能得到极大加强。
该软件简单易学,可以在几分钟内学会主要内容,且其自带的英文帮助(help 菜单)对如何使用surfer解释的很详细,其中的tutorial教程更是清晰的介绍了surfer 的简单应用,应该说surfer软件自带的帮助文件是相当完美且容易阅读的,只要学过英语的人都可以很快上手。
Surfer是具有插值功能的绘图软件,因此,即使你的数据是不等间距的,依然可以用它作图。
但依据作者的经验,最好不使用Surfer自带的插值功能,尤其是要精确确定等高线时。
由于surfer是美国的一个软件,它不提供对中文的支持,这可以算的上一个小的遗憾。
没有各种投影变化也是它的一大缺点。
尤其是在等高线领域,这不能不说是它的应用受到限制的地方。
尽量在Windows2000/XP系统下使用Surfer8,尤其是文档中使用了大量中文的情况下,在Windows9x下要把中文调整得好看简直就是一件杀人的工作,因为如果您在9x下使用宋体,字间距大到不能要,只有使用英文字体,这时候中文就互相交迭,在Windows2000/XP中就没有这个问题,直接使用宋体就可以看到漂亮的中文了。
此外,在各种报告中,如果是2000/XP系统,那么在报告编辑器中是可以直接看清中文的,但是格式就变得乱七八糟,在9x系统中格式是正确的,但是不能正常显示中文,必需另存为文本文件才可以看见报告的内容。
科苑论谈Surfer8.0在环境质量评价中的应用———以宜宾市屏山县新发乡浅层地下水水质现状评价为例肖慧英欧亚波(成都理工大学环境与土木工程学院,四川成都610059)1前言随着环境质量评价手段的不断升级,计算机技术的迅速发展,如何能够选择实用易行的模拟方法,增加可视化效果,更清晰直观的对水质进行分析受到越来越多的环境工作者的重视。
Surfer是美国Golden软件公司的产品,主要用于绘制等值线图以及相应的三维图形[Surfer8地理信息制图]。
Surfer8.0以其简单易学,绘图功能强大适用于众多科学领域,包括地质学、地球物理学、水文学、气象学、海洋研究等。
国内已有将Surfer软件应用于环境领域的研究。
如魏文飞等利用Surfer软件对福建省龙岩市的区域环境噪声进行了等值线绘制,为当地的声环境功能分区提供了简洁明了的依据[1];文雅等以中国1950~1998年160个测站年均降水量为例,利用Surfer软件的空间插值功能,方便快捷地绘制了中国的降水等值线分布[2],陆雍森将Surfer8.0应用于大气污染物扩散模拟及评价[3];杨卫华将Surfer与Mapinfo相结合,绘制大气污染物浓度等值线,不仅保证了绘图分析精度,而且还大大缩短了绘图和分析时间[4]。
多年来,国外许多不同应用领域的专家使用Surfer开发出了相当多的应用程序,大大扩展了应用领域。
但是国内的应用范围和水平与国外相比还存在不小差距。
这里比较具体地介绍Surfer8.0在环境科学领域中的应用,证明了该方法是一种科学的、简便易算且特别适合大范围多指标环境质量评价的新方法,同时也提出其应用存在的问题。
2Surfer8.0在环境质量评价中的具体应用采用改进的内梅罗水质指数方法对宜宾市屏山县新发乡浅层地下水水质现状评价结果如图1所示:应用Surfer8.0软件可以绘出八种图件分别是:等值线图(ContourMap)、基面图(BaseMap)、粘贴图(PostMap)、影像图(ImageMap)、地貌晕渲图(ShadedReliefMap)、矢量地图(VetorMap)、线框图(Wireframe)、曲面图(Surface)。
第28卷第1期2008年1月海 洋 测 绘H YD RO GRA PH IC SU RV EY I N G AND CHA R T I N GV ol 128,N o 11Jan .,2008收稿日期:2007209203;修回日期:2007210228作者简介:庞振兴(19822),男,河南商丘人,硕士研究生,主要从事物理大地测量研究。
Surfe r 810在重力异常数据格网化中的应用庞振兴,张传定,叶修松(解放军信息工程大学测绘学院,河南郑州 450052) 摘要:介绍了利用Surfer 810软件进行重力异常数据格网化的方法,并结合重力场数据的特点,编制了相应的数据格式转换程序,直接利用Surfer 810软件中几种数据插值方法,对陆地、山区、海洋等不同地形数据进行格网化处理,并对结果进行交叉验证和对比分析,结果表明利用Surfer 810软件进行重力异常数据格网化的方法是切实可行的。
关键词:重力异常;格网化;Surfer 810;数据转换中图分类号:P223+139 文献标识码:B 文章编号:167123044(2008)01200432031 引 言平均重力异常是地球重力场的基本数据,用途十分广泛[1]。
在垂线偏差、高程异常、扰动重力等的计算中都要用到平均重力异常数据。
高分辨率的平均重力异常数据是求解高精度、高分辨率的垂线偏差、高程异常、扰动重力数据的基本保障。
但是,对于一些地形比较复杂、重力实测点稀疏的地区,需要对实测的数据进行一定的处理以满足高分辨率重力异常数据的要求。
通常要根据测区地形特点、数据分布情况的不同,采用不同的格网化方法。
但是,总体来说这些方法不易计算且操作复杂。
Surfer 是Golden Soft w are 公司研制的基于W indo ws 操作系统的数值分析及绘图计算软件[2]。
它的功能十分强大,能支持12种内插方法对离散数据进行格网化,生成规则的格网数据。
利用Surfer 软件,采用合适的插值算法对重力异常数据进行格网化,只需编写简单的数据格式转换程序,便可以直接计算得到高分辨率的格网化重力异常数据,使用起来十分的便捷高效。
2 重力异常数据的计算、格式转换及格网化211 重力异常数据的计算(1)空间重力异常。
是指地面上一点的重力值与该点在近似地形面上投影点的正常重力值之差,其值常采用重力异常值加上空间改正的方法来计算,即:Δg 空=g -γ0+013086Hγ(1)式中,g 是地面点的观测重力值;γ0是正常椭球面上的正常重力值;H γ是该点的正常高;013086H γ是空间改正。
统计数据表明,Δg 空与高程存在线性关系,其变化与地形密切相关。
因此在地形起伏较大的山地、丘陵地区,利用拟合、插值等方法计算未测点的重力异常时,为了减小内插得到的未测点的空间重力异常值的误差,需要引入布格异常作为过渡。
布格异常是由空间重力异常加上层间改正求得的,即:Δg 布=Δg 空+δg 层(2)式中,δg 层=-010418DH,取D =2167g c m-3,得到Δg 布=Δg 空-011116H γ。
(2)重力异常数据的存储格式。
重力异常数据是以数据文件的方式存储的。
以澳大利亚内陆地区(纬度在-36°~-15°,经度在115°~130°)数据为例,其实测的重力点数据存储形式为:1969111969112617100-32114877711616664592491449794338144-1261641969111969112618100-32115044511617864592101389794381136-1611901969111969112621100-32125378511610581322811689794431169-551541969111969112622100-32125211211613031322571909794489132-431311969111969112623100-32125045111614297913321459794322143-62123……由左至右数据依次代表:测点编号,测站编号,海 洋 测 绘第28卷测点纬度,测点经度,测点高程,测点观测重力值,测点空间重力异常值。
其中空间异常是根据实测重力数据计算而得。
该数据并不是格网化的数据,而且还包含很多相关信息。
因此,欲用该数据进行计算,必须将其格式进行转换,提取其中的所需信息,进行插值、拟合等计算,最后得到格网化的重力异常数据。
212 重力异常数据格式的转换(1)Surfer 格网化计算所需要的数据格式Surfer 软件具有强大的数据分析功能,支持多达12种内插方法对离散点数据格网化,生成规则的格网数据。
Surfer在外部数据转换为内部识别的1grd 数据格式之前,只接受“X 、Y 、Z ”型数据格式,即“位置坐标(X,Y )+Z 属性值”,其中属性项可以是多项,但位置坐标却是唯一的[3]。
(2)数据格式转换首先打开数据文件,读取其中的经纬度、高程和空间重力异常值信息,然后按公式(2)将空间重力异常值转换为布格异常值,按照弧度与度之间的转换公式将用弧度表示的经纬度转换为度,最后将(B ,L )和布格异常值按照Surfer 所需要的格式存储。
得到Surfer 能够直接读取的数据格式为:-321168777 1161356432 1521479996-3211804451161506443731519997-32128211211617031342591660004-32135045111618297932451279999-32137378511619582343091779999…213 重力异常数据格网化离散数据格网化就是根据参考点上的属性值求出其他待定格网点上的属性值,在数学上属于插值问题。
影响离散数据格网化精度的一个因素就是内插方法的选取。
在Surfer 8中提供的内插方法达12种,本文选取了三种典型地形的已测布格异常数据作为试验数据,利用Surfer 中的12种插值方法对其进行1′×1′的格网化,并且选取其中一定数目的已知点作为未知点进行交叉验证。
现将插值结果统计如下。
(1)陆地数据格网化结果本文主要选取了澳大利亚内陆已测的1°×015°范围内的已知重力异常数据进行格网化。
该地区的三维地形图(L,B ,H )如图1所示。
选取其中的400个点进行1′×1′的格网化,然后利用其余30个点数据进行验证,并选取其中的6种精度较好的插值结果进行统计分析。
统计结果见表1。
图1 陆地三维地形图表1陆地地形格网化结果精度统计表插值算法检验点数绝对值最大偏差(mGal )绝对值最小偏差(mGal )差值的平均值(mGal )标准偏差标准方差反距离插值法3017.90.87.70.93.9克里金插值法305.54.92.20.62.6最小曲率法3010.03.13.70.73.1改进谢别德法304.63.0-0.20.52.2径向基函数插值法305.15.01.70.62.6局部多形式插值法309.03.03.90.62.8(2)山区数据格网化结果本文选取了澳大利亚1°×015°范围内的1050个已测重力点进行数据格网化。
该地区的三维地形图(L,B ,H )如图2所示。
图2 山区三维地形图选取其中的1000个点对该地区进行格网化,并利用其中的50个点的数据作为已知数据进行交叉验证。
选取其中6种精度较好的插值方法进行统计分析。
统计结果见表2。
44第1期庞振兴,等 Surfer 810在重力异常数据格网化中的应用表2山区格网化结果精度统计表插值算法检验点数绝对值最大偏差(mGal )绝对值最小偏差(mGal )差值的平均值(mGal )标准偏差标准方差反距离加权插值5014.70.96.70.84.1克里金插值法5051.03.316.711.08.7最小曲率法50361.44.811.24.7改进谢别德法5038.63.612.712.44.3移动曲面法5034.74.913.616.611.5多项式衰减法5047.43.016.417.025.8(3)海洋地区数据格网化结果本文选取了澳大利亚海域1°×015°范围内的978个已测重力点进行数据格网化。
该海域的三维地形图(L,B ,H )如图3所示。
图3 三维海底地形图由于海域的空间重力异常与布格异常相差很小,习惯上利用空间重力异常数据进行计算。
选取其中的900个点对该地区进行格网化,并利用其中的78个点的数据作为已知数据进行交叉验证。
选取其中6种精度较好的插值方法进行统计分析。
统计结果见表3。
表3海洋地形格网化结果精度统计表插值算法检验点数绝对值最大偏差(mGal )绝对值最小偏差(mGal )差值的平均值(mGal )标准偏差标准方差反距离插值法7812.30.96.71.24.1克里金插值法786.50.63.20.91.4最小曲率法7813.21.45.84.23.7改进谢别德法7814.53.85.62.04.8移动曲面法7820.34.911.69.68.5多项式衰减法7816.83.49.48.06.83 插值结果精度分析(1)插值结果与实测结果差比较对于陆地,改进谢别德法计算得到的相对精度最好,偏差的最大值为416mGal,最小仅为310mGal 。
对于山区,由于地形起伏较大,数据不充足,内插结果相对较差,但反距离插值法最大偏差为1417mGal,最小偏差仅为019mGal,能满足实际应用的需要。
而在海域,反距离加权插值法、克里金插值法、最小曲率法与改进谢别德法的插值结果均不错,特别是克里金插值法最大偏差仅615mGal,最小偏差019mGal,精度相对较高。
(2)插值的标准偏差与标准方差[4]标准偏差是指插值结果与实测结果之差的绝对值之和除以插值个数,即:(|α1|+|α2|+…+|αn |)/n;标准方差是指插值结果与实测结果之差的平方之和除以插值个数的平方根,即:((a 1-α)2+(a 2-α)2+…+(a n -α)2)/n 。
由表1~3可知:陆地地区的标准偏差结果在小数量级,标准方差在个位量级;山区的反距离插值法的标准偏差为018mGal,标准方差为411mGal;而海域克里金插值法的标准偏差为019mGal,标准方差为114mGal,能满足精度要求。
(3)不同地形对插值方法的要求由表1~3统计结果可知:对于陆地地形的改进谢别德法,对于山区的反距离加权插值法和对于海域的克里金插值法均能得到精度较好的插值结果,实验表明利用这些插值方法计算得到的1′×1′的重力异常格网化数据精度能满足实际需求。