3.Mapgis大地坐标转经纬度和经纬度换大地坐标
地质工作中常要对进行大地坐标转经纬度和经纬度换大地坐标,以下步骤请大家熟记:
一、大地坐标→经纬度(地理坐标)
1、在文本文件中输入大地坐标数据,格式为Y空格X。
如下,原始的大地坐标由一个8位的Y和一个7位的X组成,“新建文本文档.txt -记事本”显示如下:
31560000 4503500
31565000 4503500
31565000 4507500
31568500 4507500
这组坐标数据中的Y的前两位为31,是分带号,一般使用的分带有三分带,六分带,这里的坐标是三分带的,记下Y前的这两位数,在原始数据中去除掉, 现在数据变为:Y—6位,X—7位。“新建文本文档.txt -记事本”显示如下: 560000 4503500
565000 4503500
565000 4507500
568500 4507500
保存这个TXT的文本文件。
2、打开MAPGIS,启动坐标投影变形程序
如果是MAPGIS6.7版,请选择“实用服务→投影变换系统→用户文件投影转换”,点击打开文件,打开刚才的大地坐标的文本文件。
“指定数据起始位置”中出现刚才的的文本文档,显示如下:
在设置用户文件选项中,一般选:按行读取数据,X→Y顺序,生成点。 接着设置输入数据的格式,点击用户投影参数,并完成设置:
坐标系类型——大地坐标系
投影类型——5:高斯克吕格投影
比例尺分母——1
椭球面高程——0
投影面高程——0
投影带类型——3度带或6度带
投影带序号——31
X,Y的平移均设0
这里我们的大地坐标为3度带的第31带,注意填好,坐标单位为米。
填好以后设置输出的格式,我们要求输出的是经纬度,点击“结果投影参数”,完成以下设置。
坐标系类型——地理坐标系
我们输出的经纬度的单位应该是DDDMMMSS.SS(度.分.秒),注意点击投影变换、写到文件,保存就大功告成了,注意:保存的文件要写上.TXT的后缀 最后,在文本文件中计算出的结果如下:
x= 560000 y= 4503500 xp=934234.357 yp=403950.255
x= 565000 y=4503500 xp=934607.202 yp=403948.893
x=565000 y=4507500 xp=934608.692 yp=404158.555
X=568500 y=4507500 xp=934837.762 yp=403157.535
xp为经度,934234.357就是93度42分34.357秒,
yp为纬度,403950.225就是40度39分50.255秒(纬度最大为90,所以没有三位数)
二、经纬度→大地座标
同样,输入文本文件格式如下,
1234234 403950
1234607 403948
1234608 404158
1234837 403157
这里面的数据前面的为经度,格式为DDDMMMSS.SS,后面的为纬度,格式为DDDMMMSS.SS,接下来的转换过程和大地坐标转换一样,只要将刚才的用户转换参数和结果转换参数交换即可,要注意分带号的确定,如果你不知道分带号,就应该先计算分带号,算法是
经度/3得到的整数为三度带的分带号
(经度+3)/6得到的整数为六度带的分带号
计算所得的结果格式如下
x=1234234 y=403950 xp=559991.687 yp=4503492.058
。。。。。。。。。。。。。。。。。。。。。。。
。。。。。。。。。。。。。。。。。。。。。。。
。。。。。。。。。。。。。。。。。。。
其中的xp为地图上的y坐标,记得在前面加上带号,
其中的yp为地图上的x坐标