我正在处理包含地理空间数据的 .tif 图像文件,重点关注地形的海拔。我可以通过使用 dem1_band 的数组索引来获取点的海拔(如下所示),但我希望根据纬度和经度获取海拔。有没有办法使用 rasterio 库来做到这一点?
dem1 = rasterio.open('/content/DEMs/dem1.tif')
dem1_band = dem1.read(1).astype('float64')
print(dem1_band[3000,3000])
>>451.13458251953125
最佳答案
是的。来自 https://gis.stackexchange.com/a/346458/55948 :
import numpy as np
import rasterio
xs = np.array([130.5, 146.0])
ys = np.array([-25.5, -42.0])
with rasterio.open("my.tif") as src:
rows, cols = rasterio.transform.rowcol(src.transform, xs, ys)
关于raster - 使用经度和纬度通过 rasterio 获取值,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/63122557/