python读取grd数据_从等值线图形获取DAT和GRD格式的数据文件_王宜芳

中图分类号

:

P

208

文献标识码

:

B

文章编号

:

1672-

4097

(

2004

)

02-

0045-

02

1

引

言

在盆地模拟研究工作中

,

需要研究、

处理大量的

地质数据和图形

,

研究工作需要我们随时将各种数

据转化为可视的图形

,

这是我们经常在做的

,

也是驾

轻就熟的工作

;

而另一方面

,

也需要随时能够从图形

中得到数据

,

以供研究使用

,

例如

,

从构造图获取研

究区域的深度数据文件

,

这就是本文要讨论的问题。

以

从构造图获

取研究区域的

深度数据文件为

例

,

介绍这种方法。其工作流程如图

1

:

图

1

从构造

图获取深度数据文件的工作流程

2

主要工作流程

2

.

1

输入构造图

在

M

APGIS

的编辑子系统中

,

输入并编辑完成

构造图

,

输入和编辑方法显然不必多做介绍。现在

我们的计算机中有了这张构造图的

MAPGIS

矢量

文件

(

图

2

)

。

图

2

构造图

2

.

2

编辑线属性结构

在

MAPGIS

的编辑子系统中

,

用

/

线编辑

0

下的

/

编辑线属性结构

0

命令

,

在原有的

/

ID

0

、

/

线长度

0

的属性之后

,

增加一个

/

高程

0

属性。

2

.

3

编辑线属性

还是在

MAPGIS

的编辑子系统中

,

用

/

线编辑

0

下的

/

修改线属性

0

命令

,

在每

条深度等值线的

/

高

程

0

属性中输入该段等值线的高程值。

2

.

4

等值线抽稀

在现有的线文件中

,

删除所有其

它的断层线、

边框线等线段

,

只留下深度等值线

,

并另存为

/

gut

-

3

.

w

l

0

(

图

3

)

。这里需要

注意的是

,

数据点不

能过

密

,

用

/

设置参数设置坐标点可见

0

命令

,

看一下剩

下的等值线

,

如果坐标

点过密的

话

,

需要用

/

线编

辑抽稀线

0

命令

,

将等值线的数据点适当抽稀。

2

.

5

浏览等值线的高程

在

MAPGIS

的属性管理子系统中

,

打开上述线

文件

/

gut

-

3

.

w

l

0

,

浏览各等值线的高程值

,

看是否还

有线段没有输入高程值

,

如果有的话

,

可以在这里进

行编辑

,

输入其高程值

;

还要注意的是

,

不能有两条

以上的线段具有相同的

ID

号

,

否则在以后生成数