地图投影与高中斯投影

1、第七章地图投影与高斯投影本章提要本章介绍从椭球面上大地坐标系到平面上直角坐标系的正形投影过程。研究如何将大地坐标、大地线长度和方向以及大地方位角等向平面转化的问题。重点讲述高斯投影的原理和方法，解决由球面到平面的换算问题，解决相邻带的坐标坐标换算。讨论在工程应用中，工程测量投影面与投影带选择。知识点及学习要求1高斯投影的基本概念；2正形投影的一般条件；3高斯平面直角坐标与大地坐标的相互转换高斯投影的正算与反算4椭球面上观测成果（水平方向、距离）归化到高斯平面上的计算；5高斯投影的邻带换算；6工程测量投影面与投影带的选择。难点在对本章的学习中，首先要理解和掌握高斯投影的概念。高斯正算和反算计算；

2、方向改化和距离改化计算；高斯投影带的换算与应用；工程测量中投影面与投影带的选择。7.1高斯投影概述1投影与变形地图投影：就是将椭球面各元素（包括坐标、方向和长度）按一定的数学法则投影到平面上。研究这个问题的专门学科叫地图投影学。可用下面两个方程式（坐标投影公式）表示：xF(L,B)1yF2(L,B)式中L,B是椭球面上某点的大地坐标，而x,y是该点投影后的平面直角坐标。投影变形：椭球面是一个凸起的、不可展平的曲面。将这个曲面上的元素（距离、角度、图形）投影到平面上，就会和原来的距离、角度、图形呈现差异，这一差异称为投影变形。投影变形的形式：角度变形、长度变形和面积变形。地图投影的方式：（1）等

3、角投影投影前后的角度相等，但长度和面积有变形；（2）等距投影投影前后的长度相等，但角度和面积有变形；（3）等积投影投影前后的面积相等，但角度和长度有变形。1/172控制测量对地图投影的要求（1）应当采用等角投影（又称为正形投影）采用正形投影时，在三角测量中大量的角度观测元素在投影前后保持不变；在测制的地图时，采用等角投影可以保证在有限的范围内使得地图上图形同椭球上原形保持相似。（2）在采用的正形投影中，要求长度和面积变形不大，并能够应用简单公式计算由于这些变形而带来的改正数。（3）能按分带投影3高斯投影的基本概念（1）基本概念：如图1所示，假想有一个椭圆柱面横套在地球椭球体外面，并与某一条子午

4、线（此子午线称为中央子午线或轴子午线）相切，椭圆柱的中心轴通过椭球体中心，然后用一定投影方法，将中央子午线两侧各一定经差范围内的地区投影到椭圆柱面上，再将此柱面展开即成为投影面，如图2所示，此投影为高斯投影。高斯投影是正形投影的一种。图1图2（2）分带投影高斯投影6带：自0子午线起每隔经差6自西向东分带，依次编号1,2,3,。我国6带中央子午线的经度，由75起每隔6而至135,共计11带（1323带），带号用n表示，中央子午线的经度用L表示，它们的关系是L6n3,如图所示。00高斯投影3带：它的中央子午线一部分同6带中央子午线重合，一部分同6带的分界子午线重合，如用n表示3带的带号，L表示3带

5、中央子午线经度，它们的关系L3n图8-4所示。我国3带共计22带（2445带）。2/17（3）高斯平面直角坐标系在投影面上，中央子午线和赤道的投影都是直线，并且以中央子午线和赤道的交点O作为坐标原点，以中央子午线的投影为纵坐标x轴，以赤道的投影为横坐标y轴。xx500KmBxByByAAxAyBxByByAAxAy在我国x坐标都是正的，y坐标的最大值（在赤道上）约为330km。为了避免出现负的横坐标，可在横坐标上加上500OOOm。此外还应在坐标前面再冠以带号。这种坐标称为国家统一坐标。例如，有一点Y=19123456.789m，该点位在19带内，其相对于中央子午线而言的横坐标则是：首先去掉带

6、号，再减去500000m,最后得y=-376543.211m。（4）高斯平面投影的特点：5椭球面三角系化算到高斯投影面3/17将椭球面三角系归算到高斯投影面的主要内容是：（1）将起始点P的大地坐标(L.B)归算为高斯平面直角坐标x,y；为了检核还应进行反算，亦即根据x,y反算L,B。（2）通过计算该点的子午线收敛角及方向改正，将椭球面上起算边大地方位角APK归算到高斯平面上相应边PK的坐标方位角pK。（3）通过计算各方向的曲率改正和方向改正，将椭球面上各三角形内角归算到高斯平面上的由相应直线组成的三角形内角。（4）通过计算距离改正s，将椭球面上起算边PK的长度S归算到高斯平面上的直线长度s。（

7、5）当控制网跨越两个相邻投影带，需要进行平面坐标的邻带换算。7.2正形投影的一般条件高斯投影首先必须满足正形投影的一般条件。图1为椭球面，图2为它在平面上的投影。在椭球面上有无限接近的两点P和P，投影后为p和p，其坐标均已注在图上，dS为大1212地线的微分弧长，其方位角为A。在投影面上，建立如图2所示的坐标系，dS的投影弧长为ds。qlxyxyql图2图37椭球面到平面的正形投影一般公式称柯西-黎曼条件：xylq平面正形投影到椭球面上的一般条件：qlxy4/17sinBcos5B(6158t2t4)l66720NNcos5B(518t2t4142582t2)l5xXNN23224sinBls

8、inBcosB(5t9)l2224NNNxXyNcosBlcos3B(1t22)l37.3高斯平面直角坐标系与大地坐标系1高斯投影坐标正算公式（1）高斯投影正算：已知椭球面上某点的大地坐标L,B，求该点在高斯投影平面上的直角坐标x,y，即L,B(x,y)的坐标变换。（2）投影变换必须满足的条件中央子午线投影后为直线；中央子午线投影后长度不变；投影具有正形性质，即正形投影条件。（3）投影过程在椭球面上有对称于中央子午线的两点P和P,它们的大地坐标分别为（L,B）及12（l,B），式中l为椭球面上P点的经度与中央子午线(L)的经度差：lLL,P点在00中央子午线之东,l为正，在西则为负，则投影后的

9、平面坐标一定为P(x,y)和P(x,y)。12（4）计算公式ycosBlB(1t22)l3cos5B(518t2t4)l5631205当要求转换精度精确至0.OOlm时，用下式计算：NNsinBl2sinBcos3B(5t29244)l422244N6372052高斯投影坐标反算公式（1）高斯投影反算：已知某点的高斯投影平面上直角坐标x,y，求该点在椭球面上的大地坐标L,B，即x,y(L,B)的坐标变换。（2）投影变换必须满足的条件x坐标轴投影成中央子午线，是投影的对称轴；5/17l，最后得到BB(BB)、LLl。BBty2t(53t3292t2)y4f24MN32MNt(6190t245t4

10、)y6720MN5f11ly(12t22)y3NcosBf6N3cosBf1ffffffBBtyt(53t292t2)y42MN24MN3fly(12t22)y3NcosB6N3cosBff1120N5cosBx轴上的长度投影保持不变；投影具有正形性质，即正形投影条件。（3）投影过程根据x计算纵坐标在椭球面上的投影的底点纬度B，接着按B计算（BB）及经差fffff0（4）计算公式ffffffffffffffffff(528t224t46282t2)y5120N5cosBf当要求转换精度至0.01时，可简化为下式：f2f2ffffffff11ffff(528t224t4)y5ffff3高斯投影相

11、邻带的坐标换算（1）产生换带的原因高斯投影为了限制高斯投影的长度变形，以中央子午线进行分带，把投影范围限制在中央子午线东、西两侧一定的范围内。因而，使得统一的坐标系分割成各带的独立坐标系。在工程应用中，往往要用到相邻带中的点坐标，有时工程测量中要求采用3带、1.5带或任意带，而国家控制点通常只有6带坐标，这时就产生了6带同3带（或1.5带、任意带）之间的相互坐标换算问题，如图所示：6/17（2）应用高斯投影正、反算公式间接进行换带计算计算过程把椭球面上的大地坐标作为过渡坐标。首先把某投影带（比如带）内有关点的平面坐标(x,y)，利用高斯投影反算公式换算成椭球面上的大地坐标(l,B)，进而得到L

12、LIl；I0然后再由大地坐标(B,l)，利用投影正算公式换算成相邻带的（第带）的平面坐标(x,y)。II在这一步计算时，要根据第带的中央子午线LII来计算经差l，亦即此时lLLII。00算例在中央子午线LI123的带中，有某一点的平面直角坐标x5728374.726m，01y210198.193m，现要求计算该点在中央子午线LII129的第带的平面直角坐标。10计算步骤根据x,y利用高斯反算公计算换算B,L，得到B513843.9042，11111L1260213.1362。1采用已求得的B,L，并顾及到第带的中央子午线LII129，求得110l25746.864，利用高斯正算公式计算第带的直

13、角坐标x,yII为了检核计算的正确性，要求每步都应进行往返计算4子午线收敛角公式（1）子午线收敛角的概念如图所示，p、pN及pQ分别为椭球面p点、过p点的子午线pN及平行圈pQ在高斯平面上的描写。由图可知，所谓点p子午线收敛角就是pN在p上的切线pn与pt坐标北之间的夹角，用表示。在椭球面上，因为子午线同平行圈正交，又由于投影具有正形性质，因此它们的描写线pN及pQ也必正交，由图可见，平面子午线收敛角也就是等于pQ在p点上的切线pq同平面坐标系横轴y的倾角。（2）由大地坐标(L,B)计算平面子午线收敛角公式IIsinBl1sinBcos2Bl3(13224)3115sinBcos4Bl5(2t

14、2)（3）由平面坐标(x,y)计算平面子午线收敛角的公式7/17y23N3NfytanB1(1t22)ffff上式计算精度可达1。如果要达到0.001计算精度，可用下式计算：3N3Nfyty2fft(1t22)fffy515N5ft(25t23t4)fff（4）实用公式已知大地坐标(L,B)计算子午线收敛角1(0.333330.00674cos2B)(0.2cos2B0.0067)l2l2cos2BlsinB已知平面坐标(x,y)计算子午线收敛角1(0.333330.00225cos4B)(0.20.067cos2B)Z2Z2ZsinBfff7.4椭球面上观测成果归化到高斯平面上计算1概述由于

15、高斯投影是正形投影，椭球面上大地线间的夹角与它们在高斯平面上的投影曲线之间的夹角相等。为了在平面上利用平面三角学公式进行计算，须把大地线的投影曲线用其弦线来代替。控制网归算到高斯平面上的内容有：（1）起算点大地坐标的归算将起算点大地坐标(L,B)归算为高斯平面直角坐标(x,y)。（2）起算方向角的归算。（3）距离改化计算椭球面上已知的大地线边长（或观测的大地线边长）归算至平面上相应的弦线长度。（4）方向改计算椭球面上各大地线的方向值归算为平面上相应的弦线方向值。2方向改化（1）概念如图所示，若将椭球面上的大地线AB方向改化为平面上的弦线ab方向，其相差一个角值ab,即称为方向改化值。（2）方向

16、改化的过程如图所示，若将大地线AB方向改化为弦线ab方向。过A,B点，在球面上各作一大圆弧与轴子午线正交，其交点分别为D,E，它们在投影面上的投影分别为ad和be。由于是把地球近似看成球，故ad和be都是垂直于x轴的直线。在a,b8/17点上的方向改化分别为和。当大地线长度不大于10km，y坐标不大于l00km时，二abba者之差不大于0.05，因而可近似认为=。abba（3）计算公式球面角超公式为：R2(xy)ab(yy)ab2适用于三、四等三角测量的方向改正的计算公式：aby(xx)2R2may(xx)2R2mabbab2式中y1(yy)，为a、b两点的y坐标的自然的平均值。mab3距离改

17、化（1）概念如图所示，设椭球体上有两点P,P及其大地线S，在高斯投影面上的投影为P,P及s。1212s是一条曲线，而连接PP两点的直线为D如前所述由S化至D所加的改正，即为距离改正12S。（2）长度比和长度变形长度比m:指椭球面上某点的一微分元素dS，其投影面上的相应微分元素ds，则mdsdS称为该点的长度比。2长度变形:由于长度比m恒大于1，故称(m1)为长度变形。（3）长度比m的计算公式m1y22R2m式中：R表示按大地线始末两端点的平均纬度计算的椭球的平均曲率半径。my1(yy)为投影线两端点的平均横坐标值。mab9/17（4）长度比和长度变形的特点长度比m随点的位置而异，但在同一点上与

18、方向无关；当y=0(或l=0)时，m=1，即中央子午线投影后长度不变；当y0(或l0)时，即离开中央子午线时，长度设形（m-1）恒为正，离开中央子午线的边长经投影后变长。长度变形（m1）与y2（或l2）成比例地增大，对于在椭球面上等长的子午线来说，离开中央子午线愈远的那条，其长度变形愈大。（5）距离改化计算公式DS1y2m22Rm或DS12R224R2y2y2mmm7.5工程测量投影面与投影带选择1概述对于工程测量，其中包括城市测量，既有测绘大比例尺图的任务，又有满足各种工程建设和市政建设施工放样工作的要求。如何根据这些目的和要求合适地选择投影面和投影带，经济合理地确立工程平面控制网的坐标系，

19、在工程测量是一个重要的课题。2工程测量中选择投影面和投影带的原因（1）有关投影变形的基本概念平面控制测量投影面和投影带的选择，主要是解决长度变形问题。这种投影变形主要是由于以下两种因素引起的：实测边长归算到参考椭球面上的变形影响，其值为s：1RssH1m式中：H为归算边高出参考椭球面的平均高程，s为归算边的长度，R为归算边方向m参考椭球法截弧的曲率半径。归算边长的相对变形：10/17s值是负值，表明将地面实量长度归算到参考椭球面上，总是缩短的；s值与H,sH1msR11m成正比，随H增大而增大。m将参考椭球面上的边长归算到高斯投影面上的变形影响，其值为s：22Rs1ysm2m20式中：sss，

20、即s为投影归算边长,y为归算边两端点横坐标平均值，R为参考椭球010mm面平均曲率半径。投影边长的相对投影变形为2Rss1y2m0m2s值总是正值，表明将椭球面上长度投影到高斯面上，总是增大的；s值随着y平22m方成正比而增大，离中央子午线愈远，其变形愈大。（2）工程测量平面控制网的精度要求工程测量控制网不但应作为测绘大比例尺图的控制基础，还应作为城市建设和各种工程建设施工放样测设数据的依据。为了便于施工放样工作的顺利进行，要求由控制点坐标直接反算的边长与实地量得的边长，在长度上应该相等，这就是说由上述两项归算投影改正而带来的长度变形或者改正数，不得大于施工放样的精度要求。一般来说，施工放样的

21、方格网和建筑轴线的测量精度为1/50001/20000。因此，由投影归算引起的控制网长度变形应小于施工放样允许误差的1/2，即相对误差为1/100001/40000，也就是说，每公里的长度改正数不应该大于102.5cm。3投影变形的处理方法（1）通过改变H从而选择合适的高程参考面，将抵偿分带投影变形，这种方法通常m称为抵偿投影面的高斯正形投影；（2）通过改变y，从而对中央子午线作适当移动，来抵偿由高程面的边长归算到参考m椭球面上的投影变形，这就是通常所说的任意带高斯正形投影；（3）通过既改变H（选择高程参考面），又改变y（移动中央子午线），来共同抵偿mm两项归算改正变形，这就是所谓的具有高程抵

22、偿面的任意带高斯正形投影。4工程测量中几种可能采用的直角坐标系（1）国家3带高斯正形投影平面直角坐标系当测区平均高程在l00m以下，且y值不大于40km时，其投影变形值s及s均小于m122.5cm，可以满足大比例尺测图和工程放样的精度要求。，在偏离中央子午线不远和地面平均高程不大的地区，不需考虑投影变形问题，直接采用国家统一的3带高斯正形投影平面直角坐标系作为工程测量的坐标系。（2）抵偿投影面的3带高斯正形投影平面直角坐标系在这种坐标系中，依然采用国家3带高斯投影，但投影的高程面不是参考椭球面而是依11/17据补偿高斯投影长度变形而选择的高程参考面。在这个高程参考面上，长度变形为零。令smsR

23、y2m22RmH2ss0于是，当y一定时，可求得：1m则投影面高为：H投HHmHy2m2R算例：某测区海拔H=2000（m），最边缘中央子午线100（km），当s=1000（m）时，m则有22R2ms0.313(m)sms0.123(m)Rs1H1y22mm而ss0.19(m)12超过允许值（102.5cm）。这时为不改变中央子午线位置，而选择一个合适的高程参考面，经计算得高差：H780(m)将地面实测距离归算到：20007801220(m)（3）任意带高斯正形投影平面直角坐标系在这种坐标系中，仍把地面观测结果归算到参考椭球面上，但投影带的中央子午线不按国家3带的划分方法，而是依据补偿高程面归

24、算长度变形而选择的某一条子午线作为中央子午线。这就是说，在（8-173）式中，保持H不变，于是求得my2RHmm算例：某测区相对参考椭球面的高程H=500m，为抵偿地面观测值向参考椭球面上归m算的改正值，依上式算得y263700.580(km)即选择与该测区相距80km处的子午线。此时在y=80km处，两项改正项得到完全补偿。m但在实际应用这种坐标系时，往往是选取过测区边缘，或测区中央，或测区内某一点的子午线作为中央子午线，而不经过上述的计算。（4）具有高程抵偿面的任意带高斯正形投影平面直角坐标系在这种坐标系中，往往是指投影的中央子午线选在测区的中央，地面观测值归算到测区平均高程面上，按高斯正

25、形投影计算平面直角坐标。由此可见，这是综合第二、三两种坐标系长处的一种任意高斯直角坐标系。显然，这种坐标系更能有效地实现两种长度变形改正的补偿。（5）假定平面直角坐标系当测区控制面积小于100km2时，可不进行方向和距离改正，直接把局部地球表面作为平面建立独立的平面直角坐标系。这时，起算点坐标及起算方位角，最好能与国家网联系，如果联系有困难，可自行测定边长和方位，而起始点坐标可假定。这种假定平面直角坐标系只限于某种工程建筑施工之用。12/17本章作业1为什么要研究投影？我国目前采用的是何种投影？2控制测量对投影提出什么样的基本要求？为什么要提出这种要求？3椭球是一个不可展曲面，将此曲面上的测量

26、要素转换到平面上去，必然会产生变形，此种变形一般可分为哪几类？我们可采取什么原则对变形加以控制和运用?4高斯投影应满足哪些条件？60带和30带的分带方法是什么？如何计算中央子午线的经度？5为什么在高斯投影带上，某点的y坐标值有规定值与自然值之分，而x坐标值却没有这种区分?在哪些情况下应采用规定值？在哪些情况下应采用自然值？6正形投影有哪些特征？何谓长度比？7投影长度比公式的导出有何意义？导出该公式的基本思路是什么？8写出正形投影的一般公式，为什么说凡是满足此式的函数，皆能满足正形投影的条件？9学习了正形投影的充要条件和一般公式之后，你对高斯投影的实质是怎样理解的？10设ABC为椭球面上三等三角

27、网的一个三角形，试问：(1)依正形投影A、B、C三点处投影至平面后的长度比是否相等?(2)如若不等，还能保持投影的等角性质和图形相似吗？如若相等，岂不是长度比和点的位置无关吗?11写出按高斯平面坐标计算长度比m的公式，并依公式阐述高斯投影的特点和规律。12已知投影公式xf（B、L），yf（B、L），求椭球面上一点附近任意方向上1213/17长度比的计算公式，并写出主方向的长度比(提示：tgrdl)。MdB13在讨论高斯投影时提出了正形投影的充要条件(又称柯西黎曼条件)，它对问题的研究有什么作用？这个条件是如何导出的？14高斯投影坐标计算公式包括正算公式和反算公式两部分，各解决什么问题？15试述

28、建立高斯投影坐标正算公式的基本思路及主要过程。16高斯投影正算是已知求，由于值不大，故此公式可以认为是在点上展开的幂级数；反算公式中底点纬度Bf是指，由于值不大，故此公式可认为是在点上展开的幂级数。，17试证明高斯投影所求得的经线投影影像向中央子午线弯曲（凹向中央于午线）平行圈投影像向两极弯曲（凸向赤道）。18某点的平面直角坐标x、y是否等于椭球面上该点至赤道和中央子午线的距离?为什么?19什么是平面子午线收敛角？试用图表示平面子午线收敛角之下列特性：(1)点在中央子午线以东时，为正，反之为负；(2)点与中央子午线的经差愈大，值愈大；(3)点所处的纬度愈高，值愈大。20高斯投影既然是正形投影，

29、为什么还要引进方向改正？21高斯投影既然是一种等角投影，而引入方向改正后，岂不破坏了投影的等角性质吗？22试推导方向改正计算公式并论证不同等级的三角网应使用不同的方向改正计算公式。14/1723怎样检验方向改正数计算的正确性？其实质是什么？24椭球面上的三角网投影至高斯平面，应进行哪几项计算？并图示说明为什么？25试推导城市三、四等三角网计算方向改正值的计算公式，并分析所用概略坐标的精度。26已知距离改化计算公式为：Dsy2m2R2s若要求改正数的精度为107，问坐标的精度为多少（已知R=6370km，ym300km）？27回答下列问题：(1)试述高斯正形投影的定义；(2)绘图说明平面子午线收

30、敛角，方向改化和距离改化的几何意义；(3)写出大地方位角和坐标方位角的关系式；(4)估算（用最简公式和两位有效数字）高斯投影六度带边缘一条边长50KM的最大长度变形(m1)，己知y330Km。m28在高斯投影中，为什么要分带？我国规定小于一万分之一的测图采用60投影带，一万分之一或大于一万分之一的测图采用30投影带，其根据何在？29如果不论测区的具体位置如何，仅为了限制投影变形，统称采用30带投影优于60带投影，你认为这个结论正确吗？为什么？30高斯投影的分带会带来什么问题？31高斯投影的换带计算共有几种方法？各有什么特点？32利用高斯投影正、反算公式间接进行换带计算的实质是什么？已知某点在6带内的坐标为，x1944359.608m、y20740455.456m求该点在3带内第40带的坐标x、y。112215/1733在推导坐标换带表的换带公式中，对于对称点的选择有什么要求？对辅助点的选择又有什么要求？各起什么作用？34若已知高斯投