如何读取GDAL中像素的地理坐标？

Posted 2023-02-23

【中文标题】如何读取GDAL中像素的地理坐标？

【英文标题】：How to read the geographical coordinates of a pixel in GDAL?

【发布时间】：2013-10-08 10:19:23

【问题描述】：

我在 UTM 投影中有一个图像，给定一个像素值 (px, py)，我需要获取其对应的 utm 坐标 (utx, utmy)。 是否可以？我试过了

gdallocationinfo  input.tif 256 256

但我没有得到预期的结果。

谢谢

M.

【参考方案1】：

gdallocationinfo 用于从栅格中获取像素值。对于转换，您应该使用gdaltransform。如果您提供目标图像作为参数，您可以输入像素坐标并获得投影坐标作为回报。

所以：

gdaltransform input.tif
256 256

【讨论】：

感谢您的回复。我做了gdaltransform -t_srs EPSG:32629 input_red.tif 256 256，但我还需要指定源空间参考集（gdaltransform -s-srs ? -t_srs EPSG:32629 input_red.tif 256 256）。我不知道如何引用像素空间作为源引用。

源参考是你的图像，目标参考是图像的投影。您无需指定s_srs 或t_srs。

谢谢，成功了！我确实尝试了您的解决方案，但我执行了整行，我没有在第一行之后输入像素值。

是的，工作与其他 GDAL cmd 工具略有不同。有关编程 python 解决方案，请参见例如：osgeo-org.1560.n6.nabble.com/…

【参考方案2】：

这是一道数学题，您只需要fetch the affine transformation coefficients，您可以通过 Python 绑定轻松完成。此外，它取决于您的像素坐标系（例如，它是以 [1,1] 还是 [0,0] 开头？以及，哪个角是原点？）

【讨论】：

它以 [1,1] 开头。谢谢你的链接。我需要从 groovy 脚本中的像素坐标中获取 utm 坐标。我会检查你建议的 python 绑定。由于我对这些东西很陌生，我想知道是否有命令行指令进行翻译而无需包装 python 代码。