聊聊GIS中的坐标系|再版 识别各种数据的坐标系及代码中的坐标系
本篇讲讲在GIS桌面软件和实际数据中，以及各路GIS有关API的编程中，如何寻找坐标系信息。惯例：
本文约2000字，建议阅读时间10分钟。
作者：博客园/B站/知乎/csdn/小专栏 @秋意正寒
版权：转载请告知，并在转载文上附上转载声明与原文链接（https://www.cnblogs.com/onsummer/p/12082359.html）
有关“未定义坐标系”的数据如何粗略判别是什么坐标系，已经在上篇的第5.1节写明了，这里想说的是，已经定义好坐标系统的程序/数据，在哪里找。
目录：
0. 预备知识（WKT的概念）
1. 硬盘上的数据/数据库（shp、geojson、tif栅格、几个单文件数据库、kml/gml/gpx）
2. GIS数据服务种的坐标系（2.1 ArcGIS家的/2.2 geoserver的）
3. 代码种的坐标系（3.1 AO/AE 3.2 ArcGIS JsAPI 3.3 openlayers 3.4 ceisum）
0. 预备知识
0.1. WKT
和json在网络传输的地位一样，是一个标记语言，全称Well-Known Text，和上文提及的Well-Known ID出自同一个组织——OGC。
WKT是用文本形式记录地理信息的一个标记语言。什么是地理信息？除开我们常见的真·地理信息点线面之类，还包括坐标系统的定义。我们可以在epsg.io这个网站查询想知道的坐标系的WKID，也可以查到想知道坐标系的WKT。
例如，WKID=4326的坐标系的WKT如下：
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG",""]],
AUTHORITY["EPSG",""]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG",""]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG",""]],
AUTHORITY["EPSG",""]]
此处不展开WKT的每一项定义。因为WKT的博客实在太多了，本文略过不作为主要内容。
通常，如果拿到一个WKT，我会观察最外层定义是GEOGCS还是PROJCS，这能判断是地理坐标系或投影坐标系统。
然后，若为GEOGCS，我会在AUTHORITY处判断它被epsg定义的WKID，上面这个wkt即4326.
WKT还有一种二进制存储格式，叫WKB。有关这两个东西的文档在ogc官网是可查的。
1. 硬盘上的数据/数据库
1.1. shapefiles
简称shp文件。shp文件至少有*.shp、*.shx、*.dbf三个同名文件组成，若想为shp文件定义一个坐标系，告诉它“你的空间信息落于哪个坐标系下”，那就要再多一个同名的*.prj文件。
一个shp文件通常要有以上4个文件，若无prj文件只会在加载时不知道它的坐标系信息，实际数据并无影响，所以说shp文件的最低下限是shp、shx、dbf三大文件。
有关shp文件的构成格式不再补充，也能在ogc官网查到对应标准。
此prj文件内坐标系的信息，就是WKT文本。
PROJCS["Xian_1980_3_Degree_GK_CM_114E",
GEOGCS["GCS_Xian_1980",
DATUM["D_Xian_1980",
PHEROID["Xian_1980",6378140.0,298.257]],
PRIMEM["Greenwich",0.0],
UNIT["Degree",0.0174532925199433]],
PROJECTION["Gauss_Kruger"],
PARAMETER["False_Easting",500000.0],
PARAMETER["False_Northing",0.0],
PARAMETER["Central_Meridian",114.0],
PARAMETER["Scale_Factor",1.0],
PARAMETER["Latitude_Of_Origin",0.0],
UNIT["Meter",1.0]]
此处若无AUTHORITY信息，可以用PROJCS内第一个字符串"Xian_1980_3_Degree_GK_CM_114E"去epsg.io查询（epsg.io不支持模糊查询，用空格替换下划线，尽量删减一些关键词）
1.2. GeoJson
geojson是一种用json标记的地理数据文件，记录二维信息多一些，通常用于webgis传递。与turf、arcgis jsapi、openlayers等结合较好。
geojson也是一种ogc规范，开源。
打开一个geojson文件：
其中，头部“crs”字段就是坐标系信息，我们可以访问其properties.name获取有关EPSG信息，此处WKID为32650。
但是，在大多数时候geojson是没有crs字段的，例如：
这就要问数据的提供者了。
1.3. 影像数据（tif文件为例）
tif文件为栅格图像文件，后缀为tif或tiff，是ogc规范的一种，全称GeoTiff。
通常不能在资源管理器中查看tif栅格影像数据的坐标系信息，需要用GIS软件查看，因为它的坐标系信息写在数据文件内部。
如下图：
此为ArcMap中Catalog查看数据的属性。
1.4. postgis/geopackage/esri geodatabase等数据库
这些地理数据库的数据均可在客户端查看，以QGIS为例，打开这些数据库里的要素查看坐标系信息。
①postgis
postgresql的一个插件。
②geopackage
geopackage是一种轻量的单文件数据库，是sqlite的特例应用。
geopackage在QGIS里的支持还不算很强，我的版本是3.8.1，还不支持右键导入，拖拽的方式倒是可以。
读取坐标系的方法一样。
③gdb
QGIS只能读取gdb里的表格和要素数据，复杂的数据因为esri商用问题不能读取。
1.5. kml(kmz)/gml/gpx
kmz用压缩软件解压后可获取kml文件。kml是xml格式的一种特殊化，原来用于Google Earth的数据交互，它与geojson的区别在于它可以带符号、颜色、样式信息。
kml由于早期为Google Earth使用，由上篇可知默认使用WGS84（WKID=4326）坐标系，所以大多数kml文件是没有写入坐标系信息的（至少我在QGIS里导出的kml没有）。
但是，在ogc文档中有定义坐标系：http://docs.opengeospatial.org/is/12-007r2/12-007r2.html#1274
此处不展开。
gml也是类似的一种xml扩展标记语言，同样为ogc规范的一种。
其坐标系信息写在这里，比较显眼：
最后一个gpx文件，是GPS卫星的信息记录，也是一种xml的应用。因为是GPS（早期GPS指美国GPS），所以用的是WGS84坐标系。
2. GIS数据服务
2.1. ArcGIS Server
①地图服务
访问MapServerURL，可查看地图服务的基本信息：
https://sampleserver6.arcgisonline.com/arcgis/rest/services/Census/MapServer
圈出的Spatial Reference即为此地图服务（包括所有图层）的WKID。
②影像服务
同地图服务，但是影像服务并没有子图层。
例子：https://landsat2.arcgis.com/arcgis/rest/services/Landsat8_Views/ImageServer
③场景服务
场景服务，URL结尾是SceneServer，访问后返回json，在控制台可以通过这样来访问：
即可获取其坐标系。
注意，fetch函数在各主流浏览器中较新版本均支持，IE不清楚未测试。
本例URL为：http://scene.arcgis.com/arcgis/rest/services/Hosted/Building_Hamburg/SceneServer
④要素服务
基本方法同地图服务，但是要注意的是要素服务既可单独作为FeatureServer，也可以存在于MapServer中（作为要素图层）。
URL例子：
https://services.arcgis.com/V6ZHFr6zdgNZuVG0/arcgis/rest/services/ks_earthquakes_since_2000/FeatureServer
https://sampleserver6.arcgisonline.com/arcgis/rest/services/USA/MapServer/0（若为MapServer中的要素图层，请访问上一级的MapServer查看坐标系）
⑤网络分析服务
网络分析服务，URL以NAServer结尾，因博主暂无网络分析服务的URL，暂时不能给例子，用法应该与MapServer类似。
⑥高程服务
高程服务，以ImageServer为载体。方法同影像服务。
例子：https://sampleserver6.arcgisonline.com/arcgis/rest/services/Elevation/MtBaldy_Elevation/ImageServer
2.2. GeoServer
wms/wfs/wcs
访问其描述xml，然后查看layer对应的"CRS"属性即可。
例如：
wms的URL：http://<domain>:<port>/geoserver/wms?service=wfs&version=1.1.0&request=GetCapabilities
wfs的URL：http://<domain>:<port>/geoserver/wfs?service=wfs&version=1.1.0&request=GetCapabilities
wcs的URL：http://<domain>:<port>/geoserver/wfs?service=wcs&version=1.1.0&request=GetCapabilities
访问此URL，浏览器会下载一个xml文件，通过寻找以下标签的值，即可获取对应图层的坐标系：
<WMT_MS_Capabilities>
今天在我的springMVC项目--图书管理系统中,希望在登录时将登录的Users存在session中,开始是准备在controller中使用Servlet API中的对象,可是一直无法引用,不知道为 ... 在Android中获取网页里表单中的数据具体实现代码如下,感兴趣的各位可以参考过下哈,希望对大家有所帮助 MainActivity如下: 复制代码 代码如下: package cn.testjavas ... package cn.jy.demo; import java.sql.Connection;import java.sql.PreparedStatement;import java.sql.Res ... 最近研究将K3系统与生产管理系统结合起来,减少工作量,但如何确定其各自后台数据库的构成,其对应数据各自位于那张表内,总结了一下: 1.从百度搜索,查看表结构,然后找到目标表 另:K3数据库中单独 ... <!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ... A B C D 1 10 3 有 2 6 e 无 3 3 6 有 判断c列的值在A列中是否存在(假定C列为需要判断列,A列为目标列) 在D1中输入以下公式,然后下拉公式即可 =IF(C ... <!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ... 链表插入示意图:(图是个人所画)因为链表指针指来指去,难以理解,所以辅助画图更加方便. 插入某个学号后面图: 定义的结构体: struct student { ]; //学生学号 ]; //学生姓名 ... insert into tableA select * from tableB b where not exists(select 1 from tableA a where a.id = b.id) ... WPF窗口充满了各种元素,但这些元素中只有一部分是控件.在WPF领域,控件通常被描述为与用户交互的元素——能接收焦点并接受键盘或鼠标输入的元素.明显的例子包括文本框和按钮.然而,这个区别有时有些模糊. ... 这是 HTTP 系列的第三篇文章,此篇文章为 HTTP 的进阶文章. 在前面两篇文章中我们讲述了 HTTP 的入门,HTTP 所有常用标头的概述,这篇文章我们来聊一下 HTTP 的一些 黑科技. HT ... 啦啦啦啦啦 恶搞别人吗? 把下面代码做成html文件发给别人,用浏览器打开就可以看见效果了 <!DOCTYPE html> <html><head><meta ... 一.Spring Boot与缓存. JSR-107.Spring缓存抽象.整合Redis 一.JSR107 Java Caching定义了5个核心接口,分别是CachingProvider, Cach ... 代码写不对队:Beta冲刺版本 这个作业属于哪个课程 http://edu.cnblogs.com/campus/xnsy/GeographicInformationScience/homework ... 感觉这题还可以 因为总空间比输入数量 不知高到哪里去了 ,所以完全不需要考虑放不下的问题 从贪心的角度考虑,如果要使相差数量巨大的\(b\)和\(w\)能够成功放下来,应该使这些方块尽量分散(似乎有点 ... <!DOCTYPE html> <html> <head> <meta charset="UTF-8"> <title> ... Servlet 梳理 概述 Web 技术成为当今主流的互联网 Web 应用技术之一,而 Servlet 是 Java Web 技术的核心基础. 要介绍 Servlet 必须要先把 Servlet 容器 ... Java基于OpenCV实现走迷宫(图片+路线展示) 由于疫情,待在家中,太过无聊.同学发了我张迷宫图片,让我走迷宫来缓解暴躁,于是乎就码了一个程序出来.特此记录. 原图: 这张图,由于不是非常清晰, ... 1. Two Sum Given an array of integers, return indices of the two numbers such that they add up to a ...聊聊GIS中的坐标系|再版 识别各种数据的坐标系及代码中的坐标系的更多相关文章
随机推荐