资源描述：
《在MAPGIS下快速实现柱状剖面图的绘制》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、第26卷第2期2004年5月物探化探计算技术Vol126No.2May2004COMPUTINGTECHNIQUESFORGEOPHYSICALANDGEOCHEMICALEXPLORATION文章编号:1001—1749(2004)02—0173—04在MAPGIS下快速实现柱状剖面图的绘制李沙园,陈昕华,戈永怡(中国石化勘探开发研究院合肥化探所,安徽合肥230022)摘要:将油气化探数据用FORTRAN短程序生成MAPGIS点、线、面明码文件后,再利用MAPGIS图形编辑软件对图形进行编辑加
2、工,达到快速制作油气化探柱状剖面图的目的。使用该方法能充分利用MAPGIS系统资源,降低编制绘图程序的难度。关键词:MAPGIS;制作柱状剖面图;油气化探中图分类号:P632文献标识码:AAFASTDRAWINGMETHODFORCOLUMNARSECTIONUSINGMAPGISLISha2yuan,CHENXin2hua,GEYong2yi(PetroleumGeochemicalProspectingInstitute,SINOPEC,HefeiAnhui230022,China)Abst
3、ract:ThisarticleintroducesamethodtoconvertdataofgeochemicalexplorationforoilandgasintoMAPGISpoint,lineandsurfacecodefilesusingashortFORTRANprogramandthenthecolumnarsectionofgeochemicalexplorationcanbedrawnbyMAPGISeditsoftware.ThemethodwithMAPGISresou
4、rcescanreducethedifficultiesencounteredinconventionalprogramdesignforcolumnarsectiondrawing.Keywords:MAPGIS;columnarmapdrawing;geochemicalexplorationforoilandgas0引言获2002年度国家科技进步二等奖的MAPGIS软件,作为国产工具型GIS软件,其强大的功能,良好的技术支持,已经获得普遍的认可,使用该系统软件编辑的图件具有方便、速度快、精
5、度高的特点,在地学领域具有广大的用户群。油气化探柱状剖面图是将岩性剖面和多个化探指标浓度曲线在同一个图上按其取样深度表示出来,并在图上标注地层岩性说明、化探采样点及样号。笔者利用自编程序,把已有的地层岩性数据及化探数据文件转换成MAPGIS点、线、面明码文件,再利用MAPGIS图形编辑软件对图形进行编辑加工,最后形成一幅理想的柱状剖面图。该其方法简便、易行且速度快,在实际应用中取得了较好的效果。1作图的基本思想为了表示不同的岩性,需要编程形成很多图库,以便绘图时调用。其工作量较大且有一定难度,又
6、不直收稿日期:2003-09-02·174·物探化探计算技术26卷观。而MAPGIS本身带有近千种子图库(符号库、线型库、图案库),要编辑添加子图也很方便、直观。因此利用MAPGIS本身的子图库,再根据需要编辑添加一些子图,就可作为岩性图库,而使绘图程序简单化,稍有编程能力即可实现。2实现化探柱状剖面图的绘制由上所述,在做图前先要编辑添加一些子图,作为岩性图库时,可利用MAPGIS本身的子图库就能方便地制作添加。图1是我们制作的部分图案。2.1化探柱状剖面图的绘制过程用此方法绘制油气化探柱状剖面
7、图的过程见图2所示。图1部分岩性图例图2简单绘图流程Fig.1ThepartsoflithologylegendsFig.2Thegeneraldrawingflowchart从上述流程图可以看出,所绘制的剖面图是将化探数据用程序生成MAPGIS明码文件后再转换为点、线、面文件,以图形叠加的方式制作的。212编程实现将化探数据转换成MAPGIS明码文件的FORTRAN程序主要部分如下。参数说明:nscale:比例尺ix0=0ix1=1y(i):每层厚度x3=0.3zl(i):累积厚度x0=0.0
8、outfile2:.wal文件名iy0=0name1:.wap文件名dd=nscaleö1000.0name2:.wat文件名y0=höddh:图的长度xd=wdwd:图的宽度n1:岩性分层数C打开MAPGIS点、线、面文件并写文件头n2:化探采样数open(5,file=outfile2)a(i):第i列宽度open(6,file=name1)ic(i):第i层岩性代号,和区参数open(7,file=name2)中的图案对应write(7,`(1x,a)`)`WMAP9022`⋯⋯略(程序头
