Je ne peux pas obtenir une image de grille de sol avec Python - OWSLib

Question

J'essaye d'obtenir un GTIFF imagen avec OWSLib en Python mais j'obtiens une erreur. Voulez-vous me aider?

Code:

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()

Production:

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 (): erreur de gestion d'image. La tentative d'allocation de l'image brute a échoué, mémoire insuffisante. \ N \ n \ n '

Answer 1

À votre demande, vous travaillez avec EPSG:4326un système de coordonnées géographiques utilisant des degrés. La demande resx=250, resy=250demande une grille avec une résolution de 250 degrés, et le serveur essaiera de lire toutes les données et de créer une grille agrégée et de signaler un manque de mémoire.

Si vous essayez quelque chose comme 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')

Vous aurez plus de succès à obtenir une grille (résolution d'environ 111 mètres, selon l'endroit où vous vous trouvez sur la planète).