椭球面元素归算至高斯平面高斯投影幻灯片.ppt
《椭球面元素归算至高斯平面高斯投影幻灯片.ppt》由会员分享,可在线阅读,更多相关《椭球面元素归算至高斯平面高斯投影幻灯片.ppt(92页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、椭球面元素归算至高斯平面高斯投影第1页,共92页,编辑于2022年,星期一一、长度比一、长度比5.1 5.1 高斯投影概述高斯投影概述或者或者长度比不仅随点的位置,而且随线段的方向而发生变化。长度比不仅随点的位置,而且随线段的方向而发生变化。第2页,共92页,编辑于2022年,星期一n高斯投影是高斯投影是等角横轴切椭圆柱投影。等角横轴切椭圆柱投影。n高斯投影是一种等角投影。它是由德国数学家高斯高斯投影是一种等角投影。它是由德国数学家高斯(Gauss(Gauss,1777 1777 1855)1855)提出,后经德国大地测量学家克吕格提出,后经德国大地测量学家克吕格(Kruger(Kruger,
2、185718571923)1923)加以补充完善,故又称加以补充完善,故又称“高斯高斯克吕格投影克吕格投影”,简称,简称“高斯投影高斯投影”。二、高斯投影的基本概念二、高斯投影的基本概念第3页,共92页,编辑于2022年,星期一NSc中中中中央央央央子子子子午午午午线线线线赤道赤道高斯投影平面高斯投影平面赤道赤道中中央央子子午午线线 高斯投影高斯投影采用分带投影。将椭球面按一定经差分带,分采用分带投影。将椭球面按一定经差分带,分别进行投影。别进行投影。1、高斯投影的原理、高斯投影的原理第4页,共92页,编辑于2022年,星期一(1)(1)中央子午线投影后为直线,且长度不变。中央子午线投影后为直
3、线,且长度不变。(2)(2)除中央子午线外,其余子午线的投影均为凹除中央子午线外,其余子午线的投影均为凹向中央子午线的曲线,并以中央子午线为对向中央子午线的曲线,并以中央子午线为对称轴。投影后有长度变形。称轴。投影后有长度变形。(3)(3)赤道线投影后为直线,有长度变形赤道线投影后为直线,有长度变形(4)(4)除赤道外的其余纬线,投影后为凸向赤道除赤道外的其余纬线,投影后为凸向赤道的曲线,并以赤道为对称轴。的曲线,并以赤道为对称轴。(5)(5)经线与纬线投影后仍然保持正交。经线与纬线投影后仍然保持正交。(6)(6)所有长度变形的线段,其长度变形比均大于所有长度变形的线段,其长度变形比均大于l
4、l。(7)(7)离中央子午线愈远,长度变形愈大。离中央子午线愈远,长度变形愈大。赤道赤道中央子午线中央子午线平行圈平行圈子午线子午线O Ox xy y2、高斯投影的特点、高斯投影的特点第5页,共92页,编辑于2022年,星期一 我国规定按经差我国规定按经差66和和33进进行投影分带。行投影分带。6 6带自首子午线开始,按带自首子午线开始,按66的经差自西向东分成的经差自西向东分成6060个带。个带。3 3带自带自1.5 1.5 开始,按开始,按33的的经差自西向东分成经差自西向东分成120120个带。个带。高斯投影带划分高斯投影带划分3、投影带的划分、投影带的划分第6页,共92页,编辑于202
5、2年,星期一 66带与带与33带中央子午线之间的关系如图带中央子午线之间的关系如图:33带的中央子午线与带的中央子午线与66带中央子午线及分带子午线重带中央子午线及分带子午线重合,减少了换带计算。合,减少了换带计算。工程测量采用工程测量采用3 3 带带特殊工程可采用特殊工程可采用1.5 1.5 带或任意带带或任意带第7页,共92页,编辑于2022年,星期一 按照按照66带划分的规定,第带划分的规定,第1 1带中央子午线的经度为带中央子午线的经度为33,其余各带中央子午线经度与带号的关系是:,其余各带中央子午线经度与带号的关系是:L L。=6N=6N3 3(N N为为66带的带号)带的带号)例:
6、例:2020带中央子午线的经度为:带中央子午线的经度为:L L。6 206 2033117 117 按照按照33带划分的规定,第带划分的规定,第1 1带中央子午线的经度为带中央子午线的经度为33,其余各带中央子午线经度与带号的关系是:,其余各带中央子午线经度与带号的关系是:L L。=3n =3n (n n为为33带的带号)带的带号)例:例:120120带中央子午线的经度为带中央子午线的经度为 L L。3 1203 120360 360 第8页,共92页,编辑于2022年,星期一 若已知某点的经度为若已知某点的经度为L L,则该点的,则该点的66带的带号带的带号N N由下由下式计算:式计算:若已
7、知某点的经度为若已知某点的经度为L L,则该点所在,则该点所在33带的带号带的带号按下式计算:按下式计算:(四舍五入)(四舍五入)第9页,共92页,编辑于2022年,星期一x x轴轴 中央子午线的投影中央子午线的投影y y轴轴 赤道的投影赤道的投影原点原点O O 两轴的交点两轴的交点OxyP(X,Y)高斯自然高斯自然坐标坐标注:注:X X轴向北为正,轴向北为正,y y轴向东为正。轴向东为正。赤道赤道中央子午线中央子午线4、高斯平面直角坐标系的建立、高斯平面直角坐标系的建立第10页,共92页,编辑于2022年,星期一由于我国的位于北半球,由于我国的位于北半球,东西横跨东西横跨1212个个66带,
8、各带带,各带又独自构成直角坐标系。又独自构成直角坐标系。故:故:X X值均为正,值均为正,而而Y Y值则有正有负。值则有正有负。u为了免出现负的横坐标,在横坐标上加上为了免出现负的横坐标,在横坐标上加上500 000m500 000m。此外还在坐标前面再冠以。此外还在坐标前面再冠以带号。这种坐标称为国家统一坐标。带号。这种坐标称为国家统一坐标。例如:例如:有一点有一点Y=19 123 456.789mY=19 123 456.789m,该点位在,该点位在1919带内,其相对于中央子午线的坐带内,其相对于中央子午线的坐标标y=376 543.211m。u为了把各带联成整体,一般规定各投影带要有一
9、定的重叠度,其中每一为了把各带联成整体,一般规定各投影带要有一定的重叠度,其中每一66带向东加宽带向东加宽3030,向西加宽,向西加宽1515或或7.5.7.5.这样在带边缘,控制点将有两套相邻带的坐标值,地形图将有两套公里格网。这样在带边缘,控制点将有两套相邻带的坐标值,地形图将有两套公里格网。保证了控制点间的互相应用,地图的顺利拼接和使用。保证了控制点间的互相应用,地图的顺利拼接和使用。第11页,共92页,编辑于2022年,星期一例:例:有一国家控制点的坐标有一国家控制点的坐标:x=3102467.280m,y=19367622 x=3102467.280m,y=19367622380m3
10、80m,(1 1)该点位于)该点位于6 6 带的第几带?带的第几带?(2 2)该带中央子午线经度是多少?)该带中央子午线经度是多少?(3 3)该点在中央子午线的哪一侧?该点距中央子午线和)该点在中央子午线的哪一侧?该点距中央子午线和赤道的距离为多少?赤道的距离为多少?(1)第)第19带带(2)L。=619-19-3=111(3)y367622.380-500000-132377.620m,在西侧),在西侧)(距中央子午线(距中央子午线132377.620m,距赤道,距赤道3102467.280m)第12页,共92页,编辑于2022年,星期一三、椭球面三角系化算到高斯平面三、椭球面三角系化算到高
11、斯平面椭球面内(椭球面内(中央子午线中央子午线ON,赤道,赤道OE)三角网三角网PKTMQ:起始点起始点P大地坐标(大地坐标(B,l),),l=L-L0,L、L0分别为分别为P和轴子午线的大地经度;和轴子午线的大地经度;起始边起始边PK=S;起始边的大地方位角起始边的大地方位角APK;PC为垂直于中央子午线的大地线,为垂直于中央子午线的大地线,C点大地坐标(点大地坐标(B0,l=0););PP1为过为过P点平行圈,点平行圈,P1点的大地坐标(点的大地坐标(B,l=0););X为赤道至纬度为赤道至纬度B的平行圈子午弧长。的平行圈子午弧长。第13页,共92页,编辑于2022年,星期一三、椭球面三角
12、系化算到高斯平面三、椭球面三角系化算到高斯平面高斯投影面上高斯投影面上:中央子午线和赤道分别为直线中央子午线和赤道分别为直线ON 及及OE ,其他子午线和平行圈均变为曲线。其他子午线和平行圈均变为曲线。PN是是PN的投影,的投影,P P 1是是PP1的投影;的投影;P的直角坐标为(的直角坐标为(x,y););因是等角投影,大地方位角因是等角投影,大地方位角APK投影后没有变化。投影后没有变化。三角形投影后变为边长三角形投影后变为边长si的的曲线曲线三角形三角形(长度大于椭球面上的边长长度大于椭球面上的边长),且,且曲线都凹向纵坐标轴曲线都凹向纵坐标轴;第14页,共92页,编辑于2022年,星期
13、一1、椭球面三角系化算到高斯投影、椭球面三角系化算到高斯投影面问题分析面问题分析(1)投影后需用连接各点间的弦线来代替曲线。为此,必须在每个方)投影后需用连接各点间的弦线来代替曲线。为此,必须在每个方向上引进曲改直的水平方向改正;向上引进曲改直的水平方向改正;(2)根据始点)根据始点P的大地坐标的大地坐标B,L计算其平面坐标的坐标正算公式;计算其平面坐标的坐标正算公式;(3)反算公式;)反算公式;第15页,共92页,编辑于2022年,星期一1、椭球面三角系化算到高斯投影面、椭球面三角系化算到高斯投影面问题分析问题分析(4)确定平面三角形各边坐标方位角)确定平面三角形各边坐标方位角a。(5)确定
14、平面三角形各边长。)确定平面三角形各边长。第16页,共92页,编辑于2022年,星期一(1)高斯投影坐标计算高斯投影坐标计算 将起始点的大地坐标将起始点的大地坐标B,L归算为高斯平面直角坐标归算为高斯平面直角坐标x,y;根据(;根据(x,y)反算(反算(B,L)。)。(2)通过计算该点的子午线收敛角及方向改正,将通过计算该点的子午线收敛角及方向改正,将椭球面上起算边大地椭球面上起算边大地方位角归算到高斯平面上方位角归算到高斯平面上相应边的相应边的坐标方位角坐标方位角。(3)通过计算各方向的曲率改正和方向改正,将通过计算各方向的曲率改正和方向改正,将椭球面上各三角形内角椭球面上各三角形内角归算到
15、高斯平面上归算到高斯平面上的由相应的由相应直线组成的三角形内角直线组成的三角形内角。(4)通过计算距离改正,将通过计算距离改正,将椭球面上起算边的长度归算到高斯平面上椭球面上起算边的长度归算到高斯平面上的直线长度的直线长度。(5)控制网跨越两投影带时,需要进行控制网跨越两投影带时,需要进行平面坐标的邻带换算平面坐标的邻带换算。2、将椭球面三角系化算到高斯投影面、将椭球面三角系化算到高斯投影面的主要内容的主要内容n将椭球面三角系归算到平面上,包括坐标、曲率改正、距离改正和子午将椭球面三角系归算到平面上,包括坐标、曲率改正、距离改正和子午线收敛角等项计算工作。线收敛角等项计算工作。第17页,共92
16、页,编辑于2022年,星期一第一类称高斯投影正算公式,亦即由(第一类称高斯投影正算公式,亦即由(B,L)求()求(x、y);第二类称高斯投影反算公式,亦即由(第二类称高斯投影反算公式,亦即由(x、y)求()求(B,L)。)。5.2 5.2 高斯投影坐标正反算高斯投影坐标正反算一、高斯投影坐标正算公式一、高斯投影坐标正算公式高斯投影必须满足以下三个条件:高斯投影必须满足以下三个条件:中央子午线投影后为直线;中央子午线投影后为直线;中央子午线投影后长度不变;中央子午线投影后长度不变;投影具有正形性质,即正形投影条件。投影具有正形性质,即正形投影条件。第18页,共92页,编辑于2022年,星期一(1
17、)在在经经差小于差小于3.5时时,精度,精度为为0.1m第19页,共92页,编辑于2022年,星期一(2)扩展到精度展到精度0.001m的公式:的公式:自赤道量起的到所求点的子午线弧长自赤道量起的到所求点的子午线弧长所求点的大地经度与该所求点的大地经度与该点所在带的中央子午线点所在带的中央子午线的大地经度之差的大地经度之差第20页,共92页,编辑于2022年,星期一二、高斯投影坐标反算公式二、高斯投影坐标反算公式投影方程投影方程x坐标轴投影后为中央子午线坐标轴投影后为中央子午线,是投影的对称轴;是投影的对称轴;x坐标轴投影后长度不变;坐标轴投影后长度不变;投影具有正形性质。即高斯面上的角度投影
18、到椭球面上投影具有正形性质。即高斯面上的角度投影到椭球面上后角度没有变形,仍然相等。后角度没有变形,仍然相等。投影函数投影函数1 1、2 2应应满足以下三个条件:满足以下三个条件:首先根据首先根据x计算纵坐标在椭球面上的垂足纬度计算纵坐标在椭球面上的垂足纬度Bf,接着按,接着按Bf计计算算(Bf-B)及经差及经差l,最后得到最后得到反算公式的推导方法的基本思想:反算公式的推导方法的基本思想:1.1.高斯投影坐标反算基本思想高斯投影坐标反算基本思想第21页,共92页,编辑于2022年,星期一2.2.精度为精度为0.00010.0001的高斯投影坐标反算公式的高斯投影坐标反算公式垂足纬度。垂足纬度
19、。其值由子午线其值由子午线弧长计算公式反弧长计算公式反算求得算求得第22页,共92页,编辑于2022年,星期一1.1.正算实用公式正算实用公式5.3 5.3 高斯投影坐标计算的实用公式高斯投影坐标计算的实用公式第23页,共92页,编辑于2022年,星期一第24页,共92页,编辑于2022年,星期一2.2.反算实用公式反算实用公式第25页,共92页,编辑于2022年,星期一第26页,共92页,编辑于2022年,星期一5.4 5.4 方向改化公式方向改化公式u方向改正数方向改正数就是指大地线的就是指大地线的投影曲线和连接大地线两点投影曲线和连接大地线两点的弦之夹角。的弦之夹角。第27页,共92页,
20、编辑于2022年,星期一设地球椭球为一圆球;设地球椭球为一圆球;OD为轴子午线;为轴子午线;AB为一条大地线为一条大地线(是球面上一条是球面上一条大圆弧大圆弧),投影为曲线,投影为曲线ab;AD、B E是与轴子午线正交点大圆是与轴子午线正交点大圆弧,投影分别为垂直于弧,投影分别为垂直于x轴的直线轴的直线ad和和be。一、方向改化近似公式的推导一、方向改化近似公式的推导第28页,共92页,编辑于2022年,星期一方向改化概略数值方向改化概略数值误差小于误差小于0.1,可适用于三、四,可适用于三、四等三角测量的计算等三角测量的计算由表可见,对于各等三角测量计算,方向改正都不能忽略。由表可见,对于各
21、等三角测量计算,方向改正都不能忽略。第29页,共92页,编辑于2022年,星期一1)用勒让德尔定理解算球面三角形用勒让德尔定理解算球面三角形 1 1、准备知识、准备知识二、方向改化较精密公式的推导二、方向改化较精密公式的推导 勒让德尔定理:勒让德尔定理:如果平面三角形和球面三角形对应边相等,则平面如果平面三角形和球面三角形对应边相等,则平面角等于对应球面角减去三分之一球面角超。角等于对应球面角减去三分之一球面角超。第30页,共92页,编辑于2022年,星期一P为平面三角形的面积。为平面三角形的面积。可直接用球面角代替平面角计算球面角超,虽然带有误差,但研究可直接用球面角代替平面角计算球面角超,
22、虽然带有误差,但研究表明:当边长不大于表明:当边长不大于90km时,这种误差小于时,这种误差小于0.0005,可忽略。,可忽略。2)球面角超的计算球面角超的计算第31页,共92页,编辑于2022年,星期一设地球椭球为一圆球;设地球椭球为一圆球;AB为轴子午线为轴子午线;小圆弧小圆弧P1Q与轴子午线平行,垂直于与轴子午线平行,垂直于BQ、AP1,投影为,投影为PQ;大圆弧大圆弧P1CQ的投影为曲线的投影为曲线P1 C Q。2 2、方向改化较精密公式的推导、方向改化较精密公式的推导第32页,共92页,编辑于2022年,星期一勒让德尔定理勒让德尔定理第33页,共92页,编辑于2022年,星期一代入代
23、入我国二等三角网平均边长为我国二等三角网平均边长为13KM13KM,当,当y ym m250km250km时,上时,上式精确至式精确至0.010.01,故通常用,故通常用于于二等二等三角测量计算。三角测量计算。整理,以Rm代替R第34页,共92页,编辑于2022年,星期一该式精确至该式精确至0.0010.001,适用于一等三角测量计算。,适用于一等三角测量计算。第35页,共92页,编辑于2022年,星期一5.5 5.5 距离改化公式距离改化公式由由S化至化至D所加的所加的S改正称为改正称为距离改正距离改正.1)研究平面曲线长度研究平面曲线长度s与其弦线长度与其弦线长度D的关系的关系;2)研究用
24、大地坐标研究用大地坐标B、L和平面坐标和平面坐标x、y计算长度比计算长度比m的公式的公式;3)导出距离改化的计算公式。导出距离改化的计算公式。m1S:大地线长;:大地线长;s:大地线:大地线S在高斯平面的投影;在高斯平面的投影;D:直线长;:直线长;思路:第36页,共92页,编辑于2022年,星期一1 1、平面曲线长度、平面曲线长度s s与其弦线长度与其弦线长度D D的关系的关系 v是一个小角,最大不会超过方向改化值是一个小角,最大不会超过方向改化值,因此可,因此可把把cosv展开为级数:展开为级数:式中用式中用v的最大值的最大值代替代替 v第37页,共92页,编辑于2022年,星期一 目前,
25、最高的距离测量精度约为目前,最高的距离测量精度约为10-8,弧线与直线的长度差异,弧线与直线的长度差异完全可以忽略完全可以忽略.完全可以认为:完全可以认为:大地线的平面投影曲线长度大地线的平面投影曲线长度s等于其弦线长度等于其弦线长度D。第38页,共92页,编辑于2022年,星期一2 2、长度比和长度变形、长度比和长度变形 长度比长度比m是指椭球面上某一点的微分元素是指椭球面上某一点的微分元素dS,与其投影面,与其投影面上的相应的微分元素上的相应的微分元素ds之比,即:之比,即:为为长度变形长度变形。1)用大地坐标表示的长度比公式)用大地坐标表示的长度比公式第39页,共92页,编辑于2022年
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 椭球 元素 至高 平面 投影 幻灯片
限制150内