如何获取 geotiff 中单元格的坐标？

Question

我有一个包含地理信息的 tif。使用 gdal，我可以将光栅文件转换为数组 (numpy)。

如何获取该数组中一项的坐标？

Answer 1

使用仿射变换矩阵，将像素坐标映射到世界坐标。因此，例如，使用 affine 包。（使用简单的数学运算，还有其他方法可以做到这一点。）

from affine import Affine
fname = '/path/to/raster.tif'

这里有两种方法可以得到仿射变换矩阵，T0。例如，使用 GDAL/Python:

from osgeo import gdal
ds = gdal.Open(path, gdal.GA_ReadOnly)
T0 = Affine.from_gdal(*ds.GetGeoTransform())
ds = None  # close

例如，使用 rasterio:

import rasterio
with rasterio.open(fname, 'r') as r:
    T0 = r.affine

GDAL (T0) 使用的变换数组的约定是引用像素角。您可能想要引用像素中心，因此需要将其平移 50%：

T1 = T0 * Affine.translation(0.5, 0.5)

现在要从像素坐标转换为世界坐标，将坐标乘以矩阵，这可以用一个简单的函数来完成：

rc2xy = lambda r, c: T1 * (c, r)

现在，获取第一行第二列（索引 [0, 1]）中光栅的坐标：

print(rc2xy(0, 1))

另请注意，如果需要从世界坐标中获取像素坐标，可以使用倒仿射变换矩阵，~T0。

How to I get the coordinates of a cell in a geotif?