Arcgis中地理坐标系和投影坐标区别及操作

GIS基础教程之坐标系 在ArcGIS中对数据定义坐标系,可以使用Catalog的数据管理功能,也可以使用ArcToolbox的定义投影(Define Projection)工具。在Catalog中查看数据的坐标系属性并修改定义投影工具注意,上面这两种方法都是直接修改数据的坐标系属性,适用条件是:数据坐标系信息缺失或坐标系信息错误如何走好GIS数据处理的第一步显示坐标系会直接影响地图出图的形状。专题要素到标准地图的整合——以中国5A级景区分布图为例。 阅读详情
AI摘要

本文深入解析GIS中的地理坐标系与投影坐标系,阐述两者的关系及应用,并探讨数据的真实坐标系、属性坐标系与地图坐标系的区别。

(一)两种坐标系

坐标系(Coordinate System)的概念为:“In geometry, a coordinate system is a system which uses one or morenumbers, or coordinates, to uniquely determine the position of a point or othergeometric element on a manifold such as Euclidean space”(https://en.wikipedia.org/wiki/Coordinate_system)。简单的说,有了坐标系,我们才能够用一个或多个“坐标值”来表达和确定空间位置。没有坐标系,坐标值就无从谈起,也就无法描述空间位置。 

在ArcGIS中,或者说在GIS中,我们遇到的坐标系一般有两种:

1)地理坐标系(Geographic Coordinate System);

2)投影坐标系(Projected Coordinate System)。

地理坐标系进行地图投影后就变成了投影坐标系。地图投影(Map Projection)是按照一定的数学法则将地球椭球面上点的经维度坐标转换到平面上的直角坐标。地图投影的理论知识请参考其他资料,此处不做叙述。需要说明的是,也有将“坐标系(CoordinateSystem)”称为“空间参考(Spatial Reference)”的情况,例如在ArcGIS中栅格数据的属性里面。


重要的事情说三遍:

这里要讲的不是“投影(Projection)”而是“坐标系(Coordinate System)”!

这里要讲的不是“投影(Projection)”而是“坐标系(Coordinate System)”!

这里要讲的不是“投影(Projection)”而是“坐标系(Coordinate System)”!

尽管投影是介绍坐标系的一个绕不开的重要内容。但是,首先,此文是围绕坐标系展开的。其次,说三遍是为了强调投影和坐标系的本质区别。坐标系是数据或地图的属性,而投影是坐标系的属性。一个数据或一张地图一定有坐标系,而一个坐标系可以有投影也可以没投影。只有投影坐标系才有投影,地理坐标系是没有投影的。因此,一个数据或一张地图亦是可以有投影也可以没投影的。当然,非要较真,把具有地理坐标系的数据显示在平面地图上肯定也有一个投影的过程。严格来讲:我们只能说“数据或地图的坐标系”和“坐标系的投影”,而不能说“数据或地图的投影”。也许是大家平时都比较随意,尽管都是知道二者的区别的,但是却在很多想说坐标系的时候就随口说成了投影。因此,当你说“数据的投影”和“投影转换”时,可以考虑下你是不是想说“数据的坐标系”和“坐标系转换”。


先抬出重要的总结:地理坐标系经过投影后变成投影坐标系,投影坐标系因此由地理坐标系和投影组成,投影坐标系必然包括有一个地理坐标系。图1概括了两种坐标系的联系:


图1 ArcGIS 中“地理坐标系(GCS)”与“投影坐标系(PCS)”的联系


下面以一个具体示例来初识ArcGIS中的坐标系,其全部参数拷贝在下面。这一示例是一个“投影坐标系(Projected Coordinate System)”,其名称是“WGS_1984_UTM_Zone_50N”。“WKID”是该坐标系的编号,“ESPG”是“European Petroleum Survey Group”的缩写,表示其由“欧洲石油调查组织”发布。可知,“WGS_1984_UTM_Zone_50N”这个投影坐标系由两部分组成:名为“Transverse_Mercator”的“投影(Projection)”和名为“GCS_WGS_1984”的“地理坐标系(GeographicCoordinate System)”。


WGS_1984_UTM_Zone_50N

WKID:32650 Authority: EPSG

 

Projection:Transverse_Mercator

False_Easting:500000.0

False_Northing:0.0

Central_Meridian:117.0

Scale_Factor:0.9996

Latitude_Of_Origin:0.0

LinearUnit: Meter (1.0)

 

GeographicCoordinate System: GCS_WGS_1984

AngularUnit: Degree (0.0174532925199433)

PrimeMeridian: Greenwich (0.0)

Datum: D_WGS_1984

Spheroid:WGS_1984

       Semimajor Axis: 6378137.0

Semiminor Axis: 6356752.314245179

Inverse Flattening:298.257223563


    地理坐标系由三个参数来定义:角度单位(Angular Unit)、本初子午线(Prime Meridian)和大地测量系统(Datum)。地理坐标系“GCS_WGS_1984”使用的角度单位为“度(Degree)”,0.0174532925199433这个数字等于“π/180”,使用的本初子午线为0.0度经线,即格林威治皇家天文台(Greenwich)所在位置的经线,使用的大地测量系统则为“D_WGS_1984”。


地理坐标系的最重要的参数是“大地测量系统(Datum)”,而大地测量系统的最重要的参数是“椭球(Spheroid)”。椭球相同,大地测量系统不一定相同,因为原点(origin)和方位(orientation)可以不同。想象一下,同一个椭球,首先可以固定在三维空间中的任意一个点,并且在固定于某点后还能以三个自由度任意地旋转其方位(朝向)。当然,具体国家或地区在选择大地测量系统时,总是选择与这一国家或地区的地面最吻合的大地测量系统,而不是拍脑袋随便选的。我们拿到的境内的许多数据使用的都是“D_Xian_1980”大地测量系统,因为“D_Xian_1980”是我们依据我国疆域的地面自己定义出来的,因而较“D_WGS_1984”与我国疆域的地面更吻合。“D_WGS_1984”大地测量系统使用的椭球为“WGS_1984”,而“WGS_1984”椭球的“长半轴(Semimajor Axis)”和“短半轴(Semiminor Axis)”分别为6378137.0和6356752.314245179,其“反扁率(Inverse Flattening)”为298.257223563,等于Semimajor Axis/( Semimajor Axis - Semiminor Axis)。


投影的参数对不同的投影方法有一定差别,在此也不详述各投影的具体参数。投影坐标系“WGS_1984_UTM_Zone_50N”使用的“投影(Projection)”名为“横轴墨卡托(Transverse_Mercator)”,然而这个名称并不能完全准确概括其投影。事实上,投影坐标系“WGS_1984_UTM_Zone_50N”这个名称中的“WGS_1984”指出了其地理坐标系为“GCS_WGS_1984”,而“UTM_Zone_50N”则指出了其投影。“UTM_Zone_50N”这个名称指出,其投影方法是“通用横轴墨卡托(Universal Transverse Mercator,UTM)”,其投影带为北半球第50带,这个“Zone_50N”的“中央经线(Central Meridian)”正是117.0度,在“Transverse_Mercator”的参数中得到了体现。举一反三,“Xian_1980_GK_CM_117E”这个坐标系使用的地理坐标系为“GCS_Xian_1980”,而投影名称“GK_CM_117E”指出其使用以东经117度为中央经线的“高斯-克吕格(Gauss-Kruger,GK)”投影。投影的另一个重要参数是“东偏(False Easting)”。有些投影会在X坐标值前加上投影带号,比如:“Xian_1980_GK_Zone_20”的“false_easting”参数为20500000.0,其中20为投影带号,而“Xian_1980_GK_CM_117E”的“false_easting”参数为500000.0,尽管它们的中央经线都为东经117度。


(二)三个半概念

在ArcGIS中,有三个概念容易混淆(另外半个最后揭晓),需要特别进行区分:

1)数据的真实坐标系,简称为“真实坐标系”;

2)数据属性所标称的坐标系,简称为“属性坐标系”;

3)ArcMap/ArcScene中Layers的坐标系,简称为“地图坐标系”。


数据的真实坐标系是指数据记录本身所对应的坐标系。比如,国科大雁栖湖校区图书馆在地理坐标系“GCS_WGS_1984”下的经纬坐标大概为(116.679267°E,40.408265°N),在投影坐标系“WGS_1984_UTM_Zone_50N”下的平面坐标为(472786.066803m,4473121.59882m)。假如用一个Point Shapefile数据来记录国科大雁栖湖校区图书馆的位置。如果使用经纬坐标(116.679267,40.408265)来记录此位置,那么数据的真实坐标系就是地理坐标系GCS_WGS_1984,如果使用平面坐标(472786.066803,4473121.59882)来记录此位置,那么数据的真实坐标系就是投影坐标系WGS_1984_UTM_Zone_50N。


数据的真实坐标系是什么,可以通过以下方式进行验证。在ArcMap中加载这一个Point Shapefile数据,打开其属性表(Attribute Table),右击任意已有或新建的数值类型为Double(或Float)的Filed进行“Calculate Geometry”运算,如果数据的真实坐标系为地理坐标系GCS_WGS_1984,那么计算的“X Coordinate of Point”就为116.679267(图2),如果数据的真实坐标系为投影坐标系WGS_1984_UTM_Zone_50N,那么X坐标值就为472786.066803。


图2 使用“Calculate Geometry”检查Point Shapefile 数据的“真实坐标系”


对于栅格数据,也有方法判断其真实坐标系。通常,栅格分辨率数值很小的是地理坐标系,栅格分辨率的数值很大则为投影坐标系。以经常使用的SRTM DEM栅格数据为例。从网上下载的SRTM DEM数据是地理坐标系,其分辨率数值为0.00083333333(图3),表示其每个栅格的长宽都为0.00083333333度。而0.00083333333度这个距离在国科大的纬度位置大致相当于地面距离81.4520173米。因此,如果使用投影坐标系,让81.4520173这个数值作为这个纬度位置的SRTM DEM数据的分辨率是合理的选择。81.4520173是0.00083333333的接近10万倍(图3),因此分辨率(Cell Size)是判断栅格数据的坐标系是地理坐标系还是投影坐标系的重要依据。当然,大尺度的全球栅格数据在地理坐标系下其分辨率数值可以很大,比如1度,而小区域的栅格数据在投影坐标系其分辨率数值也可以很小,比如1米。


图3 同一个栅格数据在“地理坐标系”和“投影坐标系”下的分辨率数值相差巨大


数据属性所标称的坐标系是指数据文件的属性所标称的坐标系。一般来说,我们有两个入口来查看数据属性所标称的坐标系:一种是在Catalog里面右击该数据文件打开Properties,Shapefile文件可直接在Properties中查看或者修改“XY Coordinate System”,Raster文件可以通过Edit(编辑)“SpatialReference”来打开“XY Coordinate System”进行查看或者修改;另一入口是当把数据加载进ArcMap或ArcScene后,在此数据的Layer Properties的“Source”标签中查看。第二个入口只能查看而不能修改属性坐标系。图4为分别在Shapefile Properties和Layer Properties中查看国科大雁栖湖校区图书馆PointShapefile数据的属性坐标系。


图4 在Shapefile Properties 和Layer Properties 中查看Point Shapefile 数据的“属性坐标系”


这里要特别强调的是:数据的真实坐标系和属性坐标系可以不同,当二者不同时就出现错误。例如:数据的真实坐标系为地理坐标系GCS_WGS_1984,而属性坐标系为投影坐标系WGS_1984_UTM_Zone_50N,或真实坐标系为投影坐标系WGS_1984_UTM_Zone_50N,而属性坐标系为地理坐标系GCS_WGS_1984,都是错误的。类似的,数据的真实坐标系为地理坐标系GCS_WGS_1984,而属性坐标系为地理坐标系GCS_Xian_1980,或真实坐标系为投影坐标系WGS_1984_UTM_Zone_50N,而属性坐标系为投影坐标系WGS_1984_UTM_Zone_49N,也都是错误的。只要二者不统一,就是错误的。在数据处理过程中,误将属性坐标系改动,造成属性坐标系与真实坐标系不符合,是ArcGIS操作的常见错误。使数据的属性坐标系和真实坐标系吻合,是进行所有数据处理和分析的必要前提。


同样以国科大雁栖湖校区图书馆的Point Shapefile数据为示例说明真实坐标系和属性坐标系不吻合的情况。如果Point Shapefile数据的真实坐标系为地理坐标系GCS_WGS_1984,而将属性坐标系设置为投影坐标系WGS_1984_UTM_Zone_50N,图书馆点就会跑到赤道附近的印度尼西亚去(图5)。此时,系统会认为此点在投影坐标系WGS_1984_UTM_Zone_50N中坐标值为(116.679267,40.408265),而在投影坐标系WGS_1984_UTM_Zone_50N中正确的坐标值应为(472786.066803,4473121.59882)。相反,如果Point Shapefile数据的真实坐标系为投影坐标系WGS_1984_UTM_Zone_50N,而属性坐标系为地理坐标系GCS_WGS_1984,图书馆点则会“跑出地球”。因为系统会认为此点的经纬度为(472786.066803,4473121.59882),地球上的点的经纬度值最大为180度,最小为-180度,因而这两个经纬度值远远超过了180度的最大值。


图5 数据的“属性坐标系”与“真实坐标系”不吻合导致的点位错误


为何在实际操作中我们一般都不必检查数据的属性坐标系是否与真实坐标系吻合呢?这是因为在大多数时候,我们拿到的数据不但有属性坐标系并且还与真实坐标系是吻合的。一般来说,有两个明显的迹象可以判断数据的属性坐标系是不对的。一个是将数据加载到ArcMap里面后报错;另外一个是数据加载到ArcMap里面后尽管没报错,但是位置明显不对。比如上面的例子,位于国科大的点跑到印度尼西亚去了。图6是将真实坐标系为投影坐标系“Asia_Lambert_Conformal_Conic”的数据的属性坐标系改为地理坐标系“GCS_WGS_1984”并加载到ArcMap里面后的报错信息。这个信息是说“数据的范围(extent)与其空间参考的信息不吻合”。具体来讲,投影坐标系“Asia_Lambert_Conformal_Conic”中的坐标值有负值,而负值超出了地理坐标系“GCS_WGS_1984”的范围(正值)。以下叙述,如无特别说明,都隐含数据属性坐标系与真实坐标系吻合的假设,并以“数据坐标系”统称。


图6 数据的“属性坐标系”与“真实坐标系”不吻合导致的“inconsistent extent”错误信息


ArcMap/ArcScene中Layers的坐标系是指当我们把若干矢量数据或栅格数据加载进ArcMap或ArcScene里面组成“Layers”的时候,这个Layers的坐标系。“地图坐标系”在ArcMap中可以通过右击“Layers”打开“Data Frame Properties”后在“Coordinate System”标签中查看或者修改(图7),在ArcScene中可以通过右击“Scene Layers”打开“Scene Properties”后在“Coordinate System”标签中查看或者修改。



图7 在ArcMap 中查看或者修改Layers 的“地图坐标系”


可以想象,一个Layers可以包括多个数据Layer,这些数据Layer的坐标系也可以各不相同。因此,Layers的地图坐标系可以不同于各Layer的数据坐标系。可以选择Layers中的某一Layer的数据坐标系作为地图坐标系,也可以选择其他任意坐标系,只要这个坐标系的覆盖范围能覆盖所有Layer数据的范围。当在ArcMap中加载多个数据Layer时,系统会自动将第一个加载进来的有坐标系的数据的坐标系作为Layers的地图坐标系。为了使所有具有不同坐标系的数据Layer都在同一地图坐标系下进行显示等操作,当某个Layer的数据坐标系与Layers的地图坐标系不同时,系统会自动用一定的算法将数据坐标系(的坐标值)临时转换为地图坐标系(的坐标值)。这种临时的坐标系转换,并不改变每个数据本身的坐标系。


如图8所示:“library”、“library_utm”、“国科大20.tif”和“World Physical Map”等四个Layer的数据坐标系分别为:地理坐标系“GCS_WGS_1984”、投影坐标系“WGS_1984_UTM_Zone_50N”、投影坐标系“WGS_1984_Web_Mercator”和投影坐标系“WGS_1984_Web_Mercator_Auxiliary_Sphere”;Layers的地图坐标系则使用了“国科大20.tif”这一栅格Layer的数据坐标系(投影坐标系“WGS_1984_Web_Mercator”)。


当数据在与其坐标系不同的地图中显示时,会出现“变形”。这种变形是由“临时的坐标系转换”引起的。图9为使用地理坐标系“GCS_WGS_1984”作为Layers的地图坐标系的情况。可以看出,相较使用“国科大20.tif”的数据坐标系(投影坐标系“WGS_1984_Web_Mercator”)作为地图坐标系,使用地理坐标系“GCS_WGS_1984”作为地图坐标系使得“国科大20.tif”这一栅格Layer的变形明显,其轮廓(覆盖范围)在东西经度方向上被拉长了,在南北维度方向则被压缩了。如果不希望某个Layer出现变形,那么就使用该Layer的数据坐标系作为Layers的地图坐标系。


                                         图8 使用投影坐标系“WGS_1984_Web_Mercator”作为Layers 的“地图坐标系”


图9 使用地理坐标系“GCS_WGS_1984”作为Layers 的“地图坐标系”造成显示“变形”


如果地图坐标系与数据坐标系差别很大,数据显示的变形也会很大。图10为使用投影坐标系“WGS_1984_UTM_Zone_49N”作为地图坐标系的情况,图11为将“国科大20.tif”Layer的右上角放大后的情况。可以看出,“国科大20.tif”这一栅格Layer的变形明显,整体轮廓和每个栅格都“倾斜”了一个角度。尽管投影坐标系“WGS_1984_UTM_Zone_49N”和 “WGS_1984_UTM_Zone_50N”二者都与“国科大20.tif”这一栅格Layer的数据坐标系“WGS_1984_Web_Mercator”不同。但是,试验可以发现,使用“WGS_1984_UTM_Zone_50N”作为地图坐标系的时候,“国科大20.tif”的变形不会这么明显。原因之一是,坐标系“WGS_1984_UTM_Zone_49N”中投影带的中央经线为111度,而坐标系“WGS_1984_UTM_Zone_50N”中投影带的中央经线为117度,后者的投影带与“国科大20.tif”的范围更符合。选择“WGS_1984_Web_Mercator”和“WGS_1984_UTM_Zone_50N”,都会使“国科大20.tif”在投影中的变形更小。


                             图10 使用投影坐标系“WGS_1984_UTM_Zone_49N”作为“地图坐标系”造成栅格Layer 轮廓“倾斜”


图11 使用投影坐标系“WGS_1984_UTM_Zone_49N”作为“地图坐标系”造成每个栅格“倾斜”变形


当然,有时在地图显示中,出于合理或者美观的需要,这种“变形”是需要的。例如,我们拿到的行政边界矢量数据的坐标系是地理坐标系“GCS_Krasovsky_1940”。我们可以就使用地理坐标系“GCS_Krasovsky_1940”作为Layers的地图坐标系(图12),但是这样的显示和图9类似地使数据在东西方向被“拉长”。使用地理坐标系作为地图坐标系,显示出来的数据与我们在多数时候看到的全国地图的“形状”有差异。当我们换为使用投影坐标系“Asia_Lambert_Conformal_Conic”作为Layers的地图坐标系时,显示出来的数据就与我们在多数时候看到的全国地图的“形状”相仿了(图13)。如果数据经常需要在不同于自身坐标系的地图中显示,可使用“坐标系转换”让这种在显示中的临时性坐标系转换变成永久性的,也即改变数据本身的坐标系,这是下一节中将提到的。对于栅格数据,如图10和图11这种情况,坐标系转换后的数据在目标坐标系中显示时轮廓可能还是“倾斜”的,但数据坐标系变成目标坐标系后的新数据的每个栅格在目标坐标系中将不再“倾斜”(见3.3节)。


                                       图12 使用地理坐标系“GCS_Krasovsky_1940”作为“地图坐标系”显示行政边界数据


                                

图13 使用投影坐标系“Asia_Lambert_Conformal_Conic”作为“地图坐标系”显示行政边界数据


地图坐标系的重要性还体现在:在地图中编辑要素时,如果需要输入坐标值,那么此坐标值是相对于地图坐标系的。例如,在坐标系为“Asia_Lambert_Conformal_Conic”的地图中编辑坐标系为“GCS_WGS_1984”的Point Shapefile,增加点所需输入的坐标值就不能再是该点的经纬度,而是该经纬度在坐标系“Asia_Lambert_Conformal_Conic”中对应的坐标值(图14)。


也许你已经注意到,图13的地图坐标系是投影坐标系,但窗口右下角显示的坐标单位是度而不是米。类似的,图9的地图坐标系是地理坐标系,但窗口右下角显示的坐标单位是米而不是度。这个疑问就引出了那“半个概念”:“地图显示单位”。之所以称为“半个”是因为严格来说这不算一个可以与真实坐标系、属性坐标系和地图坐标系等并列的“概念”。地图显示单位可以在Layers的Properties中的“General”标签中进行设置(图15)。地图单位由地图坐标系决定,地图显示单位可以与地图单位相同,也可以与地图单位不同。


图14 在地图坐标系为“Asia_Lambert_Conformal_Conic”的地图中编辑数据坐标系不同的点数据


图15 使用与“地图单位”相同的单位(Meters)作为“地图显示单位”


本节小结:只要能做到以下“三个重要的区分”,就理解了以上“三个半概念”:

1)区分数据的“真实坐标系”和“属性坐标系”;

2)区分“数据坐标系”和“地图坐标系”;

3)区分“地图单位”和“地图显示单位”。


(三)四类常用操作

在理解清楚了以上“三个半容易混淆的概念”的基础上,以下四类常用的坐标系操作就很好理解与熟悉了。
3.1坐标系选择、新建与编辑坐标系的其他操作都会涉及到选择、新建与编辑等,这些操作可在打开的数据或地图等的“Properties”窗口中的“Coordinate System”或“XY Coordinate System”标签中进行(图16)。

图16 坐标系选择、新建与编辑等操作的入口


用户可以从“Favorites”、“Geographic Coordinate System”、“ProjectedCoordinate System”和“Layers”等文件夹中选择坐标系,或使用“Import”导入其他数据的坐标系。“Layers”文件夹只有当ArcMap或ArcScene中的Layers加载有数据(Layer)时才会显示出来,而在ArcCatalog中操作则不会出现。另外,使用“Clear”可以清除数据或地图的坐标系。

用户也可以根据需要新建坐标系。图17为新建投影坐标系(Projected Coordinate System)的操作窗口。在这一操作窗口中,用户需要设置投影坐标系的名称(Name)、投影(Projection)、单位(Linear Unit)和地理坐标系(Geographic CoordinateSystem)等参数。


图17 新建投影坐标系(Projected Coordinate System)的操作窗口


右击任意已有坐标系,点击“Copy and Modify”,可对已有坐标系进行编辑。坐标系编辑的操作窗口与坐标系新建的操作窗口类似。对于投影坐标系的编辑,用户也需要编辑投影坐标系的名称、投影、单位和地理坐标系等参数。
3.2坐标系定义

坐标系定义是指定义数据的属性坐标系。也即,将与数据的真实坐标系相同的坐标系赋予给数据的属性。坐标系定义可以通过两种方式来实现:一是如图2所示,在Catalog里面访问该数据的Properties进行修改(定义);二是使用ArcToolbox -> Projections and Transformations -> DefineProject工具(图18)。如果该数据已有属性坐标系,数据输入处会出现警告符号,警告用户已有属性坐标系将被覆盖。这两种坐标系定义方式的区别是:如果数据在ArcMap或ArcScene中打开(被锁定)后,那么就不能在Catalog中改变数据的坐标系了,但是仍然可以利用ArcToolbox中的Define Project工具来定义数据的坐标系,除非这一数据同时也在另外的ArcMap、ArcScene或ArcCatalog中被打开使用。


图18 使用ArcToolbox 中“Define Project”工具进行坐标系定义


另外,个人以为,“Define Project”工具的名称使用“Define Coordinate System”更为贴切。因为这个工具是定义坐标系(Coordinate System),而不是定义投影(Project)。有意思的是,“Define Project”工具窗口中“坐标系的输入框”的名称是“CoordinateSystem”,而不是“Project”(图18)。因此,为何ESRI会使用“Define Project”这个工具名,颇令人费解。 3.3坐标系转换

坐标系转换是指转换数据的真实坐标系。坐标系转换的前提是数据的属性坐标系与数据的真实坐标系吻合。矢量数据的坐标系转换通过ArcToolbox -> Projections and Transformations -> Project工具来实现(图19),栅格数据的坐标系转换通过ArcToolbox -> Projectionsand Transformations -> Raster -> Project Raster工具来实现(图20)。(注:不同ArcGIS版本工具在ArcToolbox中的位置可能略有不同,此处为ArcGIS 10.3中的路径)。


                                                     图19 使用ArcToolbox 中“Project”工具进行矢量数据的坐标系转换


                                               图20 使用ArcToolbox 中“Project Raster”工具进行栅格数据的坐标系转换

在选择好输入数据,并定义好输出数据及其坐标系后,用户有可能被要求定义“Geographic Transformation”。这个“Geographic Transformation”被用于两个不同地理坐标系(大地测量系统)之间的转换。只有当输入坐标系和输出坐标系拥有不同的大地测量系统时,这个输入才被要求。在某些情况下,系统会从已有的“Geographic Transformation”中根据输入和输出坐标系自动选择出合理的“Geographic Transformation”(图20),而某些情况下则需要用户从已有的“Geographic Transformation”中进行选择(图19)。
坐标系转换理论上可以在任意两个覆盖了数据范围的坐标系之间进行。坐标系转换可以在地理坐标系与投影坐标系之间,可以在地理坐标系与地理坐标系之间,也可以在投影坐标系与投影坐标系之间。
坐标系转换包括了两种过程:1)大地测量系统(地理坐标系)转换;2)投影(或反投影)。例如,将地理坐标系“GCS_WGS_1984”转换为投影坐标系“Xian_1980_GK_CM_117E”包括了两个过程:分别是一个将大地测量系统“D_WGS_1984”转换为大地测量系统“D_Xian_1980”的过程和一个将地理坐标系“GCS_Xian_1980”投影为投影坐标系“Xian_1980_GK_CM_117E”的过程。再例如,将投影坐标系“WGS_1984_Web_Mercator”转换为投影坐标系“WGS_1984_UTM_Zone_50N”包括了三个过程:分别是一个将投影坐标系“WGS_1984_Web_Mercator”转为地理坐标系“GCS_WGS_1984_Major_Auxiliary_Sphere”的过程(反投影),一个将大地测量系统“D_WGS_1984_Major_Auxiliary_Sphere”转换为大地测量系统“GCS_WGS_1984”的过程和一个将地理坐标系“GCS_WGS_1984”投影为投影坐标系“WGS_1984_UTM_Zone_50N”的过程。

栅格数据的坐标系转换与矢量数据的坐标系转换的一个重要不同是:栅格数据的坐标系转换需要设置栅格分辨率(图21)。如果是投影坐标系之间或地理坐标系之间的转换,系统一般会默认分辨率数值不变,如果是地理坐标系与投影坐标系之间的转换,系统会自动估计出一个合理的分辨率数值。例如第2节中提到的81.4520173米对应于0.00083333333度,也即地理坐标系中的0.00083333333度对应于投影坐标系中的81.4520173米。由于这里的转换是对数据的真实坐标系的永久性转换,和在地图显示中的“临时性坐标转换”不同,栅格数据的坐标系转换不但使栅格数据的分辨率数值改变,也使新数据在目标坐标系中显示时其每个栅格不再“倾斜”,尽管其整体轮廓可能还是“倾斜”的(图21)。


图21 使用“Project Raster”工具进行栅格数据的坐标系转换时需要设置栅格分辨率


同样,个人以为,“Project”和“ProjectRaster”的工具名用“Transfer ”替代“Project”也许更为贴切。因为这两个工具做的工作是坐标系转换(Transfer),而不是投影(Project)。另外,ESRI将用于栅格数据坐标系转换的工具命名为“Project Raster”,却不将用于矢量数据坐标系转换的工具命名为类似对应的“Project Feature”,不知有何考虑。
因此,ArcGIS中坐标系定义和转换的工具名称(包括Define Project、Project和Project Raster)可能对用户使用工具和理解坐标系的有关概念有一定误导。
3.4坐标系猜测坐标系猜测是指猜测数据的真实坐标系,并且把数据的属性坐标系设置为真实坐标系。这在不知道数据的坐标系或数据的属性坐标系是错误的时候是需要的。比如,给你若干点的坐标值,但是不告诉你这些坐标值是在什么坐标系下的坐标值。在这种情况下,最好的方式是向数据的提供方询问数据的真实坐标系是什么。 如果不得不猜测数据的真实坐标系,其基本原理是根据数据的坐标值特征来猜测。经验将会很重要。如果数据的坐标值看起来像是经纬度数值,那么认为其真实坐标系为地理坐标系“GCS_WGS_1984”将不会带来很大误差。如果数据的坐标值是地面长度(通常数值很大),那么可以基本认为其真实坐标系为投影坐标系。如果知道该数据大概的空间范围,可以找出该空间范围内的坐标系正确的任一“参考数据”,将其与属性坐标系未知的“目标数据”加载于同一地图中,并将地图坐标系设置为你怀疑的投影坐标系。当目标数据与参考数据范围吻合时,可以猜测此怀疑的投影坐标系即为目标数据的真实坐标系,并可辅以其他手段进行验证。怀疑的投影坐标系可以有很多个,选择哪些投影坐标系作为怀疑对象,需要基于经验根据数据的坐标值特征进行判断。比如,加投影带号与不加投影带号的投影坐标系中的数据X坐标值的大小在很多时候是有明显差异的。(四)总结

在数据分析之前,进行以下有关坐标系的准备工作是一个良好的习惯:1)使所有数据都有坐标系定义,也即使所有数据都有属性坐标系;2)保证所有数据的坐标系都是正确的,也即使所有数据的属性坐标系都与其真实坐标系吻合;3)统一所有数据的坐标系,利用坐标系转换实现。并且,这种统一的坐标系一般宜为投影坐标系,因为只有在投影坐标系下才能正确地进行距离、长度、面积、坡度等度量的计算。当然,有时我们做大尺度的工作,比如进行全球尺度的显示或分析时,一般就用地理坐标系,不用投影坐标系。

数据处理(二):wgs84坐标的UTM投影 前言 wgs84 坐标常用于我们身边的各种gps模块。从前面的文章(点这里查看文章)我们可以从gps模块中得到wgs坐标下的经纬度值,然而经纬度值有时候并不是以距离为单位的坐标,有时候更多的是用到以距离为单位的二维坐标,所以就需要从wgs84(三维坐标)投影转换为二维坐标。 python有一个模块叫做pyproj可以帮助我们完成坐标的转换 安装pip install pyproj 1 代码 imp... 阅读详情

相关推荐

地图(Gis)中的坐标系

地图(Gis)中的坐标系 文章目录地图(Gis)中的坐标系前言1.什么是坐标系(地图投影)1.1坐标系介绍1.2地理坐标系1.3投影坐标系1.4为什么要有这么多坐标系1.5什么是EPSG1.6不同坐标系的数据要如何转换2.坐标系在项目中的应用2.1地图框架和坐标系的确定2.2 Gis数据坐标系2.3 其他可能遇到的问题3.参考文献 前言 最近暂时离开的三维Gis的开发,转身投入到了二维Gis的开发,研究二维地图最好的一点终于不用考虑模型建模,贴图,着色器,粒子特效等三维方面甚至是计算机图形学的东西了,可

weixin_42066016的博客 4011

地理坐标系与投影坐标系

地理坐标系和投影坐标系

AIGIS的博客 1万+

ArcGIS中地理坐标系与投影坐标系以及投影转换

简单介绍arcgis软件中坐标系定义、设置和常见问题。

weixin_41399650的博客 8万+

常用坐标系汇总(更新)

前言     随着接触的地图种类越来越多,每种产品对地图服务的坐标系的要求不同,今天遇到了整理的好文,整理记录分享。 投影坐标系:墨卡托坐标系地理坐标系:经纬度坐标系 常用坐标系 (一)WGS84坐标系     WGS-84坐标系(World G...

Learn_Fighting的博客 5593

arcgis中dem坐标定义_重新认识ArcGIS中的坐标系

重新认识ArcGIS中的坐标系ArcGIS中的坐标系:基本概念和常用操作本文转载仅供学习分享,如有侵权联系删除:原文链接:http://blog.sciencenet.cn/blog-290812-1016263.html李郎平,Email: lilp@lreis.ac.cn(中国科学院地理科学与资源研究所,资源与环境信息系统国家重点实验室)缘由:介绍GIS(地理信息系统)中的坐标系的想...

weixin_35975197的博客 1178

ArcGIS坐标系实战:从基础概念到投影变换全解析

本文全面解析ArcGIS坐标系的核心概念与实战操作。从地理坐标系与投影坐标系的基础差异讲起,深入剖析基准面、椭球体以及中国广泛使用的高斯投影及其分带原理。重点介绍了在ArcGIS中进行坐标系定义、动态投影以及投影变换的实战步骤与常见错误排查方法,帮助用户夯实GIS数据空间分析的基石。

weixin_29284201的博客 504

arcgis 投影坐标系转地理坐标系_解决ArcGIS坐标之惑:地理坐标系和投影坐标区别及常用操作...

点击标题下方水土保持之点滴,关注后可查阅所有文章两种坐标系坐标系(Coordinate System)的概念为:“In geometry, a coordinatesystem is a system which uses one or more numbers, or coordinates, to uniquelydetermine the position of a point o...

weixin_39525097的博客 2408

从Excel坐标到GIS地图:一次讲清ArcGIS中地理坐标系与投影的那些事儿(以WGS84生成面数据为例)

本文深入解析ArcGIS中地理坐标系与投影坐标系的核心区别,重点讲解WGS84作为全球标准坐标系的优势及应用场景。通过实战案例演示如何将Excel经纬度坐标转换为GIS面数据,并揭示坐标系选择对面积计算和空间分析的关键影响,帮助用户避免常见陷阱,提升地理数据处理精度。

weixin_30363509的博客 635

坐标的_解决ArcGIS坐标之惑:地理坐标系和投影坐标区别及常用操作

点击标题下方水土保持之点滴,关注后可查阅所有文章两种坐标系坐标系(Coordinate System)的概念为:“In geometry, a coordinatesystem is a system which uses one or more numbers, or coordinates, to uniquelydetermine the position of a point o...

weixin_30568317的博客 2179

别再乱用DefineProjection工具了!ArcGIS中定义、投影、转换坐标系的区别与正确操作指南

本文详细解析ArcGIS中定义坐标系、投影转换和坐标转换的核心区别,重点纠正DefineProjection工具的常见误用。通过三层架构解析、操作场景说明和实战案例,帮助GIS工作者避免数据偏移问题,提升国土调查、环境监测等项目的坐标精度与工作效率。

weixin_30872337的博客 453

ArcGIS坐标系实战:数据框、图层属性中的投影与地理坐标系到底怎么选?

本文针对ArcGIS用户在数据框和图层属性中如何选择坐标系的核心难题,深入解析了地理坐标系与投影坐标系的本质区别与应用场景。文章通过实战案例,清晰阐述了数据框坐标系作为“显示舞台”的动态投影原理,以及图层坐标系作为数据“身份证”的定义与投影操作,旨在帮助用户避免常见陷阱,实现精确测量与高效制图。

weixin_29228781的博客 82

ArcGIS中的坐标系

今天将解决以下问题: 1.分清地理坐标系和投影坐标系 2.什么地区选什么带号 3.如何使用有关“坐标系”的三个工具 ArcGIS 无疑是应用最广的GIS 软件,上手的第一个“拦路虎”就是坐标系,很多人对坐标系的概念理解方面和实际操作方面不太清楚,对于有些学规划的小伙伴简直就是噩梦啊,会出现cad文件导入ArcGIS里面咋就没坐标了呢?图层要素为啥不在一个平面上?等问题。 所以呢,我就以ArcGIS为工具介绍坐标系的基本概念以及常用操作作为开始(大佬请绕行)。 两种坐标系 在ArcGIS中,

潘达小新 1万+

【数据库】银行/信贷业务测试工作中涉及的数据库表格

银行/信贷业务测试工作中涉及的数据库表格,金融信贷测试人必看!

孤之鹜的博客 1169

ArcGIS中的坐标系统定义与投影转换【转】

ArcGIS中的坐标系统定义与投影转换 坐标系统是GIS数据重要的数学基础,用于表示地理要素、图像和观测结果的参照系统,坐标系统的定义能够保证地理数据在软件中正确的显示其位置、方向和距离,缺少坐标系统的GIS数据是不完善的,因此在ArcGIS软件中正确的定义坐标系统以及进行投影转换的操作非常重要。 1.ArcGIS中的坐标系统 ArcGIS中预定义了两套坐标系统,地理坐标系(G...

weixin_30569033的博客 1万+

例题讲解拉格朗日乘子法、线性可分支持向量机(SVM)的推导

支持向量机(Support Vector Machine, SVM)的简单推导及例题

zzzssszzzzsz的博客 5681

别再傻傻分不清了!用Design Compiler和PrimeTime实操对比PBA与GBA模式

本文深度解析静态时序分析(STA)中Design Compiler与PrimeTime工具的PBA与GBA模式差异,通过实际案例演示两种模式在转换时间计算、时钟路径分析等方面的对比,并提供工程实践中模式选择的策略与调试技巧,帮助工程师在芯片设计流程中平衡精度与效率。

weixin_42552315的博客 414

arcgis根据7参转坐标_在ArcGIS Desktop中进行三参数或七参数精确投影转换

假设原投影坐标系统为Xian80坐标系统,本例选择为系统预设的Projected Coordinate Systems\Gauss Kruger\Xian 1980\Xian 1980 GK Zone 20投影,中央经线为117度,要转换成Beijing 1954\Beijing 1954 GK Zone 20N。在ArcMap中加载了图层之后,打开View-Data Frame Properti...

weixin_40001245的博客 1505

三相逆变器LCL滤波设计实战:从αβ转换到dq解耦控制的完整建模流程

本文详细解析了三相逆变器LCL滤波设计中的关键步骤,从αβ转换到dq解耦控制的完整建模流程。重点介绍了dq坐标系下的解耦控制策略,通过Simulink仿真验证模型准确性,并提供谐振峰偏移、解耦不彻底等实际问题的解决方案,助力电力电子工程师提升系统稳定性。

weixin_28678151的博客 229

WGS84坐标系

WGS84与WGS84 Web Mercator 1. WGS84与WGS84 Web Mercator 1.1 关于WGS1984投影坐标系 UTM (Universal Transverse Mercator)坐标系是由美国军方在1947提出的。虽然我们仍然将其看作与“高斯-克吕格”相似的坐标系统,但实际上UTM采用了网格的分带(或分块)。除在美国...

hhhhhhhhh的博客 4万+

一维光子晶体态密度计算与应用解析

态密度(DOS)是凝聚态物理和光子学中的核心概念,描述了量子态在能量空间的分布规律。基于麦克斯韦方程组的周期性介电结构分析表明,一维光子晶体通过布拉格散射形成光子带隙,其态密度曲线呈现典型的范霍夫奇点和禁带特征。这种特性使光子晶体在光学工程中具有重要价值,特别是在滤波器设计和传感器开发领域。通过转移矩阵法等数值计算手段,可以精确预测带隙位置和宽度,其中折射率对比度与周期数是关键设计参数。实际应用中,1550nm通信波段滤波器常采用Si/SiO2交替结构,而引入缺陷层则能产生高Q值缺陷模。计算优化时需注意数值

weixin_36487018的博客 193

LC-3 机器语言 指令集

生活本就沉闷,但跑起来就有风

笔记 1万+

gis投影中未定义的地理转换_ArcGIS中5分钟搞懂坐标系相关知识

原标题:ArcGIS中5分钟搞懂坐标系相关知识坐标系定义投影说明1、坐标系基础知识坐标系分为地理坐标系与投影坐标系。 1.1、地理坐标系常见地理坐标系:国家2000(CGCS2000);西安80(Xian_1980).地理坐标系坐标值为经纬度格式,如下:118 13.862 34 58.789度分1.2、投影坐标系常见投影坐标系:高斯-克吕格投影(Gauss_Kruger);UTM投影.投影坐标系...

weixin_39766109的博客 2074

【WorkBuddy专栏58】腾讯为什么要「自己打自己」——7款AI智能体产品全景对比与「赛马」逻辑深度拆解

【WorkBuddy专栏57】你的左侧空间不是「日任务清单」——工作空间才是 WB 变聪明的核心秘密

AI人人懂 670

GIS基础-坐标系01

GIS基础教程之坐标系 1、坐标系基本概念(下文1,2,3节) 2、如何选择投影坐标系(第4节) 3、各种坐标系之间的转换(第5节) 4、经纬度格式转换方法(第6节) 01 坐标系的分类 首先介绍坐标系的类别,坐标系分为两大类: 1、地理坐标系(Geographic Coordinate System, GCS) 2、投影坐标系(Projected Coordinate System, PCS) 地图投影——理解为拨开一个橘子皮,从球面变换到平面。为了便于计算啥的 简单来说,把地图投影的过程就是一个剥橘

qq_32230783的博客 5165

14、利用Helm和GitOps实现Kubernetes资源管理与CI/CD流程优化

本文介绍了如何结合Helm、GitOps和Terraform等工具,优化Kubernetes资源管理和CI/CD流程。内容涵盖了Helm提供程序的配置、GitOps的核心理念、CI/CD管道的构建与测试策略、GitHub源代码管理模型以及流程优化建议,旨在提升软件开发和部署的效率与质量。

y7z8a9的博客 154

文档:亲测可用解决vnc灰屏,黑屏,鼠标大黑叉情况

主要解决vnc-viewer灰屏,黑屏,鼠标大黑叉等现象,详细见文档。

上一篇: Cesium中的Entity API和Primitive API的区别(2)
下一篇: 将excel中的若干列数据生成json格式文件
龙之吻number7
博客等级 码龄10年 99粉丝 7原创
评论 8
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付:10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值