如何在geotif中获得单元格的坐标？

Question

我有地理信息的tif。使用gdal，我可以将光栅文件转换为数组（numpy）。在

如何获得数组中一个条目的坐标？在

Answer 1

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

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

这里有两种方法可以得到仿射变换矩阵T0。E、 g.使用GDAL/Python：

^{pr2}$

例如，使用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: (c, r) * T1

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

print(rc2xy(0, 1))

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