Surfer8_0在重力异常数据格网化中的应用
Surfer8.0线框图与曲面图绘图实例
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线，当线框图前⾯遮挡住后⾯的部分时，如果该选框被选中，则图形被遮挡住的部分将不画出来，否则，线框图将以透明的⽅式显⽰。
利用Surfer8.0软件绘制煤底板隔水层厚度等值线图
图ｌ ４ １ 。简要 步骤 如下 ：
Ｓ ｕ ｒ ｆ ｅ ｒ ８ ． ０ 软件 基 于 Ｗｉ ｎ ｄ ｏ ｗｓ 操作 系 统 的三 维绘 图 软 件 Ｉ 。Ｓ ｕ ｒ ｆ ｅ ｒ ８ ． ０ 软 件 囊 括 了插 值 、 绘图 、 数 据 分 析 等
（ １ ） 数据准备 。通过野外测量 、 采集或者既有 的格 式（ 如： ． ｄ ｘ ｆ 、 ． ｍｉ ｆ 等） ， 将 空 间离 散 的数 据 （ ｘ、 ｙ、 Ｚ ） 转 化为 Ｓ ｕ ｒ ｆ ｅ ｒ ８ ． ０软 件 可 以 接 受 的形 式 （ ． ｄ ａ ｔ 、 ． ｔ ｘ ｔ 等） 。 选择合理的网格方法进行插值 ． 得到所需要 的格 网文
更 改标题 、 图例 。 ２ Ｓ ｕ ｒ ｆ ｅ ｒ ８ ． ０ 软 件绘 制等 值线 图 的原理 与步 骤
收稿 日期 ： ２ ０ １ ７ － ０ ６ — ０ ７ 修回 日期 ： ２ ０ １ ７ — ０ ６ — ０ ８
对话框 ， 通过 “ 常规” 、 “ 等级” 、 “ 背景” 等操作 ， 来完成等 值线 图 的调整 修饰 。
隔水层厚度等值线 图是煤矿生产过程中的重要 图 件 之一 ， 尤 其 在 矿 井 防治 水 工 作 中广 泛 使用 。 图 中 内 容一般包括等值线 、 钻孔 、 隔水层厚度 、 矿界。传统 的 手 工 做法 是 先做 好 地 形 地 质 底 图 ， 然 后将 对应 钻 孔 隔 水层厚度标注在底网上 ， 再根据三角网插值法插值 ， 最 后连接相 同值绘制成网。这种方法制 图速度慢 、 效率
如何使用Surfer8_0画等值线图
图 4 经过填充处理的等值线图
本文简单介绍了 Surfer 8. 0的基本作图方法 , Surfer 8. 0 还有很多辅助功能 ,如应用 Surfer 8. 0 给出数据文件的统计性质 ;在基面图中添加台站 坐标和名字等等 。这些功能可以在使用中不断摸 索。
参考文献 :
[ 1 ] 梁 亮. Surfer软件在气象主分量分析中的应用 [ J ]. 浙江气 象. 2005, 26 ( 3) : 29 - 33.
第 3期 2006年 8月
文章编号 : 1007 - 6190 (2006) 03 - 0064 - 02
广东气象 Guangdong M eteorology
No. 3 Aug, 2006
如何使用 Surfer 8 . 0画等值线图
曾志雄 , 陈慧娴
(广州市番禺区气象局 , 广东广州 511400)
2. 1 把. dat文件转换成. grd文件
打开 Surfer 8. 0,在菜单栏单击“网格 ( G) ”| “数据 (D ) ”,选择要网格化的数据文件 gdr. dat并 确定 。此时弹出“网格化数据 ”对话框 ,要求确定 作图的 XYZ3 列 。由于我们的数据只有 3 列 ,而 且按照相应的顺序排列 ,所以不用选择 。如果这 些. dat文件是多列的 ,那么要进行选择 [ 2 ] 。然后 在“格点化方法 ”里选择一种插值方法 。在网格 化过程中 Surfer会自动进行插值 。在“输出网格 文件 ”里填入要输出的文件名和路径 ,然后在“网 格线索几何学 ”里设置网格点数以及作图范围的 上下 、左右边界值 。最后点击保存 ,所需用的 gdr.
1 资料处理
首先把相应站点的春季降雨量值资料 ,处理
grd文件就生成了 。用 Surfer 8. 0 打开该文件 ,如 图 1所示 。
美国Surfer 8.0软件在地下水环境质量评价中的应用
Ｓ ｒ ｒ ． 提供 了反距离加权 插值法 、 ｒ 一 ｕｆ ０ ｅ８ Ｋｉ ｇ
收稿 日期 ：００— ２ ５修 改 日期 ：０ ０ ５— ５ ２ １ ０ —２ ； ２ １ —０ ０
有强大的绘图、 插值和分析计算等功能 ， 已广泛应
一
作者简介 ： 王敏 （９ ９ ， ， １ ７一） 男 工程 师 ， 硕士 ， 事油 气 从
Ｓ ｒ ｒ ． 具有强大的绘图能力 ， ｕｆ ０ ｅ８ 不仅能将 离 散数据经过插值 生成规则 的格 网数据 ， 而且还 可
将 其 转换 成 丰富多 彩 的基 面 图 、 据 点 位 图 、 类 数 分 数 据 图 、 值线 图 、 形 地 貌 图 、 势 图 、 量 图 以 等 地 趋 矢 及 三维 表 面 图 ， 在 各种 图上 很 方 便 地 标 注 、 改 并 修 和制作 可 视化 图件 。 １３ 计 算分析 功 能 ．
．
美 国 Ｓｒｅ ． ｕｆｒ ０软件在地下水环境质量评价 中的应用 ８
王 敏 孙 ， 敏 曾 睿。 ， ６ ０８ ； ．昆明理工大学信息工程与 自动 １０ １ ２
昆明 ６０ ５） ５０ １
（．中国石化西南油气分公司勘探开发研究院 ， １ 四川 成都
化学院， 云南
摘
要： 通过 Ｓ ｒ ｒ ． 软件在某地地下水环境评价中的应用实例 ， 阐述 了其在浅层地下水污染 物浓度空 间分 布图绘 ｕｆ ０ ｅ８ 重点
田环境监测与环保科技研 发工作。
３ 一 ０
第 ３ 卷
第 ２ 期
王
敏等． 国 Ｓ ｒ ｒ ． 美 ｕｆ ０软件在地下 水环境质量评价 中的应用 ｅ８
２１ ０ １年 ４月
ｉ 插值法、 ｎ ｇ 最小 曲率法 、 改进 Ｓ ｅａｄ ｈ ｐｒ 插值法 、 自 然邻点插值法 、 最近邻 点插值法 、 多元 回归法 、 径 向基 函数 法 、 性 插 值 三 角 网法 、 动 平 均 法 、 线 移 局
Surfer8.0绘图软件的使用教程
第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的初步使用
实习三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）可以使用布尔表达式对原始数据点进筛选。
surfer8中文说明书
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作等值线图图心得
一天学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格式也行。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第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′的重力异常格网化数据精度能满足实际需求。