WFS avec Python

Bonjour,

Je souhaite ouvrir le lien WFS sur la DBparcellaire via Python et les bibliothèques ogr ou Owslib.
Mais j'obtiens une erreur 403 sur la première et de certificat sur la seconde...

En fait j'ai une licence Professionnels.ign.fr, pour avoir mes fonds dans Qgis ou Arcgis, est-ce que mon soucis peut venir de là? Est-ce qu'en passant par Python, le webservice IGN attend une cle API?

Merci

Bonjour,

Effectivement, l'utilisation des webservices est conditionnée à la fourniture d'une clé API et des informations d'authentification choisies.
Si vous obtenez une 403, c'est que vous n'avez pas fourni de clé, ou que les informations envoyées avec la clé sont incorrectes.
Pouvez-vous m'envoyer votre clé, ainsi que l'URL et une copie de l'erreur associée à votre 403 svp ? (en MP)

Bonjour,

Comme vous me l'avez indiqué (en message privé), que ma clé est de type *.GIS et dans sans authentification nécessaire, je viens de retenter ma connexion en python. En simple http et non https, elle fonctionne bien;

Code:

---

1 2 3 4 5 6 7 8 9 10

from osgeo import ogr driver = ogr.GetDriverByName('WFS') # le serveur WFS wfs = ogr.Open("WFS:http://wxs.ign.fr/MA_CLE/geoportail/wfs?") couche = wfs.GetLayerByName("BDPARCELLAIRE-VECTEUR_WLD_BDD_WGS84G:parcelle") import json for elem in couche: print json.loads(elem.ExportToJson())

---

Cela m'affiche bien tous les éléments;

Citation:

---

...{u'geometry': {u'type': u'MultiPolygon', u'coordinates': [[[[4.7800638, 47.2974763], ..., u'type': u'Feature', u'properties': {u'feuille': 1, u'section': u'ZA', u'code_com': u'002', u'numero': u'0061'...}
{u'geometry': {u'type': u'MultiPolygon', u'coordinates': [[[[4.776876, 47.2966752], ..., u'type': u'Feature', u'properties': {u'feuille': 1, u'section': u'ZA', u'code_com': u'002', u'numero': u'0060', ...}...

---

Mon soucis est désormais d'apporter une requête, dans ma connexion.
Dans Qgis, j'ai noté que la requête est intégré à la source;

Code:

---

http://wxs.ign.fr/MA_CLE/geoportail/wfs?SERVICE=WFS&VERSION=1.0.0&REQUEST=GetFeature&TYPENAME=BDPARCELLAIRE-VECTEUR_WLD_BDD_WGS84G:parcelle&SRSNAME=EPSG:4326&FILTER= code_dep 39 code_com 240 section ZD

---

J'ai donc tenté dans la variable qui définit ma couche d'y insérer la requête:

Code:

---

couche = wfs.GetLayerByName("BDPARCELLAIRE-VECTEUR_WLD_BDD_WGS84G:parcelle&FILTER= code_dep:'39',code_com:'240',section:'ZD'")

---

Mais là je n'ai plus rien...
A quel moment dois-je ajouter ma requête? Connaissez-vous la commande, svp?

Merci

Bonjour,

Je suppose que vous devez mettre seulement le nom de la couche.
Le filtre doit probablement être fourni avec une autre fonction.