实用经纬度换算工具
- 格式:xls
- 大小:57.50 KB
- 文档页数:5


经纬度转换公式一、用EXCEL进行高斯投影换算从经纬度BL换算到高斯平面直角坐标XY(高斯投影正算),或从XY换算成BL(高斯投影反算),一般需要专用计算机软件完成,在目前流行的换算软件中,存在一个共同的不足之处,就是灵活性较差,大都需要一个点一个点地进行,不能成批量地完成,给实际工作带来许多不便。
笔者发现,用EXCEL可以很直观、方便地完成坐标换算工作,不需要编制任何软件,只需要在EXCEL的相应单元格中输入相应的公式即可。
下面以54系为例,介绍具体的计算方法。
完成经纬度BL到平面直角坐标XY的换算,在EXCEL中大约需要占用21列,当然读者可以通过简化计算公式或考虑直观性,适当增加或减少所占列数。
在EXCEL中,输入公式的起始单元格不同,则反映出来的公式不同,以公式从第2行第1列(A2格)为起始单元格为例,各单元格的公式如下:单元格单元格内容说明A2输入中央子午线,以度.分秒形式输入,如115度30分则输入115.30起算数据L0B2=INT(A2)+(INT(A2*100)-INT(A2)*100)/60+(A2*10000-INT(A2*100)*100)/3600 把L0化成度C2以度小数形式输入纬度值,如38°14′20〃则输入38.1420起算数据BD2以度小数形式输入经度值起算数据LE2=INT(C2)+(INT(C2*100)-INT(C2)*100)/60+(C2*10000-INT(C2*100)*100)/3600 把B化成度F2=INT(D2)+(INT(D2*100)-INT(D2)*100)/60+(D2*10000-INT(D2*100)*100)/3600 把L化成度G2=F2-B2L-L0H2=G2/57.2957795130823化作弧度I2=TAN(RADIANS(E2))Tan(B)J2=COS(RADIANS(E2)) COS(B)K2=0.006738525415*J2*J2L2=I2*I2M2=1+K2N2=6399698.9018/SQRT(M2)O2=H2*H2*J2*J2P2=I2*J2Q2=P2*P2R2=(32005.78006+Q2*(133.92133+Q2*0.7031))S2=6367558.49686*E2/57.29577951308-P2*J2*R2+((((L2-58)*L2+61)*O2/30+(4*K2+5)*M2-L2)*O2/12+1)*N2*I2*O2/2计算结果XT2=((((L2-18)*L2-(58*L2-14)*K2+5)*O2/20+M2-L2)*O2/6+1)*N2*(H2*J2)计算结果Y表中公式的来源及EXCEL软件的操作方法,请参阅有关资料,这里不再赘述。
利用MAPGIS制图软件换算大地坐标和经纬度地质工作中常要对进行大地坐标转经纬度和经纬度换大地坐标,以下步骤请大家熟记:一、大地座标→经纬度(地理坐标)1、在文本文件中输入大地坐标数据,格式为Y空格X。
如下,原始的大地坐标由一个8位的Y和一个7位的X组成,“新建文本文档.txt -记事本”显示如下:31560000 450350031565000 450350031565000 450750031568500 4507500这组坐标数据中的Y的前两位为31,是分带号,一般使用的分带有三分带,六分带,这里的坐标是三分带的,记下Y前的这两位数,在原始数据中去除掉,现在数据变为:Y—6位,X—7位。
“新建文本文档.txt -记事本”显示如下:560000 4503500565000 4503500565000 4507500568500 4507500保存这个TXT的文本文件。
2、打开MAPGIS,启动坐标投影变形程序如果是MAPGIS6.7版,请选择“实用服务→投影变换系统→用户文件投影转换”→点击打开文件,打开刚才的大地坐标的文本文件。
“指定数据起始位置”中出现刚才的的文本文档,显示如下:560000 4503500565000 4503500565000 4507500568500 4507500在设置用户文件选项中,一般选:按行读取数据,X→Y顺序,生成点。
最后点击确定。
3、设置输入数据的格式,点击用户投影参数,并完成设置。
坐标系类型——大地坐标系投影类型——5:高斯克吕格投影比例尺分母——1椭球面高程——0投影面高程——0投影带类型——3度带或6度带投影带序号——31X,Y的平移均设0这里我们的大地座标为3度带的第31带,注意填好,坐标单位为米接着为:设置输出的格式,我们要求输出的是经纬度,点结果转换参数,完成设置。
4、输入投影参数坐标系类型——地理坐标系我们输出的经纬度的单位应该是DDDMMMSS。