中国核动力研究设计院建筑设计所,四川成都, 610213
摘要:在设计工作中经常需要将地质勘察报告中的边界拐点坐标导入遥感卫星地图中显示,为后期设计作图提供基础资料,但大批量坐标转换、单点输入费时费力且缺少准确度,本文提供了大批量坐标一次转化生成遥感卫图的操作方法、并提供相关操作代码。
关键词1:大批量坐标 关键词2:遥感卫图 关键词3:操作代码
在设计工作中经常需要将地质勘察报告中的井田拐点坐标导入卫星地图(如奥维地图)中显示,从而为矿区总体规划提供基础资料。转换过程比较简单,现已某矿山边界拐点坐标为例,说明转换过程,具体过程如下:
一、整理原始数据
将地勘报告中的井田拐点坐标拷入新建的word文档中(出于对建设项目隐私,本文对坐标做了简单处理),内容如表1所示:
表1 某矿拐点坐标
北 纬 | 东 经 |
x1°x1′x1″ | y1°y1′y1″ ′11″ |
x2°x2′x2″ | y2°y2′y2″ ′11″ |
x1°x1′x1″ | y1°y1′y1″ ′11″ |
转换小程序中纬度在左,经度在右,所以需要将表中的经纬度左右互换位置,并在表格末尾加一行起始坐标使井田范围闭合,如表2所示:
表2 某矿拐点坐标
东 经 | 北 纬 |
y1°y1′y1″ ′11″ | x1°x1′x1″ |
y2°y2′y2″ ′11″ | x2°x2′x2″ |
y1°y1′y1″ ′11″ | x1°x1′x1″ |
二、数据转换
将表2中的数据在excel中按照下列公式转换为10进制并保留6位有效数字,另外还需增加高度一列,如表3所示:
10进制(保留6位有效数字)经纬度=LEFT(EXX,FIND("°",EXX)-1)+MID(EXX,FIND("°",EXX)+1,FIND("′",EXX)-FIND("°",EXX)-1)/60+MID(EXX,FIND("′",EXX)+1,FIND("″",EXX)-FIND("′",EXX)-1)/3600
注:上式中XX即为上表3中经度(或维度)所在excel表格的序号,如第四行则XX=4
表2 某矿拐点坐标转换
东 经 | 北 纬 |
111.111111° ′11″ | 222.222222° |
333.333333° ′11″ | 444.444444° |
….. | ….. |
注:上述坐标为示意处理 |
三、整理数据格式
将转换好后的坐标拷入新建的word文档中,修改成卫星地图(奥维地图)中认可的kml格式,如表4所示:
表4 转换后坐标
东 经 | 北 纬 | 高度 |
111.111111° ′11″ | 222.222222° | H |
333.333333° ′11″ | 444.444444° | H2 |
….. | ….. | .. |
注:上述坐标为示意处理 |
选择表格,在“布局”中选择“转化为文本”,“文件分隔符”选择逗号,如下所示:
111.111111 , 222.222222 , H1
333.333333 , 444.444444 , H2
…
…
选择以上数据,在“替换”中将带有空格的逗号“ , ”替换为不带空格的逗号“,”。在“查找和替换”的“特殊格式”中选择“段落标记”,将其替换为空格,整理好的数据格式如下所示:
111.111111 , 222.222222 , H1 333.333333 , 444.444444 , H2 …
四、生成卫星地图(如奥维地图)认可的程序
将第三步中的坐标拷入卫星地图(如奥维地图)认可的命令流中,并将命令流拷贝在一个新建的txt格式文件中,另存文件名修改为“某矿井.kml”,然后保存,命令流如下:操作示意图见图1:
<?xml version="1.0" encoding="UTF-8" standalone="no"?>
某矿.kmlnormal#sn_ylw-pushpinhighlight#sh_ylw-pushpin某矿#msn_ylw-pushpin1111.111111 , 222.222222 , H1 333.333333 , 444.444444 , H2 …
图1 命令流程序操作示意
五、在卫星地图(如奥维地图)中生成井田边界
电脑上安装卫星地图(如奥维地图)程序后,右键点击打开文件“某矿井.kml”文本文件后,即可直接在该软件中生成某井田边界。