PROJ坐标系转换库使用

基本工具使用

proj：用于同一参考椭球下，不同投影平面坐标系的转换；cs2cs：支持proj的功能，支持可自定义坐标系对的转换；支持不同参考椭球的转换；geod：提供测地线计算功能；gie：提供回归测试功能；cct：提供4D坐标框架变换，比如ITRF93 到 ITRF2000；projinfo：地理对象和坐标操作查询；

常用操作

wkt转换PROJ字符串：projinfo -o PROJ <WKT字符串>

不同椭球框架坐标转换1：cs2cs <PROJ坐标定义字符串> +to <PROJ坐标定义字符串>

不同椭球框架坐标转换2：cs2cs <PROJ坐标定义字符串> +to +init= epsg:id

查看PROJ库支持的投影方法：cs2cs -lp

查看PROJ库支持的椭球框架：cs2cs -le

查看PROJ库支持的单位：cs2cs -lu

PROJ库逆变换选项：+inv

例子:

正变换：echo 17.756132 45.3935192042 133.12 .8 | cct +proj=cart +ellps=GRS80

逆变换：echo 4272922.1553 1368283.0597 4518261.3501 .8000 | cct +proj=cart +inv +ellps=GRS80

PROJ字符串示例定义（以CGCS2000 100度4分 带为例）

+proj=tmerc +a=6378137 +rf=298.257222101 +pm=greenwich +x_0=500000 +y_0=0 +lon_0=100d4’E +lat_0=0d +k_0=1.0 +units=m +type=crs +no_defs

PROJ库变换坐标说明

大地坐标：longitude, latitude, elevation above ellipsoid, Time

平面坐标（投影坐标或者地心笛卡尔积坐标）：X, Y, Z, Time

例子：

17.756132 45.3935192042 133.12 .8