zoukankan      html  css  js  c++  java
  • GIS中的空间参考

    GIS一般都是研究的基于地球的某个区域,例如一个国家、省或市的情况,既然地球上都有经纬度来标识,那么直接用经纬度来标识物体的位置不可以吗?但如果应用不一样,在实际中选择的空间参考系也是不一样的。例如我们使用GIS系统在做一个房间的布置的时候,就应该不会经纬度和高程数据来标识物体的位置。

    但其实我们在GIS行业中主要研究的还是比较大的一片区域,可以都是相对于地球可以用经纬度表示的。但实际上地球是一个两级稍扁,赤道略鼓的不规则球体。并且在地球表面上既有高山也有海洋、还有低于海平面的盆地等,这么复杂的形状是没法用数学公式定义的。那么可以把地球当成一个圆球的形状吗?这样很容易用数学公式描述。但这种情况下做直径几厘米或几十厘米的地球仪是可以的,但绘制大比例尺地图或者进行精密的计算时,这样的计算精度就太差了,根本不能满足需求。所以一般情况下会把地球模拟成一个椭球体,但如果有更高精度要求话应该模拟的更复杂。

    地球大致的形状

    地球可能实际的形状

    • 1979年,国际大地测量及地球物理协会给出的地球椭球体的参数数据如下:
    • 赤道半径(长轴半径): 6378.14 km
    • 极半径(短轴半径): 6356.76 km
    • 赤道周长:40075.7 km
    • 地球表面积:510100934 km2
    • 地球平均半径:6371 km

    地理坐标系

    椭球体

    如果要定义一个空间参考,首先要确定模拟地球的椭球体的一些参数,在历史上有很多椭球体模型的,例如我国在1952以前采用海福特(Hayford)椭球体,从1953年-1980年采用克拉索夫斯基椭球体,随着人造卫星的发射和测量技术的发展,我们开始使用16届国际大地测量及地球物理联合会公布的椭球体,名称为GRS(1975)。但实际上我们拿到的数据和做的系统使用的椭球体是很多的。可以说是各种各样。

    基准面

    有了椭球体来模拟地球,我们就可以比较准确对应地球表面的点和椭球体表面的点,然后用经纬度描述。那基准面是干嘛的呢?我们前面也提到,地球是一个表面高低起伏不平的椭球体,大部分是海洋,最低处和最高处相差近20km。椭球体我们确定了,那椭球体的怎么和地球体去贴和?基于整个地球的还好,但如果是基于亚洲,基于中国的一个空间参考呢?例如一个基于中国的空间参考,椭球体已经定义好了,在没有基准面的情况下,直接把椭球体贴合到地球上,地球表面是起伏不平的,因为要考虑整体贴合,所以对中国而言就贴合的不是很好。但如果我们只研究中国这块区域,是不是可以把这个椭球体调整一下,让椭球体的表面更加贴近中国的区域?答案是可以的,调整的参数其实就是基准面了。所以这就不难理解每个州、国家甚至区域都有自己的基准面。例如我们通常说的北京54和西安80指的就是基准面。

    当地球和椭球体整体贴合时的情况

    当椭球体要贴合某一区域时的情况(贴合左边的某一区域)

    我们一旦说基准面的时候,这个基准面肯定是包含一个椭球体的,在一个椭球体基础上定义出不同的基准面。椭球体和基准面是1对多的关系。像我们国家的北京54基准面使用的椭球体是克拉索夫斯基椭球体,而西安80基准面使用的椭球体是我们上面所说的GRS(1975)椭球体。

    基准面分为地心基准面和大地基准面。如果是基于整个地球的基准面呢,应该是是为地心为原点的,例如WGS84就是以地球质心为原点的地心基准面,是目前使用最广放的全球的基准面。WGS84基准面使用的椭球体就叫WGS84椭球体,这些椭球体的参数都可以查到。

    那么基准面是怎么定义的呢?

    在测量学中,大地基准面(Geodetic datum),设计用为最密合部份或全部大地水准面的数学模式。它由椭球体本身及椭球体和地表上一点视为原点间之关系来定义。此关系能以6个量来定义,通常(但非必然)是大地纬度、大地经度、原点高度、原点垂线偏差之两分量及原点至某点的大地方位角。

    GIS中的基准面通过当地基准面向WGS1984的转换7参数来定义,转换通过相似变换方法实现,假设 Xg、Yg、Zg表示WGS84 地心坐标系的三坐标轴,Xt、Yt、Zt 表示当地坐标系的三坐标轴,那么自定义基准面的 7参数分别为:三个平移参数ΔX、ΔY、ΔZ表示两坐标原点的平移值;三个旋转参数 εx、εy、εz表示当地坐标系旋转至与地心坐标系平行时,分别绕 Xt、Yt、Zt 的旋转角;最后是比例校正因子,用于调整椭球大小。(以上两段从网上摘抄,我自己也不怎么明白)

    我国主要采用的地理坐标系

    ●1954年北京坐标系(Beijing Geodetic Coordinate System,l954)

    该坐标系是通过与原苏联1942年坐标系联测而建立的,其原点不在北京,而是在苏联普尔科沃。该坐标系采用克拉索夫斯基椭球(Krasovsky-1940)作为参考椭球,高程系统采用正常高,以1956年黄海平均海水面为基准。

    ●1980年西安坐标系

    其大地原点设在西安西北的永乐镇,简称西安原点。椭球参数选用1975年国际大地测量与地球物理联合会第16届大会的推荐值。简称IUGG-75地球椭球参数或IAG-75地球椭球。2000年后的空间数据常采用该坐标系。

    ●WGS84坐标系(WGS一84 Coordinate System)

    在GPS定位中,定位结果属于WGS-84坐标系。该坐标系是使用了更高精度的VLBL、SLR等成果而建立的。坐标系原点位于地球质心,Z轴指向BIH1984.0协议地极(CTP)。用于GPS定位系统的空间数据采用该坐标系。

    ArcGIS对地理坐标系的定义

    ArcGIS地理坐标系中的类是GeographicCoordinateSystemClass类,看下该类的主要接口IGeographicCoordinateSystem,主要的定义代码如下:

    //角度单位  
    IAngularUnit CoordinateUnit { get; }  
    //基准面  
    IDatum Datum { get; }  
    //本初子午线定义  
    IPrimeMeridian PrimeMeridian { get; }  

    基准面接口为:IDatum,代码定义如下:

    ISpheroid Spheroid { get; }  

    Spheroid是椭球体,定义如下:

    //椭球体扁率  
    double Flattening { get; }  
    //长半轴值  
    double SemiMajorAxis { get; }  
    //短半轴值  
    double SemiMinorAxis { get; }  

    在ArcGIS中定义地理坐标系的界面如下:

    从界面上我们可以看出该地理坐标系是WGS84地理坐标系,使用的基准面名称为D_WGS_1984基准面,基准面使用的椭球体为WGS_1984,椭球体的长半轴长度为6378137,椭球体的短半轴长度为6356752.3142451793,扁率为198.25722356300003,因为长半轴、短半轴和扁率只要知道两个值就能知道第三个,所以只有两个可设置。下面角度单位设置的时度,每单位的弧度值是0.017453292519943295,因为π=3.1415925...=180度,所以设置的这个值0.017453292519943295*180应该等于π。

    下面就是设置本初子午线,设置为格林尼治,也就是经度0度0分0秒。

    下面截图是地理坐标系beijing54的定义:

    该坐标系的基准面是D_beijing_1954,该基准面使用的椭球体为krasovsky_1940,也就是克拉索夫斯基_1940。

    下面是xian80的定义:

    椭球体使用的时西安80,西安80的长半轴和短半轴基本上和GRS(1975)差不多。

    投影坐标系

    为什么会有投影的概念呢?

    因为地球是一个立体的球,而我们在二维地图上展示或者打印出图的都是平面的,怎么把一个球上面的曲面的东西展示到一个二维平面上。如果按照把地球沿着本初子午线切开,然后平铺,因为地球是从赤道到两级地球的沿着纬线的周长是越来越短的,如果沿赤道展开,两级肯定是展不开的,并且有裂缝。或者赤道有褶皱才行。我们在看二维地图的时候,有的经纬网是横平竖直的,有的却是向两极弯曲的,出现这种情况就是投影的不同导致的。

    按变形性质地图投影可以分三类:等角投影、等积投影和任意投影。

    国内外常用的投影有:

    ●通用横轴墨卡托投影(UNIVERSAL-TRANSVRSE-MERCATOL,UTM) 

    该投影为横轴等角割圆柱投影,该投影将地球分为60个投影带,每带经度差为6度,已被许多国家作为地形图的数字基础。一般能从南纬80度到北纬84度的范围内使用该投影,对于两极地区则采用UPS(通用球面极)投影。

    涉及的参数为:向东平移值、向北平移值、中线子午线经度、投影原点纬度、坐标单位、比例尺(若用实地坐标,该值为1)、椭球参数等。

     ●亚尔伯斯等积圆锥投影(Albers)

    即为双标准纬线投影,也即正轴等面积割圆锥投影。该投影经纬网的经线为辐射直线,纬线为同心圆圆弧。Albers投影的应用在编制一些行政区划图、人口地图、地势图等方面应用较广。如中国地势图,即是以第一标准纬度Q1=25度,第二标准纬度Q2=45度的该投影;水利部水利信息中心建立的全国1:25万水利空间数据库采用的即为该投影(中央子午线经度110度,Q1=25度,Q2=47度)。

    涉及的参数为:向东平移值、向北平移值、中线子午线经度、第一标准纬度、第二标准纬度、投影原点纬度、坐标单位、比例尺(若用实地坐标,该值为1)、椭球参数等。

    ●兰伯特等角圆锥投影(Lambert-Conformal-Conic)

    也称兰伯特正形圆锥投影,该投影的微分圆投影后仍为圆形。经线为辐射直线,纬线为同心圆圆弧。指定两条标准纬度线Q1、Q2,在这两条纬度线上没有长度变形,此种投影也叫做等角割圆锥投影,常用来编制中、小比例尺地图。等角圆锥投影有广泛的应用,特别适宜于作为中纬度处沿纬线伸展的制图区域之投影,我国的分省地图即为Q1=25度,Q2=45度的兰伯特等角圆锥投影。1962年以后,1:100地图采用了等角圆锥投影(南纬度80度,北纬度84度),极区附近,采用等角方位(极球面投影)。

    涉及的参数为:向东平移值、向北平移值、中线子午线经度、第一标准纬度、第二标准纬度、投影原点纬度、坐标单位、比例尺(若用实地坐标,该值为1)、椭球参数等。

    ●墨卡托投影(Mercator)

    等角正圆柱投影也称墨卡托投影,经纬线投影为互相正交的平行直线。该投影在航海、航空应用很广。使用该投影,等角航线在地图上是一条直线。值得注意的是,等角航线是球面上两点间对所有经线保持等方位角的特殊曲线,不是两点间的最近路线,是一条以极点为渐近点的螺旋曲线。

    涉及的参数为:向东平移值、向北平移值、中线子午线经度、无变形纬度、坐标单位、比例尺(若用实地坐标,该值为1)、椭球参数等。

    ●高斯-克吕格投影(Gauss-Kruger)

    该投影是等角横切椭圆柱投影,由德国数学家高斯提出,后经克吕格扩充并推算出计算公式,故称为高斯-克吕格投影,简称高斯投影。该投影以中央经线和赤道投影后为坐标轴,中央经线和赤道交点为坐标原点,纵坐标由坐标原点向北为正,向南为负,规定为X轴,横坐标从中央经线起算,向东为正,向西为负,规定为Y轴。所以,高斯-克吕格坐标系的X、Y轴正好对应一般GIS软件坐标系中的Y和X。

    为了控制变形,本投影采用分带的办法,我国1:2.5万——1:50万地形图均采用6度分带,1:1万及更大比例尺地形图采用3度分带,以保证必要的精度。6度分带从格林威治零度经线起,每6度分为一个投影带,该投影将地区划分为60个投影带,已被许多国家作为地形图的数字基础。一般从南纬度80到北纬度84度的范围内使用该投影。3度分带法从东经1度30分算起,每3度为一带。这样分带的方法在于使6度带的中央经线均为3度带的中央经线;在高斯克吕格6度分带中中国处于第13带到23带共12个带之间;在3度分带中,中国处于24带到45带共22带之间。

    对于两极地区则采用UPS投影(通用球面极投影)。高斯-克吕格投影通常投影带为6度范围或3度,超过了6度后变形会增大。一般常用来制作中、大比例尺的地图投影,如1:50万、1:10万、1:5万、1:1万等。

    由于高斯-克吕格投影每一个投影带的坐标都是对本带坐标原点的相对值,所以各带的坐标完全相同,使用时只需要一个带号即可。

    在高斯坐标系中,为了避免横坐标Y有负值,将其起算原点向西移动500公里,即对横坐标Y值按代数法加上500000米。此外,可以在计算出米的和数前面加上带号,以便识别该点位于何带。高斯-克吕格坐标系的横坐标最多为6位,纵向最多为7位。若横向为8位,则前2位为带号。

    涉及的参数为:向东平移值、向北平移值、中线子午线经度、投影原点纬度、坐标单位、比例尺(若用实地坐标,该值则为1)、椭球参数等。

    ●普通多圆锥投影(Polyconic)

    普通多圆锥投影的坐标为对称于中央经线和赤道的曲线,纬线投影为同轴圆圆弧,弧心位于中央直径线上,中央经线是直线,M=1。该投影适用于沿中央经线延伸的区域(15度范围内)。常用于编制中、小比例尺的数字基础。该投影在美国被广泛应用,是1:100万地图投影的基础。

    涉及的参数为:向东平移值、向北平移值、中线子午线经度、投影原点纬度、坐标单位、比例尺(若用实地坐标,该值则为1)、椭球参数等。

    国内外常采用的投影还有正轴等距离割圆锥投影、横向墨卡托投影(Transverse-Mercator)、等角方位投影(Azimuthal-Conformal)、兰伯特等积方位投影(Lambert-Azimuthal-Equal-Area)、波斯托等距离方位投影(Azimuthal-Equidistant)、心射切面投影或球心投影、正射方位投影、垂直近距离方位投影或外心投影(Vertical-Near-Side-Perspective)、正弦曲线等面积伪圆柱投影(Simusoidal)、等矩形圆柱投影(Equidistant-Cylindrical)、米勒圆柱投影(Miller-Cylindrical)、斜轴墨卡托投影、两极球极平面投影、范·德·格林顿投影(Van-Der-Grinten-I)。

    一个地理坐标系可以使用不用的投影,例如现在中国常用xian80坐标系可以有Gauss_Kruger投影、Albers投影。

    例如全球地图,在WGS84下,使用Web墨卡托的投影显示如下:

    web墨卡托的显示参数

    web墨卡托地图显示效果

    xian80 Albers投影坐标系参数

  • 相关阅读:
    【PAT甲级】1043 Is It a Binary Search Tree (25 分)(判断是否为BST的先序遍历并输出后序遍历)
    Educational Codeforces Round 73 (Rated for Div. 2)F(线段树,扫描线)
    【PAT甲级】1042 Shuffling Machine (20 分)
    【PAT甲级】1041 Be Unique (20 分)(多重集)
    【PAT甲级】1040 Longest Symmetric String (25 分)(cin.getline(s,1007))
    【PAT甲级】1039 Course List for Student (25 分)(vector嵌套于map,段错误原因未知)
    Codeforces Round #588 (Div. 2)E(DFS,思维,__gcd,树)
    2017-3-9 SQL server 数据库
    2017-3-8 学生信息展示习题
    2017-3-5 C#基础 函数--递归
  • 原文地址:https://www.cnblogs.com/rainbow70626/p/7989732.html
Copyright © 2011-2022 走看看