利用Surfer绘制双色平面剖面图

自动绘制的可视化处理技术， 提出了彩色渐变填充的 ２ ０ ０ ６ ） 实现了用 处理方 法； 马 培 仙、 李 百 祥、 陈卫东（ Ｓ ｕ ｒ ｆ ｅ ｒ 绘制单色的平剖图； 李文杰、 李军峰、 孟庆 敏 等 （ ２ ０ ０ ７ ） 研究了将航空物探数据转换为 Ｇ Ｒ Ｄ网格化数 据绘制平剖图的方法。本文在研究前人资料基础上， 介绍采用图形分解的方式， 运用 Ｓ ｕ ｒ ｆ ｅ ｒ 基面图分别绘 制 “ 正值区域” 、“ 负值区域” 的方式进行平剖图的绘 制方法， 作者经过多次试验应用， 取得了满意的效果。
Ｘ Ｔ×ｓ ｉ ｎ （ ） φ ｔ ＝Ｘ－Δ Ｙ Ｔ×ｃ ｏ ｓ （ ） φ ｔ ＝Ｙ＋Δ
Ｏ ｐ ｅ ｎＤ ｉ ｒ １ Ｐ ａ ｔ ｈ＋″ ＼ ″ ＋″ 正剖面 ｂ ｌ ｎ ″ Ｆ ｏ ｒＯ ｕ ｔ ｐ ｕ ｔＡ ｓ＃ ２ ′ 正剖面文件 Ｏ ｐ ｅ ｎＤ ｉ ｒ １ Ｐ ａ ｔ ｈ＋″ ＼ ″＋″ 负剖面 ｂ ｌ ｎ ″ Ｆ ｏ ｒＯ ｕ ｔ ｐ ｕ ｔＡ ｓ＃ ６ ′ 负剖面文件 Ｏ ｐ ｅ ｎＤ ｉ ｒ １ Ｐ ａ ｔ ｈ＋″ ＼ ″ ＋Ｆ ｉ ｌ ｅ １ Ｌ ｉ ｓ ｔ （ ｉ ）Ｆ ｏ ｒＩ ｎ ｐ ｕ ｔＡ ｓ＃ １ ′ 打开数据文件 Ｏ ｐ ｅ ｎＤ ｉ ｒ １ Ｐ ａ ｔ ｈ＋″ ＼ ″ ＋″ ｌ １ ｔ ｘ ｔ ″ Ｆ ｏ ｒＯ ｕ ｔ ｐ ｕ ｔＡ ｓ＃ ３ ′ 临时文件用来存储一条线 Ｄ ｏＷ ｈ ｉ ｌ ｅＮ ｏ ｔＥ Ｏ Ｆ （ １ ） １ Ｉ ｎ ｐ ｕ ｔ＃ １ ，ａ （ ２ ） ，ｂ （ ２ ） ，ｃ （ ２ ） ，ｄ （ ２ ） ，ｅ （ ２ ） Ｉ ｆｅ （ ２ ） ＞ｘ ｚ ｚＴ ｈ ｅ ｎ ′ 加入数据限制 ｅ （ ２ ）＝ｘ ｚ ｚ Ｅ ｎ ｄＩ ｆ Ｉ ｆｅ （ ２ ）＜ｘ ｚ ｆＴ ｈ ｅ ｎ （ ２ ）＝ｘ ｚ ｆ ｅ Ｅ ｎ ｄＩ ｆ Ｉ ｆａ （ ２ ）＜＞ ａ （ １ ）Ｔ ｈ ｅ ｎ ′ 线号改变 Ｃ ｌ ｏ ｓ ｅ＃ ３ Ｉ ｆＦ ｉ ｌ ｅ Ｌ ｅ ｎ （ Ｄ ｉ ｒ １ Ｐ ａ ｔ ｈ＋″ ＼ ″ ＋″ ｌ １ ｔ ｘ ｔ ″ ）＝ ０Ｔ ｈ ｅ ｎＧ ｏ Ｔ ｏ２ ′ 判断 ｌ １是不是空文件 Ｃ ａ ｌ ｌｃ ｈ ａ ｚ ｈ ｉ ｔ ｏ ｕ ｙ ｉ ｎ ｇ ′ 测线进行插值投影 ２ Ｏ ｐ ｅ ｎＤ ｉ ｒ １ Ｐ ａ ｔ ｈ＋″ ＼ ″ ＋″ ｌ １ ｔ ｘ ｔ ″ Ｆ ｏ ｒＯ ｕ ｔ ｐ ｕ ｔＡ ｓ＃ ３ ′ 重新建立临时文件存储线 Ｗ ｒ ｉ ｔ ｅ＃ ３ ，ａ （ ２ ） ，ｂ （ ２ ） ，ｃ （ ２ ） ，ｄ （ ２ ） ，ｅ （ ２ ） ａ （ １ ）＝ａ （ ２ ） ｂ （ １ ）＝ｂ （ ２ ） ｃ （ １ ）＝ｃ （ ２ ） ｄ （ １ ）＝ｄ （ ２ ） ｅ （ １ ）＝ｅ （ ２ ） Ｇ ｏ Ｔ ｏ１
其中程序调用过程“ ｃ ｈ ａ ｚ ｈ ｉ ｔ ｏ ｕ ｙ ｉ ｎ ｇ ” 代码如下： ｂ ｌ＝Ｔ ｅ ｘ ｔ １ Ｔ ｅ ｘ ｔ′ 倍率控制投影幅度 ｊ ｄ＝Ｔ ｅ ｘ ｔ ２ Ｔ ｅ ｘ ｔ／ １ ８ ０ １ ４ １ ５ ９ ２ ６ ５ ３ ５ ８ ９ ７ ９ ′ 测线角度 ３ ｚ ｌ＝Ｔ ｅ ｘ ｔ ３ Ｔ ｅ ｘ ｔ ′ 点号增量， 可以控制弃点忽略的程度 Ｏ ｐ ｅ ｎＤ ｉ ｒ １ Ｐ ａ ｔ ｈ＋″ ＼ ″ ＋″ ｌ １ ｔ ｘ ｔ ″ Ｆ ｏ ｒＩ ｎ ｐ ｕ ｔＡ ｓ＃ ４ ′ 开始插入剖面与测线交点坐标 Ｏ ｐ ｅ ｎＤ ｉ ｒ １ Ｐ ａ ｔ ｈ＋″ ＼ ″ ＋″ ｌ ２ ｔ ｘ ｔ ″ Ｆ ｏ ｒＯ ｕ ｔ ｐ ｕ ｔＡ ｓ＃ ５ Ｉ ｎ ｐ ｕ ｔ＃ ４ ，ｆ （ １ ） ，ｇ （ １ ） ，ｈ （ １ ） ，ｊ （ １ ） ，ｋ （ １ ） ｆ （ ２ ）＝ｆ （ １ ） ｇ （ ２ ）＝ｇ （ １ ） ｈ （ ２ ）＝ｈ （ １ ） ｊ （ ２ ）＝ｊ （ １ ） ｋ （ ２ ）＝ｋ （ １ ） Ｗ ｒ ｉ ｔ ｅ＃ ５ ，ｆ （ １ ） ，ｇ （ １ ） ，ｈ （ １ ） ，ｊ （ １ ） ，ｋ （ １ ） Ｄ ｏＷ ｈ ｉ ｌ ｅＮ ｏ ｔＥ Ｏ Ｆ （ ４ ） Ｉ ｎ ｐ ｕ ｔ＃ ４ ，ｆ （ １ ） ，ｇ （ １ ） ，ｈ （ １ ） ，ｊ （ １ ） ，ｋ （ １ ） Ｉ ｆＡ ｂ ｓ （ ｇ （ １ ）－ｇ （ ２ ） ） ＞ｚ ｌＴ ｈ ｅ ｎ ′ 将弃点归为 ０值 Ｗ ｒ ｉ ｔ ｅ＃ ５ ，ｆ （ ２ ） ，ｇ （ ２ ）＋ ０ ５ ，ｈ （ ２ ） ，ｊ （ ２ ） ，０ ，ｆ （ １ ） ，ｇ （ １ ）－ ０ ５ ，ｈ （ １ ） ，ｊ （ １ ） ，０ Ｗ ｒ ｉ ｔ ｅ＃ ５ Ｗ ｒ ｉ ｔ ｅ＃ ５ ，ｆ （ １ ） ，ｇ （ １ ） ，ｈ （ １ ） ，ｊ （ １ ） ，ｋ （ １ ） Ｇ ｏ Ｔ ｏ１ ２ Ｅ ｎ ｄＩ ｆ Ｉ ｆｋ （ １ ） （ ２ ）＜ ０Ｔ ｈ ｅ ｎ ｋ ｇ （ ２ ）＝ｇ （ １ ）＋ ０ ５ ｈ （ ２ ）＝ｈ （ ２ ）＋（ ｈ （ １ ）－ｈ （ ２ ） ） Ａ ｂ ｓ （ ｋ （ ２ ） ） ） ／ （ Ａ ｂ ｓ （ ｋ （ １ ） ）＋Ａ ｂ ｓ （ ｋ （ ２ ） ） ） （ ｊ （ ２ ）＝ｊ （ ２ ）＋（ ｊ （ １ ）－ｊ （ ２ ） ） Ａ ｂ ｓ （ ｋ （ ２ ） ） ） ／ （ Ａ ｂ ｓ （ ｋ （ １ ） ）＋Ａ ｂ ｓ （ ｋ （ ２ ） ） ） （ Ｗ ｒ ｉ ｔ ｅ＃ ５ ，ｆ （ ２ ） ，ｇ （ ２ ） ，ｈ （ ２ ） ，ｊ （ ２ ） ，０ ′ 插入交点， 幅值设为 ０ ， 投影时直接投影到测线上 Ｅ ｎ ｄＩ ｆ Ｗ ｒ ｉ ｔ ｅ＃ ５ ，ｆ （ １ ） ，ｇ （ １ ） ，ｈ （ １ ） ，ｊ （ １ ） ，ｋ （ １ ） １ ２ ｆ （ ２ ）＝ｆ （ １ ） ｇ （ ２ ）＝ｇ （ １ ） ｈ （ ２ ）＝ｈ （ １ ） ２ ５ ６
２ 绘制平面剖面图的程序实现
先将 一 个 工 作 区 所 有 数 据 制 作 成 一 个 文 件， 文 件格式是以 ｄ ａ ｔ 或ｃ ｓ ｖ 为扩展名的 Ａ Ｓ Ｃ Ｉ Ｉ 明码文 件，
图１ 剖面拆解闭合示意图
Ｘ 坐 标、 Ｙ坐 标、 Ｔ值 等 ５列。 数 包括线号、 点 号、 Δ 据先按照线号、 点号进行排序， 然后以同一测线的测 点为 循 环 单 元 进 行 插 值 投 影 计 算， 形成 ｂ ｌ ｎ文 件。 通过 程 序 可 以 很 容 易 完 成 ｂ ｌ ｎ文 件 的 创 建， 以下是 作者用 Ｖ Ｂ编写的程序原代码：
２ ５ ５
矿 产 勘 查 ２ ０ １ ２年 Ｅ ｎ ｄＩ ｆ ａ （ １ ）＝ａ （ ２ ） ｂ （ １ ）＝ｂ （ ２ ） ｃ （ １ ）＝ｃ （ ２ ） ｄ （ １ ）＝ｄ （ ２ ） ｅ （ １ ）＝ｅ （ ２ ） Ｗ ｒ ｉ ｔ ｅ＃ ３ ，ａ （ １ ） ，ｂ （ １ ） ，ｃ （ １ ） ，ｄ （ １ ） ，ｅ （ １ ） ′ 向临时文件写入剖面数据 Ｌ ｏ ｏ ｐ Ｅ ｎ ｄＩ ｆ Ｃ ｌ ｏ ｓ ｅ＃ ３ Ｃ ａ ｌ ｌｃ ｈ ａ ｚ ｈ ｉ ｔ ｏ ｕ ｙ ｉ ｎ ｇ ′ 最后一条线因为跳出了循环需要单独进行插值投影
［ 收稿日期］ ２ ０ １ １－ ０ ７－ １ ５ ［ 第一作者简介］ 王身龙， 男， １ ９ ７ ７年生， 学士， 工程师， 现主要从事矿产勘查和物探数据处理工作。
２ ５ ４
第 ３卷 第 ２期 王身龙等： 利用 Ｓ ｕ ｒ ｆ ｅ ｒ 绘制双色平面剖面图
矿 产 勘 ２ ＭＩ Ｎ Ｅ Ｒ Ａ ＬＥ Ｘ Ｐ Ｌ Ｏ Ｒ Ａ Ｔ Ｉ Ｏ Ｎ ２ ０ １ ２年 ３月 ， ２ ０ １ ２ Ｍａ ｒ ｃ ｈ
利用 Ｓ ｕ ｒ ｆ ｅ ｒ 绘制双色平面剖面图
王身龙， 夏训银， 王洪生， 张进国
（ 华北地质勘查局五一九大队， 保定 ０ ７ １ ０ ５ １ ）
摘 要 文章介绍了使用 Ｓ ｕ ｒ ｆ ｅ ｒ 制作双色平面剖面图的基本思路， 阐述了绘图数据格式 和 创 建 正、 负异常的 ｂ ｌ ｎ文 件方法， 并用 Ｖ Ｂ编写了相应的程序， 经实际应用证明效果比较明显。 关键词 平面剖面图 Ｓ ｕ ｒ ｆ ｅ ｒ 程序 中图分类号： Ｐ ６ ２ ３ ． ６ 文献标识码： Ａ 文章编号：１ ６ ７ ４－ ７ ８ ０ １ （ ２ ０ １ ２ ） ０ ２－ ０ ０ ０ １－ ０ ６
曲线与测线的 交 点 则 以 利 用 正、 负相的值利用 线性求解的方式确定： Ｘ Ｘ ）× ｊ ＝Ｘ １ ＋（ ２ －Ｘ １ Ｙ Ｙ ）× ｊ ＝Ｙ １ ＋（ ２ －Ｙ １ Ｚ １ ＋ Ｚ Ｚ ２ １ Ｚ １ ＋ Ｚ Ｚ ２ １
在建立 ｂ ｌ ｎ文件时要准确地把正、 负 异 常 分 开， 建立 正 异 常 文 件 时 把 负 异 常 点 当 作 ０ ， 同样在建立 负异常文件时把正异常点当作 ０ ， 并 按 测 线、 测点顺 序求出节点坐标， 写入相应的 ｂ ｌ ｎ文件。
［ １－ ４ ］
１ Ｓ ｕ ｒ ｆ ｅ ｒ 绘制平面剖面图的基本思路
用Ｓ ｕ ｒ ｆ ｅ ｒ 绘制平 剖 图 的 基 本 思 路 是 利 用 Ｓ ｕ ｒ ｆ ｅ ｒ 绘 来 基 面 图 的 功 能。 Ｓ ｕ ｒ ｆ ｅ ｒ支 持 白 化 文 件 “ ｂ ｌ ｎ ” ， 它是一个 Ａ Ｓ Ｃ Ｉ Ｉ 明 码 文 件， ｂ ｌ ｎ文 件 格 式 ． 为： Ｎ ，０（ Ｎ为 折 线 线 上 节 点 的 个 数， ０为 白 化 特 征值） Ｘ １ ， Ｙ １（ 第一个节点的坐标） Ｘ ２ ， Ｙ ２（ 第二个节点的坐标） Ｘ ３ ， Ｙ ３ ………… Ｘ ｎ ， Ｙ ｎ（ 第 ｎ节点的坐标） 当首尾两个节点 坐 标 相 同 时， 即 Ｘ １＝Ｘ ｎ且 Ｙ １ ＝Ｙ ｎ ， 那么 表 示 这 是 一 个 封 闭 曲 线， 可以进行颜色 充填。如果 Ｘ １ ｎ或 Ｙ １ ｎ时， 是一条折线。 ≠Ｘ ≠Ｙ 我们在绘制平剖图时先将曲线拆解为正负两个 ） ， 这样就需要建立两个分别表 示正异 常、 区域（ 图１ 负异常的 ｂ ｌ ｎ文件， 在建立 ｂ ｌ ｎ文件的关键是要求出 曲线中每个节点 的 坐 标。 因 此， 我们首先根据绘图 比例尺计算出表示各观测点场值大小的节点坐标和 曲线与测线的交点坐标。假设测线与北方向的夹角 为φ ， 那 么 结 合 测 点 的 坐 标， 任意观测点场值为 Δ Ｔ 的节点坐标为： 。刘浩军（ ２ ０ ０ ４ ） 研究了航空物探平剖图
０ 引言
在规则网磁法勘 探 中， Ｔ平 面 剖 面 图 （ 以下简 Δ 称平剖图） 是 一 种 重 要 的 成 果 图 件， 它在确定异常 形态、 追索异常走向、 确定数据网格化效果的研究工 作中起着不可 或 缺 的 作 用。 通 常 在 绘 制 平 剖 图 时， 把Δ Ｔ的正值部 分 赋 予 暖 色 调 （ 红色系） ， 而负值部 分则赋予冷色调（ 蓝色系 ） ， 这样在图面中能起到突 出异常的作用。 现阶段还没有一种用来制作平剖图的成熟软 件， 通常是用 ｇ ｒ ａ ｐ ｈ ｅ ｒ 来绘制， 但ｇ ｒ ａ ｐ ｈ ｅ ｒ 在 绘 制单 方 向测线还可以实现， 且只能是一条一条的绘制， 效率 太低， 要绘制多方 向 测 线 等 较 复 杂 的 平 剖 图 时 就 非 常困难。 我国物探工作者在平剖图的制作方面作了较 多 的研究