MAPGIS的坐标系为数学坐标系,与投影平面直角坐标系中的X、Y坐标相反,即横坐标为X,纵坐标为Y,未经投影变化之前均为毫米表示。MAPGIS的用户坐标系是指由用户指定的相对二维坐标系,一般与实际地物定位无关;地理坐标系是以经纬度表示的,经度的起点在格林威治,向东为正,纬度自赤道起,向北为正,常用来坐标定位;投影平面直角坐标系是将地球球面投影到平面后所设定的坐标系。我们常说的1954年北京坐标系,1980年西安坐标系均为高斯投影的投影平面直角坐标系,只不过它们采用了不同的椭球参数;北京坐标系使用克拉索夫斯基椭球,西安坐标系采用IAG1975年推荐椭球。
一、大地坐标→经纬度(地理坐标)
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版,请选择“实用服务→投影变换系统→用户文件投影转换”→点击打开文件,打开刚才的大地坐标的文本文件。
“指定数据起始位置”中出现刚才的的文本文档,显示如下:
560000 4503500
565000 4503500
565000 4507500
568500 4507500
在设置用户文件选项中,一般选:按行读取数据,X→Y顺序,生成点。最后点击确定。
3、设置输入数据的格式,点击用户投影参数,并完成设置。
坐标系类型——大地坐标系
投影类型——5:高斯克吕格投影
比例尺分母——1
椭球面高程——0
投影面高程——0
投影带类型——3度带或6度带
投影带序号——31
X,Y的平移均设0
这里我们的大地座标为3度带的第31带,注意填好,坐标单位为米
接着为:设置输出的格式,我们要求输出的是经纬度,点结果转换参数,完成设置。
4、输入投影参数
坐标系类型——地理坐标系
我们输出的经纬度的单位应该是DDDMMMSS.SS注意点写到文件,保存就大功告成了,注意:保存的文件要写上.TXT的后缀
最后,在文本文件中计算出的结果如下:
x= 560000.000 y=4503500.000 xp= 934234.357 yp= 403950.255
x= 565000.000 y=4503500.000 xp= 934607.202 yp= 403948.893
x= 565000.000 y=4507500.000 xp= 934608.692 yp= 404158.555
x= 568500.000 y=4507500.000 xp= 934837.762 yp= 404157.535
xp为经度,934234.357就是93度42分34.357秒,
YP为纬度,403950.225就是40度39分50.255秒(纬度没有最多90,所以没有三位数)
二、经纬度(地理坐标)→大地坐标
同样,输入文本文件格式如下,
934234 403950
934607 403948
934608 404158
934837 403157
这里面的数据前面的为经度,格式为DDDMMSS,后面的为纬度,格式为DDMMSS
接下来的转换过程和大地坐标转换一样,只要将刚才的用户转换参数和结果转换参数交换即可,
要注意分带号的确定,如果你不知道分带号,就应该先计算分带号,算法是
经度/3得到的整数为三度带的分带号
经度/6得到的整数为六度带的分带号
计算所得的结果格式如下
x= 934234.000 y= 403950.000 xp= 559992.662 yp=4503569.573
x= 934607.000 y= 403948.000 xp= 564996.549 yp=4503549.933
x= 934608.000 y= 404158.000 xp= 564984.958 yp=4507560.321
x= 934837.000 y= 404157.000 xp= 568483.383 yp=4507560.911
其中的xp为地图上的y坐标,记得在前面加上带号,
其中的yp为地图上的x坐标