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”
- <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