使用经度和纬度获取光栅值答案

【问题标题】：Using longitude and latitude to get values with rasterio使用经度和纬度获取光栅值
【发布时间】：2020-07-28 07:18:10
【问题描述】：

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

【问题讨论】：

标签： raster rasterio

【解决方案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)

【讨论】：