3374519744利用Surfer Scripter自动生成高质量高密度电法成果图

1、利用surfer scripter自动生成高质量高密度电法成果图侯继荣1，2 、张会强1，2、黄小年1，2（1，中交第一公路勘察设计研究院有限公司 西安市 710075；2，西安中交公路岩土工程有限公司 西安市 710075）摘要：本文针对目前高密度电法工作中成图质量和和制图效率中存在的问题，提出了一种利用surfer scripter自动生成高质量的成果图的方法，给出了方法的具体实现过程和详细代码，并通过在青藏高原某高速公路实际工程中的应用，验证了该方法有效性和实用性。关键词：surfer；高密度电法；自动成图中图分类号：u412.22 文献标识码：a1引言高密度电法是近几十年发展起来的的一

2、种工程物探方法，因其勘探点密度大、精度高、效果明显，在工程和环境勘察中应用极为广泛。但在实际勘探中，其解释精度除与野外勘探点密度相关外，与最终成果图的质量也密切相关。目前国内高密度电法成果图大都利用瑞典反演软件res2dinv直接生成，生成的成果图为位图格式，横向最大像素点个数与成图计算机显示器最大像素点个数一致，当勘探剖面较长或者为提高勘探精度而增加电极数量时，受制于显示器的分辨率，成果图质量无法相应提高，当电极数量大于一定时成果图甚至无法反映全部勘探点的数据特征，使解释精度受到影响。解决这个问题的一种方法是将res2dinv反演结果保存为surfer数据格式的文件，在surfer下手工绘制

3、能反映全面数据特征且不受的计算机显示器影响的高质量成果图，但这种成图方法操作步骤繁琐，成图效率极低，限制了其在实际生产中应用。本文通过在surfer自带的开发工具surfer scriper上进行二次开发，将保存后的数据在surfer下自动、快速成图，能极大地提高制作高质量成图的效率。2 surfer及surfer scripter 简介surfer是美国golden software公司推出的一套二维和三维图形绘制软件，可以利用其绘制等值线图、表面图、阴影地貌图等多种地质图形。除了强大的手动绘图功能外，surfer还提供多种activex自动化对象（如图1所示），每一个自动化对象都有其自己的

4、属性和方法，其中属性描述的是对象的状态，方法是对象的执行功能。通过编程给属性赋予不同的参数值和利用对象的方法执行各种操作，就可以实现surfer 软件提供的所有绘图功能。surfer scripter是surfer自带的二次开发工具，利用其提供编程语言scripter basic language可以在不依赖外部编程工具的情况下操作控制surfer自动化对象，实现图形绘制自动化。图1：surfer对象的层序关系3数据保存为surfer格式的方法、保存后的文件及其用途为了实现surfer scripter自动绘图，必须清楚绘图所使用的数据文件的形式和内容。将数据保存为surfer格式目前在res

5、2dinv软件上已经实现，但在保存时，选项较多，如图2所示，除第1项迭代次数外，其它选项选择不同，最后保存的文件名和内容也不同，为了方便后面的自动化成图，首先需统一这些选项的选择如下：第2项选择保存为电阻率对数值；第3项选择底图覆盖方法；第4项需要保存的截面只勾选含地形的模型电阻率值。图2：将反演数据保存为surfer 格式的选项按照上面的选择点击确定，将会自动根据原数据文件产生4个文件，下面以原数据文件abc.dat为例，说明保存为surfer 格式后产生的文件及其用途：abc _topreslog.dat：反演模型电阻率对数值的surfer格式数据文件，用于制作电阻率等值线图。abc _t

6、opreslog.bln：surfer边界文件，用于绘制基面图，以显示剖面图地形和隐藏绘图边界外的图形元素。abc _topreslog.lvl：surfer等级数据文件，用于用色谱方式表示电阻率的大小。abc _topreslog_post.dat：张贴图数据文件，用于表示勘探点的位置、高程等信息，该文件数据主要有三列：第一列用于用平距方式表示勘探点位置；第二列为勘探点的高程；第三列用于用地面斜距方式表示勘探点位置。4 自动生成成果图的实现利用scripter生成成果图就是利用scripter basic language编制代码操作各种surfer自动化对象，在图形文档的同一位置依次添加等

7、值线图、基面图、张贴图等图形对象，并通过调整图形，使各图形叠加组合成最终成果图。程序流程为：创建surfer程序对象建立图形文档读取电法数据数据网格化利用网格化文件绘制等值线图调用等级数据文件充填等值线，生成用色谱方式表示电阻率的彩色等值线图调用边界文件绘制基面图，绘制张贴图，标注勘探点的位置完成绘图。在绘图过程中，数据网格化是控制成图质量的关键，网格化的行数和列数选择较小时，成果图将无法反映全部数据的特征，影响精度，太大是又会造成计算浪费，因此设置成根据用户输入动态确定网格化的行、列数；网格化插值方法种类较多，克里格法或自然邻近点法比较适合高密度电法数据，但自然邻近点法容易使图形边界产生锯齿

8、状条纹，影响图形美观，所以这里网格方法选择克里格法。surfer下的图形为矢量图格式，可以任意缩放，但如需要将图形输出为其他格式（如位图）的图形或者打印，就需要设置合适的绘图比例，这里将绘图比例也设置成由用户根据需要输入。下面是生成成果图的主要代码：创建surfer对象并建立图形文档set srfapp = createobject(surfer.application)srfapp.visible = true 使surfer应用程序可见set plotdoc = srfapp.documents.add(srfdocplot)提示输入网格化使用的行数和列数cols =inputbox(输入

9、网格化使用的行数)rows= inputbox(输入网格化使用的列数)数据网格化infile=ysname + _topreslog.datgridfile = ysname + .grdsrfapp.griddata datafile:=infile, numcols:= cols,numrows:= rows, _algorithm:= srfkriging, showreport:=false, outgrid:=gridfile绘制等值线set mapframe = plotdoc.shapes.addcontourmap(gridfile)调用等级数据文件充填等值线set conto

10、urmap = mapframe.overlays(1)infile=ysname + _topreslog.lvl绘制基面图，确定绘图边界和显示地形轮廓infile=ysname+_topreslog.blnset mapframe= plotdoc.shapes.addbasemap(infile)set basemap = mapframe.overlays(1)basemap.fill.pattern= diagonal crossbasemap.fill.backcolor=srfcolorwhitebasemap.fill.forecolor=srfcolorwhite用地面斜距方

11、式标注勘探点位置infile=ysname+_topreslog_post.datset mapframe= plotdoc.shapes.addpostmap(infile)set postmap=mapframe.overlays(1)postmap.labcol=3 用斜距方式标注勘探点位置这样即初步完成了成果图的绘制，如果需要将图形进一步调整或美化，只需再添加改变相关对象属性的代码即可。5应用实例为验证自动化成图方法的实际效果，我们选择了青藏高原某高速公路隧道勘察中的高密度电法资料进行测试，该电法剖面起止桩号为k123+010k127+110，长度4100米，首先在res2dinv软件

12、下直接成图，成果图如图3所示，然后用surfer scripter程序自动成图，输入网格化的行数为2000，列数为200，自动成图过程耗时约两分钟左右，其中大部分时间主要耗费在数据网格化过程中，真正成图过程仅需几秒种，生成的成果图见图4。通过比较可以发现，图4 和图3很相似，这也证明了自动生成的成果图是合理的，但图3为位图格式，而图4为矢量图格式，可以任意缩放而不影响显示效果，将两张图同时放大，就可以发现利用surfer scipter自动生成的成果图质量有明显提高，对局部异常细节和异常边界反应更为清晰，如图5所示。如果用传统手工方法在surfer下制作与图4完全相同的一张图，因为制图步骤较多

13、、操作繁琐，一般需要20分钟以上，利用surfer scipter自动成图技术较传统成图方法具有明显的优势。图3：res2dinv直接生成的高密度电法成果图图4：利用surfer scripter自动生成的高密度电法成果图 a. res2dinv生成的成果图 b.利用surfer scripter自动生成的成果图图5：两种成果图的局部放大图6结语目前surfer软件在工程物探成图中上应用较为广泛，利用其附带surfer scipter进行二次开发，将反演后的电法数据保存为surfer格式，然后自动生成成果图，解决传统手工绘制高质量成果图中操作繁琐，效率低下的问题，能够简单快速地提高高密度电法成果图的质量，节省大量时间和人力成本。因此，这种自动成图技术可以实际生产中推广，具有较高的实际应有价值。参 考 文 献（1）闫昊明，surfer绘图指南m，武汉：中国科学院