Para obter um único valor de pixel em um ponto de uma varredura usando rasterio, há um exemplo aqui: https://github.com/mapbox/rasterio/pull/275
No entanto, existe uma API direta no rasterio (e não no cli) que pode ser usada para extrair valor em um único ponto em um raster?
- EDITAR
with rasterio.drivers():
# Read raster bands directly to Numpy arrays.
#
with rasterio.open('C:\\Users\\rit\\38ERP.tif') as src:
x = (src.bounds.left + src.bounds.right) / 2.0
y = (src.bounds.bottom + src.bounds.top) / 2.0
vals = src.sample((x, y))
for val in vals:
print list(val)
vals = src.sample((x, y))
ouvals = src.sample(x, y)
? Não parece funcionar