Non riesco a ottenere l'immagine del suolo con Python - OWSLib

Question

Sto cercando di ottenere un'immagine GTIFF con OWSLib in Python ma ottengo un errore. Mi aiuteresti?

Codice:

from owslib.wcs import WebCoverageService
wcs = WebCoverageService('http://maps.isric.org/mapserv?map=/map/phh2o.map', version='1.0.0')

bbox = (-54.7968690628577, -28.35600162342199, -54.77831702365789, -28.33600341642725)

response = wcs.getCoverage(
    identifier='phh2o_0-5cm_mean', 
    crs='urn:ogc:def:crs:EPSG::4326',
    bbox=bbox, 
    resx=250, resy=250, 
    format='GEOTIFF_INT16')

response.read()

Produzione:

b '\ n <ServiceExceptionReport version = "1.2.0" \ nxmlns = "http://www.opengis.net/ogc" xmlns: xsi = "http://www.w3.org/2001/XMLSchema-instance" xsi: schemaLocation = "http://www.opengis.net/ogc http://schemas.opengis.net/wcs/1.0.0/OGC-exception.xsd">\nmsImageCreate (): errore di gestione dell'immagine. Tentativo di allocare immagine raw non riuscito, memoria esaurita. \ N \ n \ n '

Answer 1

Su tua richiesta, stai lavorando con EPSG:4326un sistema di coordinate geografiche che utilizza i gradi. La richiesta resx=250, resy=250richiede una griglia con una risoluzione di 250 gradi, e il server proverà a leggere tutti i dati e farà una griglia aggregata e segnalerà un esaurimento della memoria.

Se provi qualcosa di simile 0.001

response = wcs.getCoverage(
    identifier='phh2o_0-5cm_mean', 
    crs='urn:ogc:def:crs:EPSG::4326',
    bbox=bbox, 
    resx=0.001 , resy=0.001, 
    format='GEOTIFF_INT16')

Avrai più successo ottenendo una griglia (risoluzione intorno ai 111 metri, a seconda di dove ti trovi sul pianeta).