pyproj.Proj也支持直接通过epsg代码设置投影坐标系。
注:通过pyproj.Proj可以将球面上的点投影到平面上,具体用法可参考 利用pyproj将经纬度投影为平面坐标以及地理坐标系背景知识.
EPSP的英文全称是European Petroleum Survey Group,中文名称为欧洲石油调查组织。它负责维护并发布坐标参照系统的数据集参数,以及坐标转换描述,该数据集被广泛接受并使用,通过一个Web发布平台进行分发,同时提供了微软Acess数据库的存储文件,通过SQL脚本文件,mySQL, Oracle和PostgreSQL等数据库也可使用。
目前已有的参考椭球体,大地坐标系,投影坐标系等不同组合都对应着不同的epsg ID号,在EPSG中被称为EPSG code,它代表特定的椭球体、单位、地理坐标系或投影坐标系等信息。EPSG对世界的每一个地方都制定了地图,但是由于座标系不同,所以地图也各不相同。具体不同的EPSG code可在Coordinate Systems Worldwide上查询。
比较常用的几个EPSG CODE如下表:
EPSG CODE | 含义 |
---|---|
4479 | 中国大地2000坐标系,中国专用 |
4326 | WGS84坐标系,GPS使用 |
3857 | 基于WGS84的伪墨卡托投影坐标系,互联网常用 |
pyproj.CRS表示一个坐标参数系统(Coordinate Reference System),pyproj.CRS.from_epsg()支持epsg代码。比如:
In [34]: crs=pyproj.CRS.from_epsg(4479)
In [35]: crs
Out[35]:
<Geocentric CRS: EPSG:4479>
Name: China Geodetic Coordinate System 2000
Axis Info [cartesian]:
- X[geocentricX]: Geocentric X (metre)
- Y[geocentricY]: Geocentric Y (metre)
- Z[geocentricZ]: Geocentric Z (metre)
Area of Use:
- name: China
- bounds: (73.62, 16.7, 134.77, 53.56)
Datum: China 2000
- Ellipsoid: CGCS2000
- Prime Meridian: Greenwich
然而,epsg代码既可以表示投影坐标系(如EPSG:3070表示NAD83/Wisconsin Transverse Mercator),
也可以表示大地坐标系(如EPSG:4479表示China Geodetic Coordinate System 2000 ),
还可以表示大地基准面(如EPSG:1043表示China 2000 ),
还可以表示参考椭球(如EPSG:1024表示CGCS2000 )。
注:上述四种坐标系的定义和区别可参考 利用pyproj将经纬度投影为平面坐标以及地理坐标系背景知识.
但pyproj.CRS.from_epsg()只支持表示投影坐标系和大地坐标系的EPSG代码,比如:
In [40]: crs=pyproj.CRS.from_epsg(3070)
In [41]: crs
Out[41]:
<Projected CRS: EPSG:3070>
Name: NAD83 / Wisconsin Transverse Mercator
Axis Info [cartesian]:
- X[east]: E