Obteniendo valores de píxel en un punto único usando rasterio

Question

Para obtener un valor de píxel único en un punto en un ráster usando rasterio, hay un ejemplo aquí: https://github.com/mapbox/rasterio/pull/275

Sin embargo, ¿hay una API directa dentro de rasterio (y no la interfaz de línea de comandos) que se pueda utilizar para extraer el valor en un solo punto en un ráster?

-- EDIT

with rasterio.drivers():

    # Leer bandas de ráster directamente en matrices Numpy.
    #
    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)

Answer 1

El método de la API de Python que admite el comando rio-sample está documentado aquí: https://rasterio.readthedocs.io/en/latest/api/rasterio._io.html#rasterio._io.DatasetReaderBase.sample

src.sample() toma un iterador sobre tuplas x, y, así que haz lo siguiente:

for val in src.sample([(x, y)]): 
    print(val)