Nous utilisons des cookies pour vous garantir la meilleure expérience sur notre site. Si vous continuez à utiliser ce dernier, nous considèrerons que vous acceptez l'utilisation des cookies. J'ai compris ! ou En savoir plus !.
banniere

Le portail francophone de la géomatique


Toujours pas inscrit ? Mot de passe oublié ?
Nom d'utilisateur    Mot de passe              Toujours pas inscrit ?   Mot de passe oublié ?

Annonce

Printemps des cartes 2024

#1 Fri 11 December 2020 15:02

Samuel Branchu
Juste Inscrit !
Date d'inscription: 11 Dec 2020
Messages: 1

Fond de carte Geoportail avec matplotlib : erreur "bbox invalide"

Bonjour,
J'espère être au bon endroit pour poser ma question.

J'ai installé python sous Windows (avec python-3.9.0-amd64.exe) puis les modules suivants avec pip : matplotlib, pyproj, OWSLib.
Enuite, j'ai installé basemap avec un fichier wheel trouvé sur le site
https://www.lfd.uci.edu/~gohlke/pythonlibs/#basemap (Unofficial Windows Binaries for Python Extension Packages)
pip install basemap-1.2.2-cp39-cp39-win_amd64.whl

Je veux réaliser un fond de carte issu du Géoportail identique à celui obtenu avec la requête suivante qui marche (EPSG:4326 étant WGS84 2D) :
https://wxs.ign.fr/choisirgeoportail/ge … T=256&

Avec le code python suivant, en choisissant EPSG:2154 (RGF93 / Lambert-93) l'image produite est correcte, avec un fond blanc hors des limites du RGF93 ; mais avec EPSG:4326, comme dans la requête web ci-dessus, j'obtiens le message "bbox invalide" (détail ci-dessous).
Le code minimal est le suivant (pyplot n'est utile que pour faire l'image, il n'est pas en cause).

Code:

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
plt.figure(figsize=(15,22))
map = Basemap(projection='merc', lat_0=25, lon_0=25, resolution='i', area_thresh=0.1,
              llcrnrlon=0, llcrnrlat=0, urcrnrlon=50, urcrnrlat=50, 
              # epsg=4326)
              epsg=2154)
map.wmsimage("https://wxs.ign.fr/choisirgeoportail/geoportail/r/wms?SERVICE=WMS&VERSION=1.3.0&CRS=EPSG:4326", layers=["GEOGRAPHICALGRIDSYSTEMS.PLANIGNV2"], verbose=True, format='jpeg')
plt.show()

Ma question : comment comprendre le message "bbox invalide" ?
Eventuellement, y'a-t-il une autre layer permettant d'avoir un fond de carte à toute échelle en tout point du globe, avec la clé choisirgeoportail, accessible par wmsimage (qui supporte WMS 1.1.1 et 1.3.0) ?
Je vous remercie.

Retour complet :

https://wxs.ign.fr/choisirgeoportail/ge … =EPSG:4326
id: WMS, version: 1.1.1
title: Service WMS Raster Geoportail Publics, abstract: Ce service permet la visualisation de couches de données raster IGN au travers d'un flux WMS
available layers:
['GEOGRAPHICALGRIDSYSTEMS.PLANIGNV2']
projection options:
['EPSG:2978', 'EPSG:3305', 'IGNF:RGPFUTM6S', 'EPSG:4642', 'EPSG:4622', 'EPSG:4641', 'EPSG:3171', 'IGNF:LAMB1C', 'IGNF:LAMBE', 'EPSG:2977', 'IGNF:RGM04UTM38S', 'EPSG:2990', 'EPSG:4633', 'IGNF:RGF93CC42', 'IGNF:TAHI51UTM06S', 'EPSG:27571', 'EPSG:32701', 'EPSG:27583', 'EPSG:23032', 'EPSG:2989', 'EPSG:3172', 'IGNF:WGS84UTM20', 'EPSG:32707', 'EPSG:3296', 'EPSG:3043', 'EPSG:3306', 'IGNF:RGF93CC46', 'IGNF:REUN47GEO', 'EPSG:32739', 'EPSG:27564', 'IGNF:LAMB1', 'IGNF:RGF93CC48', 'EPSG:3042', 'EPSG:3942', 'EPSG:3044', 'EPSG:32620', 'EPSG:4171', 'CRS:84', 'EPSG:27592', 'EPSG:3857', 'EPSG:4621', 'EPSG:2970', 'IGNF:UTM32W84', 'EPSG:4624', 'EPSG:2973', 'IGNF:TAHI79UTM6S', 'EPSG:4687', 'IGNF:RGF93G', 'IGNF:LAMB93', 'EPSG:3034', 'EPSG:4630', 'EPSG:4623', 'EPSG:27594', 'EPSG:4627', 'EPSG:4749', 'EPSG:4639', 'EPSG:27572', 'EPSG:3035', 'EPSG:3949', 'IGNF:UTM20W84GUAD', 'EPSG:27563', 'EPSG:4692', 'EPSG:4629', 'IGNF:GUAD48UTM20', 'EPSG:32632', 'EPSG:3950', 'EPSG:4559', 'EPSG:4690', 'EPSG:4636', 'IGNF:RGSPM06U21', 'EPSG:3947', 'EPSG:23030', 'EPSG:4638', 'EPSG:4643', 'IGNF:REUN47GAUSSL', 'IGNF:RGFG95GEO', 'IGNF:GUADFM49U20', 'IGNF:UTM20W84', 'IGNF:CSG67UTM22', 'EPSG:4275', 'IGNF:LAMB4C', 'EPSG:4258', 'EPSG:2975', 'IGNF:WGS84G', 'IGNF:UTM39SW84', 'EPSG:4470', 'EPSG:27591', 'EPSG:32621', 'EPSG:4626', 'EPSG:2972', 'EPSG:3944', 'EPSG:3298', 'EPSG:3170', 'EPSG:2981', 'IGNF:UTM22RGFG95', 'IGNF:CSG67GEO', 'EPSG:4326', 'IGNF:STPM50GEO', 'EPSG:4558', 'EPSG:21781', 'EPSG:4632', 'IGNF:LAMB2C', 'IGNF:RGF93CC43', 'EPSG:2971', 'IGNF:LAMB2', 'EPSG:4691', 'IGNF:UTM30W84', 'EPSG:32740', 'EPSG:27581', 'IGNF:RGF93CC45', 'EPSG:2987', 'EPSG:3304', 'IGNF:LAMB3C', 'EPSG:2988', 'EPSG:2976', 'EPSG:4628', 'EPSG:2980', 'IGNF:MAYO50UTM38S', 'EPSG:3302', 'EPSG:27593', 'IGNF:MAYO50GEO', 'EPSG:27584', 'IGNF:GUAD48GEO', 'EPSG:3943', 'IGNF:RGR92UTM40S', 'IGNF:LAMB3', 'EPSG:32705', 'IGNF:RGF93CC44', 'EPSG:2969', 'EPSG:4471', 'EPSG:4637', 'EPSG:4689', 'EPSG:4807', 'IGNF:UTM20W84MART', 'IGNF:MART38GEO', 'IGNF:STPM50UTM21', 'EPSG:3303', 'EPSG:27573', 'EPSG:23031', 'EPSG:32738', 'EPSG:4467', 'EPSG:3945', 'IGNF:UTM31W84', 'EPSG:4644', 'IGNF:MART38UTM20', 'IGNF:RGPFUTM7S', 'EPSG:32631', 'IGNF:RGPFUTM5S', 'IGNF:RGF93CC49', 'EPSG:27582', 'EPSG:2154', 'IGNF:UTM01SW84', 'EPSG:3727', 'EPSG:3948', 'IGNF:LAMB4', 'IGNF:GUADFM49GEO', 'IGNF:RGF93CC47', 'EPSG:4688', 'EPSG:3946', 'EPSG:27561', 'EPSG:32622', 'EPSG:32706', 'IGNF:RGF93CC50', 'EPSG:4625', 'EPSG:4463', 'EPSG:27562', 'IGNF:IGN63UTM7S', 'EPSG:32630', 'EPSG:27574', 'EPSG:3297', 'IGNF:WALL78UTM1S']
Traceback (most recent call last):
  File "xxx.py", line 21, in
    map.wmsimage("https://wxs.ign.fr/choisirgeoportail/geoportail/r/wms?SERVICE=WMS&VERSION=1.3.0&CRS=EPSG:4326", layers=["GEOGRAPHICALGRIDSYSTEMS.PLANIGNV2"], verbose=True, format='jpeg')
  File "PythonPython39libsite-packagesmpl_toolkitsbasemap__init__.py", line 4350, in wmsimage
    img = wms.getmap(service='wms',bbox=(xmin,ymin,xmax,ymax),
  File "PythonPython39libsite-packagesowslibmapwms111.py", line 262, in getmap
    u = openURL(base_url, data, method, timeout=timeout or self.timeout, auth=self.auth)
  File "PythonPython39libsite-packagesowslibutil.py", line 206, in openURL
    raise ServiceException(req.text)
owslib.util.ServiceException:

  bbox invalide

Hors ligne

 

Pied de page des forums

Powered by FluxBB