2018年各幅图按椭球面法线投影的计算公式和结果.doc

各幅图按椭球面法线投影的计算公式和结果 　　/week114 　　　　各幅图按椭球面法线投影的计算公式和结果 　　　　施一民1,2, 陈月梅1,孙伟伟1 　　(1.同济大学 测量与国土信息工程系,上海 200092;2.现代工程测量国家测绘局重点实验室 　　上海,200092; 　　摘 要:为更大幅度地减小地图投影变形,提出了按椭球面法线投影的一种新的地图投影方法。它与 　　楔形投影一样也是以规则的经纬网格为单元，以每 4 个网格角点所构成的平面来切割并逼近椭球面 　　[1] ，有所不同的是,椭球面上各点投影到相应的切割平面是按照椭球面的法线方向。藉助于严密的几 何关系,推导出由椭球面上大地经纬度到投影平面上平面直角坐标的投影公式, 由数据验证充分说 　　明,由此所产生的各类投影变形都非常小。 　　关键词: 地图投影;投影变形;切割平面;椭球面法线 　　　　1 引言 　　　　为了减小长度投影变形,作者在文献[1]中提出了一种新的地图投影方法: 用多个连续吻 合的平面图幅(等腰梯形)来切割椭球面, 将各投影区域分别限定在单幅图内,,即按某一比 例尺地图图幅规定的经差和纬差构成规则的经纬网格, 分别对每一网格按某种纯几何方式 [1] 来进行投影。文献[2] 采用了四种投影方式 中的一种: 楔形投影,亦即椭球面上同纬度的点 　　子均以该纬线所在平面与椭球短轴交点为投影中心来投影到所在的切割平面。与现有的地图 投影方法相比,楔形投影的优点是能使纬线族投影为平行直线, 经线族投影为向两极收敛的 斜直线,且减小了最大的长度投影变形,但它并非是等角投影,角度投影变形还是略大些。 　　投影范围既已限制在单幅图内,能否由此寻找出一种能使各类投影变形均很小的多面体 　　[1] 投影,为此,本文再探索按各点的椭球面法线方向的地图投影 。 　　2 对按椭球面法线方向投影的初探 　　　　将椭球面上各点沿着其法线方向投影到相应的切割平面上, 它近似于中心投影，因对圆 球面而言，这正是以球心为投影中心的严格的中心投影。可证每个网格的东西边界经线投影 为等腰梯形的两腰，一族经线映射为该平面图幅中的一族向北收敛(在北半球)的斜直线, 并对称于梯形上下底中点的连线;但纬线族的投影则是接近于平行直线的向下凸的曲线。 　　作为一种纯几何的投影方式,在此并不加上任何限制某类投影变形的条件,只是基于与 此投影方式相对应的几何关系来推求由椭球面到各投影平面的投影公式。 　　　　3 按椭球面法线的投影公式 　　　　3.1 投影平面方程和法向量的表示式 　　　　若在椭球面上用规则的经纬格网覆盖整个制图区域，例如对应于 1:10000 图分幅，网 格的经纬度间隔 , B , , L 分别为 3′45″和 2′30″，文献[1]已证得 4 个经纬网格角点必 　　　　定共面,且构成等腰梯形,将其作为各网格相应的投影平面。4 个网格角点作为该平面图幅的 4 个图廓点,由其已知的大地经纬度就能求得该投影平面的法线向量。设西南图廓点 (WS) 的 ???????????????????????????? 　　基金项目:国家自然科学基金资助项目 ,地球空间环境与大地测量教育部重点实验室开放基金 资助项目(03-04-02) 　　作者简介:施一民(1942- ), 男 ,浙江宁波人, 教授，博士生导师,E-mail: yimshi@; 　　), 女 ,江苏镇江人, 硕士生。E-mail:chym204@; 陈月梅(1982- 　　1 　　/week114 　　　　大地经纬度为 ，BS , LW ， ,其相应的三维直角坐标为 　　　　X WS , N S cos BS cos LW 　　(1) YWS , N S cos BS sin LW 　　2 WSZ N , sin ， SS1eB 　　　　该投影平面的法线向量即为平面上两个相交的向量(例如西图廓线和南图廓线向量)的 矢量积,在此,这三个向量都表示在三维大地直角坐标系之中。于是即可推导出投影平面法线 向量的三个分量为 　　V X , ( N N sin B N N S sin BS )(1 e 2 )(sin LW sin LE ) N S cos BS 　　　　(2) VY , ( N N sin B N N S sin BS )(1 e 2 )(sin LE sin LW ) N S cos BS 　　　　VZ , ( N N sin BN N S sin BS ) sin，LE LW ，N S cos BS 　　　　式(2)中, BS , BN 分别为该经纬网格南北边纬线的