WebGIS系列(一):坐标系

# 1. 坐标系

1. 大地坐标系与地理坐标系有何不同?

        大地坐标系是以地球质心为原点建立的三维坐标系,用于描述地球上任意点的空间位置。而地理坐标系是以经纬度作为坐标值,用于描述地球表面上任意点的位置。

        换句话说,大地坐标系是基于物理概念和数学模型所建立的坐标系,而地理坐标系则是基于地球表面上的地理位置所建立的坐标系。

  • 地理坐标:天文坐标和大地坐标

    ⭐ 摘录自 https://www.jianshu.com/p/8ef5ebf5559a 

    地面点在球面上的位置用经纬度表示时,称为地理坐标。地理坐标按坐标所依据的基准线、基准面及求坐标方法的不同,可分为天文坐标大地坐标两种。

    天文坐标又称天文地理坐标,是用天文经度和天文纬度表示地面点在大地水准面上的位置。

    大地坐标又称大地地理坐标,用大地经度和大地纬度表示地面点在旋转椭球面上的位置。

    天文经纬度是用天文测量方法测定的。大地经纬度是根据一个起始的大地点(大地原点,这点的大地经、纬度与天文经度一致)的大地坐标,按大地测量所得的数据推算而得的。在国家各等级平面控制测量中,实际测定和应用的是大地坐标。

2. 通常所说的西安80,北京54是指什么?

    椭球体基准面之间是一对多的关系(我国的北京54、西安80、国家2000坐标系都使用了不同的椭球体)。大地基准面用于控制参考椭球和地球的相对位置。

    2.1 地心基准面→地心坐标系:由卫星数据得到,使用地球的质心作为原点,使用最广泛的是 WGS 1984,我国的国家2000也是使用地心基准面(地心基准比参心基准更适用于卫星定位系统,如GPS

    2.2 参心基准面→参心坐标系:特定区域内与地球表面吻合,大地原点是参考椭球与大地水准面相切的点,例如Beijing54(局部平差)、Xian80(整体平差)。我们通常说的Beijing54、Xian80坐标系实际上指的是我国的两个大地基准面。当确定参考椭球体和大地基准面后,就可以建立地理坐标系,所以通常称谓的Beijing54、Xian80在没有说明有没有投影及投影信息下,也可以被视为是地理坐标系。

    👉 更严谨的说法是 [来源]:地理坐标系 (GCS) 由基准面Datum、角度测量单位(通常是度)和本初子午线(用于指定 0° 经度的位置)定义。投影坐标系 (PCS) 由线性测量单位(通常是米或英尺)、地图投影、地图投影所使用的特定参数和地理坐标系组成。

  地理坐标系 = 基准面 + 本初子午线 + 角度单位

  投影坐标系 = 地理坐标系 + 投影方法 + 线性单位

  • WGS84和CGCS2000

    CGCS2000(EPSG:4490)与WGS84(EPSG:4326)这两种坐标系都是地心坐标系,采用的参考椭球非常接近。由于扁率差异引起椭球面上的纬度和高度变化最大达 0.1mm,在当前测量精度范围内是可以忽略这点差异的,可以说两者相容至cm级水平,但若是一个点的坐标精度达不到cm水平时,就不认为CGCS2000和WGS84的坐标是相容的,反之两者皆可 [链接]。

    在现有的测绘精度和非专业情况下,可以将CGCS2000坐标系视为等同于WGS84。

    CGCS2000椭球和WGS84椭球极为相似,偏差仅有0.11mm,完全可以兼容使用 [链接]。

  • GCJ-02和BD-09

    谷歌地图OSM地图,都是WGS84坐标系。

    天地图是CGCS2000坐标系。

    高德腾讯是GCJ02坐标系,GCJ02坐标系是在WGS84基础上加密而成的(火星坐标系)。

    百度是BD09坐标系,是在GCJ02上再次加密而成的。

    ① GCJ-02:火星坐标系,国测局02年发布的坐标体系,在WGS84基础上加密而成,高德、腾讯、谷歌地图中国大陆板块等使用。(G表示Guojia国家,C表示Cehui测绘,J表示Ju局  

   如果再以投影为基础作细分,则可以分为GCJ-02 经纬度投影GCJ-02 Web 墨卡托投影

   高德坐标转换Web API(非高德坐标→高德坐标):坐标转换-API文档 | 高德地图API

  📍 WGS84-GCJ02坐标转换的公式由一个关于经纬度的线性多项式(次数从0.5到2)加上经纬度的正弦函数组成。如果都是线性多项式的话,可以很容易推到反函数,但是后面加了一个非线性的函数(正弦函数应该是为了周期性的增加误差用的),这样反函数就非常不容易推导出来。所以关于从GCJ02坐标系转到WGS84只能用迭代法来进行求解了 [来源]

   ② BD-09:百度坐标系,在GCJ-02基础上再次加密,百度地图使用

    百度坐标转换Web API(非百度坐标→百度坐标):坐标转换服务-API文档 | 百度地图API

    BD09LL表示百度经纬度坐标,BD09MC表示百度墨卡托米制坐标。


    WGS84的坐标转化为GCJ02的坐标是单向的,即WGS84的坐标能够准确地变换为GCJ02坐标;但GCJ02坐标转换为WGS84时会存在精度损失。 GCJ02的坐标和BD09的坐标转换是双向的 [来源] [转换代码]

   🧰 工具:QGIS火星坐标转换插件

   🧰 资源:什么是GCJ02火星坐标系   

   🧰 历史:中国地图坐标(GCJ-02)偏移算法

3. 度分秒(DMS)和十进制度(DD)是两种地图单位吗? 

   是,它们是用不同格式表示同一位置经纬度的单位 [转换工具]。

    度分秒表示法中,1度等于60分、1分等于60秒( 1°=60′ ;1′=60″ ),比如经纬度 45°52′30″,那我们转换为十进制来表示就是 45.875°(45+52/60+30/3600)

   地图单位由地图坐标系决定,地图显示单位可以与地图单位相同,也可以与地图单位不同(如:WGS84的地图单位是度,但地图显示单位可以设置成 “度” 或 “米”等)。

   地理坐标系 (GCS) 的单位是角度(angular),通常是度(degree)。投影坐标系 (PCS) 的单位是线性的(linear),最常见的是米(metre)。

  关于地图单位的查阅可以从 EPSG.io: Coordinate Systems Worldwide中的unit参数获取,例如EPSG:4326 Unit: degree ;EPSG:4294 Unit: degree minute second hemisphere;EPSG:3857 Unit: metre;EPSG:3426  Unit: US survey foot...

ArcGIS中地图单位、地图显示单位的设置

4. ArcGIS工具箱中的【定义投影工具】和【投影工具】有什么不同?

     定义投影工具主要用于数据具有未知不正确坐标系的情况。它重新分配元数据中的坐标系信息。它不会改变实际的坐标值。

     投影工具(以及【栅格】-【投影栅格】工具)将所有坐标从一个坐标系转换为另一个坐标系,方法是将数据取消投影为地理坐标。它将原投影的地理坐标转换为新的地理坐标系统,然后将这些地理坐标投影到新的投影坐标系统中。它还会更新元数据信息。

   📌 投影、投影栅格工具若是在两个属于不同基准面的坐标系间转换,需要先进行地理(坐标)变换。(👉 即当要将A坐标系转为B坐标系时,要先判断A、B所属的地理坐标系是否相同,不同则要进行地理变换。

   📌 地理基准面变换是一种用于在两个不同的地理坐标系之间进行变换以确保数据正确对齐的计算 [链接]。

  >>> 同一基准面 / 不同基准面的转换:链接1链接2链接3

  👍 来源:Coordinate Systems: What’s the Difference?

  • 地图投影类型

  👉 摘录自 《地理信息系统导论》Kan-tsung Chang 著

    地图投影可以根据所保留的性质或投影面进行分组。制图者通常根据地图投影所保留的性质将其分成四类:正形(等角)等积等距等方向或等方位

    正形投影保留了局部角度及其形状;等积投影以正确的相对大小显示面积;等距投影保持沿确定路线的比例尺不变;等方位投影保持确定的准确方向。地图投影的名称通常包含它所能保留的性质,例如兰勃特正形圆锥投影或阿伯斯等积圆锥投影。

    正形等积两种性质是相互排斥的,否则一个地图投影所能保留的性质就不止一种,例如会同时保留正形和等方向。正形等积的性质是全局性质,即可应用于整幅地图投影。等距等方位性质是局部性质,只能在距地图投影中心较近的地方实现。

  • 带号

  📍 并不是所有投影都有分带,分带是高斯克吕格投影自带的

  在ArcGIS中,区分有无带号,从投影坐标系统的名称中可以看出:

  • CGCS2000_3_Degree_GK_CM_111E:名称中有CM,在Y坐标(经线方向)不带带号
  • CGCS2000_3_Degree_GK_Zone_30:名称中有Zone,则在Y坐标(经线方向)带带号

  👉 带带号就是在Y坐标值前多加投影带数字,例如,某第38投影带上的不带投影带的坐标值为(494832,2546025),则带投影带的坐标为(38494832,2546025)

  👉 是否带带号和原始数据有关,不带带号的数据是X6位、Y7位,带带号的数据是X8位、Y7位(其中X的前两位表示带号) [链接]。

  👉 我国区域内 六度带带号范围 13-23(共11个),三度带带号范围24-45(共22个)

  👉 在比例尺1:2.5万—1:50万图上采用6°分带,对比例尺为1:1万及大于1:1万的图采用3°分带。

  📌 对于有带号的数据,可以直接根据带号判断是3度分带(≥24)还是6度分带(≤23);对于没有带号的数据,根据比例尺判断是几度分带。

  >>> 通过坐标推断投影带类型

  >>> 判定影像或矢量是否带带号并做投影变换

  • ArcGIS投影命名方式

在ArcGIS中,有四种不同的命名方式。

(1)Beijing 1954 3 Degree GK CM 75E:北京1954坐标系,三度带法,中央经线在东经75°,横坐标前不加带号

(2)Beijing 1954 3 Degree GK Zone 25:北京1954坐标系,三度带法,分带号为25,横坐标前加带号

(3)Beijing 1954 GK Zone 13N:北京1954坐标系,六度带法,分带号为13,横坐标前不加带号

(4)Beijing 1954 GK Zone 13:北京1954坐标系,六度带法,分带号为13,横坐标前加带号

5.  什么是动态投影

    数据框中含有不同坐标体系的数据时,数据会动态的进行投影变换。新的伪投影仅仅对查询和显示有意义,而实际的数据不会被改变。

    例如,如果添加的第一个图层使用的是2000国家大地坐标投影坐标系统,其它的图层也会按照该投影系统动态地进行变换。同样,如果添加的第一层图如果使用的是WGS84地理坐标系统,其它的图层也会自动调整为与之一致的坐标系统。即使添加的是一个使用了投影坐标系统的数据,也会动态地取消投影 [链接] 。

6. Google、Microsoft、ArcGIS提供的地图服务分别采用了什么坐标系?

    现在都统一为Web Mercator,即正轴墨卡托投影,和墨卡托投影的主要区别是为了实现上的方便,把地球模拟为球体而非椭球体,精度理论上差别0.33%,比例尺大时基本可以忽略。同时纬度范围变成了(-85,85) [计算过程] ,南北极显示不了,但不影响正常使用,这样也减少了切图数量 [来源]。

    由于该投影仅用于地图显示,而不用于显示数字坐标,因此我们不需要椭球投影的额外精度。球面投影会导致 Y 方向上大约 0.33% 的比例失真,这在视觉上是不明显的 [来源]

7. Web Mercator(EPSG:3857)

    Web Mercator 是一个投影坐标系统,其基准面是 WGS 1984 Web Mercator 坐标系使用的投影方法不是严格意义的墨卡托投影,而是一个被 EPSG称为伪墨卡托的投影方法,这个伪墨卡托投影方法的大名是 Popular Visualization Pseudo Mercator,PVPM。 看起来就觉得这个投影方法不是很严谨的样子,大众化的?受欢迎的?可视化伪墨卡托投影……因为这个坐标系统是 Google Map 最先使用的,或者更确切地说,是Google 最先发明的。在投影过程中,将表示地球的参考椭球体近似的作为正球体处理(使得全世界可以在一个正方形里面刚好放下,极大简化了Web地图的涉及和计算的效率 [来源])。这也是为什么在 ArcGIS 中我们经常看到这个坐标系叫 WGS 1984 Web Mercator (Auxiliary Sphere 辅助、备用球形)。Auxiliary Sphere 就是在告知你,这个坐标在投影过程中,将椭球体近似为正球体做投影变换(球体半径等于 WGS 1984 长半轴的长度),虽然基准面是WGS 1984 椭球面。后来,Web Mercator 在 Web 地图领域被广泛使用,这个坐标系就名声大噪。尽管这个坐标系由于精度问题一度不被GIS专业人士接受,但最终 EPSG 还是给了 WKID:3857

    Web Mercator 最早是由 Google 提出的,当前已经成为 Web Map 的事实标准。但是也许是由于上面”伪“的原因,最初 Web Mercator 被拒绝分配EPSG代码。于是大家普遍使用 EPSG:900913(Google的数字变形) 的非官方代码来代表它,多用于一些开源项目中。直到2008年,才被分配了EPSG:3785的代码,但在同一年没多久,又被弃用,重新分配了 EPSG:3857 的正式代码,使用至今。至今已是Google Maps和几乎所有其他Web制图应用程序使用的坐标系

  📍 EPSG当初拒绝分配代码,并不是冷酷无情无理取闹,从技术角度看是有原因的。简而言之,主要原因在于基准面被篡改后,本来是等角投影的Mercator坐标变换算法,不再等角了,而是近似等角,也就是出现角度变形。这种变形势必影响了坐标的精度 

  📍 EPSG:3785 这是 EPSG 在 2008 年给 Web Mercator 设立的WKID,但是这个坐标系的基准面是正圆球,不是WGS 1984。 存在了一段时间后被弃用。

  📍 有两种用于仿真 Web 服务所用墨卡托投影的方法

    链接:EPSG:900913与EPSG:3857等效的

    链接:Web Mercator 公开的小秘密_慕名ArcGIS-CSDN博客

8. 为什么选择墨卡托投影?

    等角正轴圆柱投影,等角保证了对象的形状不变形,也保证了方向和相互位置的正确性(在航海、航空中应用)。

     谷歌地图Virtual EarthBing Maps百度地图MapabcArcGIS Online等采用Web Mercator(或称Spherical Mercator,球面墨卡托)坐标系(EPSG:3857),天地图采用CGCS2000国家大地坐标系(可选择的投影方式:Web墨卡托投影 或 经纬度投影) [来源]

9. 墨卡托投影的缺点

Web墨卡托投影下每个国家的大小和实际大小的差异

    等角的代价是面积的巨大变形,特别是两极地区; 墨卡托投影的地图最大的缺点就是和现实差别太大,面积变形非常严重

    在墨卡托投影的地图上,变形最严重的就是格陵兰岛了。地图上非洲的大小和格陵兰岛差不多大,但如果计算一下的话:非洲面积约是3020万平方千米,格陵兰岛面积约是217万平方千米,也就是说,非洲的面积是格陵兰岛的14倍!再比如墨卡托投影的地图上的加拿大,看起来是个瘦瘦的长方形,但实际上,加拿大是个类似正方形的形状。它不仅扭曲得不像样,而且还被放大了好几倍。在地图上加拿大的面积大约是美国的三倍,实际上加拿大比美国大不了多少 [来源]

10. UTM 与 高斯克吕格

  UTM投影全称为“通用横轴墨卡托投影”,英文名称为Universal Transverse Mercator,是一种等角横轴割圆柱投影,圆柱割地球于南纬80度、北纬84度两条等高圈,被许多国家用作地形图的数学基础,如中国采用的高斯-克吕格投影就是UTM投影的一种变形,很多遥感数据,如Landsat和Aster数据都应用UTM投影发布的 [来源] 我国规定1:1万、1:2.5万、1:5万、1:10万、1:25万、1:50万比例尺的地形图均采用高斯克吕格投影。该投影在英美等国家被称为横轴墨卡托投影 [来源]

  我国大于或等于1 : 50万的大中比例尺地形图多采用六度带的高斯-克吕格投影,三度带的高斯-克吕格投影多用于大比例尺1 : 1万测图。我国小于1 : 50万比例尺的地形图采用正轴等角割圆锥投影,又称为兰勃特正形圆锥投影。 (📍 并不是所有投影都有分带,分带是高斯克吕格投影自带的

高斯克吕格投影:中央经线和带号的计算

📌 以下哪一比例尺的地形图不适合采用高斯-克吕格投影的是(A)

A. 1:100万     B. 1:5万       C. 1:1万         D. 1:10万

📌 以下说法不正确的是(D)
A、1∶5万的地形图采用6度分带投影,即经差为6度,从零度子午线开始,自西向东每个经差6度为一投影带,全球共分60个带,用1,2,3,.....表示。
B、1∶1万的地形图采用3度分带投影,从东经1.5度的经线开始,每隔3度为一带,用1,2,3,......表示,全球共划分120个投影带。
C、地形图上公里网横坐标前2位就是带号,例如:1∶5万地形图上的横坐标为21210000,
其中21即为带号,210000为横坐标值。
D、6度带中央经线经度的计算方法:中央经线经度=6°×当地带号 (应为L=6N-3

📌 下列关于各坐标系描述中,正确的是(C)
A、基准面是在椭球的基础上建立的,即基准面和椭球体之间存在一对一的关系。
B、正高表示地面点到似大地水准面的垂直距离,正常高是指空间任一点沿该点处垂线至大地水准面的距离。
C、在高斯-克吕格投影中,由于每一个投影带的坐标都是此点对本带坐标原点的相对值,所以如果不在横轴坐标前加上带号,就会造成各带的坐标相同的现象。
D、我国基本比例尺地形图中,小于1:50万的地形图采用正轴等角圆柱投影,又叫墨卡托投影;航海图小于1:50万的地形图多用正轴等角割圆锥投影,又叫兰伯特投影。

11. 墨卡托投影、横轴墨卡托投影、通用横轴墨卡托投影(UTM)的区别和联系

  相同点是:三个投影都是圆柱投影

  区别如下:

  • 墨卡托投影是正轴等角切圆柱投影(web墨卡托则是把地球模拟成球形而不是椭球
  • 横轴墨卡托投影是横轴等角切圆柱投影(在中国称为高斯-克吕格投影,常用于地形图
  • 通用横轴墨卡托投影,也称通用墨卡托投影(UTM),是横轴等角圆柱投影(6°分带,遥感影像图大多用的就是UTM

12. 经纬度直投

  👉 摘录自 聊聊GIS中的坐标系|再版 - 四季留歌 - 博客园

  👉 摘录自 https://segmentfault.com/a/1190000041452965?sort=newest

  在我们使用二维方式展示地图,而坐标系为地理坐标系时,用到了是一种特殊的投影方式,经纬度等间隔直投

  经纬度等间隔直投:英文叫法是Platte Carre projection,是等距矩形投影Equirectangular projection)基准点纬度取0°(赤道)时的特殊情况。它的特点是相同的经纬度间隔在屏幕上的间距相等,没有复杂的坐标变换。我们可简单的理解为,在笛卡尔坐标系中,将赤道作为X轴,子午线作为Y轴,然后把本来应该在南北两极相交的经线一根一根屡直了,成为了互相平行的经线,而每条纬线的长度也在这个过程中都变为与赤道等长。

  在经纬度等间隔直投中,经度范围是-180180,纬度范围是-9090,因此投影出的地图是长方形,且长宽比是2:1

  把经纬度等同直角坐标来绘制,这样就有一个后果:

  当纬度=90度时,理论上有无限个点(即经度有无限个),若经纬度=直角坐标,那么赤道有多长,极点就展开拉的有多长。

  这种拉扯是非常恶心的,天地图就这么干,在我国高纬度地区,使用经纬度直投的绘制方法,会使得本来直角相交的两条道路斜变

上图为经纬度直投

  这也是ArcMap的默认“投影”,通常比较丑,尤其是高纬度地区变形非常严重,但是胜在渲染速度快。

  📍 网页上现在采用的主流是Web墨卡托,其次以经纬度直投为辅助。

   左图表示地球球面上大小相同的圆形,右上为墨卡托投影,投影后仍然是圆形,但是在高纬度时物体被严重放大了。右下为经纬度等间隔直投,圆的大小变化相对较小,但是高纬度时的图像明显被拉长了 [链接]。 

  

13. 罗宾森投影

    罗宾森投影是伪圆柱投影。罗宾森投影既不等角也不等积。它通常会造成形状、面积、距离、方向和角度畸变。罗宾森投影主要适用于一般世界地图。“国家地理”将罗宾森投影用于其世界地图约十年,直到 1998 年。罗宾森投影可能是一种最常用于绘制世界地图的折衷伪圆柱地图投影 [来源]。罗宾森投影也称为矫形法投影,是用于世界地图的折衷投影,用于一般及专题世界地图 [来源]。

    GIS 常用的地图投影与我们平时在教室或杂志所见投影未必相同。例如,罗宾逊投影在全球尺度的一般制图当中应用广泛,因其有美观效果。但是罗宾逊投影在 GIS 中不适用。 GIS 中使用的地图投影通常都有前述提及的一种保存性质,特别是正形性质。因为它保存了局部形状和角度,正形投影可以使相邻地图在角落处正确连接。这对开发地图系列十分重要,比如美国地质调查局( USGS )的标准地形图图幅地图。

14. 贝尔曼等积圆柱投影

    这种投影是一种适用于绘制世界地图的等积圆柱投影。 贝尔曼投影是圆柱等积地图投影的一种,其标准纬线设置为南北纬 30°。由于其等积的属性,它可以高度压缩极地地区。[来源1] [来源2]

15. 坐标系转换

  👉 参考 [GIS]地理坐标系与投影坐标系 - 知乎

  👉 参考 GDAL影像投影转换

  • 同一椭球的坐标系转换

    >>> 同一椭球间的转换是严密的。

    1️⃣ 同一椭球体投影换带

    2️⃣ 同一椭球不同投影坐标系转换

    3️⃣ 同一椭球体地理坐标系与投影坐标系的转换

  • 不同椭球的坐标系转换

    >>> 不同椭球间的转换是不严密的。

    1、平面坐标系之间的转换:四参数(至少要两个已知点)

  适用于已知的范围较小时,不考虑高程的影响,不同椭球之间高斯平面坐标的转换可以采用四参数。将两个坐标系都视为平面坐标。转换需要至少两个已知公共点坐标。

    2、椭球体之间的转换:三参数(最少要一个已知点)、七参数(至少要三个已知点)

  七参数:当范围较广阔时,根据已知公共点的多少可以采用七参数来计算,两个坐标系均为空间直角坐标系。至少需要三个已知点。

  三参数(特殊的七参数法):如果区域范围不大,最远点间的距离不大于 30km(经验值),也可以仅用三参数,即 X 平移、Y 平移、Z 平移,而将X 旋转、Y 旋转、Z 旋转、尺度变化四个参数取值为0,此方法需要在两个坐标系下至少提供一个共同已知点。

     🧰 算法:坐标转换四参数解算

  同一地理坐标系的投影转换是严密的,但不同地理坐标系之间需要先转换到地心立体坐标系,然后通过七参数转换 [链接]。 

16.投影转换

  从一种地图投影变为另一种地图投影的理论和方法。

  投影转换的方式:1)正解变换 x,y→X,Y  2)反解变换 x,y→B,L→X,Y

  根据转换的方法不同,投影转换可分为:1)解析变换;2)数值变换;3)解析-数值变换

17.逻辑坐标与窗口坐标转换

  地图发布到Web网页上时,涉及地理空间位置在网页容器中的表示。逻辑坐标指实际的地理坐标,表示真实的地理空间位置;窗口坐标指在Web网页中逻辑坐标对应的屏幕坐标,是根据Web网页中地图容器布局(大小与位置),对逻辑坐标进行转换而得到的。当在客户端实现图形交互绘制、地图查询、编辑等功能时,鼠标交互获取到的是窗口坐标,通常要将其转换为对应的逻辑坐标,进而实现具体功能。

  转换方法:将屏幕的某点和地理空间实际的某点进行配对,然后给定一个地理坐标根据该坐标和地理空间那点的偏移和比例尺来进行计算就可以获得屏幕坐标。

📃参考资料 WGS84与WGS84 Web Mercator

📃参考资料【GIS】地图与投影

# 2. SRID

    空间参考标识符 (SRID) 是与特定坐标系、容差和分辨率关联的唯一标识符。目前有多种公认的标准 SRID,例如欧洲石油标准组织EPSG(European Petroleum Standards Group,现已解散)或ESRI定义的 SRID [来源]

    WKID是分配给一个坐标系的唯一编号。每个WKID号码的权限将是Esri或EPSG的组织,这些来源的数字并不重叠,所以不必担心哪个机构定义了WKID [来源]

WKID ≤ 32767 权限:EPSG
WKID > 32767 权限:ESRI

    在ArcGIS中可以看到坐标系权限:

    (1)某些数据库和空间类型(如 PostgreSQL 中的 PostGIS 几何或 SQL Server 中的地理类型)使用预定义的 EPSG 代码子集,只可使用具有这些 SRID 的空间参考。

    (2)其他情况下(如 Oracle、PostgreSQL、IBM DB2 和 Informix 中的 ST_Geometry 类型)可使用具有 Esri 定义的 SRID 的空间参考。

  EPSG 查询网站 :EPSG.io: Coordinate Systems Worldwide ( GitHub:maptiler/epsg.io

  WKID 查询网站 :ArcGIS WKID

  EPSG维护着空间参照对象的数据集,OGC标准中空间参照系统的SRID与EPSG的空间参照系统ID相一致 [来源]

# 3. WKT

    此处的WKT 指 Well-Known Text Representation of Spatial Reference Systems,是OpenGIS定义地理参照坐标系的方式(表达坐标系的方式有好几种,有OpenGIS的WKT;有Proj4的表达方式;有EPSG权威编码方式等等)。👉 还有一种WKT指Well-known Text Representation of Geometry,是对几何对象进行文本描述

  如果保存任何坐标系统的投影文件(.prj),并在文本编辑器中打开它,将看到它的WKT。

  最新标准是 ISO 19162:2019,该标准取代了 ISO 19162:2015

  📍 版本历史:ISO 19125-1:2004WKT1)→ ISO 19162:2015WKT2)→ ISO 19162:2019(WKT2 version1)

  📍 Esri也参与了WKT 2规范的编写,在OGC的WKT 1规范上有一些变化,使之更严格,包含更少的歧义。因此,一些数据库区分“OGC WKT”“ESRI WKT”表示。这个问题在WKT 2中得到了很大程度的解决,因为它的定义更好了。一个区别是它没有TOWGS84,更像WKT2  [维基百科]

  📍 ArcGIS主要依据WKT1坐标系定义标准,即ISO 19125-1:2004;QGIS则符合WKT2坐标系定义标准,即ISO 19162:2019 [来源]

  链接:WKT介绍 - 百度文库

# 4. 软件中的坐标系

  ## 4.0 EPSG:40400

     EPSG:404000被认为是一个二维的通配符坐标参考系(单位为米)[官方链接],是由GeoTools定义并由GeoServer使用的自定义代码 [官方链接]。

    所有的投影坐标系都是2D平面,它们之间唯一的区别是,它们有不同的方程,用于将2D平面映射到球体上,404000是一个特例,因为它没有这个方程。只要你不使用这个方程(通过重新投影你的数据),普通投影坐标系和“真正的”2D平面之间就没有区别 [Using QGIS with EPSG:404000]。

    如果图层没有定义投影,OpenLayers就不能对它进行重投影,而且必须将它用作视图投影 [Add EPSG:404000 WMS to OpenLayers]。

  ## 4.1 ArcGIS

    在ArcGIS 10中Web Mercator有三种编号。他们分别是3857 、102100、102113(已弃用)。其实三者表示同一个投影 [来源]  

    👉 在ArcGIS Server的REST服务页面中,Web Mercator的空间参考会记做:102100(3857)

   4.1.1 ArcGIS中的WKID:3395WKID:3857

      两者在ArcGIS中的名称分别为WGS_1984_Web_Mercator_Auxiliary_Sphere(3857)WGS_1984_World_Mercator(3395),即一个是Web墨卡托投影,一个是(常规)墨卡托投影。比较两者的定义文件,差异如图所示。

  “辅助球体类型”参数可以是 0(使用地理坐标系的长半轴或半径)、1(使用短半轴或半径)、2(计算并使用等积半径)或 3(使用等积半径并将大地纬度转换为等积纬度)。 

    在ArcGIS Desktop的帮助文档“墨卡托投影”中有这样一句话:有两种用于仿真 Web 服务所用墨卡托投影的方法。如果墨卡托实现支持椭球体(椭圆体),则投影坐标系必须以基于球体的地理坐标系为基础。这要求必须使用球体方程。墨卡托辅助球体的实现仅具有球体方程。此外,如果地理坐标系是基于椭圆体的,它还具有一个投影参数,用于标识球体半径所使用的内容。默认值为零 (0) 时,将使用长半轴。

    在 EPSG.io 中写到,EPSG:3857相对于EPSG:3395,比例尺误差为0.7%,地图上的北角差异可达43公里(相当于地面上的21公里)。EPSG:3857的纬度范围为 [ 85.06°S , 85.06°N ],EPSG:3395的纬度范围为 [ 80°S , 84°N ]

  📃参考资料 投影知识点 - 太一吾鱼水 - 博客园

  ## 4.2 GeoServer 提供的坐标系切片方案

    注:GeoServer新版本也提供了瓦片大小为512*512的切片方案 


评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值