使用经度和纬度通过rasterio获取值

Question

我正在处理包含地理空间数据的.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

Answer 1

是的。来自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)