《空间数据坐标变换.pptx》由会员分享,可在线阅读,更多相关《空间数据坐标变换.pptx(63页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、地球与地理参数 地球是一个近似球体,其自然表面是一个极其复杂的不规则曲面为描述和表达地球表面,必须选择一个与地球形状、大小相接近的椭球体来近似代替它地球椭球的两个主要参数为长轴半径a,短轴半径b,以及三个派生参数:扁率、第一偏心率和第二偏心率 第1页/共63页梨形球体WEPNPSAabO地球椭球体第2页/共63页扁率:第一偏心率:第二偏心率:第3页/共63页常用的地球椭球体数据椭球体名称提出年代长 轴 半 径(m m)短 轴 半 径(m m)扁率EverestEverest1830183063772766377276635607563560751 1:300.8300.8BasselBasse
2、l1841184163373796337379635607963560791 1:299.15299.15ClarkeClarke1880188063782496378249635651563565151 1:293.5293.5ClarkeClarke1886188663782066378206635658463565841 1:295.0295.0HayfordHayford1910191063783886378388635691263569121 1:297297克拉索夫斯基 1940194063782456378245635686363568631 1:298.3298.3IUGGIU
3、GG1976197663781606378160635677563567751 1:298.25298.25WGS84WGS841984198463781376378137635675263567521 1:298.26298.26第4页/共63页不同的地方,由于变形规律的特点,因此根据自己国家的具体位置和采用的投影选择适合自己的椭球体。例如我国。第5页/共63页中国使用的地球椭球体 1)海福特椭球(1910)我国52年以前基准椭球 a=6378388m b=6356911.9461279m =0.336700336702)克拉索夫斯基椭球(1940 Krassovsky)北京54坐标系基准椭
4、球 a=6378245m b=6356863.018773m =0.33523298692第6页/共63页中国使用的地球椭球体3)1975年I.U.G.G推荐椭球(国际大地测量协会1975)西安80坐标系基准椭球 a=6378140m b=6356755.2881575m =0.00335281317784)WGS-84椭球(GPS全球定位系统椭球、17届国际大地测量协会)WGS-84 GPS 基准椭球 a=6378137m b=6356752.3142451m =0.00335281006247第7页/共63页地图投影建立地球表面上点与投影平面上点之间的一一对应关系 通过将不可展的球面投影到
5、一个可展曲面上,然后将该曲面展开成为一个平面实质是建立地球椭球面上的地理坐标(经纬度)和平面上直角坐标之间的函数关系第8页/共63页地图投影地球表面上一点:地理坐标、(经纬度)地图采用直角坐标x、y 第9页/共63页地图投影的基本分类 基于投影面与球面相互位置的分类 通常地图投影过程所采用的可展曲面有圆锥面、圆柱面和平面,对应的地图投影有圆锥投影、圆柱投影和方位投影 第10页/共63页方位投影:以平面作为投影面,使平面与球面相切或相割,将球面上的经纬线投影到平面上而成。圆柱投影:以圆柱面作为投影面,使圆柱面与球面相切或相割,将球面上的经纬线投影到圆柱面上,然后将圆柱面展为平面而成。圆锥投影:以
6、圆锥面作为投影面,使圆锥面与球面相切或相割,将球面上的经纬线投影到圆锥面上,然后将圆锥面展为平面而成。第11页/共63页投影面与地球的相对位置:正轴投影(投影面的中心轴与地轴重合)斜轴投影(投影面的中心轴与地轴斜交)横轴投影(投影面的中心轴与地轴垂直)第12页/共63页投影面与地球椭球体的相交位置切投影(投影面与椭球体相切)割投影(投影面与椭球体相割)第13页/共63页第14页/共63页第15页/共63页基于投影方法的分类 透视-几何投影:依据透视原理,根据视点、物点与像点之间的几何关系,来建立投影方程 几何-解析投影:首先根据经纬线形状确定投影方程的基本形式,然后再依据某种约束条件解析地推求
7、出特定投影的基本方程 解析投影:依据人们给出的约束条件逐步推求经纬线的形状与投影方程 第16页/共63页基于投影变形的分类等角投影:投影后任意点处由任意两条微分线段构成的角度不产生变形。该种投影可以使得投影前后的形状保持不变,。等面积投影:无论是微分面积还是区域面积在投影前后均保持不变。任意投影:既不保持角度不变,也不保持面积不变,可能还存在长度变形。但当投影方向选择合适时,可以保证某一方向上的长度不变。第17页/共63页长度变形、面积变形和角度变形的关系1 1 1 1)等积投影上不能保持等角特性,)等积投影上不能保持等角特性,2 2 2 2)等角投影上不能保持等积特性;)等角投影上不能保持等
8、积特性;3 3 3 3)任意投影上不能保持等角和等积的特性;)任意投影上不能保持等角和等积的特性;4 4 4 4)等积投影的形状变形比较大;)等积投影的形状变形比较大;5 5 5 5)等角投影的面积变形比较大;)等角投影的面积变形比较大;6 6 6 6)等距投影的面积变形小于等角投影,角度)等距投影的面积变形小于等角投影,角度变变 形小于等积投影。形小于等积投影。第18页/共63页地图投影变形 地球上同纬度带经差相同的网格必具有相同的大小和形状。但是它们在地球上同纬度带经差相同的网格必具有相同的大小和形状。但是它们在投影中不一定能保持原来的大小和形状,甚至彼此间有很明显的差异投影中不一定能保持
9、原来的大小和形状,甚至彼此间有很明显的差异.第19页/共63页7.2.1 地图投影变形种类 可以发现变形表现在长度、面积和角度三个方面。分别用长度比、面积比的变化显示投影中长度变形和面积变形。如果长度变形或面积变形为零,则没有长度变形或没有面积变形。角度变形即某一角度投影后角值与它在地球表面上固有角值之差。第20页/共63页1)长度变形地图上的经纬线长度与地球仪上的经纬线长度特点并不完全相同,地图上的经纬线长度并非都是按照统一比例缩小的。地球仪上经纬线的长度具有下列特点:纬线长度不等,所有的经线长度都相等同一条纬线上,经差相同的纬线弧长相等 第21页/共63页2)面积变形地图上经纬线网格面积与
10、地球仪经纬线网格面积的特点不同,在地图上经纬线网格面积不是按照统一比例缩小的。地球仪上经纬线网格的面积具有下列特点:同一纬度带内,经差相同的网格面积相等;同一经度带内,纬线越高,网格面积越小。第22页/共63页3)角度变形 指地图上两条所夹的角度不等于球面上相应的角度。地球仪上经线和纬线处处都呈直角相交。第23页/共63页第24页/共63页 变形椭圆变形椭圆 变形椭圆变形椭圆 :论述和显示投影变形。论述和显示投影变形。变形椭圆:地面一点上的一个无穷小圆微分圆(也称单位圆),在投影后一般地形成的一个微分椭圆。利用这个微分椭圆能较恰当地、直观地显示变形的特征。第25页/共63页 实地上半径为单位值
11、(r1)的微分圆,在不同投影中具有不同的形状和大小。第26页/共63页 不同位置的变形椭圆形状差异很大,但面积大小差不多等积投影第27页/共63页 变形椭圆保持为圆形,但在不同位置上面积差异很大等角投影第28页/共63页椭圆的形状与大小都有着不同的变化 任意投影 第29页/共63页GISGIS中地图投影选择 数据获取(数据源地图的投影)数据标准化预处理(按某一参照系数字化)数据存储(统一的坐标基础)数据输出(具有相应投影的地图)数据应用(检索查询、覆盖分析)数据处理(投影转换)地图基础(地图投影)第30页/共63页地图投影选择的一般原则 GIS所采用的投影系统应与本国的基本地图系统所采用的投影
12、系统一致;各比例尺GIS中的投影系统应与相应比例尺主要信息源地图的投影一致;各地区的GIS投影系统应与该地区所使用的投影系统一致;一般选择1-3种投影系统。第31页/共63页中国的地图投影选择 1 1:5050万、1 1:2525万、1 1:1010万、1 1:5 5万、1 1:2.52.5万、1 1:1 1万、1 1:50005000采用高斯克吕格投影;1 1:100100万地形图采用LambertLambert投影,大部分省区图以及大多数这一比例尺的地图也多采用LambertLambert投影或属于同一投影系统的AlbersAlbers投影;第32页/共63页 德国数学家、物理学家、天文学
13、家高斯于19世纪20年代拟定,后经德国大地测量学家克吕格于1912年对投影公式加以补充,故称为高斯克吕格投影。7 7种国家基本比例尺地形图中,规定1 1:1 1万、1 1:2.52.5万、1 1:5 5万、1 1:1010万、1 1:2525万、1 1:5050万比例尺地形图,均采用高斯克吕格投影。高斯克吕格投影 第33页/共63页高斯克吕格投影将椭圆柱横切于地球椭球体上,其切线为某一经线(称为中央经线),然后根据一定的约束条件(即投影条件)将中央经线两侧规定范围内的点投影到椭球面上,从而得到各点的高斯投影。第34页/共63页高斯克吕格投影图示第35页/共63页投影条件 中央经线和地球赤道线投
14、影成直线且作为该投影带的对称轴;等角投影;中央经线上没有长度变形。第36页/共63页投影特点 没有角度变形,沿任意方向的长度比相等;中央经线上无长度变形;同一纬线上,离中央经线越远,变形越大;同一经线上,纬度越低,变形越大;等变形线为平行于中央经线的直线;最大变形位于各投影带的赤道边缘处。第37页/共63页1)6度带是从0o子午线起,自西向东每隔经差6为一投影带,全球分为60带,带号用自然序数1,2,3,60表示。2)3度带,是从东经1度30分的经线开始,每隔3度为一带,全球划分为120个投影带。高斯克吕格投影分带第38页/共63页1 1)高斯克吕格投影分带图示第39页/共63页我国 1:2.
15、5万1:50万地形图采用6分带方案,我国领土位于72-136之间,共覆盖11个投影带(13带-23带)。1:1万及更大比例尺地形图采用3分带方案,全球共分为120个投影带,我国领土共覆盖23个投影带(24-46带)。第40页/共63页LambertLambert投影 用一标准圆锥面相割地套在地球球体上,按某种规则将地球表面上的要素投影到圆锥面上;然后将圆锥面沿某一母线(通常选择某一经线)展开,来得到地球表面的平面投影结果。第41页/共63页LambertLambert投影 PNPSAA第42页/共63页LambertLambert投影特点:投影后经线交于圆锥顶点;纬线表现为同心圆,圆心为圆锥顶
16、点。第43页/共63页坐标系经线和纬线是地球表面上两组正交的曲线,由其构成的坐标系称为地理坐标系坐标系是一组参数,它告诉你如何判读对象的定位坐标,其中一个参数就是投影 第44页/共63页1)X坐标值在赤道以北为正,以南为负;Y坐标值在中央经线以东为正,以西为负。我国在北半球,X坐标皆为正值。Y坐标在中央经线以西为负值,为此将各带的坐标纵轴西移500公里,即将所有Y值都加500公里。2)由于采用了分带方法,各带的投影完全相同,某一坐标值(x,y),在每一投影带中均有一个,在全球则有60个同样的坐标值。因此,在Y值前,需冠以带号,这样的坐标称为通用坐标。高斯克吕格投影坐标 第45页/共63页第46
17、页/共63页高斯克吕格投影计算 s是由赤道到纬度的经线弧长,可根据上面公式估算,也可查高斯克吕格投影计算表;e为第二偏心率,e为第一偏心率;为某经线与中央经线的经度之差(弧度表示);a为地球长半轴。第47页/共63页坐标变换 由于数字化仪桌面坐标系与地图直角坐标系的不一致以及图纸变形等原因,由数字化仪读出的坐标与地图直角坐标系中的坐标不可能一致,为获得一致的数据,必须进行坐标转换以使坐标系统一。第48页/共63页相似变换 假设地图的实际比例尺在X、Y方向是一致的。YXyxOoA0B0第49页/共63页若令 第50页/共63页仿射变换 图纸数字化时两方向的比例尺存在差异,假设分别为mx、my。第
18、51页/共63页投影变换 不同来源的地图还可能存在地图投影与地理坐标的差异,以及地图比例尺之间、地图比例尺与数字化仪的长度单位之间的不一致。投影变换是地图制图的理论基础。第52页/共63页正解变换:通过建立两种投影之间的严密数学解析关系,直接由一种投影的数字化坐标精确变换为另一种投影的坐标 反解变换:以地理坐标为中间媒介,先将一种投影坐标反解出其地理坐标,然后再将其地理坐标带入另一种投影的坐标计算公式中,计算出新的投影坐标。第53页/共63页 中国地形图分幅与编号第54页/共63页7.4.2 中国地形图分幅与编号旧标准我国基本比例尺地形图分幅与编号,以1:100万地形图为基础延伸出1:50万、
19、1:25万、1:10万;再以1:10万为基础,延伸出1:5万、1:2.5万、1:1万三种比例尺;1:100万从赤道起向两极每纬差4为一行,至88,南北半球各分为22横列,依次编号A、B、.V;由经度180西向东每6一列,全球60列,以1-60表示,如北京所在1:100万图在第10行,第50列,其编号为 J-50 第55页/共63页7.4.2 中国地形图分幅与编号旧标准(续1)在1:100万图上,按经差3纬差2分成四幅1:50万地形图,编为A、B、C、D,如 J-50-A 按经差130纬差1分成16幅1:25万地形图,编为1、.16,如 J-50-1。按经差30纬差20分成144幅1:10万地形
20、图,编为1、.144,如 J-50-1。这三种比例尺各自独立地与1:100万地图的图号联系。第56页/共63页1:10万图的基础上每经差15纬差10分成四幅1:5万地形图,编为A、B、C、D,如 J-50-1-A1:5万图上每经差730纬差5分成四幅1:2.5万,编为1、2、3、4,如 J-50-1-A-11:10万图上每经差345纬差230分成64幅1:1万地形图,编为(1)、.(64),如 J-50-1-(1)1:1万图上每经差152纬差115分成四幅1:5000地形图,编为a、b、c、d,如 J-50-1-(1)-a7.4.2 中国地形图分幅与编号旧标准(续2)第57页/共63页7.4.
21、2 中国地形图分幅与编号旧标准(续3)第58页/共63页7.4.3 中国地形图分幅与编号新标准在100万的基础上划分,比例尺有明确的代号。行号数字码列号数字码比例尺编号100万图幅列号数字码100万图幅行号字符码ABCDEFGH,分别代表100万、50万、25万、10万、5万、2.5万、1万、5000。例如:J-50-E-013020第59页/共63页7.4.3 中国地形图分幅与编号新标准(续)第60页/共63页7.4.4 新旧标准的转换旧标准图幅的序号/(百万分幅该比例尺行向上图幅数)Row=商1Col 余数第61页/共63页例如 1:50万 的图幅 旧:J-50-C 新:J-50-B-0020011:25万 的图幅 旧:J-50-(11)新:J-50-C-0030031:10万的图幅 旧:J-50-13 新:J-50-D-0020017.4.4 新旧标准的转换(续)1:50万1:25万第62页/共63页感谢您的观看!第63页/共63页