如何从地图坐标 [经度，纬度] 获取波段计数 1 图像的像素值？

Question

所以我正在使用波段 4 的 Landsat-8 图像。我想获得给定经度和纬度值的 n 数组索引。我试过下面的代码

xs=np.array([92.8])
ys=np.array([26.6])
row,col=rasterio.transform.rowcol(src.transform,xs,ys)
print(row,col)
>>[99709] [-11627]     #output from the code 

获得的输出是错误的，因为我的图像大小是 (7791, 7651) 并且这些值不在索引范围内。 谁能建议我做一些改变。

Answer 1

你能解决这个问题吗:)，我有同样的问题！