利用mapgis软件实现经纬度和大地坐标转换.docVIP

利用mapgis软件实现经纬度（地理坐标）与大地坐标间的转换 MAPGIS的坐标系为数学坐标系，与投影平面直角坐标系中的X、Y坐标相反，即横坐标为X，纵坐标为Y，未经投影变化之前均为毫米表示。 MAPGIS的用户坐标系是指由用户指定的相对二维坐标系，一般与实际地物定位无关；地理坐标系是以经纬度表示的，经度的起点在格林威治，向东为正，纬度自赤道起，向北为正，常用来坐标定位；投影平面直角坐标系是将地球球面投影到平面后所设定的坐标系。我们常说的1954年北京坐标系，1980年西安坐标系均为高斯投影的投影平面直角坐标系，只不过它们采用了不同的椭球参数；北京坐标系使用克拉索夫斯基椭球，西安坐标系采用IAG1975年推荐椭球。 一、大地坐标→经纬度(地理坐标) 1、在文本文件中输入大地坐标数据，格式为Y空格X。 如下，原始的大地坐标由一个8位的Y和一个7位的X组成，“新建文本文档.txt －记事本”显示如下： 4503500 4503500 4507500 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 这里面的数据前面的为经度，格式为ＤＤＤＭＭＳＳ，后面的为纬度，格式为ＤＤＭＭＳＳ 接下来的转换过程和大地坐标转换一样，只要将刚才的用户转换参数和结果转换参数交换即可， 要注意分带号的确定，如果你不知道分带号，就应该先计算分带号，算法是 经度／３得到的整数为三度带的分带号 经度／６得到的整数为六度带的分带号 计算所得的结果格式如下 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.0