hello~~今天yogurt要和大家分享的是如何玩转地图投影的公式,由正解公式求解:(Lat,Lon)-->(E,N),或者由反解公式求解:(E,N)-->(Lat,Lon)。关于怎么进行投影?大家可以参考我之前写过的一篇《.gen地图文件的投影编程实现(以墨卡托投影和兰伯特投影为例)》http://www.cnblogs.com/to-sunshine/p/6048438.html。接下来,yogurt先给大家简单介绍一下地球椭球及其参数,还有兰伯特投影和墨卡托投影公式,让大家对于投影这个概念和过程有个大致地了解。然后,通过两个例子分别对兰伯特投影和墨卡托投影进行正反算求解。
============================yogurt小课堂开课了===========================
什么叫做地球椭球?我们都知道地球是一个两边宽上下扁的不规则椭球体。那么为了简化计算,我们把这个椭球体近似看作一个标准椭球体,它拥有以下参数:
长半径:a ; 短半径:b ;
扁率:
第一偏心率: 第二偏心率:
扁率与偏心率的关系为:e ² = 2 X f - f ² 。
当然,地球椭球体与普通标准椭球体还是有一定区别的,下面的参数则是地球椭球体特有的: 首先,我们要来了解这两个概念:
法截面:设过椭球表面上任意一点A做法线AL,通过法线的平面所截成的截面;
子午圈截面:包含子午圈的椭球体截面; --> 对应曲率半径 (该半径对应的圆弧是与所有经线相平行的)
卯酉圈截面:垂直于子午圈截面的法截面。 --> 对应曲率半径(该半径cos之后得到的半径是对应的纬线圈上的半径,对应的圆弧是与所有纬线相平行的&