01第一章 地球椭球的基本公式

1地图投影MapProjections蒲英霞南京大学地理信息科学系2009年9月2第一章地球椭球体的基本公式3地球椭球体的形状和大小地球椭球面上的基本点、线、面和地理坐标子午圈、卯酉圈曲率半径及平均曲率半径地球球半径经线弧长纬线弧长地球椭球面上的梯形面积4大地水准面与大地体（Geoid）设想当海水面完全处于静止状态下，并延伸到大陆内部，使它成为一个处处与铅垂线（重力线）正交的连续的闭合曲面，这个曲面叫做大地水准面。由它所包围的球体，叫做大地体。GeoidNPSP§1.1地球椭球体的形状和大小大地水准面仍然不是一个规则的曲面，而是一个起伏不平的重力等位面，即地球物理表面。5高程基准是推算国家统一高程控制网中所有水准高程的起算依据，它包括一个水准基面和一个永久性水准原点。其中水准基面在理论上通常采用大地水准面，实际上是取验潮站长期观测结果计算出来的平均海面来确定的。中国以青岛港验潮站的长期观测资料推算出的黄海平均海面作为中国的水准基面，即零高程面。中国水准原点建立在青岛验潮站附近，并构成原点网。用精密水准测量测定水准原点相对于黄海平均海面的高差，即水准原点的高程，定为全国高程控制网的起算高程。大地水准面的作用——高程基准6地球椭球面与地球椭球体(Ellipsoid)选择一个大小和形状同大地水准面极为接近的，以椭圆短轴为旋转轴的旋转椭球面。这个旋转椭球面可代表地球的形状，又称为地球椭球面或参考椭球面（原面）。由它所围成的球体，称为地球椭球体或地球椭球。GeoidNPSPEllipsoid7地球椭球体的形状和大小扁率(FlatteningorCompression)eeeeabfa第一偏心率(FirstEccentricity)22212eeeabea第二偏心率(SecondEccentricity)22222eeeabebPOE1EP1Aaebe8第一偏心率和第二偏心率之间的关系：2221221eee2212211eeePOE1EP1Aaebe9椭球名称年代长半径（m）短半径（m）扁率附注埃弗斯特(Everest)1830637727663560751:300.8英国白塞尔(Bessel)1841637739763560791:299.15德国克拉克(ClarkeⅠ)1866637820663565841:295.0英国克拉克(ClarkeⅡ)1880637824963565151:293.5英国海福特(Hayford)1909637838863569121:297.01942年国际第一个推荐值克拉索夫斯基(Krassovsky)1948637824563568631:298.3苏联1967年大地坐标系1967637816063567551:298.2471971年国际第一个推荐值1975年大地坐标系1975637814063567551:298.2571975年国际第一个推荐值1980年大地坐标系WGS84197919846378137637813763567551:298.2571:298.2572235631979年国际第一个推荐值美国世界各国常用的地球椭球体数据101952年以前采用海福特椭球体从1953年起采用克拉索夫斯基椭球体1980年国家大地坐标系采用的地球椭球体（1975年大地坐标系）的参数：ae=6378140，扁率为1：298.2572008年7月1日起启用2000国家大地坐标系，具体参数为：ae=6378137，be=6356752.31414,扁率为1：298.257222101我国采用的椭球体：111980西安坐标系采用国际地理联合会（IGU）第十六届大会推荐的椭球参数，大地坐标原点在陕西省泾阳县永乐镇的大地坐标系，又称西安坐标系。1954年北京坐标系采用克拉索夫斯基椭球参数，又称北京坐标系。2000国家大地坐标系采用地心坐标系。12§1.2地球椭球面的基本点、线、面和地理坐标平行圈(parallel)子午圈(meridian):长半径为ae，短半径为be的椭圆纬线(parallel)地理经度(longitude)地理纬度PP1EE1起始经线FGaebe赤道r平行圈子午圈地理经度点两极(pole)线地理纬度(latitude)经线(meridian)面地理坐标13子午圈卯酉圈P1PEE1ADrQLW子午圈(PEP1E1)和卯酉圈(AQW)§1.3子午圈、卯酉圈曲率半径及平均曲率半径通过地面任一点的法线可以有无数法截弧，它们与椭球面相交则形成无数法截弧，其中有一对互相垂直的法截弧，称为主法截弧。主法截弧都是椭圆，其中一个是子午圈。与子午圈垂直的另一个圈称为卯酉圈。地球椭球面上的子午圈始终代表南北方向；卯酉圈除了两个极点外，代表东西方向。14子午圈曲率半径(radiusofcurvatureinmeridian)(M)：地球椭球体表面上某点法截弧曲率半径中最小的曲率半径。卯酉圈曲率半径(radiusofcurvatureinprimevertical)：地球椭球体表面上某点法截弧曲率半径中最大的曲率半径(N)。2122321(1)(1sin)eaeMe22121(1sin)eaNeP1PEE1ADrQLW子午圈(PEP1E1)和卯酉圈(AQW)15子午圈曲率半径(M)和卯酉圈曲率半径（N）之间的关系在赤道上：201(1)eMae0eNa在极点上：9090211eaMNe子午圈曲率半径(M)和卯酉圈曲率半径（N）除在两极处相等外，在其它纬度相同的情况下，同一点上卯酉圈曲率半径均大于子午圈曲率半径。2122321(1)(1sin)eaeMe22121(1sin)eaNe≤16曲率和曲率半径设一曲线的直角坐标方程是y=f(x)，且f(x)具有二阶导数。则该曲线的曲率为：232(1)yKy曲线上一点处的曲率半径与曲线在该点处的曲率互为倒数。xyy=f(x)MD17三轴平均球半径：椭球体三个半轴取平均数等面积球半径：球体的面积等于地球椭球体的面积等体积球半径：球体的体积等于地球椭球体的体积)(63711183mabaRe)(6371116)36017611(22meeaRa)(637111032mbaRv§1.4地球球半径18§1.5经线弧长设在经线上有一点A，其纬度为1，当A点沿经线移动一无限小距离达到A点时，纬度增加d，设A点的经线曲率半径为M，则经线上这一微分弧AA为MddsAAm'如果求1至2一段经线弧长，可将上式积分，即21MdsmPEP1BAA12dE119将子午圈曲率半径M带入上式，经整理，最后得：2121212121(1)()(sin2sin2)(sin4sin4)(sin6sin6)246meBCDsaeA24611124611146116134517511.00505177394642563155250.0056623776416512151050.000010624564256350.0000000208512AeeeBeeeCeeDe其中同纬差的子午线弧长由赤道向两极逐渐增长纬差1的子午线弧长在赤道为110576m,在两极为111695m。20§1.6纬线弧长在地球椭球体上，纬圈是圆。若求纬线上两点间经差为（2-1）的一段弧长，可由下式求得：cosnsrNNrP1PE1Esn同经差的纬线弧长由赤道向两极缩短经差1的纬线弧长在赤道为111321m，在纬度45处为78848m，在两极为0。21§1.7地球椭球面上的梯形面积设在地球椭球面上，有两条无限靠近的经线和两条无限靠近的纬线，其经度为和+d，纬度为和d。它们构成了一个微分梯形ABCD（如图所示）。这个微分梯形的边长为经线和纬线的微分弧长。因而有：MdBCADdNCDABcos于是微分梯形的面积为：ddMNdTcos经度为1和2的经线和纬度为1和2所围成的球面梯形的面积，可由下式积分求得：222211112212221coscos(1)(1sin)edTMNddaedeDACBdd22上式经过积分并加以整理，得：)7cos27sin5cos25sin3cos23sincos2sin(mmmmDCBAKT2212211212461112461114611612(1)()213511.003363605728161330.001124027261616310.0000016989801610.0000000027112meKaeAeeeBeeeCeeDe其中23实习作业1.计算纬度为3230N的纬圈上某一点的子午圈曲率半径（M）、卯酉圈曲率半径（N）和纬圈半径（r）。2.计算从赤道（0）至3230N的经线弧长。3.计算在3230N纬圈上从本初子午线（0）至11930E的纬线弧长。4.计算球面梯形（30N~35N，115E~120E)的面积。几点注意事项：1.比较克拉索夫斯基椭球体与2000国家大地坐标系中所采用的椭球体之间的区别；2.变量名称的定义统一；MajorRadius(a),MinorRadius(b),Eccentricity(e)MeridianRadius(M),PrimeVerticalRadius(N),ParallelRadius(r)MeridianLength(Sm),ParallelLength(Sn),Area(T)3.度分秒与弧度的换算；（DegreeMinuteSecondToRadian)4.基本函数BasicFun.c5.参数所在的头函数Parametr.h