实测剖面图的数字化试验

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

实测剖面图的数字化试验
杨学俊郑先纬
(四川地勘局科研所,成都,610081)
摘要:用电子表格软件,实现音4面数据计算电算化;用MAPGIS的自动转换功能,将电算化的系统数据转换为MAPGIS的图形数据,再用MAPGIS的编辑手段,完成岩
性花纹、素描图、野外照片等剖面要素的输入。

关键词:区嘎实测剖面数字嗨MAPGl8
一、目的、意义
实测剖面是区调工作的必修项目,目前基本上还是用手工完成。

由于地质认识不可能一蹴而就,在复杂区的地质剖面从开始测制到野外资料验收,一般会几易其稿。

这给剖面的制图工作带来许多麻烦,有时为了图面整洁,一经改动就要将剖面图全部重新绘制一遍,作了许多重复工作,效率低下。

当报告中需要引用这些剖面资料时,又不得不再次复制,许多时间就这样浪费在简单的重复劳动中。

如果剖面图也能全部数字化存储在计算机中,那么对剖面资料的查询、修改、调用就会方便很多,许多手工作图时的重复劳动就可避免,工作效率会成倍提高。

由于剖面观测数据客观、准确、系统,已经具备了数字化的良好基础。

用电子表格软件,首先可以使剖面的数据计算电算化,从而减少了人工计算的负担,也减少了人工计算出错的几率。

其次,有了电算化的系统数据,也容易被MAPGIS利用,可直接由计算数据文件转换为图形数据文件。

二、试验内容
(一)剖面数据计算的程序化
剖面数据计算的程序化不需要复杂的编程就可实现。

有现成的电子表格应用软件(如:WPS、Excel、金算盘电子表格)可以利用。

在电子表格中输入原始的剖面测量数据(方位、斜距、坡度、产状)和相应的计算公式,就可得出其它数据(平距、高差、倾剖夹角、走剖夹角、平面导线图上测点坐标、地形图上测点坐标、视倾角、分层厚度等)。

数据的输入工作只包括一次原始数据的输入,和一次公式输入,中间过程自动计算完成,不需要人工干预。

比手工计算工作量减少5.10倍,出错的概率也大大减少。

在剖面基本数据计算表的基础上,我们就可以编制用代数坐标法取代尺、规作图的计算表格,MAPGIs制图数据表格,以及地层厚度计算表格等。

[作者简介]杨学俊(1968一),男,工程师,区域地质与旷产普查专业,从事互调工作。

=;
・衍
表一1剖面基本数据计算表
,ABCDEF
誓。

导线方位坡度角斜距平距(m)=高差(m产
罄导线号
(角度)(角度)(m)D+COS(RADlANS(C))D+SIN(RADIANS(C))20一l100.51009962.872
31—27i.810099.03.13.92
42—362.2323l98.1.12
(--)用代数坐标法取代尺、规作图
1。

作图原理
常规的作图方法是按基本数据计算表中的数据,用量角器量出方位、用三角板确定距离把线条和点绘制在米厘纸上。

它的致命弱点是,由于每一次的量画都是以前一次为基础进行,所以会将每一次的量度误差累计起来,后画的部分精度越来越差,如果哪一导线画错,从它往后的导线就完全报废。

而且,这一方法没有充分利用米厘纸可以直接读取坐标值的特性,作图时用三角板的刻度来读数,读图时又常用米厘纸本身的刻度来读数.标准不一致,致使在读图时产生系统误差。

代数坐标法制图.就是利用米厘纸可以直接读取坐标值的特性,按三角函数关系计算出各导线终点以剖面起点为原点的直角坐标值,然后用该坐标值按米厘纸自身的刻度标绘在图纸上。

三角关系如下图所示:
YJL
———寺剖面制图
总方位
2号导线终点(x2,y2
,一觚u,,蒜////
音蛹旋E点彳
(0。

0)
———‘—巳二~
l号±|i溅(xbz森i)r—~之里2号地形

’(x2,z2)说明:L=导线平距,Q=导线与剖面制图总方位的夹角,B=导线坡度角
(x,Y):绘导线平面圈坐标,(X,z):绘地形剖面图坐标
图l代数坐标法绘制导线平面圈和剖面地形线示意图
≮h

2.代数坐标法绘制导线平面图
设导线起点坐标为(0,0),选取一定的剖面总方位,设定作图比例尺、作图单位,则各导线的起点与终点坐标可由F表2算出。

按表中的坐标值,在米厘纸上直接点绘出各导线的起点与终点,以及它们之间的连线。

这样就取代了直尺和量角器的作图方法.提高了精度和效率。

表2导线平面图测点坐标计算表
ABCDEFGHIJ
平面导线图t测点坐标(按115000作图,单位:mm)
导剖面导坡
斜距平距XY
l线总方
线度
方(m)
(m)
X=F*COS(RADIANS(B.c))
Y=PSIN(RADIANS(B..C)1号位角∑X艺Y
位‘1000/LI:例尺倒数+1000/比倒尺倒致
起点900000000030~l90100-51009962196196-35.3.54I~29071.81009903187383642.9
代数坐标法绘制剖面地形线
设地形线起点坐标为(0,0).注意正负地形的拐点位置,则地形线可由F表3计算出的各坐标点连接而成。

如果剖面上地形线的起点需要在垂向上平移一段距离11,计算时则可设±也形线起点坐标为(0,11)就楚r。

表3削曲地形线测点坐标计算表
ABCDEFGHlJK忖削…"削l血I地膨战洲点坐标(拽1/5000佧蹦.单位・mm)
l地电}lIt垤利
’}’蚺I向肚XZ
Ji旭蚺I
X_Z_
(¨1)(111)U似『『I“¨1
”‘’‘)S(I(AI)IANtq(I|-(’11XXIc.NIN(H^I)I^N川11))ZZ
怖‘_…l●
’1000/lL驯Jt恻数+10001IE例八f“数2起点0000000000
1r—,●1——
JIJIU¨IOO-5IUO99:1'2.H7219b19.6-I.7・I.7-_“^_--
4l~29071-81009903.13918.738.3—2.8-4.54.确定剖面测点在地形图上的位置
n:地形I鳘|I:青出剖而起点的商斯坐标值(x。

Y’)。

必须采用900制l鳘I总方位,按表4可iI绅…^l,战^.地形…¨Iq岛埘小h、似首。

JIlllLh‘法II,以比较力。

使地腮绘、lI盘址安测削lhIb地形I玺I的吻合程度。

‰引
表4剖面测点在地形图上的位置计算表
ABDEIGI14卜1lJ
地形圉上测点高斯坐标(单位。

km)
导坡平距
高差纬向坐标x.经向坐标Y。

I线导线

斜距
(m相对上一点的相对上一点的方位(m)(m)
廿角)增量=
£X’
增量=£r
E’SIN(RADIANSE+COS(RADIANS
(90・C))/1000(90・C))/1000
2起点00000O3436.84017767.5630-l100.51009962-8.72.0.023436.8201017767.664l、27I.810099.03.13.90033436850.0917767.75(三)用MAPGIS自动转换功能绘制剖面基本图形和符号在电子表格计算出导线平面测点坐标、剖面地形线测点坐标、地形图上测点高斯坐标的基础上,MAPGIS就可以利用这些系统化的数据.将数值信息转换为可视化的圈形信息。

1.用MAPGIS制作常规剖面导线平面圈
MAPGIS绘制导线平面图,可将其分为两个基本要素:导线起始和终结的点位与编号;起始与终结点位间的连线,前者由点图元(一般是小圆点和数字)实现,后者用线图元(折线)实现。

采H:|明码格式文件,将导线的点位和编号组成点明码格式文件(・.wat)。

将导线的连线组成线明码格式文件(+.wal),再用MAPGIS的数据转换系统将其转换为点文件(+.wt)和线文件f+w0。

编写导线起始和终结点位明码格式文件:
1.按前表2《导线平面图测点坐标计算表》的计算方法,以剖面总方位为制图方向.按作幽比例尺,以毫米(mm)为单位,计算出导线平面图上测点坐标(X11,Yn)。

2.用字处理软件(WPS、WORD等)。

按MAPGIS明码格式编排出符合要求的}.wat文本文件。

格式如F:
氐乳
文r卜说明:
第一行:文件头,8个字节,统一为:wMAP9022(5.32版之前的文件为:WMAP8022或WMAP7022或wMAP6022;5.32版后可能还有所不同。


第二行:点陶元的数目,按真实个数填写。

从第三行起,每个点图元的数据占一行,各项数据之间用逗号分隔。

依次为:
点图元的计算机坐标(以mm为单位),Xn,Yn,
点图元的序号,n,
点图元的类型.(从0,1,2,3,4.5中取值。

0为字符串:1为子图:2为圆:3为弧;4为幽像:5为文本),
点图元的图形参数和图层号.(点图元的类型不同,参数设置也不同,应视具体情况而定。

如果对此参数设置不明确,也可以反其道而行之,先在MAPGlS中编辑输入一个此类点。

调整好它的参数和图层,保存后在MAPGIS数据接口转换子系统中把它转成+.war文件,用文本编辑手段把它的相关数据粘贴过来即可。


3.字处理软件(WPS、WORD等)保存的文本文件,~般以+.txt为后缀。

因此,用计算机操作系统的文件重命名工具,将文件的后缀・.txt改为MAPGIS要求的+.wat后缀。

4.用MAPGIS数据接口转换子系统将+.war文件转换成MAPGIS点文件。

编写导线编号的1VIAPGIS明码格式文件:
把计算山的测点坐标,加上一个平移量(如:Xn+0.5,Yn+0.5.即把编号放在点位的右上角),作为MAPG[S明码格式文件中的“点图元的计算机坐标”数据项,参照“导线起始和终结点位明码格式文件编辑”所用的步骤,“点图元的类型”数据项取…0’值(字符串),按字符串的要求设置图形参数和图层,完成导线编号的MAPGIS点图形文件。

编写导线的连线明码格式文件:
第1步:按前表2《导线平面圈测点坐标计算表》的计算方法,以剖面总方位为制图方向.按作圈比例尺,以毫米(mm)为单位,计算出导线平面图上测点坐标(xn,Yn)。

第2步:把计算出的测点坐标,按顺序作为MAPGIS明码格式文件中的“线图元的点的坐标”数据项,_};fj字处理软件(WPS、WORD等),按MAPGIS明码格式编排出符合要求的}.waI文本文件。

格式如F:
文件说明:
第一行:文件头,8个字节。

统一为:WMAP9021(5.32版之前的文件为:mIAP8021或wMAP702i或WMAP6021,5.32版后可能还有所不同。


第二行:线图元的条数+1(参考手册中写明是“线数”,以我们的应用来看,5.32
版实际t}=Ij的是“线数+l”即:如果有3条线,就要弓成4.否则转换厉会少
一条线或没有线。


第三行:第】条线圈元的图形参数、幽偿号及是否透明输出(如果对此参数设置不明确,也可以反其道而行之,先在MAPGIS中编辑输入一条线,调鹩好它的参数羽I
图层及是否透明输出,保存后在MAPGIS数据接口转换子系统中把它转成
+.wal文件,用文本编辑手段把它的相关数据粘贴过米即可。

)第四行:第1条线的点数
第五行:第1条线第1个点的坐标
第六行:第1条线第2个点的坐标
第n行:第1条线最后一个点的坐标
第n+l行:第】条线的fD号+逗号(参考手册中写明灶“ID号,线I(:J变”。

以我们的应用米看,绘制剖面时,在这里不需要计算线艮度。

内此我们试着把“线艮度”
一项省略或取任意值,结果.532版MAPGIS数据接口转换子系统接受这种
修改,但ID号后面的逗号必须保留,否则要出现非法操作。

)第n+2行:如果有第2条线,重复上面第三行及其后备行的编辑。

最后,用MAPGIS数据接[J转换子系统将以上方法编辑五,的+.wal文{;.{=转换成MAPGIS剖面导线平面图+.wI线文件。

2.MAPGIS剖面导线平面图与数字地形图的套合
作常规剖面平面幽时,考虑到作幽与心fⅨ:fl"J方便,要选择一个剖面总体方位.以此方位向图纸的最大延伸方向作圈。

因此,如袭2所示,要将剖面总体方位引入剑制图坐标的计算中米。

而需要与数字地形幽套合的剖面导线平面目的MAPGIS制剧坐标的计算,与上述竹!常规剖面平面圈时的计算略有不同。

它是以处于理想套合位置的数字地形图(即方里网摆平的地形图)为参考对缘,选择的剖面总体方位必须是90。

.在表日《剖面测点在地形图上的位置计算表》的基础上,按方里网高斯坐标与MAPGIS直角坐标之间的转换公式,就能计算出与数字地形圈套合的剖面导线平面图的MAPGIS制图坐标(表5)。

设:数字地形图已处在理想套合位置上(即方里网是拦平的)
x。

;为标准削框左。

r角数值最小的方里网交义点的南北向离斯直角座标(单位:km)Y。

为其对应的MAPGIS纵-J:E)/,J,(单位:mill);
Y:;为标准幽框左r角数值最小的方里1169交义,:i的k,i,ti向高斯直角库标(单位:km)x.。

.为其对应的MAPGIS横坐标(单位:ltlm):
X为任意点的MAPGIS横坐标(单位:mm):
Y为任意点的MAPGIS纵坐标(单位:mm);
Y’为该任意点f1"3尔pq向高斯直角废标(单位:kin):
x’为该任意点的南北向赢斯直角庳标(单位:km)
则:x=(Y’一k){i000000/比例尺倒数+x。


Y=(x’一x驺)+1000000/比例尺倒数+Y瑚
≮W
表5与数字地形图套合的剖面导线平面MAPGl8制图坐标计算表ABCDE
地形圈上测点商斯MAPGIS制图坐标(比例尺:1150000,单位:mill)导
坐标(单位:km)方里网交叉点(343500,1776300)对应MAPGIS坐标原点(0,0)
1线
纬向坐标经向坐标X=Y=

X。

(C一17763){lOOOOOO/LL例尺倒数(B-3435){1000000/比例尺倒数起点34368417767569l2036.80
30一l343682l776766931636.45
4l一234368517767.75950337.10
52—33436.87177677895.603740
参照前述《用MAPGIS制作常规剖面导线平面图》的一系列步骤,就可得到与MAPGIS数字地形图坐标系和比例尺完全一致的数字剖面导线平面图,用它来观察实测剖面与地形圈的吻合程度,既方便又直观。

3.MAPGIS剖面图地形线的制作
地形线也有准确的测量记录,注意正负地形的拐点位置,按茎卫《剖面地形线测点坐标计算表》可以得到地形的起点、终点和拐点的制图坐标。

MAPGIS绘制剖面图地形线蚓,将其作为线型来表达:由起始点、正负地形的拐点与终结点位间的连线构成,刚线倒元实现(MAPGIS5.32版增加了光滑曲线)。

采j{{明码格式文件,将地形线的连线组成线明码格式文件(+.wal),再用MAPGIS的数据转换系统将其转换为线文件(+.w1)。

具体步骤参j!{i《用MAPGtS制作常规剖面导线平面图》第4步“编写导线起始与终结点位间连线的MAPGIS明码格式文件”的作法。

4.有严格位置标识意义地质内容的制作
地质内容按其位置标识意义严格程度的不同,可分为:有严格位置标识意义的内容,和I位置标识意义不很严格的内容这两类。

前者如地层界线,采样位置、小岩脉的产出位置等,后者如一段岩层的总体产状,一段地层的代号,岩脉的代号等。

有严格位置格式意义的内容,一定也有准确的测量位置记录,可以由此按制图比例用电子表格软件计算出其在导线上的交点坐标,然后用MAPGIS明码格式文件将这些交点转换成醒目的指示性文件,再进入MAPGIS编辑子系统。

先按与导线平面图交点处实际应标注的内容,输入地质现象(或按走向输入地质界线,或按位避输入采样位置,等等),然后按晴盥投影,在剖面圈的相应位置上输入有关的地质现象。

(四)用MAPGIS的编辑功能输入剖面图形和符号
MAPGIS的自动转换功能以很高的效率替代我们的手工操作,绘制山了剖面基本图形(导线平面圈、剖蕊地形线图、地层界限等)和一些有严格位置标识意义的符号(导线编号、采样位置等)。

地质剖面是一项比较复杂的专业性较强的工作,目前。

还不能把剖面中的所
=--rf
有内容通过MAPGIS的自动转换功能米实现。

那些不能通过自动转换功能来实现的剖面地质内容,一般都是位置标识意义不很严格的内容。

可以通过MAPGIS的图形编辑系统,把它们分解成点、线、面等图形要素,按地质意义的合理性与制图的美观性等原则,分别输入和编辑。

应Ⅲ起米有tt'J'T.如臼动转换那样快捷,也有一些编辑功能不那么尽如人意,这里存在一个软件的使用者与软件相适应的问题,也有软件臼身不断改进、完善以满足用户需求的问题。

对我们地质工作者来说,思想上要有不断提高各项工作数字化程度的意识,注意学习和适应软件的现有功能和操作方式,尽快改变手工时代养成的操作习惯,总结经验、提高应用技能,创造性地发挥软件的功能,使我们的数字化水平不断提高,以适应信息时代飞速发展的要求。

1.点型剖面地质内容的编辑输入
那些可以用计算机常规字库和符号库输入的文字、符号,如:图名、比例尺、图例、简单的地层代号、矿产代号、变质代号、岩脉代号、走向倾角、样品编号等各种注释内容,用MAPGIS点圈元中的注释来表达。

那些在计算机常规字库和符号库中不存在的文字、符号,如:复杂的地层代号、产状符号、构造符号、化石符号、岩性花纹等专业性的特殊符号,用MAPGlS点图元中的子圈来表达。

另外,在MAPGIS点图元输入功能中还包括:圆、口&、版面、幽像的输入。

可以由此输入其它软件编辑好的文本文件、素描圈或野外照片等。

2.线型剖面地质内容的编辑输入
线型剖面地质内容,除已经由自动转换功能实现的剖面基本图形(导线平面幽、剖面地形线圈、地层界限等)以外,还包括剖面地质界限(地层、断层、特殊岩性等)和豳面整饰内容(水平标尺、垂直标尺、图例、柱状图等)。

都可以用线图元米输入、修改。

另外,一些看似点状,但具备一向延伸规律的符号,如岩性花纹中的混合岩化:~~~石英岩:・’・-。

・...二云片岩:一一一等都可以造成线型来使_El{|,这比用点图元来实现要快捷得多。

3.面型剖面地质内容的编辑输入
在MAPGIS系统中,用丰富的色彩和花纹来表达地质内容比手jf:作幽要容易。

加之彩色打印技术的酱及,因此,我们应有意识地加大色彩的使用力度,更多地用彩色来表现我们的地质认识。

所以,MAPGIS制图的面型剖面地质内容,可以包括刚以区分变质带、地层、岩性等的面色,也可以是用以区分变质带、地层、岩性等的具二向延伸规律的花纹符号。

这样我们就可发挥计算机的优势,把那些过去手工制酗需要一笔一笔描绘的岩性界线和花纹.用计算机的面型输入手段,大面积快速完成。

(需要指出的是,MAPGIS5.32版图案库的设置对我们来说还不够完善。

我们希望MAPGIS的开发者能增加一个参数,来设置图案符号的倾斜度。

这样我们就可以让幽案花纹根据地层的产状而变化。


k吐
总之,用MAPGIS作实测地质剖面的尝试才开始,需要区调工作者不断实践,总结技巧、发现问题,以便MAPGIS的开发者有针对性地加以改进,共同提高区调地质工作的数字化水平。

主要参考文献
e1】中华人民共和国国家标准GBl2343—901:25000、1:50000地形圈编绘规范(1990-12-01实施)
[2]四川地质矿产开发局1:5万区调地质幽计算机成幽有关规定(哲行),1996年8月
[3]区域地质矿产调查工作图式图例(1:50000),地质出版社,1983
“]四J¨地质矿产局区域地质调查工作手册(1:50000),1992,
[5]国家测绘总局.地图制图.1980
[6]中国地质大学(武汉)信息工程学院.MAPGIS地理信息系统参考手册.1998
[7]中华人民共和国地质矿产行业标准数字化地质图图层及属性表格式(DZ/T0197一1997)
[8]陈克强等。

彩色地图辅助编图系统(MAPCAD)制图精度问题,中国区域地质.1997
【∞辜寄蓉等,1:5万地质图sQc,计算机制躅研究报告,四川地勘局区域地质调查队,1999
[10】地球科学一中国地质大学学报,地理信息系统专集。

1998年第23卷第4期
[11]徐绍铨、吴祖仰,大地测盘学,武汉测绘科技大学出版社,1996
[12]陈永奇,高等应用测量,武汉测绘科技大学出版社,1996
≮f7。

相关文档
最新文档