第六章_等值线自动绘制技术SURFER使用2

1、计算机地质制图计算机地质制图 讲授：陈练武讲授：陈练武第六章第六章 第六章第六章 课 时 数：8学时本章目的： 了解计算机自动生成等值线图的方法。掌握SURFER软件的使用方法。本章难点： 等值线图生成软件SURFER的使用方法教学手段：多媒体教学方法：课程讲授其它使用技巧其它使用技巧（1）空白和空白文件空白和空白文件由离散的由离散的XYZ数据文件生成网格文件时，数据文件生成网格文件时，SURFER将根据将根据原始原始X，Y的取值范围和所选用的数学模型，自动生成一的取值范围和所选用的数学模型，自动生成一个矩形网格，其网格边界为最外数据点的范围。如据表个矩形网格，其网格边界为最外数据点的范围。如

2、据表7.1数据所生成的网格如图所示，其中数据所生成的网格如图所示，其中X坐标最少值坐标最少值36 6 2 2 5 7 2 . 4 ， 最 大 值， 最 大 值 3 6 6 2 5 2 2 3 . 8 ， 网 格 间 距， 网 格 间 距26.78181818m(共共100等分等分)，其中，其中Y坐标最少值坐标最少值3 8 9 1 7 3 6 . 5 ， 最 大 值， 最 大 值 3 8 9 3 8 5 2 . 8 ， 网 格 间 距， 网 格 间 距26.78860759m(共共80等分等分)。但在实际工作中，由于某些区域缺少原始数据或由于其它但在实际工作中，由于某些区域缺少原始数据或由于其它

3、原因，有必要由规则的网格中剔除一个或多个由封闭多边原因，有必要由规则的网格中剔除一个或多个由封闭多边形定义的区域，被剔除的区域形成空白。在画图时形定义的区域，被剔除的区域形成空白。在画图时, 等值等值线图上空白部分的等值线被消除掉；另外研究区以外（如线图上空白部分的等值线被消除掉；另外研究区以外（如井田边界）的等值线也需要去掉。这时就需要空白网格。井田边界）的等值线也需要去掉。这时就需要空白网格。空白区域由空白文件定义，一个空白文件可以定义一个或空白区域由空白文件定义，一个空白文件可以定义一个或几个空白区域。几个空白区域。 格式如下：格式如下： length,flag Pname 1x1,y1

4、x2,y2.xn,ynlength,flag Pname 2x1,y1x2,y2.xn,yn length,flag Pname 1x1,y1x2,y2.xn,ynlength,flag Pname 2x1,y1x2,y2.xn,yn 空白文件为多段文件，每段一个实体（点、线、面）。每段文件空白文件为多段文件，每段一个实体（点、线、面）。每段文件由由A,B两列数据组成，首行为标志行，标志行两列数据组成，首行为标志行，标志行A列值为列值为1时，为点时，为点数据，大于数据，大于1为线数据，大于为线数据，大于3为线或多边形数据。标志行为线或多边形数据。标志行B列值列值为为1时，表示多边形内部的区域被

5、空白；标志行时，表示多边形内部的区域被空白；标志行B列值为列值为0时，表时，表示多边形外部的区域被空白。由第示多边形外部的区域被空白。由第2行开始依次为各顶点的行开始依次为各顶点的X,Y坐坐标。每段第二行与最后一行的标。每段第二行与最后一行的X,Y坐标相等时，为封闭多边形。坐标相等时，为封闭多边形。 含有空白区域的网格文件构筑方法如下：含有空白区域的网格文件构筑方法如下：（1）按常规方法先生成一个未考虑空白区域的网格文件）按常规方法先生成一个未考虑空白区域的网格文件.GRD；（2）根据需要，按以上格式用电子工作表建立一个空白文件根据需要，按以上格式用电子工作表建立一个空白文件.BLN；（3）由

6、由GridBlank 按提示先打开未考虑空白区域的网格文件按提示先打开未考虑空白区域的网格文件.GRD，再按提示输入空白文件再按提示输入空白文件.BLN 名，指定输出文件路径和名，指定输出文件路径和文件名（为包含空白文件的新网格文件文件名（为包含空白文件的新网格文件.GRD）。）。从图形菜单中选择等值线命令或线网图命令来查看这新的网格文从图形菜单中选择等值线命令或线网图命令来查看这新的网格文件，就可看到空白的效果。件，就可看到空白的效果。下面以表下面以表7.1数据为例详细说明。如果该井田的边界由数据为例详细说明。如果该井田的边界由11个拐点构成，要求井田边界以外不画等值个拐点构成，要求井田边界

7、以外不画等值线。另外，在井田内有一由线。另外，在井田内有一由5个拐点构成的多边形区域，该区域内不能画等值线。具体操作步骤如个拐点构成的多边形区域，该区域内不能画等值线。具体操作步骤如下：下：第一步，首先将两个区域的边界数据点建立一个空白文件（如第一步，首先将两个区域的边界数据点建立一个空白文件（如BJ.BLN），），每个区域的最后一个数每个区域的最后一个数据必须是第一个点的重复。这样，第一个区域共据必须是第一个点的重复。这样，第一个区域共12个点，第二个区域共个点，第二个区域共6个点，数据输入格式如下。个点，数据输入格式如下。 12,036622631.99,3893771.463662354

8、9.94,3893774.9436623809.31,3893687.3236624519.53,3893679.7836625207.12,3893832.9636625199.58,3891774.2536624189.37,3891778.8936624186.47,3891954.1336623565.03,3891958.7736623558.06,3891800.9436622622.10,3891829.9536622631.99,3893771.46 6,136623135.54,3893378.5736623685.57,3893247.8436623706.30,38927

9、74.3336623408.16,3892573.4536622988.87,3892724.9136623135.54,3893378.57第二步，用原始数据生成一个网格文件（如第二步，用原始数据生成一个网格文件（如DATA.GRD）.第三步，选择第三步，选择GridBlank ，按提示先打开未考虑空白区域的网格文件按提示先打开未考虑空白区域的网格文件DATA.GRD，再按提示输入空白文件再按提示输入空白文件BJ.BLN 名，指定输出空白后的网格文名，指定输出空白后的网格文件（默认为件（默认为OUT.GRD）。）。如果要查看空白效果可选择如果要查看空白效果可选择GridGrid Node E

10、ditor。第四步，选择第四步，选择MapContour Map选项，用选项，用OUT.GRD重新生成等值线图。重新生成等值线图。第五步，选择第五步，选择MapBase Map选项，根据提示输入空白文件选项，根据提示输入空白文件BJ.BLN，按按F2功功能键选中所有图形对象，选择能键选中所有图形对象，选择MapOverlay Maps选项。选项。 （2）粘贴数据点粘贴数据点粘贴命令用数据点位置和标记建立一粘贴图形。粘贴图能覆盖等粘贴命令用数据点位置和标记建立一粘贴图形。粘贴图能覆盖等值线图或线网图，由值线图或线网图，由MapPost MapNew Post Map或由工具栏选取被粘贴的数据文件

11、后，显示粘贴图对话框。由或由工具栏选取被粘贴的数据文件后，显示粘贴图对话框。由2个个Tab组成。组成。 General被粘贴的数据文件是一包含了被粘贴的数据文件是一包含了(X,Y) 坐标以及符号或数据标记列的数据文件。坐标以及符号或数据标记列的数据文件。X 规定粘贴点的规定粘贴点的X坐标在工作表中所用列。坐标在工作表中所用列。Y 规定粘贴点的规定粘贴点的Y坐标在工作表中所用列。坐标在工作表中所用列。Symbol指定一包含字体和符号信息的描述符号位置的数字列。如果这列是指定一包含字体和符号信息的描述符号位置的数字列。如果这列是 None, 那么所有的粘贴点使用缺省符号。例如那么所有的粘贴点使用缺

12、省符号。例如, 如果对粘贴点使用缺省符号的如果对粘贴点使用缺省符号的符号集，在符号集中第一个符号给符号位置符号集，在符号集中第一个符号给符号位置“0”, 第二个符号给符号位置第二个符号给符号位置“1”等等。在符号列内能使用符号位置和空白行的组合。当指定符号位置时等等。在符号列内能使用符号位置和空白行的组合。当指定符号位置时, 就使就使用该位置的专门符号；没有指定符号位置时，使用缺省符号。用该位置的专门符号；没有指定符号位置时，使用缺省符号。 粘贴图允许符号列使用粘贴图允许符号列使用Symbol set:index 这种格式来选择符号集。缺省符号集这种格式来选择符号集。缺省符号集为为“Defau

13、lt Symbols” 。对缺省符号索引从。对缺省符号索引从0开始。开始。例如例如, Arial:65 将使用将使用Arial字体的字母字体的字母 A.Angel 规定了包含角度值的列号，符号逆时针方向旋转为正。规定了包含角度值的列号，符号逆时针方向旋转为正。由由Default Symbol对话框点击符号键打开符号属性对话框。可选择缺省符号集对话框点击符号键打开符号属性对话框。可选择缺省符号集（缺省符号集为（缺省符号集为Default Symbols）、）、选择缺省符号，设定其颜色属性。选择缺省符号，设定其颜色属性。粘贴点间隔设置粘贴点间隔设置 为为1时时, 工作表的所有数据点都被粘贴在图上。

14、如果设置为工作表的所有数据点都被粘贴在图上。如果设置为5，则从工作表中每隔五个点取一数据用于在图上粘贴。则从工作表中每隔五个点取一数据用于在图上粘贴。符号大小设置符号大小设置设置粘贴在图上的符号大小。可选择相等大小或与标志值大小成比例。后一种设置粘贴在图上的符号大小。可选择相等大小或与标志值大小成比例。后一种情况下，打开比例对话框，可选择比例方式：线性比例或与标志值的平方根成情况下，打开比例对话框，可选择比例方式：线性比例或与标志值的平方根成比例。可以设置符号高度与标志值大小的对应关系等。比例。可以设置符号高度与标志值大小的对应关系等。 Labels设置注记内容所在列、注记位置、字体等内容。设

15、置注记内容所在列、注记位置、字体等内容。 Worksheet Column for Labels注记内容所在列。注记内容所在列。Position标志位置设置标志位置设置 设定粘贴的标志相对于粘贴点的位置。下拉式列表框可设定粘贴的标志相对于粘贴点的位置。下拉式列表框可设定标志的位置，如选择自定义选项设定标志的位置，如选择自定义选项,可指定可指定X或或Y的偏移量。的偏移量。还可设置标志的角度等。还可设置标志的角度等。在立体图中，标志平面选择。在立体图中，标志平面选择。3D标志线标志线 当图片倾斜（当图片倾斜（90除除 图图15.2.5 注记格式设置对话框注记格式设置对话框外）或迭合粘贴图到线网图时

16、，用到外）或迭合粘贴图到线网图时，用到3D标志线组对话框。标志线组对话框。3D标志线为由线网图标志线为由线网图上某点到该点标志的垂直连线。上某点到该点标志的垂直连线。 （3） 比例尺用来在图上绘制一个比例尺，比例尺按自定义参数分为相等间隔区间。它是图用来在图上绘制一个比例尺，比例尺按自定义参数分为相等间隔区间。它是图形比例的直观表示法。形比例的直观表示法。激活图形，由激活图形，由Map Scale Bar 打开比例尺对话框。打开比例尺对话框。Number of Cycles 设定比例尺使用的周期数。例如设定比例尺使用的周期数。例如, 如果你规定如果你规定6, 比例尺就分成比例尺就分成六等份六等

17、份Cycles Spacing 周期间隔，设定在每个周期内的数据单位值。例如周期间隔，设定在每个周期内的数据单位值。例如,如果设定为如果设定为5，则周期的比例尺总长将是则周期的比例尺总长将是30 单位。单位。Label Increment 标志增量标志增量, 通常与周期间隔值是相等的。也可使用与轴不同的通常与周期间隔值是相等的。也可使用与轴不同的单位，例如单位，例如, 你的图形使用厘米为轴刻度，而比例尺可以是千米表示。你的图形使用厘米为轴刻度，而比例尺可以是千米表示。 （4） 由函数关系生成由函数关系生成.GRD文件文件GridFunction 打开网格函数对话框。在函数编辑框中输入数学打开网

18、格函数对话框。在函数编辑框中输入数学函数函数Z= f(X,Y)。数学函数逐点的计算每一个数学函数逐点的计算每一个Z值，并将其写入到值，并将其写入到输出的网格文件中。其中网格结点的数目取决于所选定的输出的网格文件中。其中网格结点的数目取决于所选定的 X和和 Y的最小值、的最小值、 最大值和增量。所采用的数学函数中能使用系统最大值和增量。所采用的数学函数中能使用系统 的任的任何内部函数。何内部函数。例如例如, 按给定的函数关系，按给定的函数关系，X、Y的取值范围和增量产生一个网格的取值范围和增量产生一个网格文件。文件。缺省的输出网格文件名为缺省的输出网格文件名为OUT.GRD，可改变文件的路径和名

19、字可改变文件的路径和名字。 （5）网格代数处理）网格代数处理由具有相同由具有相同X,Y网格尺度的两个网格文件按照一定的函数关系网格尺度的两个网格文件按照一定的函数关系C= f(A,B) 生成一个新的输出网格文件。这里生成一个新的输出网格文件。这里A和和B是输入的网格文件，是输入的网格文件，C是输出网格文件。从输入网格文件的相应结点上是输出网格文件。从输入网格文件的相应结点上(具有相同具有相同X,Y坐坐标的网格结点标的网格结点) 计算函数值，计算的结果送到具有相同坐标的输出计算函数值，计算的结果送到具有相同坐标的输出网格文件内。网格文件内。网格数学也能用于单一的网格文件。例如网格数学也能用于单一

20、的网格文件。例如,使用函数使用函数C= log10(A)，则网格文件则网格文件C的每一个结点的值将为网格的每一个结点的值将为网格A相应的每一网格结点值相应的每一网格结点值的以的以10 为底的对数。也可以进行简单的数学操作为底的对数。也可以进行简单的数学操作,例如例如 C=A-100。这时输出网格文件这时输出网格文件C比输入网格文件比输入网格文件A低低100 单位。单位。由由GridMath打开网格代数对话框。输入网格文件打开网格代数对话框。输入网格文件A组指定执行组指定执行网格数学操作的第一个网格文件。输入网格文件网格数学操作的第一个网格文件。输入网格文件B组指定用于网格组指定用于网格数学操作

21、的第二个网格文件。如果只想在一个网格文件上执行操数学操作的第二个网格文件。如果只想在一个网格文件上执行操作作, 可使其为空白。第二个网格文件应包含与第一个网格文件相同可使其为空白。第二个网格文件应包含与第一个网格文件相同的网格行、列数的网格行、列数, 和相同的和相同的X和和Y列。缺省的输出网格文件名为列。缺省的输出网格文件名为MATHT.GRD，可改变文件的路径和名字。底部的编辑框用于输可改变文件的路径和名字。底部的编辑框用于输入函数入函数C= f(A,B)。 。 （6）计算体积和面积）计算体积和面积由由GridVolume选取要进行体积计算的网格文件，显示选取要进行体积计算的网格文件，显示G

22、rid网格网格体积对话框。体积对话框。Upper Surface组：可选定一网格文件作为上表面（曲面）；或设组：可选定一网格文件作为上表面（曲面）；或设定上表面的定上表面的Z为常数值（水平平面）。为常数值（水平平面）。Lower Surface组：可选定一网格文件作为下表面（曲面）；或设组：可选定一网格文件作为下表面（曲面）；或设定下表面的定下表面的Z为常数值（水平平面）。为常数值（水平平面）。在选定上表面的和下表面之后在选定上表面的和下表面之后, 点击点击OK执行体积和面积计算。结执行体积和面积计算。结果显示在编辑器窗口，可保存成果显示在编辑器窗口，可保存成ASCII文本文件文本文件, 或拷

23、贝和粘贴到或拷贝和粘贴到绘图窗口。绘图窗口。 例例1：计算层状矿体储量：计算层状矿体储量在层状矿体的勘探中，每个钻孔都可以得到矿层顶、底板的三维在层状矿体的勘探中，每个钻孔都可以得到矿层顶、底板的三维坐标，即坐标，即X,Y,Z值。由勘探网的多个钻孔资料，可以得到代表了层值。由勘探网的多个钻孔资料，可以得到代表了层状矿体上层面和下层面层面的两组数据。分别用状矿体上层面和下层面层面的两组数据。分别用Surfer进行网格化进行网格化处理，就能得到代表上、下表面的两个处理，就能得到代表上、下表面的两个.GRD文件。利用体积计文件。利用体积计算功能，可以得到该矿体的近似体积，乘以矿石的体重，便可求算功能，可以得到该矿体的近似体积，乘以矿石的体重，便可求出该矿体的矿石储量。出该矿体的矿石储量。例例2：计算湖水容量：计算湖水容量湖底是一个下凹的曲面，通过一定数量的定位测深，用湖底是一个下凹的曲面，通过