化探元素等值线图绘制技巧初探之离散数据网格化方法选择及白化方法简介
离散点插值方法、等值线的绘制及平滑技巧
离散点插值方法、等值线的绘制及平滑技巧2008-06-10 22:45由于等值线图看起来非常直观、形象,因此在天气预报、气候预测分析等方面用得非常多,已成为预报员不可缺少的工具之一。
如各等压面层的位势高度图、高空环流、温度及降水分布图等等。
目前也有一些非常好的微机用绘图软件,如SUFER、GRADS 等。
这些软件一般都只有DOS 版,在流行的WINDOWS平台上,虽然可以调用,但不能使用鼠标操作,故不如使用在WINDOWS 环境中开发的软件方便。
因此,许多希望在自己开发的应用程序中能方便地显示及打印各类等值线图的人,都想知道绘制等值线图的原理方法。
如何用格点资料绘制等值线图在文献〔1〕中已有介绍,而离散点(如气象台站) 的资料必须通过插值才能绘制等值线图。
插值的方法有几种,比如三角网插值,它是将相邻的三个点连成一个个三角形,然后用文献〔1〕介绍的追踪法或其它方法在三角形边上进行插值。
此方法的优点是需要插的值少,插值算法简单,数据处理量少。
缺点是三角网的生成随意性很大,任意四个点可生成二组不同的三角形,不同的三角网插值得出的等值线也不可能完全相同。
若人为固定三角网,当有资料缺测时,就不得不重新调整三角网,因此程序的通用性不高。
目前比较流行的是通用性好的网格化方法。
1 离散点网格化理论上,离散点网格化可采用局部曲面拟合方法,用多元回归方法建立V ( X , Y) = a + bX + c Y或二次方程V ( X , Y) = a + bX + c Y + dX2 + eY2 + f X Y所谓局部,是指采用拟合点周围一定范围的部分离散点进行拟合。
这种做法虽然在许多情况下效果不错,但通常总是有些地方与实际情况有较大出入。
因此根据人工绘制等值线时的直接内插方式,采用以下几个步骤。
111 定点(1) 根据离散点的分布范围确定整个网格的范围及网格距,网格距不宜过大或过细。
(2) 确定每个离散点属于哪一个网格,也即每个网格包含哪些离散点。
等值线图及其计算机绘制
等值线图及其计算机绘制摘要：本文介绍了等值线图的概念，网格化处理，并且给出了二维、三维等值线的计算机算法四维等值曲面绘制的基本思想。
关键词：等值线等值线图网格填充曲面绘制在各种工程及船体、汽车车身设计中大量使用等值线和等值线图。
等值线是平面区域内某些物理量数值相同的点构成的曲线，等值线图呈现平面区域内标量参数的分布状况。
等值线图是一种十分有用的图形，常常用来表示连续分布逐渐变化的绘图现象的数量特征，反映绘图对象的差异变化及其某些现象的分布规律，它能在形态和数值两个方面给人以完整概念。
计算机绘制等值线图快速、高效、图形简洁、清楚，因此得到了广泛的应用。
网格化处理是等值线图计算机绘制的基础并且等值线图可分为二维等值线、三维等值线及四维空间等值曲面，其计算机算法也不尽相同。
1、计算机绘制等值线图的基础——网格化平面区域可以由任意多边形网格构成，即可以是纯三角形网格，纯四边形网格等各种边数相同的多边形网格构成的区域，或者由各种边数不相同的多边形网格混合构成的区域。
三角形网格主要用于散乱数据的曲面拟合，可通过三角剖分及其优化后得到。
在三角形网格上，等值线在每个三角形上都是从一条边进入后到另一边出来。
对于四边形(矩形)网格，若原始数据点是散乱分布，可利用shepard 方法或最小二乘法等方法来求出网格点上的值。
网格可进行细分，具体作法是将网格横向和纵向间距缩小为原来的一半，用原来的网格值通过按距离加权平均法(局部化)等方法，求出新网格点的值。
这样，得到的网格值误差可能较大，在实际情况中，我们可利用原始数据或其他附加条件对网格值进行修改。
2、二维等值线图二维等值线图包括：平面区域的填充，颜色标尺及区域内的二维等值线。
2.1 平面区域的填充及颜色标尺用连续变化的颜色填充平面区域，根据区域内各网格的物理量数值大小确定该网格填充的颜色。
我们建立一个颜色标尺来确定它们之间的对应关系。
（颜色的RGB值与物理量数值之间为线性关系）。
绘制等值线
MAPGIS绘制等值线第一节整体思路及制图目标物化探图件中，等值线图是物化探工作中必须的作图解释环节，是对地质现象分析解译的一个很重要的方法。
这里我们介绍和要求掌握两种用MAPGIS绘制等值线图的方法：1、掌握利用网格化剖分（GRID模型）绘制等值线图的方法；2、掌握利用三角剖分（TIN模型）绘制等值线图的方法。
第二节数据准备与制图步骤一. 网格化剖分（GRD模型）1、数据准备：生成“制表符分隔的文本文档”或“逗号分隔的文本文档”（如图7.1）。
图7.1 制表符分隔的文本文档2、打开MAPGIS→空间分析→DTM分析，选择“GRD模型”→“离散数据网格化…”（图7.2），然后选择要进行网格化的文件打开。
图7.2 打开离散数据网格化步骤3、在弹出的“离散数据网格化”对话框中（图7.3），按下列顺序设置各项：1）设置数据列：按数据情况设置坐标X、Y及物性值Z所在的列。
2）网格参数设置：这里只要设置合适的网格间距就可以了。
3）网格化方法设置：一般情况下，物探、化探的数据网格化方法都选择Kring泛克里格法网格化。
单击“选择”按钮，弹出“克里格网格化配置对话框”（图7.4），可以对网格化进行一些配置（对配置内容不是很熟悉的情况下，建议选择默认的值即可）；单击“搜索”按钮，弹出“网格化点搜索配置对话框”（图7.5），可以对网格化的搜索类型、规则等进行配置。
4）设置输出网格文件名及存储路径。
设置好以上信息后，单击“确定”即可。
图7.3 离散数据网格化参数设置对话框图7.4 克里格网格化配置对话框图7.5网格化点搜索配置对话框4、选择“GRD 模型”→“平面等值线图绘制”，选择网格化数据文件（*.GRD）打开后，弹出图7.6，然后按下列步骤进行等值线图绘制参数设置。
图7.6 设置等值线参数对话框1）“设置”栏里：选择“等值线套区”（可以对等值线进行普染色，即建立区文件）、绘制色阶和保留边界线（如图7.7）。
2）“光滑”栏里：选择“等值线光滑处理”（可以对绘制的等值线进行一些光滑处理，使线条更流畅），然后在光滑度里选择合适的程度，一般选择中程度即可（如图7.7）。
MAPGIS操作-等值线图
二、绘制等值线图
处理结果：
二、绘制等值线图
点击绘制平面等值线图后设置类似于展点中设置。 全部设置完成后点击确定。保存生成的点线 面文件。 在mapgis中打开即可。能与原图套合即操作 成功。 一般按照实际情况可以使用工程裁剪改变等 值线图外轮廓形状（自学）。
展点过程即是将野外采集的点位置数据按比 例投影于图上。 Mapgis的此功能一般与excel联用。 首先对mapgis进行基本设置。
一、展点
1、将野外定位的位置点坐标及元素各数据载 入excel。 2、另存为文本文件（制表分隔符）并保存
一、展点
实用服务—投 影变换—P投影 转换—U用户文 件投影转换， 打开保存的txt 文件。设置如 右表：
MAPGIS 初级培训 讲稿
（绘制等值线图）
成都理工大学地球化学专业 张峻基、周姗姗
概述
Mapgis绘制等值线图与sufer原理大致相似，都 是以网格化数据为基础。但mapgis此功能是 建立于标准位置点基础之上。 因此mapgis绘制等值线图分为两个步骤： 1、图上展点； 2、绘制等值线图。
一、展点
一、展点
用户投影参数设置：
一、展点
结果投影参数设置：
一、展点
设置分隔符：
一、展点
生成如下窗口，并保存，展点结束。
二、绘制等值线图
1、空间分析—DTM分析—文件—打开数据文 件—点数据文件（打开展点中生成的点文件） 2、GRD模型—离散数据网格化（打开最初 生成的txt）出现如下窗口：
二、绘制等值线图
X、Y为坐标，Z为元素值。一般使用Kring泛 克立格法网格化。设置结束后点击确定。
等值线图的制作流程课件
84
点属性标注
点属性标注 标注属性选择 可选标注域名 （字段名） 添加到文件 （保存注释的文件）
等值线图的制作流程
85
点属性标注
点属性标注 设置属性标注注释
等值线图的制作流程
86
点属性标注
点属性标注 效果
如何取消 1.关闭文件 2.不保存修改 3.重新打开
等值线图的制作流程
87
等值线图的制作流程
47
纸张大小设置
等值线图的制作流程
48
编辑图片
等值线图的制作流程
49
缓冲区分析
缓冲区分析：
包括点缓冲区分析 线缓冲区分析 区缓冲区分析
要点：
半径 缓冲区类型 图元
等值线图的制作流程
50
操作流程 1输入半径
等值线图的制作流程
51
操作流程 2选择类型
等值线图的制作流程
3375006
等值线图的制作流程
2
处理导出到 txt文档
导出：选择数据复制，粘贴到文本文件 保存为SourceData.txt，如图
等值线图的制作流程
3
2. 投影变换
准备投影变换参数： 投影前WGS84坐标 获取参考GPS参数设置 地图的投影坐标 请查看地图参数：地图的说明，
如湖南省小土壤图附投影参数说明
等值线图的制作流程
99
区属性标注
区属性标注 示例
等值线图的制作流程
100
生成裁剪边界弧段区文件、区文件
在工程中添加弧段的区文件 在区编辑中选择区编辑生成区，然后保存。
等值线图的制作流程
68
保存区文件、打开区文件
等值线图的制作流程
69
生成区，保存后从工程删除区文件
Surfer中网格化方法的选取方法解析
加权反距离法（Inverse distance to a power）
首先是由气象学家和地质工作者提出的，后来由于D．Shepard的工 作被称为谢别德法方法。 基本原理是设平面上分布一系列离散点，已知其位置坐标和属性值 ，P（x，y）为任一网格点，根据周围离散点的属性值，通过距离加 权插值求P点属性值。 实质是待插值点邻域内已知散乱点属性值的加权平均，权的大小与 待插值点的邻域内散乱点之间的距离有关，是距离ｎ次方的倒数。 加权反距离插值法认为任何一个观测值都对邻近的区域有影响，且 影响的大小随距离的增大而减小。 该方法的优点是可以通过权重调整空间插值等值线的结构，但是其 计算值容易受到数据点集群的影响，计算结果中常出现孤立点数据 明显高于周围数据点的现象。
径向基函数法（radial basis function）
又称距离基函数，是多个数据插值方法组合的一种多形式网格化方函数来 定义不同的加权方法，进行不同方式的网格化。 所有径向基函数插值法都是准确的插值器，它们都能尽量适应数据 。若要生成一个更圆滑的曲面，对所有这些方法都可以引进一个圆 滑系数。 径向基函数插值方法具有很强的拟合数据点、产生光滑曲面的能力 ，其适应范围也类似克里格法。
最近邻点法（Nearest neighbor）
是荷兰气象学家A．H. Thiessen提出的一种分析方法。 最初用于从离散分布气象站的降雨量数据中计算平均降雨量 ，GIS和地理分析中多采用其进行快速赋值。 最近邻点插值的一个隐含的假设条件是任意网格点的属性值 都是用距离它最近的位置点的属性值，用每一个网格节点的 最邻点值作为待求的节点值。 该方法适合对规则分布的数据进行网格化；或者大多数数据 点位于网格节点上；或者在一个完整的数据文件中，只有少 数点无值，可以采用该方法来填充无值的数据点。 总之，最近邻点插值法是均质无变化的，更适合于均匀间隔 的数据插值，可以有效填充无值数据区域。
Mapgis等值线图绘制方法
M a p g i s等值线图绘制方法本页仅作为文档封面，使用时可以删除This document is for reference only-rar21year.MarchMapgis等值线图绘制方法1.打开Mapgis主菜单，选择“实用服务”—“投影变换”2.选择“投影变换”—“用户文件投影变换”3.点击打开文件4.将TXT文件导入到Mapgis中5.点击“用户投影参数”6.设置用户投影参数坐标系类型设置为“投影平面直角”椭球参数设置为“西安80”投影类型为“高斯克吕格”比例尺分母为“1”坐标单位为“米”投影带类型为“6度带”投影带序号为“19”设置完毕后点击确定7.设置结果投影参数坐标系类型设置为“投影平面直角”椭球参数设置为“西安80”投影类型为“高斯克吕格”比例尺分母为“50000”坐标单位为“毫米”投影带类型为“6度带”投影带序号为“19”设置完毕后点击确定8.指定数据起始位置选择属性数据下第一行数据作为数据起始位置9.设置分隔符（1）选择“按指定分隔符”（2）选择“设置分隔符”（3）分割符号中，将（Tab键）选中（4）属性名称所在行，选择第一行的属性行（5）将数据类型都改为“5.双精度型”（6）点击确定10.设置点图元参数（1）选择点图元参数（2）修改子图号及高度宽度后点击确定11.投影变换点击投影变换后，确定生成文件12.保存文件13.在Mapgis主程序中选择“空间分析”—“DTM分析”14.加载点数据，打开“文件”—“打开数据文件”—“点数据文件”15.点击“处理点线”—“点高程点提取”选择高程属性项后，点击确定16.离散数据网格化（1）点击“Grd模型”—“离散数据网格化”（2）在网格化方法中选择Kring泛克里格法网格化（3）点击确定，生成Grd文件17. 等值线绘制（1）在“Grd模型”中选择“平面等值线绘制”（2）载入Grd文件（3）选择“等值线光滑处理”（4）制图幅面选择“原始数据范围”（5）点击确定，生成等值线18.等值线文件保存。
Surfer中网格化方法的选取方法
应用条件
数量小于２５０个点数据的网格化，对于２5０～ １０００个数据点，效果也不错 方法速度快，适合于大量（１０００个以上）数 据的网格化 适应范围类似克里格法 可以通过权重调整空间插值等值线的结构，计算 值容易受到数据点集群的影响，常出现孤立 点数据明显高于周围数据点 适合规则分布、或者大多数数据点位于网格节点 上的数据，更适合于均匀间隔的数据插值， 可以有效填充无值数据区域 方法速度快，适合中等数量、均匀分布的数据网 格化 实际上是一种趋势面分析作图程序，可用来确定 数据的大规模趋势和图案。被广泛应用于地 质科学。该方法具有速度快特点、然而其去 掉了原数据中的局部细节，不利于资料的详 细分析
11.4
12.6
13.8
15
16.2
-4.8
0.6
1.8
3
4.2
5.4
6.6
7.8
9
10.2
11.4
12.6
13.8
15
16.2 -0.8 -1.6 -2.4 -3.2 -4
50 46 42 38 34 30 26 22 18 14
G— 多 项 式 回 归
-0.8 -1.6 -2.4 -3.2 -4 -4.8 0.6 1.8 3 4.2 5.4 6.6 7.8 9 10.2 11.4 12.6 13.8 15 16.2
三角剖分法（triangulation/liner interpolation ）
是一种严格的插值方法，使用最佳的Delaunay三角形，通过 直线连接各数据点形成一系列三角形，并且所有的三角形互 不相交，每个三角形内的网格节点值由该三角平面决定。 由于采用所有的数据点去构造三角形，因而原数据能得到很 好的体现，给定三角形内的全部节点都要受到该三角形的表 面限制。 该方法速度快，适合中等数量、均匀分布的数据的网格化， 地图上稀疏区域将会形成截然不同的三角面。 当数据量足够时，该方法对断线的保留具有其他方法不可比 拟的优势。
测绘技术如何进行栅格地图绘制
测绘技术如何进行栅格地图绘制近年来，测绘技术在地理信息系统（GIS）领域得到了广泛应用。
测绘师们使用各种技术手段来绘制地图，并在栅格地图的制作中发挥了重要作用。
本文将探讨测绘技术在栅格地图绘制中的应用。
首先，让我们了解栅格地图和测绘技术的基本概念。
栅格地图是一个由像素组成的图像，每个像素代表着地球表面的一个点。
而测绘技术是通过测量和记录地球表面的各种要素，获得真实世界的准确数据。
测绘技术可以使用地面测量仪器、航空摄影、卫星遥感等多种手段来进行。
在栅格地图绘制中，首先需要获取原始数据。
这些数据可以是通过现场测量获得的，也可以是通过航空摄影或卫星遥感获取的。
原始数据通常包含地表高程、地物分布、水域边界等信息。
通过各种方式获取的数据需要进行处理和整理，以便进行后续的绘制工作。
接下来，测绘师需要将原始数据转换为栅格形式。
这涉及到数据的插值和网格化处理。
插值是根据已知点的数据推算未知点的数值。
网格化是将连续的数据分割成离散的网格。
插值和网格化的过程需要运用数学、计算机科学和统计学等知识，以保证数据的准确性和连续性。
一旦数据被转换为栅格形式，接下来的工作是栅格地图的绘制。
这可以通过计算机软件来实现。
测绘师需要使用专业的绘图软件，如ArcGIS或QGIS等，来处理和显示栅格数据。
这些软件提供了丰富的绘图和编辑工具，使得绘制栅格地图变得简单而高效。
在进行栅格地图绘制之前，测绘师还需要进行一系列的数据预处理。
这包括数据清洗、边界修正、数据压缩等步骤。
数据清洗是为了除去异常值和错误数据，以保证数据的质量。
边界修正是为了解决数据边界不连续的问题，使得地图的视觉效果更加自然。
数据压缩是为了减小地图文件的大小，提高地图的加载速度。
除了栅格地图的绘制，测绘技术在栅格地图的更新和维护中也发挥了重要作用。
地球表面的特征和地貌在不断变化，因此栅格地图需要定期更新。
测绘师需要利用最新的测量技术，获取地球表面的最新数据，并将其反映在栅格地图中。
多种软件联合实现等值线快速绘制的制图模式
（ 河南省地球物理工程勘察院， 河南 郑州 ４ ５ ０ ０ ５ ３ ）
摘 要：针对地震解释成图的特点， 提出了一种多种软件联合绘制等值线的制图模式， 它包括用 Ｅ ｘ ｃ ｅ ｌ 取数据、 用 Ｓ ｕ ｒ ｆ ｅ ｒ 网格化模型自动绘制等值线、 用Ａ ｕ ｔ ｏ Ｃ Ａ Ｄ修改等值线和用 Ｍ ａ ｐ Ｇ Ｉ Ｓ 处理成图 ４个步骤。应用该制图模式减少 计算工作量， 避免人为误差， 实现了高精度高效率地勾绘有构造的等值线。通过实例对比线性内插的三角网法、 距 离倒数乘方法、 最小曲率法、 径向基本函数法和克立金法这 ５种网格化法的特点， 说明克里金网格化模型绘制等值 线的效果最佳， 以及 Ａ ｕ ｔ ｏ Ｃ Ａ Ｄ比 Ｍ ａ ｐ Ｇ Ｉ Ｓ 在编辑修改等值线图形时的优异之处。 关键词： Ｅ ｘ ｃ ｅ ｌ ； Ｓ ｕ ｒ ｆ ｅ ｒ ； 网格化模型； Ａ ｕ ｔ ｏ Ｃ Ａ Ｄ ； Ｍ ａ ｐ Ｇ Ｉ Ｓ ； 等值线； 制图模式 中图分类号：Ｐ ６ ３ １ ； Ｔ Ｐ ３ １ 文献标识码：Ａ 文章编号： １ ０ ０ ０－ ８ ９ １ ８ （ ２ ０ ０ ９ ） ０ ４－ ０ ４ ７ ２－ ０ ５
１ 绘制等值线的制图模式
煤田地震勘探的成果图件有各地层的等 ｔ ０时 间图和等深度构造图等， 其中煤层底板等值线图是 最重要的成果图件， 现在以制作煤层底板等值线图 为例将这个制图模式详细说明如下。 １ ． １ 用 Ｅ ｘ ｃ ｅ ｌ 获取三维数据 Ｅ ｘ ｃ ｅ ｌ 具有强大的公式计算功能和 Ｖ Ｂ Ａ编程功 能， 内置了丰富的函数， 如数学与三角函数、 统计函 数、 逻辑函数等， 可以进行各种数据计算和统计。针 对地震解释的数据计算创建一些含有公式的表格来 获取数据， 例如时—深转换表格、 地层倾角计算表、 空间校正表、 速度计算表格等， 在这些表格中桩号、 时间、 速度、 深度和坐标等各自有固定的列和相应的 公式。随着地震、 地质、 钻孔等资料掌握的丰富全 面， 需要增加或调整数据时， 只需对个别数据作修改 即可得到结果， 不需重复输入公式， 非常灵活实用。 实际工作中有些计算比较复杂或专业， 用Ｅ ｘ ｃ ｅ ｌ
离散点插值方法、等值线的绘制及平滑技巧
离散点插值方法、等值线的绘制及平滑技巧离散点插值方法、等值线的绘制及平滑技巧2008-06-10 22:45由于等值线图看起来非常直观、形象,因此在天气预报、气候预测分析等方面用得非常多,已成为预报员不可缺少的工具之一。
如各等压面层的位势高度图、高空环流、温度及降水分布图等等。
目前也有一些非常好的微机用绘图软件,如SUFER、GRADS 等。
这些软件一般都只有DOS 版,在流行的WINDOWS平台上,虽然可以调用,但不能使用鼠标操作,故不如使用在WINDOWS 环境中开发的软件方便。
因此,许多希望在自己开发的应用程序中能方便地显示及打印各类等值线图的人,都想知道绘制等值线图的原理方法。
如何用格点资料绘制等值线图在文献〔1〕中已有介绍,而离散点(如气象台站) 的资料必须通过插值才能绘制等值线图。
插值的方法有几种,比如三角网插值,它是将相邻的三个点连成一个个三角形,然后用文献〔1〕介绍的追踪法或其它方法在三角形边上进行插值。
此方法的优点是需要插的值少,插值算法简单,数据处理量少。
缺点是三角网的生成随意性很大,任意四个点可生成二组不同的三角形,不同的三角网插值得出的等值线也不可能完全相同。
若人为固定三角网,当有资料缺测时,就不得不重新调整三角网,因此程序的通用性不高。
目前比较流行的是通用性好的网格化方法。
1 离散点网格化理论上,离散点网格化可采用局部曲面拟合方法,用多元回归方法建立V ( X , Y) = a + bX + c Y或二次方程V ( X , Y) = a + bX + c Y + dX2 + eY2 + f X Y所谓局部,是指采用拟合点周围一定范围的部分离散点进行拟合。
这种做法虽然在许多情况下效果不错,但通常总是有些地方与实际情况有较大出入。
因此根据人工绘制等值线时的直接内插方式,采用以下几个步骤。
111 定点(1) 根据离散点的分布范围确定整个网格的范围及网格距,网格距不宜过大或过细。
MAPGIS67化探制图教程简介
地信网MapGIS67化探制图教程简介一、内容简介地信网MapGIS67化探制图视频教程由久负盛名的地信网专家组成员、中地数码MapGIS讲师鼠马象鸡主讲，教程共分为15节课，全长近七个小时。
教程除全面介绍了在MapGis67软件平台上绘制地球化学基本图件的方法外，还特别为学员们讲授了应用高级操作技巧绘制频率统计直方图、渐变色化探剖面图、平面剖面图、重要异常剖析图以及点位数据符号图的方法，增强MapGIS67软件平台的适用性，就是对非化探相关专业的人员了解和学习软件应用技术，掌握更多的高级操作技巧，进一步提高工作效率，也有着抛砖引玉、逐类旁通的的作用，是难得的学习教程。
为配合本套教程的录制，作者还专门编制了相应的化探制图数据统计模板和化探剖面数据处理模板，有了这两套模板，学员只要把相关数据直接粘贴到模板中，就可以立即获得化探制图所需要的各项参数和投影生成化探剖面所需要的数据，教程还为学员介绍了化探剖面渐变色背景色区的制作和应用方法，大大提高了工作效率。
二、视频教程内容：第01节化探制图内容简介第02节数据资料的准备和整理第03节直方图绘制方法示例第04节数据投影及点位图生成第05节离散数据网格化处理第06节地球化学图的绘制第07节单元素异常图绘制第08节组合异常图和综合异常图第09节异常剖析图的绘制方法第10节化探图件的检查整饰第11节投影绘制化探剖面图第12节给化探剖面添加渐变色区第13节平面剖面图的绘制方法第14节点位数据符号图的绘制第15节化探数据处理辅助程序简介三、相关数据资料：教学课件：化探制图示例统计模板：化探制图数据统计模板.xls剖面模板：化探剖面数据处理模板.xls实习资料：地质地理底图（文件夹）实习资料：实习数据（1万土壤）.xls实习资料：直方图绘制（文件夹）相关软件：化探数据处理辅助程序安装包.rar相关规范：（文件夹）DZ T 0011-1991 地球化学普查规范_比例尺1：50000.pdfDZ T 0075-1993 地球化学勘查图图式、图例及用色标准.pdfDZ T 0145-1994 土壤地球化学测量规范.pdfDZ T 0167-1995 区域地球化学勘查规范1：20万.pdfDZ T 0167-2006 区域地球化学勘查规范.pdf。
等值线绘制算法
1.等值线绘制的过程分为三个步骤，离散数据网格化、等值线生成
以及等值线图的填充。
项目中的数学模型将离散数据进行了网格化，本模块需要提供等值线的生成和等值线几何体的绘制。
2.等值线的生成方法。
规则化的网格数据提供了每一个网格点上的
坐标和浓度值，基本思路是从区域边界开始，利用线性插值逐网格跟踪等值线得到等值线在其所穿过的网格边上等值点的坐标，保存这些坐标并根据一定的方位判别法对这些等值点进行排序，为下一步绘制几何体做准备。
3.等值线几何体的绘制。
一个网格内的等值线确定为直线，不需要
进行平滑处理。
首先对上一步骤中得到的各条等值线按照其起点的位置进行排序，排序的顺序与等值线追踪的顺序一致，即区域左、上、右、下、内部等值线，内部等值线按照横坐标递增排序，这样排好序的等值线按照起点坐标左、上、右、下、内部的顺序保存起来。
整个算法的基本描述如下:
1) 按起点纵坐标从下至上的顺序对起点在左边界上的等值线排序;
2) 按起点横坐标从左至右的顺序对起点在上边界上的等值线排序;
3) 按起点纵坐标从上至下的顺序对起点在右边界上的等值线排序;
4) 按起点横坐标从右至左的顺序对起点在下边界上的等值线排序;
5) 按起点横坐标从左至右的顺序对内部封闭的等值线排序;
6) 填充第三种等值线与网格下边界或左边界以及起点和终点所在的边界所围的区域。
对最后一条等值线, 则还需填充与网格上边界或右边界所围的区域;
7) 填充第二种等值线与起点和终点所在的边界以及这二边界相交的顶点所围的区域;
8) 填充第一种等值线与起点和终点所在的边界的顶点所围的区域;
9) 填充内部封闭等值线所围的区域.。
实习七 应用Surfer绘制等值线图
实验七：应用SURFER绘制等值线图1．实习目的油气地质研究涉及到地表高程、油层构造、地层厚度、储层性质等参数平面分布特征的描绘，人们常等值线图等平面图件来表示。
地质数据一般为离散分布的，等值线图的绘制过程就是将离散的数据点变换为连续的图形，以建立起区域性连续的整体概念。
计算机绘图是利用计算机所具有的高速运算能力，应用插值和拟合等数据处理方法，把数字处理过程和作图过程融为一体，提高作图的速度与效率。
需要注意的是计算机平面制图尚不能主动考虑具体的地质情况，因此，应用Surfer软件生成的图形多数情况下要以矢量方式输出并进行编辑、修改。
本次实习要求掌握应用SURFER软件绘制等值线图和输出的方法。
2. 数据准备表实7-1 某油田储层孔隙度平面分布数据。
序号X YΦ序号X YΦ1418621.00531184.008.5715418120.00532037.0012.8 2418429.00531278.0016.616417697.00531913.0010.5 3418614.00531394.0010.3717417895.00532011.008.94418094.00531290.00 3.5618418255.00532243.0017.5 5418269.00531394.0014.9719417341.00531925.00176418427.00531489.008.0920417725.00532133.00 6.47418150.00531343.0010.821418065.00532324.0012.1 8417914.00531375.00 5.3222418207.00532455.0019.2 9418431.00531701.0010.323417986.00532486.0014.5 10417897.00531586.008.424417499.00532458.0016.96 11417708.00531719.0011.825417725.00532597.0018.45 12417886.00531809.00 6.326417913.00532635.0016.03 13418067.00531920.0011.527417985.00532833.0013.9 14417439.00532191.009.428418429.00530450.007.73．实习步骤3．1 数据输入启动【Grapher】进入系统的初始画面，下拉【File】菜单，点击【New】命令，弹出【New Window】 窗口，选Worksheet项，点击【OK】打开空白Worksheet数据表。
用OpenGL绘制物化探数据等值线
用OpenGL绘制物化探数据等值线叶娜;许惠平;赵彧;覃如府【期刊名称】《世界地质》【年(卷),期】2004(023)001【摘要】在Visual C++6.0环境下调用OpenGL库用网格法绘制标准物化探数据等值线,以图示形式给出数据集的二维和三维方式的等值线图,并给出相应的数字高程模型图,方便用户分析和理解数值分布情况.用线性内差法计算所有网格边上的等值点,横边上的等值点到前一个最近网格点的距离存放在二维数组X中,纵边上的等值点到下一个最近网格点的距离存放在二维数组Y中.调用自编的搜索函数在网格边上搜索等值线的线头,用追踪函数对这个线头所在线上的点进行追踪,并把追踪到的等值点的坐标按顺序存放在二维数组DW中,用OpenGL中绘制Bezier曲线的方法分段连接DW数组中的等值点,并按颜色从红到蓝递变的顺序在屏幕上绘制等值线.【总页数】6页(P79-84)【作者】叶娜;许惠平;赵彧;覃如府【作者单位】吉林大学,地球探测与信息技术学院,吉林,长春,130026;同济大学,教育部海洋地质重点实验室,上海,200092;中国国家信息安全测评认证中心,东北测评中心,吉林,长春,130061;吉林大学,地球探测与信息技术学院,吉林,长春,130026【正文语种】中文【中图分类】P631.84;TP391【相关文献】1.用Surfer软件绘制等值线图中地质数据的处理方法 [J], 陈练武;牛菲2.基于LiDAR点云数据的等值线绘制方法 [J], 何睦;李梦珍3.多种物化探数据等值线投影分层图程序 [J], 杨生;朱添宝4.基于OpenGL的三维等值线绘制实现 [J], 成诗明;杨光5.基于OpenGL的自适应网格法的等值线绘制与填充实现 [J], 邓枚;周汝良因版权原因，仅展示原文概要，查看原文内容请购买。
Surfer化探元素等值线图绘制技巧初探之离散数据网格化方法选择及白化方法简介
①简单平面（Simple planer surface）： z(x,y)=A+Bx+Cy ②双线性鞍形（Bi-linear saddle）： z(x,y)=A+Bx+Cy+Dxy ③二次表面（Quadratic surface）： z(x,y)=A+Bx+Cy+Dx2+Exy+Fy2 ④三次表面（Cubic surface）： z(x,y)=A+Bx+Cy+Dx2+Exy+Fy2+Gx3+H x2y+Ix y2+Jy3 或由用用户自定义。 选定的曲面方程相应显示在下面；右面的参数框内 则显示X、Y和总的最高项次。 用户也可以利用Parameters框自定义多项式方程。
数据度量用来提供有关的数据信息。根据 度量所得的数据资料,可以再次利用一个网 格数据网格的其他方法。
它其实不是一种插值方法,它是一种数据的 度量方法,通过这种方法可以找到比较合适 的插值方法。
最小曲率法的对话框包括以下内容： 最大残差（Max Reciduals）：单位与数据的相同，比较 合适的值是数据精度的10%。缺省的最大残差为 0.001 (Zmax - Z min) 最大重复参数（Max lterations）：通常设为网格结点数 的1到2倍。例如，对于50×50的网格，最大重复参数在 2500与5000之间。 内部和边缘张性系数（Internal and Boundary Tension）： 设定弹性薄板内部和边缘弯曲度的参数。该值愈大，弯曲 愈小。缺省值均为0。 松弛系数（Relaxation Factor）：算法参数，通常，该值 愈大，迭代算法会聚愈快。缺省值为1，一般不用另设定。
离散数据网格化参数的确定和数学模型的选择--以Sufer 7.0、Mapgis 6.0为例
离散数据网格化参数的确定和数学模型的选择--以Sufer 7.0、
Mapgis 6.0为例
高艳芳
【期刊名称】《地质与勘探》
【年(卷),期】2002()z1
【摘要】使用计算机绘制等值线首先要对离散数据进行网格化处理.针对物化探数
据的特点,详细地讨论了利用DEM模型进行网格化时如何设定所需的系列参数,同
时对采用的数学模型进行了应用比较.
【总页数】4页(P139-142)
【关键词】等值线;网格化参数;数学模型
【作者】高艳芳
【作者单位】
【正文语种】中文
【中图分类】P628
【相关文献】
1.离散的物探数据网格化方法的对比及选取——以高精度磁测数据为例 [J], 姚文;
李琼
2.MAPGIS化探图数据网格化方法选择探讨 [J], 吴现兴;李余晖;王海丰
3.MapGIS与AutoCAD、Sufer．CorelDraw之间数据转换的特点及问题 [J], 郭
丽
4.地球物理数据网格化参数的确定及模型的选择 [J], 王玉敏;冯宁;姚敏
5.MAPGIS高程模型网格化参数的确定及网格化模型的选择 [J], 伍卓鹤;赵献军;陈明
因版权原因，仅展示原文概要，查看原文内容请购买。
离散数据网格化参数的确定和数学模型的选择--以Sufer 7.0、Mapgis 6.0为例
增 刊P P P P P P P P P P P P 高艳芳： 离散数 据网格化参数的确定和数学模型的选择P P P P P P P P P P P P P P P
差， 使原始数据中包含的客观特征能无明显损失地 传递到内插计算的数字模型中去， 是对所有的数学 模型的最基本的要求。 可用来对离散数据进 行网格化的数学模型很 多。在 !"#$%&’ ( 中， 给定了 ) 种数学模型， 分别是： 距离 倒 数 （ *+,$%-$ ./-01+2$ 03 1 435$% ） 、 克 里格 （ 6%/+7/+7 ） 、 最小曲率 （ 8/+/9"9 :"%,10"%$ ） 、 最近临 近点 （ ;10"%1< ;$/7=>3% ） 、 自然临近点 （ ;10"%1< ;$/7=? >3%） 、 多项式回归 （ @3<A+39/1< B$7%$--/3+） 、 径向基本 函数 （ B1./1< C1-/- D"+20/3+ ） 、 改进谢别德 （ 83./#/$. !=$41%.E- 8$0=3. ） 、 三 角网 F 线 型插值 （ G%/1+7"<10/3+ 5/0= H/+$1% *+0$%43<10/3+ ） 。8147/-I’ ( 也给出了 J 种 之多， 分别是： 距离倒数、 （泛） 克里格法、 稠密数据 中值选取法、 稠密数据高斯距离权法。不同的数学 模型有各自的特点， 不同类型的数据可以选择相应 的数学模型。 距离倒数 、 克里格 、 径向基本函数、 最近点这 J 种数学模型可以采用 搜索方式设定参数 进行网格 化。其他的几种进行网格化插值计算时， 采用的是 全部数据进行网格插值计算。尽管所有的 KL8 数 字高程模型类的应用程序提供给用户的网格化的方 法很多， 但并不是所有都能满足物化探数据处理的 要求。 距离倒数是一种加权平均插值法， 可以进行确 切的、 圆滑的插值。多年来， 在物化探数据处理过程 中被广泛采用。它的特征之一是在网格区域内产生 围绕观测点位置的 “ 牛眼” 。 “ 牛眼” 可以通过改变 圆滑参数进行消除。 普通克里格法或泛克里格方法是一种适合于许 多领域的网格化插值方法， 目前应用也比较广。该 方法在数学上可对所研究的对象提供一种最佳线性 无偏估计， 同时也兼顾了数据空间分布特征， 其内插 的结果可信度高。它能够比较精确地表示出隐含在 原始数据中的趋势， 避免产生距离倒数中过多 “牛 眼” 点。 径向基本函数是多个数据插值方法的组合。其 中的复二次函数被许多人认为是最好的方法， 既能 够反映原始数据的特征， 又能光滑等值线。网格化 结果相似于克里格插值结果。使用它来制作地球化 学图， 效果最好。 最近临近点法， 严格地说， 它不是一 种插值 法， 它只是将距离网格最近的点的 安插到网格 的节点上去。这个 法是针对 点均匀分布而且 空值点又少的原始 来说， 这样的 不需要插
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
加权反距离插值法是一种非常快速的网格化方法，在小于500个数据点时，可 以使用No Search(使用所有点)的搜索类型来快速生成网格。
ห้องสมุดไป่ตู้
其中hij 网格节点“j”与邻近点“i”之间的有效分离距 离。
Zj 网格节点“j”的内插值。 Zi 邻近点。 dij网格节点“j”与邻近点“i”之间的距离。 β 权重系数。 δ Smoothing参数。
设置反距离权插值高级选项Power和Smoothing参数：
权重系数Power确定随着数据点到网格节点距离的增加， 其权重降低的程度。随着power逼近0，生成的表面逼近 一个水平面，该平面通过数据文件中的所有观测点的平均 值。随着权重系数的增加，生成的表面由最邻近点插值， 导致表面变成多边形。多边形表现了最接近内插节点的观 测表面。可接受的权重系数通常在1和3之间。
项目实施过程中，化探手段是综合找矿方法中的重要的研 究手段之一
化探元素异常等值线图对于揭露研究区各元素异常分布特 征具有重要的指示意义。
一般实施采样方法为不规则测网法。 绘制化探元素异常等值线图时，要根据客观环境特征
和数据本身的特点，选择合适的网格化方法 网格化方法的特征及应用条件 网格化概念——是指通过一定的插值方法，将稀疏的、
克里金法中包含了几个因子：变异图模型，漂移类型 和 矿块效应。
其中变异图模型（Variogram Model）是用来确定插值 每一个结点时所用数据点的邻域，以及在计算结点时给予 数据点的权重。 Surfer提供了多种最常用的变异图模型，它们是指数、高 斯模型、线性、对数、矿块效应、幂、二次模型、有理数 二次模型、球面模型和波（空洞效应）。如果拿不准用哪 一种变异图，可选用线性变异图，大多数情况下，效果较 好。
里格法网格化、稠密数据中指选取网格化、稠密数据高斯距离权网格 化。
Surfer中网格化方法有12种。Surfer网格化方法基本涵盖Mapgis中 网格化方法。
Surfer中网格化方法——加权反距离法（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）。
算矿产储量的方法：按照样品与待估块段的相对空间位置和相关程度 来计算块段品位及储量，并使估计误差为最小。后来，法国学者马特 隆对克里格法进行了详细的研究，使之公式化和合理化。 克里格法的基本原理是根据相邻变量的值（如若干样品元素含量值）， 利用变差函数所揭示的区域化变量的内在联系来估计空间变量数值。 该方法总是尽可能地去描述原数据所隐含的趋势特征，以区域化变量 理论为基础，以变差函数为主要工具，在保证研究对象的估计值满足 无偏性条件和最小方差条件的前提下求得估计值。 对于高值数据点会使之沿某一“脊”分布，而不围绕该点孤立插值， 不形成“公牛眼”等值线。 克里格法极为灵活，广泛地应用于各个科学领域，适于各种类型的离 散数据，网格化精度高，是极佳的网格化方法。
Surfer应用更广泛，下面详细介绍Surfer中各种网格化方法的特征。
首先由气象学家和地质工作者提出的。
加权反距离插值法是一种加权平均内插，可以是准确插值或平滑插值，通常表 现为准确插值。
基本原理是设平面上分布一系列离散点，已知其位置坐标和属性值，P（x,y） 为任一网格点，根据周围离散点的属性值，通过距离加权插值求P点属性值。
线性变异图的克里金法应用广泛，效果较好。
这是一种在地学中广泛应用的网格化方法。由 最小曲率法构成的插值表面像一个线性弹性薄 板，试图在尽可能严格地尊重数据的同时，生 成与原始数据点尽可能吻合的最平滑的曲面。
最小曲率法不是准确插值，是典型的平滑插值。
实质是待插值点领域内已知散乱点属性值的加权平均，权的大小与待插值点的 领域内散乱点之间的距离有关，是距离n次方的倒数。
加权反距离插值法认为任何一个观测值都对邻近的区域有影响，且影响的大小 随距离的增大而减小。
在计算一个网格节点的Z值时，赋给数据点的权重是分数，所有数据点的权重 和为1。权重与数据点到节点距离成反比，愈靠近节点的原始数据点，其权重 愈大。
不规则分布的数据插值加密为规则分布的数据，以适 合绘图的需要。 原始数据的不规则分布，造成缺失数据的“空洞”。 网格化则是用外推或者内插的算法填充了这些“空 洞”。
主要等值线图绘制软件Surfer及Mapgis Mapgis中4种网格化方法，距离幂函数反比加权网格化、Kring泛克
平滑参数Smoothing把“不确定性”因素与用户输入的 数据联系起来，平滑参数愈大，计算相邻网格节点Z值时 特异观察点的绝对影响就愈小。平滑参数大于0，则没有 任何一个数据点对于某个节点的权重为1，即使该数据点 正好位于网格节点上。
克里金法是一种在许多领域都很有用的地质统计格网化方法。 最初是由南非金矿地质学家克里格根据南非金矿的具体情况提出的计