第二章地理空间数学基础

第二章地理空间数学基础§2-1地球空间参考§2-2空间数据投影§2-3空间数据的采集§2-4GIS中的地理参照系§3-5GIS的数据质量§3-6空间数据标准空间数据的采集§2-1地球空间参考§2.1.1地球形状与地球椭球§2.1.2坐标系统§2.1.3高程基准一、地理空间的数学建构---如何建立地球表面的几何模型（形状与大小）1、最自然的面：包括海洋底部、高山、高原在内的固体地球表面，难以用一个简洁的数学式描述。2、相对抽象的面，即大地水准面假设一个当海水处于完全静止的平衡状态时从海平面延伸到所有大陆下部，而与地球重力方向处处正交的一个连续、闭合的水准面。§2.1.1地球形状与地球椭球一、地球形状与大小1、最自然的面测量基准—地球模型地球表面水准面大地水准面铅垂线地球椭球体地球重力也不均匀，大地水准面不是一个十分规则的表面，且不能用简单的数学公式来表达，大地水准面不能作为测量成果的计算面。一、地球的形状与大小•2、大地水准面：无数个水准面中，其中与平静的平均海水面相吻合并延伸穿过陆地、岛屿而形成的封闭曲面称。大地水准面的特性1。同一水准面上的重力位处处相等；2。同一水准面上任一点的铅垂线都与水准面相正交；3。由于地表起伏以及地球内质量分布不均匀，所以大地水准面是个复杂的曲面无法准确描述和计算；4。水准面和铅垂线是野外观测的基准面和基准线。为了测量成果计算的需要，选用一个同大地体相近的、可以用数学方法来表达的旋转椭球来代替地球---椭球体。二、椭球体模型abc一、地球的形状与大小地球椭球体模型双轴椭球体模型（旋转椭球体）x2a2y2a2z2b2++=1如：克拉索夫斯基椭球体地球椭球的基本元素常用符号a，b，α，e和e’表示。符号的名称和公式为(衡量形状和大小参数）：长半轴=a；短半轴=b；扁率α=(a-b)/a；第一偏心率=第二偏心率=222abae222bbae已知其中两个元素（包含a或b）,就可以推算其他三个元素。地球椭球椭球名称创立年代长半径akm短半径bkm扁率αWGS8419846378.1376356.7521∶298.261975年GRS椭球1980年代以后19756378.1406356.7551∶298.257苏联克拉索夫斯基椭球1953-80年代19406378.2456356.9121∶297美国海福特椭球1953年前19106378.3886356.9121∶298.3§2.1.2坐标系统球面坐标系统天文地理坐标系大地地理坐标系平面坐标系统：投影坐标系统高斯平面直角坐标系地方独立平面直角坐标系投影过程示意图§2.1.2坐标系统一、球面坐标系1、天文地理坐标系大地水准面和铅垂线是天文地理坐标系的基准面和线地面点的坐标是它沿铅垂线在大地水准面上投影点的经度和纬度（，）天文经度是首天文子午面与测站天文子午面的夹角，天文纬度为测站垂线方向与地球平均赤道面高是地面点沿铅垂线到大地水准面的距离2、大地地理坐标系（大地坐标）NSMP格林尼治天文台LBH大地经度(L)大地纬度(B)大地高(H)过地面点的子午面与起始子午面之间的夹角过地面点的法线与赤道面之间的夹角地面点沿法线至参考椭球面的距离G是依托地球椭球用定义原点和轴系以及相应基本参考面标示较大地域地理空间位置的参照系我国常用的大地坐标系大地坐标系参考椭球坐标原点北京54坐标系克拉索夫斯基椭球体原苏联西部的普尔科夫西安80坐标系IUGG（16届）椭球体1975陕西泾阳WGS84坐标系IUGG（17届）椭球体地球质心本地坐标系国家参考椭球（修改半径）§2.1.3高程系统1、基本概念：高程（绝对高程、海拔）-----地面点到大地水准面的铅垂距离。假定（相对）高程-----地面点到假定水准面的铅垂距离。高差-----两点间的高程之差。2.高程系统与高程基准我国国家高程系统：黄海高程系我国国家高程基准：1956年黄海高程系1985年国家高程基准◆我国取青岛附近黄海平均海水面为大地水准面，即零基准面；人为而定，相对稳定我国1956年取前6年的平均潮位作为大地水准面；1985年取1952年1979年共26年观测的平均潮位作为大地水准面)；◆水准原点建在青岛市内，作为我国的国家高程基准。1956年高程基准的高程为72.289米，1985年高程基准的高程为72.260米。2、我国高程系统与高程基准国家高程原点二地图比例尺P77比例尺表示法数字比例尺,1:10000或1/10000文字比例尺,万分之一图解比例尺与实地比例尺面积比例尺比例系数SF=实际比例尺/主比例尺主比例尺：在投影面上没有变形的点或线上的比例尺。局部比例尺：在投影面上有变形处的比例尺。第三章空间数据的采集和质量控制§3-2空间数据的地理参照系和控制基础地图投影mapprojections：为什么要进行投影将三维椭球体转变成二维平面地图，称地图投影地球椭球体为不可展曲面将不可展的地球椭球面展开成平面，并且不能有断裂，则图形必将在某些地方被拉伸，某些地方被压缩，故投影变形是不可避免的。§2-2空间数据投影投影过程示意图2.1.1、地图投影的内涵地图投影依据一定的数学法则，将不可展的地球曲面运用特定的数学方法展示到平面上，最终在地表面点与地图平面点之间建立一一对应的关系),(),(21LBfyLBfx(B,L)(x,y)数学解析法2.1.2、地图投影的分类按投影面分类按投影方式分类变形特征分类几何透视法按投影面分类方位投影圆锥投影圆柱投影按投影方式分类横轴正轴斜轴变形特征分类任意投影等角投影等积投影圆锥投影、圆柱投影、局部平面投影是最常用的三大类距离变形、面积变形、角度（方位）变形距离、面积、方位三个因素总是存在矛盾，如何兼顾，何者优先，决定了采用什么投影、计算方法、相关参数由地球椭球面（或球面）投影到平面引起的变形。表现在形状、大小两个方面。2.1.3投影变形地图投影变形的图解示例（摩尔维特投影－等积伪圆柱投影）长度变形角度变形地图投影变形的图解示例（UTM－横轴等角割圆柱投影）面积变形和长度变形投影变形示意图地图投影的变形示意地图投影的变形关于地图投影的几点结论：实现等角、等面积、等距离同时做到的投影不存在投影方式有多种多样，一个国家或地区依据自己所处在的经纬度、幅员大小以及图件用途选择投影方式在大于1：10万的大比例尺图件中，各种投影带来的误差可以忽略。2.1.4常用的地图投影我国1:1万-1：50万的地形图采用高斯—克吕克投影(Gauss-Krüger)我国1：500——1：5000的图，采用城市平面局域投影英美国家常用横轴墨卡托投影（UTM）高斯——克吕格投影的几何概念一、高斯-克吕格平面直角坐标系（1）.高斯投影的概念（2）.高斯投影的分带和编号（3）.高斯平面直角坐标系（4）.高斯投影的关系式高斯-克吕格投影高斯(Gauss,CarlFriedrich,1777-1855)，德国数学家，天文学家，物理学家。（1）.高斯投影的概念等角横轴椭圆柱投影特点：•1.中央子午线和地球赤道投影成为直线且为投影对称轴；•2.等角投影，经纬线投影后保持相互垂直关系；•3.中央子午线上没有长度变形。投影演示变形特点：1.中央经线上没有长度变形；2.同一纬线上，离中央经线越远，变形越大；3.同一经线上，纬度越低，变形越大；4.等变形线为平行于中央经线的直线。高斯——克吕格投影（等角横切椭圆柱投影）优点:1.等角性质，适用于系列比例尺地图的适用与编制；2.经纬网与直角坐标的偏差不大；3.计算量小，只需计算一个带，全球通用。6分带法3分带法（2）.高斯投影的分带和编号国家坐标系：高斯-克吕格直角坐标系高斯投影的分带和编号6°分带法：从格林尼治零度经线起，每6°分为一个投影带，全球共分为多少个投影带？投影带号n和中央经线经度L0的计算公式为：12313141516172930311201352426283032345759616带3带06121872849616818039157587991711773°分带法：从东经1°30′起，每3°为一带，将全球划分为多少个投影带？投影带号n和中央经线经度L0的计算公式为：L0=(6°*n—3°)L0=(3°*n)课堂思考题1、如何根据某一点的经度来计算其所在投影带的带号？EX：A点的经度为东经119°32′，则A点所在的6°带和3°带的带号分别为多少？12313141516172930311201352426283032345759616带3带06121872849616818039157587991711772、为什么3°带是东经1.5°开始的？（3）.高斯平面直角坐标系及坐标本初子午线=0°经线纬线XYZNEWPOR••赤道=0°••-经度-纬度R–地球半径O–地球球心XYo赤道中央经线高斯平面直角坐标系的特点1、坐标轴的特点o赤道中央子午线XYⅠⅡⅢⅣ△X△Y2、象限及方位角的特点3、坐标计算的特点XA=X0+△X=X0+DcosαYA=Y0+△Y=Y0+DsinαADα高斯平面直角坐标特点我国规定将各带纵坐标轴西移500公里（实质是将所有y值加上500公里）Y坐标值前再加上所在投影带带号EX:假设某点在第18带Y=743.353kmX=464.784kmY+500km加带号Y=243.353kmX=464.784km原始坐标（真实坐标）Y=18743.353kmX=464.784km高斯坐标（名义坐标）500km第X带纵坐标:表示离赤道公里数横坐标:前二位表示投影带号后几位表示离移轴纵线的公里数如:18210表示第18带,例移轴纵线210公里,即中央经线以西290公里移轴纵线赤道中央经线500km4、高斯投影的关系式高斯-克吕格投影正解公式：（B,L）→(X,Y)，原点纬度0，中央经度L0上面公式中东纬偏移FE=500000米+带号*1000000；高斯-克吕格投影比例因子k0=1高斯-克吕格投影反解公式：(X,Y)→（B,L），原点纬度0，中央经度L0补充材料2、补充资料：高斯投影的演示高斯投影1、补充资料：坐标换算程序坐标换算程序高斯-克吕格投影与UTM通用横轴墨卡托投影(universaltransverseMercatorcoordinatesystem:UTM)：等角横轴割椭圆柱投影高斯-克吕格投影的中央经线长度比等于1UTM中央经线长度比等于0.9996长度比:投影面上一微小线段ds´和球面上相应微小线段ds之比。用公式表示为：μ=ds´/ds长度比是一个变量，它不仅随着点的位置不同而变化，还随着方向的变化而变化。长度比是指某点某方向上微小线段之比。2.1.4GIS中的地理参照系地球自然表面点位坐标系的确定包括两个方面的内容：一是地面点在地球椭球体面上的投影位置，采用地理坐标系；二是地面点至大地水准面上的垂直距离，采用高程系。一．地理参照系1、经纬度坐标系（地理坐标）对空间定位有利，但难以进行距离、方向、面积量算。2、笛卡儿平面坐标系便于量算和进一步的空间数据处理和分析。3、高程系统描述空间点在垂直高度上的特性--高程——由高程基准面起算的地面点的高度。地理坐标系1.构成经度、纬度2.获取天文测量、大地测量3.我国的大地坐标系1954年北京坐标系，克拉索夫斯基椭球，原点在原苏联西部的普尔科夫80西安坐标系，1978年建立自己的坐标系，原点在陕西泾阳县，采用1975国际大地测量协会推荐的地球参数。一、地理参照系平面坐标系P79将椭球面上的点通过地图投影的方法投影到平面上，通常使用平面坐标系。1.平面极坐标系：矢径(极径)、矢量角2.平面直角坐标：笛卡尔坐标坐标平移,坐标旋转(自学)第三章空间数据的采集和质量控制§3-2空间数据的地理参照系和控制基础一、地理参照系二、GIS的地理基础--控制基础各种GIS的数据源、服务目的和各自特征可以不同，但均有自身统一的地理基础。地理基础是地理信息数据表示格式与规