利用Grapher绘制彩色渐变剖面平面图

《利用Grapher绘制彩色渐变剖面平面图》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、万方数据第32卷第2期2008年4月物探与化探GEOPHYSICAL＆GEOCHEMICALEXPLORATIONV01．32．No．2Apr．，2008利用Grapher绘制彩色渐变剖面平面图陈明，何门贵(广东省地质物探工程勘察院，广东广州510800)摘要：在物探扫面工作中，有大量的剖面数据需要绘制成剖面平面图，要求寻找一个高效快速的制图方法。通过对Grapher自动化功能的分析和实践，编制程序实现了快速绘制彩色渐变剖面平面图的需求，说明利用自动化脚本能够在获得精美图件的同时，可以极大地提高作图效率。关键词_．

2、grapher；自动绘图；剖面平面图中图分类号：P631．2文献标识码：A文章编号：10130—8918(2008)02-0196—03用计算机自动绘制物探剖面平面图有许多方法，可以自己编写专门的绘图程序，如刘浩军⋯提出将剖面划分条带区间进行填充，从而达到渐变着色效果；可以利用GIS二次开发进行绘制，如孙中任等旧1提出构造边界文件，通过Surfer自动化接口实现剖面平面图绘制和着色。但前者需要相当成熟的编程能力，不容易推广应用，后者则不能实现渐变着色。为了寻找一个既能快速获得精美图件又简便易行的方法，笔者研究了Gr

3、apher的实现方式。作为地质工作者手中的绘图利器，Grapher拥有大多数流行的数学统计功能和数十种基于二维或三维数据的图形类型，绘制剖面图是其中一种基本功能。随着新版渐变颜色填充特性的加入，Grapher绘制的剖面图和剖面平面图也能像Surfer绘制的等值线图一样拥有多样的色彩、直观地表现物探异常。Grapher提供了全面详尽的二次开发接口，能够很好地与各类编程语言对接，如VisualC++、Del．phi、VisualBasic、ASP．NET等，其本身甚至附带了基于SaxBasic引擎的二次开发环境。笔者以

4、Basic语言为例说明如何利用新版的Grapher6实现物探彩色剖面平面图的自动绘制。1实现方法图1是用Basic语言编制的自动化绘图程序结构框图，这个实现方法最核心的地方是利用Grapher的渐变填充分别对正异常和负异常区域进行着色。正是新的Grapher引入了简化的ColorMap特性，使渐变色填充得以实现。图2是一个白色一黑色渐变填充的例子。收稿日期：2007—05—25图1程序实现框图图2白色一黑色渐变填充在自动化实现过程中，要求我们利用[．CLR]文万方数据2期陈明等：利用Grapher绘制彩色渐变剖面平

5、面图·197·件格式，其定义简单明了。例如一个蓝一绿一黄颜色渐变的[Blue—Green-Yello．CLR]文件内容如下：ColorMap1025550025501002550其中第一行是标识符、版本号和插值方法，第二至第四行是各铆点位置及相应的RGB颜色分量，注意必须有第O％位置和第100％位置。然而由于刚出现的ColorMap特性是一个简化了的实现，它允许我们定义每个百分点上的渐变色，却不允许我们指定每个百分点对应的数值。这种简化的设计在绘制单个剖面图时是很便利的，但由于剖面平面图上要表现多条剖面，每个剖面上

6、量值的图3多剖面的渐变填充不符合实际要求图4彩色渐变剖面平面万方数据·198·物探与化探32卷范围是不同的，使用同一个ColorMap导致每个剖面上填色标准不统一。图3说明使用同一个ColorMap对不同剖面进行填充，导致渐变色不能真实反映异常特征。笔者经过实践，发现通过编程的方法可以将ColorMap中的颜色分量(RGB)与量值对应：假设整个工区所有剖面的极大值为z～，某个剖面量值的极大值为y一，通过y～×255／z．。这样的算式，可以得到该剖面极大值位置的RGB分量，从而得到与整个工区一致的渐变填充。最后通过B

7、asic语言编制成自动化脚本程序(见图1)，根据每个剖面量值特征，获得相应的Col—orMap，从而得到一幅合格的彩色剖面平面图。2应用实例图4是利用自动化脚本程序绘制的某磁法工区剖面平面图。图中正异常用白色一红色渐变表示，生成ColorMap的Basic语言代码如下：clrW2R=Path&”w2r．clr”OpenclrW2RForOutputAs#fncRGB=255一Int(yMax宰255／zMax)Print#fn．”ColorMapll”Print#fn．”0255255”Print#fn．”1002

8、55”&cRGB&””&cRGBClose#fn负异常用白色一蓝色渐变表示，生成ColorMap的Basic语言代码如下：clrW2R=Path&”w2r．clr”Opencl舟2RForOutputAs#fncRGB=255一Int(yMax芈255／zMax)Print#fn．”ColorMap1l”Print#fn．”0255255”Print#fn．