导读：进行盆地模拟或其它需要地质格架的数值运算，常常需要从其它数据绘图软件中导出xyz的通用数据格式。以前我们写过双狐数据（dfd2data.exe）的导出，这里我们将关注国内另外一款比较常用的地质绘图工具——Geomap数据的导出，重点是跨带图件。

1、一般的，所需要的数据体可以通过Geomap自带的功能（数据制图|导出数据）来获取，其导出的数据格式（等值线）如格式1。根据这些格式，我们可以批量处理，把它变成xy_z的格式。这种方法对于没有跨带的图件是适用的。但如果图件是跨带的，而且以“大地坐标（格式2）”的方式导出数据，由于Geomap导出时没有自动换算到同一带上，这些数据实际不能直接使用（如格式1中所示的数据，不难发现，两条等值线的坐标是完全不同的，但它们在图上确为同一条等值线上的片段）。对于这种情况，我们先以“经纬度坐标”的方式导出，然后根据[1]中提到的坐标换算方法（相同坐标系下！），把坐标统一到同一坐标带上即可。根据这个思路，我们编写了Geomap坐标换算工具，使用时根据不同的坐标系（北京54，西安80），只需输入中央子午的经度即可（注：这些信息可以在文件|编辑图件信息中查看）。附件是针对西安80坐标系编写的，用于换算Geomap的等值线数据，适用于高斯投影，正轴等角圆锥投影参看方法2。

Geomap坐标换算工具（运行环境为Python 3.x）：附件下载 备用下载

操作步骤如下：

1、首先在Geomap中导出经纬度坐标的等值线数据

2、在Geomap中查看图件信息，获取球体参数[ IAG 75 (China 1980)]及中央子午

3、把附件中的文件（*.py,*.pyc）与要转换的数据文件放在同一目录，双击运行Geomap_contour data transfer.py，根据提示输入文件名（*.*）及中央子午

4、结果保存在：文件名_results.txt

格式1：

GmLine v3.0(Contour)   //文件头
3 500   //第一个数据是数据所包含的行数，第二个数是等值线的高程，即Z值
107 1 17.307889, 41 23 21.007318  // x坐标，y坐标；这里是地理坐标（经纬度） dms
107 0 55.726430, 41 23 15.842552
107 0 45.706968, 41 23 10.887126

格式2：

GmLine v3.0(Contour)
5  800
15728110.0 4798310.0 1
15728850.0 4799269.0 1
15729180.0 4799613.9 1
15729340.0 4799823.9 1
15729600.0 4800040.0 1
4 800
16267147.1 4808259.4 1
16266993.1 4808627.4 1
16266788.6 4808989.0 1
16266631.3 4809170.7 1

ps. Geomap导出的等值线数据有些缺失高程数据，而且看来没有好的补救方法。

2、利用双狐软件（2012年11月4日更新）

如果手头上有双狐绘图软件，则下列的方法将更为简单：

（1）新建一个空白的双狐图形文件

（2）将新建的双狐图形文件的投影信息设置成与对应的Geomap图形的地图投影信息一致

（3）在Geomap中，将需要导出的图形数据以经纬度坐标的方法导出

（4）将第3步导出的数据以经纬度坐标的方式导入双狐图形中

[1] 地理信息系统坐标系及转换的一些概念