geopandas学习(二) 地理坐标参考系
文章参考:参考
1.坐标参考基础
1.1 CRS
现实世界中的地球作为一个球体,当我们想要用同样的方式利用坐标来唯一确定地球球面上的某个位置时,需要一套适应球体形状的坐标系统。
而当我们想要在纸面或电脑屏幕上绘制平面地图时,就又需要有一套将地球球面展平的方法。
上述的这些用于在不同情况下定义对象位置信息的坐标系统,就称为坐标参考系统(Coordinate Reference System,下文统称CRS)。CRS可细分为地理坐标系和投影坐标系。
1.1.1 地理坐标系
以弧度制下度数为单位的地理坐标系(Geographic Coordinate Systems)帮助我们定位物体在地球球面上的具体位置以及绘制球体地图
地理坐标系以地表上确定的某一个点为原点,创建了包裹全球的网格,譬如WGS84,将本初子午线与赤道的交点作为原点:
1.1.2 投影坐标系
地理坐标系虽然解决了我们在地球球面上定位的问题,但纬度和经度位置没有使用统一的测量单位。
为了解决上述问题,各种各样的投影坐标系(Projected Coordinate Systems)被开发出来(其中右下角为地理坐标系,其余均为投影坐标系):
投影坐标系指的是从将3D球面展平为2D平面的一套数学计算方法,利用它可以优化形状、比例/距离以及面积的失真情况。
但实际情况中没有在整个地球表面都能“三全其美”的投影坐标系,有些投影坐标系优化形状上的失真,有些投影坐标系优化距离上的失真,有些投影坐标系专门针对面积失真进行优化,而有些投影坐标系可以对局部区域进行三个方面上的优化。
常用的投影坐标系如横轴墨卡托(Universal Transverse Mercator,简称UTM),基于经度将全球等分为编号0-60的区域,且每个区域又进一步细分为南半球区域或北半球区域。
1.2 常用CRS格式
通过前文我们了解到什么是CRS,而在计算机系统中要使用CRS,需要将其文档化,下面我们来了解CRS两种常见的文档存储格式。
1.2.1 Proj4
Proj4字符串是一种识别空间或坐标参考系统的简洁方法,通过其定义的语法规则,将想要定义的CRS全部参数信息保存到一条字符串中。
Proj4语法
Proj4字符串包含了一种CRS全部元素信息,用+连接每个元素定义部分,如下面的例子记录了横轴墨卡托北11区CRS对应的Proj4字符串:
+proj=utm +zone=11 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0
它记录了如下信息:
- proj=utm:声明投影方法为墨卡托
- zone=11:声明对应北11区(因为这里是横轴墨卡托所以拥有独立分区,但并不是所有CRS都有分区,且在Proj4中区号加S才为南半球分区如11S,否则默认为北半球分区)
- datum=WGS84:声明基准面为WGS84(基准面是椭球体用来逼近某地区用的,因此各个国家都有各自的基准面。国内常用的基准面有:BEIJING1954,XIAN1980,WGS84等)
- units=m:声明坐标系单位设置为米
- ellps=WGS84:声明椭球面(如何计算