使用仿射变换矩阵,将像素坐标映射到世界坐标。例如,使用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
.