第五章地理信息系统-坐标与投影.

物流信息管理大连海事大学第五章地理信息系统实验评分标准实验任务：全球船舶轨迹图注册、描画5分导入外部数据到地图属性表5分创建统计图、专题图5分拓扑化全球航线网络5分创建路径查询程序5分创建最短路径程序5分加分项：创建流量分配程序+5-10分与外部程序集成+5-10分其他+5-10分空间参照系统和地图投影一个要素要进行定位，必须嵌入到一个空间参照系中，因为GIS所描述是位于地球表面的信息，所以根据地球椭球体建立的地理坐标（经纬网）可以作为所有要素的参照系统。因为地球是一个不规则的球体，为了能够将其表面的内容显示在平面的显示器或纸面上，必须进行坐标变换。本章讲述了地球椭球体参数、常见的投影类型。考虑到目前使用的1：100万以上地形图都是采用高斯——克吕格投影，本章最后又对该种投影类型和相关的地形图分幅标准做了简单介绍。1.地球椭球体基本要素(一)地球形状地球自然表面是一个起伏不平、十分不规则的表面，这个高低不平的表面无法用数学公式表达，也无法进行运算。所以在量测与制图时，必须找一个规则的曲面来代替地球的自然表面。当海洋静止时，它的自由水面必定与该面上各点的重力方向（铅垂线方向）成正交，我们把这个面叫做水准面。但水准面有无数多个，其中有一个与静止的平均海水面相重合。可以设想这个静止的平均海水面穿过大陆和岛屿形成一个闭合的曲面，这就是大地水准面。1.地球椭球体基本要素1.地球椭球体基本要素(二)地球大小椭球体名称年代长半轴（米）短半轴（米）扁率白塞尔(Bessel)1841637739763560791：299.15克拉克(Clarke)1880637824963565151：293.5克拉克(Clarke)1866637820663565841：295.0海福特(Hayford)1910637838863569121：297克拉索夫斯基1940637824563568631：298.3IUGG1976637816063567751：298.25WGS8419846378137.06356752.31：298.262.坐标系2.1地理坐标系地理坐标系是以地理极(北极、南极)为极点。通过A点作椭球面的垂线，称之为过A点的法线。法线与赤道面的交角，叫做A点的纬度ψ。过A点的子午面与通过英国格林尼治天文台的子午面所夹的二面角，叫做A点的经度λ。2.坐标系2.2平面坐标系地理坐标是一种球面坐标。由于地球表面是不可展开的曲面，也就是说曲面上的各点不能直接表示在平面上，因此必须运用地图投影的方法，建立地球表面和平面上点的函数关系，使地球表面上任一点由地理坐标（φ、λ）确定的点，在平面上必有一个与它相对应的点，平面上任一点的位置可以用极坐标或直角坐标表示。2.坐标系2.2平面坐标系直接建立在球体上的地理坐标，用经度和纬度表达地理对象位置建立在平面上的直角坐标系统，用（x，y）表达地理对象位置投影2.坐标系2.2平面坐标系坐标系的平移x=x’+ay=y’+bO'OXYX'Y'baP2.坐标系2.2平面坐标系坐标系的旋转x=x’cosθ+y’sinθy=y’cosθ-x’sinθOXYX'Y'Pθ2.坐标系2.2平面坐标系坐标系的平移加旋转x=x’cosθ+y’sinθ+ay”=y’cosθ-x’sinθ+bOXYX'Y'PθO'Y''X''3.地图投影的基本概念•为什么要进行投影？•地图投影实质•投影变形•投影方法•投影选择所考虑的因素3.地图投影的基本概念(一)为什么要进行投影•将地球椭球面上的点映射到平面上的方法，称为地图投影•地理坐标为球面坐标，不方便进行距离、方位、面积等参数的量算•地球椭球体为不可展曲面•地图为平面，符合视觉心理，并易于进行距离、方位、面积等量算和各种空间分析3.地图投影的基本概念(二)地图投影的实质建立地球椭球面上经纬线网和平面上相应经纬线网的数学基础，也就是建立地球椭球面上的点的地理坐标（λ，φ）与平面上对应点的平面坐标（x，y）之间的函数关系：当给定不同的具体条件时，将得到不同类型的投影方式。),(),(21fyfx3.地图投影的基本概念3.地图投影的基本概念(三)投影变形将不可展的地球椭球面展开成平面，并且不能有断裂，则图形必将在某些地方被拉伸，某些地方被压缩，故投影变形是不可避免的。•长度变形•面积变形•角度变形3.地图投影的基本概念(四)投影分类变形分类：等角投影：投影前后角度不变等面积投影：投影前后面积不变；任意投影（等长投影）：角度、面积、长度均变形，但某一方向长度不变投影面：横圆柱投影：投影面为横圆柱圆锥投影：投影面为圆锥方位投影：投影面为平面投影面位置：正轴投影：投影面中心轴与地轴相互重合斜轴投影：投影面中心轴与地轴斜向相交横轴投影：投影面中心轴与地轴相互垂直相切投影：投影面与椭球体相切相割投影：投影面与椭球体相割3.地图投影的基本概念3.地图投影的基本概念(五)投影选择因素制图区域的地理位置、形状和范围制图比例尺地图内容出版方式3.地图投影的基本概念(六)GIS中的投影GIS以地图方式显示地理信息，而地图是平面，地理信息则在地球椭球上，因此地图投影在GIS中不可缺少。GIS数据库中地理数据以地理坐标存储时，则以地图为数据源的空间数据必须通过投影变换转换成地理坐标；而输出或显示时，则要将地理坐标表示的空间数据通过投影变换变换成指定投影的平面坐标。GIS中，地理数据的显示可根据用户的需要而指定投影方式，但当所显示的地图与国家基本地图系列的比例尺一致时，一般采用国家基本系列地图所用的投影。3.地图投影的基本概念(七)常用的地图投影方式1)世界地图的投影世界地图的投影主要考虑要保证全球整体变形不大，根据不同的要求，需要具有等角或等积性质，主要包括：等差分纬线多圆锥投影、正切差分纬线多圆锥投影（1976年方案）、任意伪圆柱投影、正轴等角割圆柱投影。2)半球地图的投影东、西半球有横轴等面积方位投影、横轴等角方位投影；南、北半球有正轴等面积方位投影、正轴等角方位投影、正轴等距离方位投影。3.地图投影的基本概念3)各大洲地图投影亚洲地图的投影：斜轴等面积方位投影、彭纳投影。欧洲地图的投影：斜轴等面积方位投影、正轴等角圆锥投影。北美洲地图的投影：斜轴等面积方位投影、彭纳投影。南美洲地图的投影：斜轴等面积方位投影、桑逊投影。澳洲地图的投影：斜轴等面积方位投影、正轴等角圆锥投影。拉丁美洲地图的投影：斜轴等面积方位投影。3.地图投影的基本概念4)我国常用投影方式1：100万：兰勃投影（正轴等积割圆锥投影）大部分分省图、大多数同级比例尺也采用兰勃投影1：50万、1：25万、1：10万、1：5万、1：2.5万、1:1万、1：5000采用高斯—克吕格投影。4.高斯——克吕格投影由于这个投影是由德国数学家、物理学家、天文学家高斯于19世纪20年代拟定，后经德国大地测量学家克吕格于1912年对投影公式加以补充，故称为高斯——克吕格投影。高斯——克吕格投影在英美国家称为横轴墨卡托投影。美国编制世界各地军用地图和地球资源卫星象片所采用的全球横轴墨卡托投影（UTM）是横轴墨卡托投影的一种变型。高斯克吕格投影的中央经线长度比等于1，UTM投影规定中央经线长度比为0.9996。在6度带内最大长度变形不超过0.04%。4.高斯——克吕格投影高斯克吕格投影的中央经线和赤道为互相垂直的直线，其他经线均为凹向并对称于中央经线的曲线，其他纬线均为以赤道为对称轴的向两极弯曲的曲线，经纬线成直角相交。在这个投影上，角度没有变形。中央经线长度比等于1，没有长度变形，其余经线长度比均大于1，长度变形为正，距中央经线愈远变形愈大，最大变形在边缘经线与赤道的交点上；面积变形也是距中央经线愈远，变形愈大。为了保证地图的精度，采用分带投影方法，即将投影范围的东西界加以限制，使其变形不超过一定的限度，这样把许多带结合起来，可成为整个区域的投影（图4-12）。4.高斯——克吕格投影4.高斯——克吕格投影我国规定1：1万、1：2.5万、1：5万、1：10万、1：25万、1：50万比例尺地形图，均采用高斯克吕格投影。1：2.5至1：50万比例尺地形图采用经差6度分带，1：1万比例尺地形图采用经差3度分带。6度带是从0度子午线起，自西向东每隔经差6为一投影带，全球分为60带，各带的带号用自然序数1，2，3，…60表示。即以东经0-6为第1带，其中央经线为3E，东经6-12为第2带，其中央经线为9E，其余类推（图4-13）。3度带，是从东经1度30分的经线开始，每隔3度为一带，全球划分为120个投影带。图4-13表示出6度带与3度带的中央经线与带号的关系。4.高斯——克吕格投影4.高斯——克吕格投影X坐标值在赤道以北为正，以南为负；Y坐标值在中央经线以东为正，以西为负。我国在北半球，X坐标皆为正值。Y坐标在中央经线以西为负值，运用起来很不方便。为了避免Y坐标出现负值，将各带的坐标纵轴西移500公里，即将所有Y值都加500公里。由于采用了分带方法，各带的投影完全相同，某一坐标值（x，y），在每一投影带中均有一个，在全球则有60个同样的坐标值，不能确切表示该点的位置。因此，在Y值前，需冠以带号，这样的坐标称为通用坐标。5.地形图的分幅和编号国家基本比例尺地形图有1：1万、1：2.5万、1：5万、1：10万、1：20万、1：50万和1：100万七种。普通地图通常按比例尺分为大、中、小三种，一般以1：10万和更大比例尺的地图称为大比例尺地图；1：10万至1：100万的称为中比例尺地图；小于1：100万的称为小比例尺地图。对于一个国家或世界范围来讲，测制成套的各种比例尺地形图时，分幅编号尤其必要。通常这是由国家主管部门制定统一的图幅分幅和编号系统。5.地形图的分幅和编号5.1地形图的分幅目前，我国采用的地形图分幅方案，是以1：100万地形图为基准，按照相同的经差和纬差定义更大比例尺地形图的分幅。百万分之一地图在纬度0o—60o之间的图幅，图幅大小按经差6o，纬差4o分幅；在60o—76o之间的图幅，其经差为12o，纬差为4o；在76o—80o之间图幅的经差为24o，纬差为4o，所以各幅百万分之一地图都是经差6o，纬差4o分幅的。每幅百万分之一内各级较大比例尺地形图的划分，按规定的相应经纬差进行，其中，1：50万、1：20万、1：10万三种比例尺地形图，以百万分之一地图为基础直接划分。一幅百万分之一地形图划分四幅1：50万地形图，每幅为经差3o，纬差2o；一幅百万分之一地图划分为36幅1：20万地形图，每幅为经差1o，纬差40’；一幅百万分之一地图划分144幅1:10万地形图，每幅为经差30’，纬差20’。每幅大于1：10万比例尺的地形图，则以1：10万图为基础进行逐级划分，一幅1：10万地形图划分四幅1：5万地形图；一幅1：5万地形图划分为四幅1：2.5万地形图。在1：10万图的基础上划分为64幅1：1万地形图；一幅1：1万地形图又划分为4幅1：5000地形图，（见表4-1）。5.地形图的分幅和编号5.1地形图的分幅比例尺（万）图幅大小图幅间的数量关系经度纬度1：1006度4度11：503度2度411：201度40分36911：1030分20分14436411：515分10分57614416411：2.57.5分5分23045766416411:13分45秒2.5分92162304256641645.地形图的分幅和编号5.2分幅编号(一)1：100万地图的编号该种地形图的编号为全球统一分幅编号。列数：由赤道起向南北两极每隔纬差4o为一列，直到南北88o（南北纬88o至南北两极地区，采用极方位投影单独成图），将南北半球各划分为22列，分别用拉丁字母A、B、C、D……V表示。行数：从经度180o起向东每隔6o为一行，绕地球一周共有60行，分别以数字1、2、3、4……60表示。由于南北两半球的经度相同，规定在南半球的图号前加一个S，北半球的图号前不加任何符号。一般来讲，把列数的字母写在前，行数的数字写在后，中间