雷 恩，孔祥福

（1.中国建筑材料工业地质勘查中心青海总队，青海 西宁 810008；2.成都理工大学地球科学学院，四川 成都 610059）

在地质工作中，经常需要从某一图面(地形图或地质图)上批量读取点的X、Y坐标值用于数据重新分析或是指导野外生产工作，如化探实际材料图上的采样点(分析数据点)或是水系沉积物测量采样设计点。如果手工一个一个在纸制图上用直尺读取点的坐标并录入计算机中，则工作量十分大，而且这种方法读取的坐标数据还存在较大的误差。为此，本文探讨基于MAPGIS软件平台上简易快速获取点的坐标值的途径，通过简单编辑就可为后续工作提供可靠的点位坐标。经笔者多年实践，有两种可靠方法可实现这一目标，第一种方法是在MAPGIS软件中的“图形处理/文件转换”子系统下实现坐标输出；第二种方法是在“实用服务/投影变换”子系统下完成。无论运用哪种方法，其前提都必须是底图数据已经完成校正，并生成了含有需要读取点图元坐标的点文件。

1 底图数据准备

底图数据可分为两种，一种是已经矢量化的点、线、区文件(*.wt、*.wl、*.wp)；另一种则是扫描图件经转换为MAPGIS可识别的MSI格式文件。这两种不同格式文件在作为底图数据使用前均需要校正到图幅的标准图框位置上。其中矢量化文件的校正在“实用服务/误差校正”子系统中完成，MSI文件的校正则在“图象处理/图象分析”子系统中完成，具体的操作步骤可参考文献[1]。

2 点图元的输入

输入点图元前还需要进行点属性结构的编辑操作。在工程管理窗口中新建名如“设计点.wt”点文件，点选该文件后，单击“点编辑→编辑点属性结构”，在弹出的对话框中进行属性结构编辑。除默认的“ID”字段外，还需添加如“采样点号、经度Y、纬度X”等字段名称，其中经度Y、纬度X的字段类型选择“双精度型或浮点型”，其字段长度可多设几位，便于显示。其中“ID”为输入点时的默认顺序号，“采样点号”用于赋存点的编号，“经度Y、纬度X”用于赋存点的坐标值。

完成以上设置，在工程管理窗口中添加底图文件后，便可输入点图元了。单击“点编辑→输入点图元”，在弹出的窗口中选择“子图”输入类型，根据需要设置“子图号、子图高度、子图宽度及子图颜色”，同时点选右下角的“即时属性输入”，单击“确定”之后，输入点图元时便可即时对其进行系统编号(赋存于“采样点号”属性栏中)，所有点图元输入完成后保存文件。

MAPGIS系统中，图幅的数据单位是以“毫米”表示的，因此直接从图幅中读出的坐标值是按当前图幅的比例尺缩小表示的。为了使从MAPGIS系统中读取的以“毫米”为单位的图幅坐标数值与以“米”为单位的实际坐标值一致，必须将图幅的比例尺变换成1∶1000(毫米)，变换之后则与实际坐标值1∶1(米)一致了。以H48E017024标准图幅为例，该图幅实际比例尺为1∶5万，变换成1∶1000，则需要放大50倍。单击“其它→整图变换→键盘输入参数”，在弹出的“图形变换”窗口中，勾选“点、线、面变换及参数变换”，“比例参数”中“X比例、Y比例”各输入“50”，单击“确定”，保存项目，此时整图就变成1∶1000的比例尺，图上读取的以“毫米”为单位的坐标值与以“米”为单位的实际坐标值一致(缺省投影带号)。当批量读取坐标值后，再用“整图变换”功能将图幅变换(缩小)回原来比例尺。

3 坐标输出的实现

点位坐标值的批量输出有两种途径，第一种是在MAPGIS软件中的“图形处理/文件转换”子系统下实现坐标转出；第二种则是在“实用服务/投影变换”子系统下完成。

3.1 “文件转换”子系统中输出坐标值

进入“图形处理/文件转换”子系统，单击“文件→装入点”，选择需要读取坐标值的点文件，接着单击“输出→输出MAPGIS明码格式”，单击“确定”，保存为“*.wat”格式的数据文件。用EXCEL软件打开该文件，显示结果如下：

其中，第一行“WMAP9022”为MAPGIS下点明码文件头标志，第二行为点图元的总个数，从第三行开始为一列以逗号“，”分隔的数据，对其进行数据分列，第一列便是“横坐标Y”，第二列为“纵坐标X”，第三列为“ID”。其他的列分别为点类型(1表示子图)、子图号、子图高度、子图宽度、旋转角度、辅助色、子图颜色、线宽、图层号、透明输出等信息。把前两列数据取整数、提取保存即可得到批量点坐标。但这种方式不能把点图元的“采样点号”属性同时输出。

3.2 “投影变换”子系统中输出坐标值

进入“实用服务/投影变换”子系统，单击“文件→打开文件”，打开需要读取坐标值的点文件。单击“工具→点位置转换为属性”，选择“图元文件”，选择打开的点文件，点选“[结]点位置X、[结]点位置Y”，并在下拉框中分别选择“经度Y、纬度X”，依次点击“转换、确定”，便把坐标赋到点图元属性上了。之后单击“工具→属性生成文本文件”，选择“图元文件”并命名“生成文件”，依次点击“转换、确定”，完成属性转出文本操作。这种方法在转出点图元坐标的同时，还能转出其他如“采样点号”等属性。至此，通过MAPGIS软件系统中的简单编辑，便可快速地获取批量点的坐标值。

4 应用及效果

笔者于2017年参加了“青海省都兰县金水口地区1∶2.5万地球化学测量”投标书的编写，需要在420km2范围内进行水系沉积物测量的采样工作。开展野外采样工作前，在室内先进行了地理底图基础上的水系沉积物采样点位设计并对其按公里网方格系统的编号，再运用MAPGIS系统下的“实用服务/投影变换”子系统，即上文所述的第二种方法，快捷地获取了工作区内的7 000余个水系沉积物测量设计采样点的理论坐标值。在野外实际采样工作时，运用这些设计采样点的理论坐标值进行GPS导航或是位置参照，能快速地定位采样点并完成采样工作，同时能使采样点分布更加均匀合理，与设计的采样位置更趋一致，提高了工作效率，取得了很好地运用效果。

5 结语

前人也探讨过MAPGIS读取地形图点坐标的方法[2]，但文中所用方法读取的坐标值(横坐标3位整数，纵坐标4位整数)不能直接用于实际工作。因此，MAPGIS系统中要批量输出地图点位坐标，关键是要进行“整图变换”操作。通过“整图变换”，把底图数据文件的比例尺变换成1∶1000，这样从MAPGIS系统中读取的以“毫米”为单位的图幅坐标数值便与以“米”为单位的实际坐标值一致了。经MAPGIS系统自动读取的横坐标(经度Y)整数位为6位，是没有加上图幅投影带号的，需要统一添加，这在EXCEL软件中可轻松完成。综合比较两种方法，在“实用服务/投影变换”子系统下转出的坐标值更好利用，含有除ID号外的自己系统的编号，以便于与图上位置对应并进行管理。

【参考文献】

[1]武汉中地信息工程有限公司.MAPGIS地理信息系统实用教程[M].湖北武汉，2003.

[2]火克俭.基于MAPGIS读取地形图点坐标方法[J].中国非金属矿工业导刊，2007(1)：63-64.