高斯投影及高斯平面直角坐标.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,第三章 高斯投影及高斯平面直角坐标系,2,3.1,地图投影概述,3.1.1,地图投影的意义与实现,由椭球面投影到平面,大地经纬度,B,L,，,与平面坐标,x,y,的关系,因椭球面是不可展曲面,要建立一一对应的关系,必然会产生投影变形，控制投影变形有各种不同的方法，对应于不同的投影。,3.1.2,地图投影变形及其表述,1,、投影长度比、等量纬度及其表示式,3,长度比：,投影平面上微分长度与椭球面上相应微分长度之比。,投影平面上微分长度：,椭球面上微分长度：,3.1.2,地图投影变形及其表述,4,上式中,q,为等量纬度，计算公式为,引入等量纬度后，使相同角度量的,dq,与,dL,所对应的椭球面上的弧长相同。,3.1.2,地图投影变形及其表述,上式中,q,为等量纬度，计算公式为,引入等量纬度后，使相同角度量的,dq,与,dL,所对应的椭球面上的弧长相同。,3.1.2,地图投影变形及其表述,6,上式中,q,为等量纬度，计算公式为,引入等量纬度后，使相同角度量的,dq,与,dL,所对应的椭球面上的弧长相同。,3.1.2,地图投影变形及其表述,7,引入等量纬度后，投影公式为：,求微分，得：,其中：,l,=,L-L,0,3.1.2,地图投影变形及其表述,根据微分几何，其第一基本形式为：,8,其中：,3.1.2,地图投影变形及其表述,9,则，长度比公式为：,将 代入上式，得：,3.1.2,地图投影变形及其表述,10,当,A=0,或180，得经线方向长度比：,当,A=90,或270，得纬线方向长度比：,要使长度比与方向无关，只要：,F=0,E=G,，,则长度比可表示为：,3.1.2,地图投影变形及其表述,长度比与,1,之差，称为长度变形，即：,11,v,m,0,，,投影后长度变大，反之，投影后长度变短。,3.1.2,地图投影变形及其表述,2,、主方向和变形椭圆,12,主方向：,在椭球面上正交的两个方向投影到平面上后仍然正交，则这两个方向称为主方向。,性质：,主方向投影后具有最大和最小尺度比。,对照第一基本形式，得：,且：,3.1.2,地图投影变形及其表述,13,代入长度比公式，得：,若使：,使长度比为极值的方向：,由三角公式得：,3.1.2,地图投影变形及其表述,14,由此得，长度比极值为：,将三角展开式代入得：,因此，最大长度比,a,与最小长度比,b,可表示为：,3.1.2,地图投影变形及其表述,15,不难得出下列关系：,3.1.2,地图投影变形及其表述,若对应于最大和,最小长度,比,方向在椭球面上为,x,轴和,y,轴方向，在投影面上为,x,1,和,y,1,方向，则有：,16,椭球面上,投影面上,3.1.2,地图投影变形及其表述,17,3,、方向变形与角度变形,某方向（以主方向起始）,投影后为,1,，则有：,由三角公式，得：,显然，当,+,1,=90,或,270,时，方向变形最大,3.1.2,地图投影变形及其表述,18,若,与,1,表示最大变形方向，则最大变形量可表示为：,顾及：,解得最大变形方向为：,3.1.2,地图投影变形及其表述,19,两方向,、所夹角的变形称为角度变形，用表示。即：,显然，当,+,1,=90,、,+,1,=270,或,+,1,=270,、,+,1,=90,时，角度变形最大，最大角度变形可表示为：,3.1.2,地图投影变形及其表述,4,、面积比与面积变形,20,椭球面上单位圆面积为,，投影后的面积为,ab,，,则面积变形为：,3.1.3,地图投影的分类,1,、按投影变形的性质分类,(1).,等面积投影,a b=,1,(2).,等角投影,a=b,(3).,等距离投影,某一方向的长度比为,1,。,21,3.1.3,地图投影的分类,2,、按采用的投影面和投影方式分类,22,(1).,方位投影,投影面与椭球面相切，切点为投影中心，按一定条件将椭球面上的物投影到平面上。,3.1.3,地图投影的分类,(2).,正轴或斜、横轴圆柱投影,正轴圆柱投影：,投影圆柱面与某纬线相切,(,切圆柱投影,),、或相割,(,割圆柱投影,),切圆柱投影：,投影圆柱面与赤道相切，纬线投影成,一组平行直线，经线投影成与纬线正交,的另一组平行直线。,割圆柱投影：,投影圆柱面与两条对称纬线相割，纬线,投影成一组平行直线，经线投影成与纬,线正交的另一组平行直线。,23,3.1.3,地图投影的分类,横轴圆柱投影：,投影圆柱面与某经线相切。,斜轴圆柱投影：,用于小比例尺投影，将地球视为圆球，,投影圆柱体斜切于圆球进行投影。,(3).,圆锥投影：,圆锥面与椭球面相切或相割，将椭球面上,物投影到圆锥面上，展开圆锥面得投影平,面。,根据圆锥顶点位置不同，分正圆锥,投影、斜圆锥投影。,24,3.1.3,地图投影的分类,25,习 题,1.,给出等量纬度的定义，引入等量纬度有何作用。,2.,投影变形与长度无关时应满足哪些条件？并给出证明。,3.,变形主方向有什么性质？,4.,最大方向变形与最大角度变形的方向满足什么条件？,5.,地图投影按变形性质分哪几类？按投影方式分哪几类？,26,27,3.2,正形投影与高斯,-,克吕格投影,3.2.1,正形投影的概念和投影方程,长度比与方位角无关的投影称为正形投影，必须满足条件,E=G,F=0,，,即：,由第二式解得：,1,3.2.1,正形投影的概念和投影方程,代入第一式，得：,28,考虑到导数的方向，开方根得：,再代入 式，得：,1,2,3,3.2.1,正形投影的概念和投影方程,29,2,，式称为,Kauchi-Rimann,方程，满足该方程的,复变,函数,为解析函数，可展开成幂级数，即有,：,3,其反函数也是复变函数，可以写成：,3.2.2,高斯,-,克吕格投影的条件和性质,30,高斯,-,克吕格投影的条件：,1.,是正形投影,2.,中央子午线不变形,3.2.2,高斯,-,克吕格投影的条件和性质,高斯投影的性质：,1.,投影后角度不变,2.,长度比与点位有关，与方向无关,3.,离中央子午线越远变形越大,31,为控制投影后的长度变形，采用分带投影的方法。常用,3,度带或,6,度带分带，城市或工程控制网坐标可采用不按,3,度带中央子午线的任意带。,3.2.2,高斯,-,克吕格投影的条件和性质,32,3.2.2,高斯,-,克吕格投影的条件和性质,33,中央子午线在平面上的投影是,x,轴，赤道的投影是,y,轴，其交点是坐标原点。,x,坐标是点至赤道的垂直距离；,y,坐标是点至中央子午线的垂直距离,，,有正负。,为了避免,y,坐标出现负值，其名义坐标加上 500 公里。,为了区分不同投影带中的点，在点的,Y,坐标值上加带号,N,所以点的,横坐标的名义,值为,y=N,1000000+500000+y,3.3,高斯投影坐标正算和反算公式,34,3.2.1,高斯投影正算公式,赤 道,因正形投影的导数与方向无关，将投影点坐标在,H,点展开，得：,3.3.1,高斯投影正算公式,因此，高斯投影级数展开式可表示为：,35,其各阶导数为：,3.3.1,高斯投影正算公式,36,将导数代入展开式，虚实分开后，得到高斯投影正算公式如下：,3.3.1,高斯投影正算公式,为便于编程计算，可将正算公式改写成如下形式：,37,3.3.2,高斯投影反算公式,38,在中央子午线投影成的,x,轴上取点,X,f,=x,，该,点称为底点，用子午弧长反算公式求得底点的纬度,B,f,和相应的等量纬度,q,f,，,以,底点,为,展开,点,进行级数,展开,，得：,3.3.2,高斯投影反算公式,相应的各阶导数为：,39,3.3.2,高斯投影反算公式,代入级数展开式，虚实分开得：,40,4,3.3.2,高斯投影反算公式,41,将大地纬度展开成等量纬度的级数式,其中：,5,3.3.2,高斯投影反算公式,由 式，得：,42,4,代入 式，得：,5,3.3.2,高斯投影反算公式,将各系数代入上式，得纬度,B,的反算公式,：,43,3.3.2,高斯投影反算公式,为便于编程计算，可将反算公式改写成如下形式：,44,3.3.2,高斯投影反算公式,利用高斯投影的正反算公式，亦可进行不同投影带坐标的换带计算。其计算步骤如下：,1.,根据高斯投影坐标,x,y,，,反算得纬度,B,和经度差,l,；,2.,由中央子午线的经度,L,0,求得经度,L=L,0,+,l,；,3.,根据换带后新的中央子午线经度,L,0,，,计算相应的经差：,4.由高斯投影正算，求得新的高斯投影坐标,x,，,y,。,45,习 题,1.,高斯投影的条件是什么？,2.,简述高斯投影投影正算公式的推导；,3.,已知某点的坐标：,B=29,0405.3373,L,=121,1033.2012,计算：1).该点的3 带和6 带带号；,2).该点的3 带高斯投影坐标并反,算检核；,46,3.4,平面子午线收敛角和长度比,47,3.4.1,平面子午线收敛角的计算公式,平行圈,子午线,沿平行圈纬度不变，求微分得：,3.4.1,平面子午线收敛角的计算公式,对高斯投影公式求偏导数，得：,48,3.4.1,平面子午线收敛角的计算公式,49,代入上式，得：,将,展开成,tg,的级数，得：,3.4.1,平面子午线收敛角的计算公式,由此可见，,是经差的奇函数，在,x,轴为对称轴，东侧为正，西侧为负。,子午线收敛角在赤道为,0,，在两极等于,经差,l,，,其余点上均小于,经差,l,。,50,3.4.1,平面子午线收敛角的计算公式,子午线收敛角也可以表示成高斯平面坐标的级数展开式。,51,平行圈,L,=,常数,L+dl,=,常数,P,点沿,与,y,轴,平行方问微分,变,动,到,P,点，子午线收敛角可表示为：,沿,y,坐标的微分，得：,3.4.1,平面子午线收敛角的计算公式,52,代入子午线收敛角公式，得：,由高斯投影反算公式求出偏导数，得：,3.4.1,平面子午线收敛角的计算公式,代入上式子午线收敛角计算公式，得：,53,将,展开成,tg,的级数，得：,3.4.2,长度比计算公式,54,由高斯投影长度比的定义式，得：,将前面的偏导数代入上式，得：,开方后得出以大地坐标表示的长度比公式：,3.4.2,长度比计算公式,为给出由高斯投影坐标表示的长度比公式，反解高斯投影的,y,坐标正算公式，得：,55,对上式求平方和四次方，得：,3.4.2,长度比计算公式,56,代入用大地坐标表示的长度比公式，得：,顾及：,代入上式，得：,可见，长度比是,y,坐标的偶函数，且只与,y,坐标有关。,3.5,高斯投影距离与方向改化以及坐标方位角,57,3.5.1,高斯投影的距离改化,椭球面上的大地线投影到高斯平面上为曲线，与平面上两点相连的直线相比，其微分线段间的差异极小，可表示为：,其中：,3.5.1,高斯投影的距离改化,58,此弧线与直线间的最大偏角即为方向投影改化，本为二次小项,故此相对长度差异仅为,4,次项,相对于距离测量的最高精度亦可忽略,因此可认为：,用辛卜生公式数值积分得：,3.5.1,高斯投影的距离改化,将长度比公式,代入上式，得：,59,3.5.1,高斯投影的距离改化,距离改化,S,可表示为：,60,其中：,在城市及工程应用中测边离中央子午线不会超过,45,公里，则距离改化公式可进一步简化为：,3.5.2,高斯投影方向改化,61,1,、高斯投影曲线的形状,高斯投影曲线的形状向,x,轴弯曲，并向两极收敛。,3.5.2,高斯投影方向改化,2,、高斯投影方向改化,62,保角投影前后角度相同，即：,3.5.2,高斯投影方向改化,将球面角超计算公式代入上式，得：,63,因方向值顺时针方向增加，考虑其正负号后，方向改化公式可表示如下：,上式具有,0.1,的计算精度，适用于三、四等控制网的方向改化计算。改化公式中的曲率半径可足够近似地取,6370km,3.5.3,坐标方位角和大地方位角的关系式,64,A,12,T,12,习 题,1.,已知某点的坐标：,B=29,0405.3373,L,=121,1033.2012,计算：1).该点的3 带高斯投影后的中央子午,线收敛角；,2).该点的3 带高斯投影的长度比。,2.已知起始点坐标：,x,3,=3239387.624 m,y,3,=40446822.368m,起始平面方位角,T,31,=192,3708.51,，,距离,S,31,=7619.245m，,各方向观测值如下：,13,：,00000.00 23,：,00000.00 31,：,00000.00,12,：,2571747.71 21,：,395112.50 32,：,372636.65,将上述边长和方向归算到高斯平面上。,65,3,1,2,3.6,通用横轴墨卡托投影,3.6.1,墨卡托投影,66,墨卡托投影为等角割圆柱投影，圆柱与椭球面相割于,B,0,的两条纬线，投影后不变形。,特性：,等角航线在投影平面上为直线。因此，该投影便于在航海中应用。,3.6.2,通用横轴墨卡托投影,67,简称,为,UTM,，与高斯投影相比，仅仅是中央子午线的尺度比为,0.9996,，其投影公式如下：,3.6.2,通用横轴墨卡托投影,长度比和子午线收敛角计算公式。,68,3.6.2,通用横轴墨卡托投影,通用横轴墨卡托投影的反算步骤：,1.,先由通用横轴墨卡托投影坐标计算高斯投影坐标；,2.,再利用高斯投影反算公式，计算大地纬度和经度。,69,3.6.2,通用横轴墨卡托投影与高斯投影的比较,70,3.7,局部区域中的高斯投影及其相应的区域性椭球,局部区域中常采用地方独立坐标系，其,高斯坐标以往并非由经纬度求得,，而是直接将边长投影到边长归算的高程基准面（投影面）,再选定过测区中心附近的坐标纵轴，计算高斯投影边长和方向改正，在平面上由起始点坐标、起始方位角来平差计算各控制点坐标。,71,3.7,局部区域中的高斯投影及其相应的区域性椭球,地方独立坐标系的参数：,72,1.,投影面：一般采用区域的平均高程面；,2.,中央子午线的经度或位置：一般取用过区域中心附近一控制点的经度，或采用整分或整度的经度。,3.,起始坐标、起始方位角、起始边长。,3.7,局部区域中的高斯投影及相应的区域性椭球,城市及工程控制网采用地方独立坐标系，边长的投影面是区域的边长归算的高程基准面而并不是国家参考椭球面。其高斯坐标所对应的椭球面应是与投影面相接近的区域性椭球面，而不是国家参考椭球面。,73,习 题,1.,已知某点的坐标：,B=29,0405.3373,L,=121,1033.2012,计算：1).该点的3 带,UTM,投影坐标；,2).该点,UTM,投影的长度变形。,74,