PDF《S下快速实现 柱状剖面图的绘制》

2 .

1 化探柱状剖面 图的绘制过程用此方法 绘制油气化探柱状剖面图的过 程见图2所 示.序号图案岩性lJ固JI白云质灰岩与l99l11l藻团块 灰岩 互层 I l l l

2 0 o I I l I 灰岩 l

1 1 f

2 0 l 白云质灰岩 I

1 l l I 够l国10l202簟粒灰 岩lI母I@JfIIII203lI、:,II含生物屑灰岩陶辑绘图敌据文件 l 程序睦序转换明码文件 I 【 形成点强 、 面文件 l 编辑点 线、面文件 I 形成 工程文 件I输出图1部分岩性图例圈2简单绘图漉程Fi g.

1 Th e p a r t s o f l i t h o l o g y l e g e n d s Fi g .

2 Th e g e n e r a l d r a wi ng f lo wc h a r t 从上述 流程 图可 以看 出,所绘制 的剖 面 图是将化 探数 据用 程序 生成 MAP GI S明码文 件后 再转 换为点、线、 面文件 , 以图形 叠加 的方式制 作的.

2 .

2 编程实现将化探数 据转换成 MA P G I S明码 文件的 F O R T R AN程序 主要部分 如下 . 参数说明:nscale:比例尺y( i ) : z l ( i ) : o u t f i l e 2: na mel: na me 2: h: wd: 每层厚度累积厚度.wa l 文件名.wa p文件名.wa t 文件名图的长度图的宽度 岩性分层 数化探采样数第i列宽度 第i层岩性代号 , 和 区参 数中的图案对应略(程序头,定义数组、变量、读数据ix0:=0CCixl一1x3― 0.3 x0― 0.0 i y0=

0 d d― n s e a l e /

1 0

0 0 .

0 y

0 = h / d d x d: = wd 打开 MAP GI S点、线、 面文件并写 文件 头open( 5, f i l e ― o u t f i l e

2 ) o p e n( 6, f i l e = n a me

1 ) o p e n( 7. f i l e ― n a me

2 ) wr i t e ( 7, '

(

1 x, a ) '

) '

W M AP9

0 2

2 '

wr i t e ( 6, '

(

1 x, a ) ) W M AP9

0 2

3 wr i t e ( 5, '

(

1 x, a )'

) W M AP9

0 2 1'

读数据xl ― x O + a(

1 ) xl 1: xl一3O. 维普资讯 http://www.cqvip.com 2期 李沙 园等 :在MAP GI S下快速 实现柱状剖 面图的绘制 ・1

7 5・ x2= =xl一一20 yl= =y0 nn= =0 y d一0.nn n= =

0 d o

6 0

0 i 一1.n1 y d= y d+ y( i ) C 判 断是否是 相 同的岩性 , 如果 两层 相 同则 归并为一层,否则进一步判断后画图if((ic(i)一ic(i+1)).ne.0) t h e n C 判断 厚度 是否 ≥

2 r am,2 r am 以上可 以作 岩性图if(yd.gt.2.*d d)t h e n n n n― n n n+

1 i a( n n n) 一inn― n i l +

1 n t ― n t +

1 h d=z l ( i ) / d d y 2= y

0 - - h d C 在面文件 中画线 ( . wa p ) wr i t e(

6 ,'

(

1 x, a) )'

1 ,0 ,1 ,1 ,1

0 .0

0 0

0 0 0,1

0 .

0 0

00 00, 1, 0, 0'

wr i '

t e (

6 , ¨(

1 x, i l, '

, '

i

1 ) ¨) i x 0, i x

0 wr i t e (

6 , '

'

(

1 x, i

6 , '

, '

, i

2 ) ) n n, i x0 wr i t e (

6 ,

4 O O )

5 wnt e (

6 ,

2 0

0 ) xl , y l wr i t e (

6 ,

2 0

0 ) x2, y l wr t t e (

6 ,

2 0

0 ) x 2, y

2 wr i t e (

6 ,

2 0

0 ) xl , y

2 wr t t e(

6 ,

2 0

0 ) xl , yl wr t t e (

6 , ¨(

1 x, i 6, '

, '

, f

3 .

1 ) ¨) n n, x

3 e n d i f C 画线(.wa d wr i t e( 5,'

(

1 x, a) )

1 ,0,1 ,1,1

0 .0

0 0

0 0 0,1

0 .

00 000 0, 1, 0,

0 wr i t e ( 5,

4 0 O )

5 wr i t e ( 5,

2 0

0 ) xl , yl wr i t e (

5 ,

2 0

0 ) x 2, y l wr i t e (

5 ,

2 0

0 ) x 2, y

2 wnt e (

5 ,

2 0

0 ) xl , y

2 wnt e (

5 ,

2 0

0 ) x l, y l wr i t e ( 5, ¨(

1 x, i 6, '

, '

, f

3 .

1 ) ¨) n t , x

3 y d一0.yl= =y2 e n d i f