Para obtener un solo valor de píxel en un punto de un ráster con rasterio, hay un ejemplo aquí: https://github.com/mapbox/rasterio/pull/275
Sin embargo, ¿hay una API directa dentro de rasterio (y no el cli) que se pueda usar para extraer valor en un solo punto en un ráster?
- 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))
ovals = src.sample(x, y)
? Ninguno parece funcionarsrc.sample()
toma un iterador sobre lasx, y
tuplas, así que hazlofor val in src.sample([(x, y)]): print(val)
.