Discussion :

[little_kevin]

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

[gpain]

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)

[little_kevin]

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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;

...{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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

[gpain]

Bonjour,

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