Surfer做等值线图

从事煤炭资源勘查工作中，需要做很多等值线图，周围同事多是用传统的三角网插值方法，一个钻孔一个钻孔地量取之间的距离，然后在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乂丫22.0，如下图：
点击 Output Format :
utput Format
V Chain conhguous ：?egments. [Tolerance in dxfZxyw.in ；, Chainl
「 Entity Headers Specify -Output Data Format
Raw
Moye to (0), Draw to (1) 广 ID.Xy(Z]
* Blanking (Surfer .BLN]
Blank I riside [1] * Blank Outride (0)
口区 I Cancel
选择输出格式为Blank Outside(0)，即可输出为标准的bln 文件。

第二步：网格化数据
1、在两个数据文件建立好了后，打开Surfer ,新建图形文档。

2、依次点击网格一一数据，在出现的对话框中选择刚才建立的钻孔数据文件*.txt ,出
现如下界面：
i ,2D (ignore
1i. valued
Decimal Places
Description:
A Header line precei header line contains folloyj, and an inside Vertices are listed or Example: 4.0 400,200,12 .300,120,12 300.150,12
250.150,10
出网格文件
码沫灌序媒犷普查'煤质1等值线图而向U
这里有几点需要注意：
①数据列中，X、丫为生成等值线图的横轴、纵轴，与CAD和数学中一致，但是与测绘中相反。

Z为需要描述的目标数据列，根据需要可以选择D、E、…等列。

②网格化方法很多，网上书上多推荐克里格法（kriging）。

③“网格线素几何学” ——这是我要说的重点。

这里的X方向、Y方向的最大最小值，圈定了需要网格化区域的范围。

我用了边界拐点中的最大最小坐标，网格化的结果是一个比矿区边界还大的矩形区域。

如果不更改，最终网格化的区域是以钻孔坐标最值圈定的矩形区域，生成的等值线图不能充满矿区范围。

间距/行数一一我选择的5〜20,个人觉得这个能保证精度需求。

3、点击确定，生成网格化数据文件*.grd。

第三步：白化
注意：白化要在绘制等值线图之前进行，我后面用到的截图，有些是未先进行白化而成的，悲剧了。

点击网格一一白化一一选择刚生成的*.grd文件打开一一选择刚生成的*.bln文件打开——保存新的grd文件。

第四步：绘制等值线图
依次点击地图一一等值线图一一新建等值线图一一选择刚生成的新的grd文件一一打开一一自动绘制出了等值线图。

如下图（已白化）：
第五步：处理等值线图
双击图形，会出现属性窗口。

在等级标签下点击等级（很蛋疼的一个按钮，找了很久都 没找到），可调整相邻等值线差值和最值，如下图:
常规 等级
此 5 国
28.5 29 29.5 30 30.5
31
31.5
和
双击等值线图四周的坐标刻度，可调整坐标的显示属性。

第六步：加钻孔
1、加孔号
点击地图—一张贴图—一新建张贴图—一选择第一步生成的*.txt 文件打开—一在左侧 列表中选中Post ——双击图形——常规标签中工作表列，X 、Y 坐标对应列B 、C ——缺省 符号选择圆形—一标注标签工作表列选列A ，位置上齐一一确定。

这里，发现生成的图形出现两个坐标轴，且钻孔位置偏离实际位置。

（下图未经白化）
3&364000 3^3^6000 3&S6BOOO 3i370000 35372000 35374000 55376000 35378000 3538QC
3256000-
3262000-
326GOOO-
3259000-
馀⑤
35364000 3536600Q 35368000 35370000 35372000 35374000 35376000 35378000 3538000
先按F2键（选中所有图层），点击地图一一覆盖地图。

然后就好了。

3262000-
3260000-
33卷4收 35366000 3536so0。

35370000 353Z200Q 35374000 353砍00 353了B000 353800
第八步：输出为Dxf 文件
由于工作要求，我需要把生成的等值线图转换成CAD 能识别的Dxf 文
件。

点击文件—一输出—一选择Dxf 格式，点保存，出现如下对话框： 35364Q00 353fifiOQO 3536S&QQ 35370000 3E372000 35374。

口0 35M76。

皿 B537BGQ0
2、加钻孔数据
除了编号，钻孔还需要数据（如灰分、硫分等）。

只需重复前面的方法，再次生成张贴 图，把标注用工作表列更改为需要的列，符号相对位置改成右齐，并再次覆盖。

如下图（未 白化）：
3253000-
点击地图—一基面图—一选择刚才生成的*.bln 文件打开—一按F2键一一点击地图 覆盖地图。

如下图（已白化）：
的364。

口0 35^66000 336800。

35^70000 35372口g 35374000 55376000 353双Q0 35僦⑭匚
第七步：加边界
3262000-
326COOO-
3256000-
3266000-
28 73,
7
AutoCAD DXF 嵬出-等期图Ad,dxf
点击默认，并确定。

在CAD中打开刚生成的Dxf文件，发现已校正。

Ps：新生成的Dxf文件中，钻孔非圆，而是由多段线圈成，文字也为图案填充，不符合平
时的作图习惯，且方便更改。

反思自己的作图步骤，发现从第四步绘制等值线完成，就可以直接转换成Dxf文件，余下的在CAD中完成即可。