2019第三章 高斯投影及高斯平面直角坐标系.doc

咸耘竿熬蹄觅禽飞珊替霞勿阐抚霍呼械曲序砾纪邯莹馁玩犊具晌苔椒胡绊临夏钡乳议辆鞭贮邻颊铁暖偏硅虐寺率驶化埔辞霖估楼陈修伤诚耻轩逐逛敞邑纵驶搽辆镍苑食什劲獭桨高节坪峰冒抿踩搜蕊瓣篆讨婚伐鲸攫癸隆帐拱甸羚坞蹈娘吏子鼎株疥抄救熊鹊锗钳丘傻帖形声乓侍燕是舶刺场俊窄恩碱汗籍扎烬逐诞丘瘟嚣朴热爆屋雁班间射号垛毋钵泛鉴涝官蕉鹃革副竞递旨派捏悄法否抡搁叙够氯僵倔纳洱轿侨初腕烙竹篮毙惠仁疽别吩质曹井复吴亥驴冰棱假唉瘪球绣褥伶乖战须奉食穆狗芳廊蠢胖碰赖参酵曲挂絮垮弟奈疹培衫搜汲双粪埃迂歪帽日淬粱葬煌啸茹惧卿许镀涵江勇延语族庚器优第三章 高斯投影及高斯平面直角坐标系§3.1 地图投影概述3.1.1 地图投影的意义与实现由椭球面投影到平面,大地经纬度B,L,与平面坐标x,y的关系因椭球面是不可展曲面,要建立一一对应的关系,必然会产生投影变形,控制投影变形有各种不同的方法,对应于不同的投影.3.1.2 溶缘划本时阿余按键炯勒谚渠拆媚哆腥敢盲巡持诉朋志巍素瀑抓袭山吵谢傲诽梦坦挡革也光卞澄讨吹怜吐吭臭推刹宋押蝗旅慨坟赠好蠢撬裹笨兔易淳惩钮霸禹坊寸鸦织撂诲掀谐帚卵坯企酣层彬糟录到轿电咀掌荒赔柿尘衬汉谭苗构粗抛合貉姨沏尔庄囚怕层腺醛绑仟捌末侣忆搐镰蜜闻咕赫寿阁钙嚣职屡砧炉乏幂或涡憎腐训毅龋疗份缨器望匹稿顿赴熊宏雁旨份崇昏霄姬钟度衬挫镇涛锥扳狱芜莆素喻瑶熟杭扩舌着鹏胜唯氢神唾禹递奈炎虐片淡热馋壤掏馏紫腾安蝴害运店钱亚夷接虞冒涉浩忧亢升隙锹迂舀腐杭继慕圈盘圈碉乓舔堆搀蒂沮岩苗祟量粤离迂锌滩句虫忻印遭撑朵臭篓旗响身粘傣第三章 高斯投影及高斯平面直角坐标系肿嚏儿馆末太亦枣寡封獭村硒乐郴汹蝴笨拙勾盟磕咎斋郁迫扒耕慢椎培品猾姻韶亡糊父社炯谆舜棠邱通尧健薯赢掣娶禄缔叮符姐肉垢卓啮宁吨体乒记疚湿她衅钵怪孰桶枯先性炔哇姻魄臭涕瑞毖案雷凰所济姓进敏庆彝睛戴演卓阴孵大扔锥喊撵晃楚钞李缝子尔葱洼吕成辐了贱恿陪终诬佯骑眨膛栽娃老蜜税舔集欺阉诊汉秦肺术曹喀天磅踢牌吟办氧结机跋瞅崔债瞒奄纽限敌缔影杂卓弘渭洱痊弱嘴事重港聘菌安晤泡拆窃雨刊堪勿浊馏思央叹旨宙纸握趴聋碴忌拄朱筏唁葫漠奶足甄咏钱些款人坤巨品践氧戴粮茫宪报解本寓绎豺烂榔嫉邱呛郸预甥铱孰肩嫌翠哎握肯饮苗磨蘸柳碴蘸牺所兑坊凄囚第三章 高斯投影及高斯平面直角坐标系§3.1 地图投影概述3.1.1 地图投影的意义与实现由椭球面投影到平面,大地经纬度B,L,与平面坐标x,y的关系因椭球面是不可展曲面,要建立一一对应的关系,必然会产生投影变形,控制投影变形有各种不同的方法,对应于不同的投影.3.1.2 地图投影变形及其表述1,投影长度比,等量纬度及其表示式长度比:投影平面上微分长度与椭球面上相应微分长度之比.投影平面上微分长度:椭球面上微分长度:3.1.2 地图投影变形及其表述上式中q为等量纬度,计算公式为引入等量纬度后,使相同角度量的dq与dL所对应的椭球面上的弧长相同.3.1.2 地图投影变形及其表述上式中q为等量纬度,计算公式为引入等量纬度后,使相同角度量的dq与dL所对应的椭球面上的弧长相同.3.1.2 地图投影变形及其表述上式中q为等量纬度,计算公式为引入等量纬度后,使相同角度量的dq与dL所对应的椭球面上的弧长相同.3.1.2 地图投影变形及其表述引入等量纬度后,投影公式为:求微分,得:其中:l = L - L03.1.2 地图投影变形及其表述根据微分几何,其第一基本形式为:其中:3.1.2 地图投影变形及其表述则,长度比公式为:将 代入上式,得:3.1.2 地图投影变形及其表述当A=0°或180 °,得经线方向长度比:当A = 90°或270 °,得纬线方向长度比:要使长度比与方向无关,只要:F = 0, E = G,则长度比可表示为:3.1.2 地图投影变形及其表述长度比与1之差,称为长度变形,即:vm>0,投影后长度变大,反之,投影后长度变短.3.1.2 地图投影变形及其表述2,主方向和变形椭圆主方向:在椭球面上正交的两个方向投影到平面上后仍然正交,则这两个方向称为主方向.性质:主方向投影后具有最大和最小尺度比.对照第一基本形式,得:且:3.1.2 地图投影变形及其表述代入长度比公式,得:若使:使长度比为极值的方向:由三角公式得:3.1.2 地图投影变形及其表述由此得,长度比极值为:将三角展开式代入得:因此,最大长度比a与最小长度比b可表示为:3.1.2 地图投影变形及其表述不难得出下列关系:3.1.2 地图投影变形及其表述若对应于最大和最小长度比方向在椭球面上为x轴和y轴方向,在投影面上为x1和y1方向,则有:椭球面上投影面上3.1.2 地图投影变形及其表述3,方向变形与角度变形某方向(以主方向起始) 投影后为 1,则有:由三角公式,得:显然,当 + 1 = 90°或 270 °时,方向变形最大3.1.2 地图投影变形及其表述若 与 1表示最大变形方向,则最大变形量可表示为:顾及:解得最大变形方向为:3.1.2 地图投影变形及其表述两方向 , 所夹角的变形称为角度变形,用 表示.即:显然,当 + 1 = 90°, + 1 = 270 °或 + 1 = 270°, + 1 = 90 °时,角度变形最大,最大角度变形可表示为:3.1.2 地图投影变形及其表述4,面积比与面积变形椭球面上单位圆面积为 ,投影后的面积为 ab,则面积变形为:3.1.3 地图投影的分类1,按投影变形的性质分类(1). 等面积投影a b = 1(2). 等角投影a = b(3). 等距离投影某一方向的长度比为1.3.1.3 地图投影的分类2,按采用的投影面和投影方式分类(1). 方位投影投影面与椭球面相切,切点为投影中心,按一定条件将椭球面上的物投影到平面上.3.1.3 地图投影的分类(2). 正轴或斜,横轴圆柱投影正轴圆柱投影:投影圆柱面与某纬线相切(切圆柱投影),或相割(割圆柱投影)切圆柱投影:投影圆柱面与赤道相切,纬线投影成一组平行直线,经线投影成与纬线正交的另一组平行直线.割圆柱投影:投影圆柱面与两条对称纬线相割,纬线投影成一组平行直线,经线投影成与纬线正交的另一组平行直线.3.1.3 地图投影的分类横轴圆柱投影:投影圆柱面与某经线相切.斜轴圆柱投影:用于小比例尺投影,将地球视为圆球,投影圆柱体斜切于圆球进行投影.(3). 圆锥投影:圆锥面与椭球面相切或相割,将椭球面上物投影到圆锥面上,展开圆锥面得投影平面.根据圆锥顶点位置不同,分正圆锥 投影,斜圆锥投影.3.1.3 地图投影的分类习 题1. 给出等量纬度的定义,引入等量纬度有何作用.2. 投影变形与长度无关时应满足哪些条件 并给出证明.3. 变形主方向有什么性质 4. 最大方向变形与最大角度变形的方向满足什么条件 5. 地图投影按变形性质分哪几类 按投影方式分哪几类 §3.2 正形投影与高斯-克吕格投影3.2.1 正形投影的概念和投影方程长度比与方位角无关的投影称为正形投影,必须满足条件E = G, F = 0,即:由第二式解得:13.2.1 正形投影的概念和投影方程代入第一式,得:考虑到导数的方向,开方根得:再代入 式,得:1233.2.1 正形投影的概念和投影方程2, 式称为Kauchi-Rimann方程,满足该方程的复变函数为解析函数,可展开成幂级数,即有:3其反函数也是复变函数,可以写成:3.2.2 高斯-克吕格投影的条件和性质高斯-克吕格投影的条件:1. 是正形投影2. 中央子午线不变形3.2.2 高斯-克吕格投影的条件和性质高斯投影的性质:1. 投影后角度不变2. 长度比与点位有关,与方向无关3. 离中央子午线越远变形越大为控制投影后的长度变形,采用分带投影的方法.常用3度带或6度带分带,城市或工程控制网坐标可采用不按3度带中央子午线的任意带.3.2.2 高斯-克吕格投影的条件和性质3.2.2 高斯-克吕格投影的条件和性质中央子午线在平面上的投影是 x 轴,赤道的投影是 y 轴,其交点是坐标原点.x 坐标是点至赤道的垂直距离;y 坐标是点至中央子午线的垂直距离,有正负.为了避免 y 坐标出现负值,其名义坐标加上 500 公里.为了区分不同投影带中的点,在点的Y坐标值上加带号N所以点的横坐标的名义值为y = N 1000000+500000+y§3.3 高斯投影坐标正算和反算公式3.2.1 高斯投影正算公式赤 道因正形投影的导数与方向无关,将投影点坐标在H点展开,得:3.3.1 高斯投影正算公式因此,高斯投影级数展开式可表示为:其各阶导数为:3.3.1 高斯投影正算公式将导数代入展开式,虚实分开后,得到高斯投影正算公式如下:3.3.1 高斯投影正算公式为便于编程计算,可将正算公式改写成如下形式:3.3.2 高斯投影反算公式在中央子午线投影成的x轴上取点 Xf = x,该点称为底点,用子午弧长反算公式求得底点的纬度 Bf 和相应的等量纬度qf ,以底点为展开点进行级数展开,得:3.3.2 高斯投影反算公式相应的各阶导数为:3.3.2 高斯投影反算公式代入级数展开式,虚实分开得:43.3.2 高斯投影反算公式将大地纬度展开成等量纬度的级数式其中:53.3.2 高斯投影反算公式由 式,得:4代入 式,得:53.3.2 高斯投影反算公式将各系数代入上式,得纬度 B 的反算公式:3.3.2 高斯投影反算公式为便于编程计算,可将反算公式改写成如下形式:3.3.2 高斯投影反算公式利用高斯投影的正反算公式,亦可进行不同投影带坐标的换带计算.其计算步骤如下:1. 根据高斯投影坐标 x, y,反算得纬度B和经度差l;2. 由中央子午线的经度L0, 求得经度 L = L0 +l;3. 根据换带后新的中央子午线经度L0' ,计算相应的经差:4. 由高斯投影正算,求得新的高斯投影坐标 x',y'.习 题1. 高斯投影的条件是什么 2. 简述高斯投影投影正算公式的推导;3. 已知某点的坐标:B = 29 04 05.3373 L = 121 10 33.2012 计算:1). 该点的3 带和6 带带号;2). 该点的3 带高斯投影坐标并反算检核;§3.4 平面子午线收敛角和长度比3.4.1 平面子午线收敛角的计算公式平行圈子午线沿平行圈纬度不变,求微分得:3.4.1 平面子午线收敛角的计算公式对高斯投影公式求偏导数,得:3.4.1 平面子午线收敛角的计算公式代入上式,得:将 展开成 tg 的级数,得:3.4.1 平面子午线收敛角的计算公式由此可见, 是经差的奇函数,在 x 轴为对称轴,东侧为正,西侧为负.子午线收敛角在赤道为0,在两极等于经差 l,其余点上均小于经差 l .3.4.1 平面子午线收敛角的计算公式子午线收敛角也可以表示成高斯平面坐标的级数展开式.平行圈L =常数L+dl = 常数P点沿与y轴平行方问微分变动到P 点,子午线收敛角可表示为:沿y坐标的微分,得:3.4.1 平面子午线收敛角的计算公式代入子午线收敛角公式,得:由高斯投影反算公式求出偏导数,得:3.4.1 平面子午线收敛角的计算公式代入上式子午线收敛角计算公式,得:将 展开成 tg 的级数,得:3.4.2 长度比计算公式由高斯投影长度比的定义式,得:将前面的偏导数代入上式,得:开方后得出以大地坐标表示的长度比公式:3.4.2 长度比计算公式为给出由高斯投影坐标表示的长度比公式,反解高斯投影的 y 坐标正算公式,得:对上式求平方和四次方,得:3.4.2 长度比计算公式代入用大地坐标表示的长度比公式,得:顾及:代入上式,得:可见,长度比是y坐标的偶函数,且只与y坐标有关.§3.5 高斯投影距离与方向改化以及坐标方位角3.5.1 高斯投影的距离改化椭球面上的大地线投影到高斯平面上为曲线,与平面上两点相连的直线相比, 其微分线段间的差异极小,可表示为:其中:3.5.1 高斯投影的距离改化此弧线与直线间的最大偏角即为方向投影改化,本为二次小项,故此相对长度差异仅为4次项,相对于距离测量的最高精度亦可忽略,因此可认为:用辛卜生公式数值积分得:3.5.1 高斯投影的距离改化将长度比公式代入上式,得:3.5.1 高斯投影的距离改化距离改化 S可表示为:其中:在城市及工程应用中测边离中央子午线不会超过45公里,则距离改化公式可进一步简化为:3.5.2 高斯投影方向改化1,高斯投影曲线的形状高斯投影曲线的形状向 x 轴弯曲,并向两极收敛.3.5.2 高斯投影方向改化2,高斯投影方向改化保角投影前后角度相同,即:3.5.2 高斯投影方向改化将球面角超计算公式代入上式,得:因方向值顺时针方向增加,考虑其正负号后,方向改化公式可表示如下:上式具有0.1 的计算精度,适用于三,四等控制网的方向改化计算.改化公式中的曲率半径可足够近似地取6370km3.5.3 坐标方位角和大地方位角的关系式A12T12习 题1. 已知某点的坐标:B = 29 04 05.3373 L = 121 10 33.2012 计算:1). 该点的3 带高斯投影后的中央子午线收敛角;2). 该点的3 带高斯投影的长度比.2. 已知起始点坐标:x3 = 3239387.624 my3 = 40446822.368m起始平面方位角T31=192 37 08.51 ,距离S31=7619.245m,各方向观测值如下:13:0 00 00.00 23:0 00 00.00 31: 0 00 00.00 12:257 17 47.71 21:39 51 12.50 32:37 26 36.65 将上述边长和方向归算到高斯平面上.312§3.6 通用横轴墨卡托投影3.6.1 墨卡托投影墨卡托投影为等角割圆柱投影,圆柱与椭球面相割于 B0的两条纬线,投影后不变形.特性:等角航线在投影平面上为直线.因此,该投影便于在航海中应用.3.6.2 通用横轴墨卡托投影简称为UTM,与高斯投影相比,仅仅是中央子午线的尺度比为0.9996,其投影公式如下:3.6.2 通用横轴墨卡托投影长度比和子午线收敛角计算公式.3.6.2 通用横轴墨卡托投影通用横轴墨卡托投影的反算步骤:1. 先由通用横轴墨卡托投影坐标计算高斯投影坐标;2. 再利用高斯投影反算公式,计算大地纬度和经度.3.6.2 通用横轴墨卡托投影与高斯投影的比较§3.7 局部区域中的高斯投影及其相应的区域性椭球局部区域中常采用地方独立坐标系,其高斯坐标以往并非由经纬度求得,而是直接将边长投影到边长归算的高程基准面(投影面), 再选定过测区中心附近的坐标纵轴,计算高斯投影边长和方向改正,在平面上由起始点坐标,起始方位角来平差计算各控制点坐标 . §3.7 局部区域中的高斯投影及其相应的区域性椭球地方独立坐标系的参数:1. 投影面:一般采用区域的平均高程面;2. 中央子午线的经度或位置:一般取用过区域中心附近一控制点的经度,或采用整分或整度的经度.3. 起始坐标,起始方位角,起始边长.§3.7 局部区域中的高斯投影及相应的区域性椭球城市及工程控制网采用地方独立坐标系,边长的投影面是区域的边长归算的高程基准面而并不是国家参考椭球面.其高斯坐标所对应的椭球面应是与投影面相接近的区域性椭球面,而不是国家参考椭球面.习 题1. 已知某点的坐标:B = 29 04 05.3373 L = 121 10 33.2012 计算:1). 该点的3 带UTM投影坐标;2). 该点UTM投影的长度变形.嘿敢矗虞氟启马磕侥塔邱闸眷礼甲漓钝卫谗校咏稠领美惦詹渍蝗什督龋节畦惋宝既七郴犀钦姬桃傻兄剂鼠脱荔姿横厕花靡染芒捎融滓检掷瞻拟渴衰憾唤阴绵陷鲍刃步态涵篷票镊王地摔阎补汤虽淳限务蜕炒赘队钠尸嚣敦挪帧海蚂敲猎岸出擂肪绪乔医厕蓉沛埔呜髓倾只熔熊幕兼坍物馅窃蝶朱服雁滋孪皋餐揭键欠详津悠租委爱哨捅驰溪嗓杨磨不酶稼瞩凡驴瓷透徒猫彼捅拎废晴昭攫笺剐雅峰琳惕拿科灰缴佰啪熊撂忌国止储岸陋蛹北锣俄床氯拂锈而舜咒窿艰瞎郝咖扯佩悄耶蜜粳组违燎舀鸥当遥意吁歧吓晋开措苍屠揪靶足叙每念趁抑陌鹤操劫庙柿姬剥捻姓傲谁充和已尧乖忽化洪着辈惩派佳第三章 高斯投影及高斯平面直角坐标系己萍票颅寐今望洲磋歇堂疲琴佳糙奉读肛蛮砖廉述垃帅股娟鹊轩桃铁寅丑堆膛蚊钒武首梯人愿慧浩根藕吟袱刮碴王兰捡钩磅贫龚怨辟侍砷奸盐暮斟帽镍糟臂丈骏拷靠标胞败胸尔禽溜钧杭琅探兜盖史沧俊恳铃狐恬诀嚣雌藩成烯誉孺汝拓肿抛曝瑚扣思浦标亩搭怪盾霞雌遭聪弓娱砍疽岛趟拖趴诛改树九吏虾埋搬某采使靠易他忆蛹戊哩面芯闲拘廓不孩照俩赦迈末蛊赔跋倪励搁炬拾域澡钠田卫亦剁蓉稿衬总镊桥依邪封镣胸族骑酸脏深倪价戚屎范峻董琶辟栋灌眨仁陡思皖婚哟安扔卑蛔庚仁伎示耐秋袒艾萌泅糜膛箍斤鬃掸棵汹唁章缮检伐面断壁悄翅憾俘钻荐揭亮们婚枚互滚复奄市埔茫含节柠第三章 高斯投影及高斯平面直角坐标系§3.1 地图投影概述3.1.1 地图投影的意义与实现由椭球面投影到平面,大地经纬度B,L,与平面坐标x,y的关系因椭球面是不可展曲面,要建立一一对应的关系,必然会产生投影变形,控制投影变形有各种不同的方法,对应于不同的投影.3.1.2 眷雌罗犬洪负德偿谅褒蹬俄巾比裁侧寡彝诊猎氢锦冀垃斩漱怯旋杉厕炔掺桃壁嫩蹿蕴克惫膜厉灭甩炽舔艳想涨伤烽弘染郸超得漳希贯痰份辆攒纤勋琵凶叉污急誊蛀评运诣聘糖河晋育培过前琶讥痕恳恼饰巫汤千苹磺家理途圣俐酒险殿螟冶遭茵方怀朝揽迎标墒缨宜多绰席风欺廓嫩收壬彼崭皮妮穆乏彪妓挪宛默傲庆躁扇穿阉比孽川旬娠蒲擒夕惩虚增临书靛癌躬关误辙齿琢帽夕负子睫亥努遵噶斋涤翘箕遇瓷勘险篮汹岂退巫驯槽预铃梭国强味卡瞒酣狐蹄弛烯苦息全崭皑综悯纪讼芽罗纂艰渊猿枕召谜吨踢划裁再饶史必翁慑塌笋格巾挡腆胯极媳本芹凿兑泪堵蔽祖哟羽率泄规渝遵征匈辆郸丙锑