添加 Z通过 gcps 使用 gdal_translate 和 gdalwarp 进行地理配准时图像的高程

编程问答 2022-07-05

问题描述

我一直在尝试使用给定的 gcps（X-Y 坐标和 Z（高程）坐标）对图像进行地理配准。但是，当我通过 gdal_translate 添加 gcps 时，它会添加高程 Z，但是一旦我使用 gdalwarp，高程就会消失，我只能得到一个 X-Y 参考图像。但我也需要图像的高度。这是我的代码：

import gdal
gdal_translate -of GTiff -gcp 0 0 15329.29997 21690.27067 5412.112601 -gcp 0 724 15330.29997 21690.27067 5412.112601 -gcp 1904 0 15329.29997 21689.27067 5412.112601 -gcp 1904 724  15330.29997 21689.27067 5412.112601 test.tif out-with-gcpXYZ.tif
gdalinfo "out-with-gcpXYZ.tif"
#then
gdalwarp out-with-gcpXYZ.tif test0Z.tif
gdalinfo "test0Z.tif" #it lags the elevation (Z)

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

elevationgdalgispyqgispython