gpt4 book ai didi

python - 如何获取 geotiff 中单元格的坐标?

转载 作者:太空狗 更新时间:2023-10-30 00:40:21 32 4
gpt4 key购买 nike

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

如何获取该数组中一个条目的坐标?

最佳答案

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

关于python - 如何获取 geotiff 中单元格的坐标?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/27861197/

32 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com