gdaltransform

转换坐标。

简介

gdaltransform [--help-general]
    [-i] [-s_srs srs_def] [-t_srs srs_def] [-to "NAME=VALUE"]
    [-ct proj_string] [-order n] [-tps] [-rpc] [-geoloc]
    [-gcp pixel line easting northing [elevation]]* [-output_xy]
    [srcfile [dstfile]]

描述

gdaltransform实用程序将坐标列表重新投影到任何支持的投影中,包括基于GCP的转换。

-s_srs <srs_def>

设置源空间引用。可以传递的坐标系是OGRSpatialReference.SetFromUserInput()调用支持的任何坐标系,它包括EPSG pc和gcse(即EPSG:4296)、PROJ.4声明(如上所述)或包含已知文本的.prj文件的名称。

-t_srs <srs_def>

设置目标空间参考。可以传递的坐标系是OGRSpatialReference.SetFromUserInput()调用支持的任何坐标系,它包括EPSG pc和gcse(即EPSG:4296)、PROJ.4声明(如上所述)或包含已知文本的.prj文件的名称。

-ct <string>

PROJ字符串(以+PROJ=pipeline开头的单步操作或多步字符串)、描述CoordinateOperation的WKT2字符串或urn:ogc:def:CoordinateOperation:EPSG::XXXX urn重写从源到目标CRS的默认转换。它必须考虑源和目标CRS的轴顺序。

3.0 新版功能.

-to NAME=VALUE

设置适合传递到的变压器选项 GDALCreateGenImgProjTransformer2() .

-order <n>

用于翘曲的多项式的阶数(1到3)。默认值是根据gcp的数量选择多项式顺序。

-tps

基于现有GCPs的薄板样条变压器的受力分析。

-rpc

强制使用RPCs。

-geoloc

强制使用地理定位阵列。

-i

反变换:从目的地到来源地。

-gcp <pixel> <line> <easting> <northing> [<elevation>]

提供用于转换的GCP(通常需要三个或更多)

-output_xy

将输出限制为“x y”而不是“x y z”

<srcfile>

具有源投影定义或GCP的文件。如果未给定,则从命令行读取源投影 -s_srs-gcp 参数

<dstfile>

具有目标投影定义的文件。

坐标以成对、三元组(对于3D)或四元组(对于X、Y、Z、time)的形式从标准输入中读取,并以相同的方式转换和输出到标准输出。gdalwarp提供的所有转换都被处理,包括基于gcp的转换。

请注意,输入和输出必须始终采用十进制形式。目前不支持DMS输入或输出。

如果提供了输入图像文件,则输入以该图像上的像素/线坐标为单位。如果提供输出文件,则输出以该图像上的像素/线坐标为单位。

实例

重投影示例

从一个投影坐标系到另一个投影坐标系的简单重投影:

gdaltransform -s_srs EPSG:28992 -t_srs EPSG:31370
177502 311865

在“Belge 1972/Belgian Lambert 72”投影中以米为单位生成以下输出:

244510.77404604 166154.532871342 -1046.79270555763

图像RPC示例

以下命令使用与命名文件关联的RPC模型请求基于RPC的转换。因为使用了-i(逆)标志,所以转换是从输出地理参考(WGS84)坐标返回到图像坐标。

gdaltransform -i -rpc 06OCT20025052-P2AS-005553965230_01_P001.TIF
125.67206 39.85307 50

生成以图像上的像素和线条为单位的输出:

3499.49282422381 2910.83892848414 50

十、 Y,Z,时间变换

历元2000.0时坐标从ITRF2000到ITRF93的15项时变Helmert坐标变换

gdaltransform -ct "+proj=pipeline +step +proj=unitconvert +xy_in=deg \
+xy_out=rad +step +proj=cart +step +proj=helmert +convention=position_vector \
+x=0.0127 +dx=-0.0029 +rx=-0.00039 +drx=-0.00011 +y=0.0065 +dy=-0.0002 \
+ry=0.00080 +dry=-0.00019 +z=-0.0209 +dz=-0.0006 +rz=-0.00114 +drz=0.00007 \
+s=0.00195 +ds=0.00001 +t_epoch=1988.0 +step +proj=cart +inv +step \
+proj=unitconvert +xy_in=rad +xy_out=deg"
2 49 0 2000

生成以经度、纬度和椭球体高度(米)为单位测量的输出:

2.0000005420366 49.0000003766711 -0.0222802283242345