Discussion :

[mathildeg]

J’ai le même problème que la discussion « Accès aux flux WFS Géoportail [Résolu] »

Etant déjà classé en résolu, et ne répondant pas à mon problème, je me permets d’en ouvrir une autre.

Mon entreprise vient d’acheter une licence géoservice de l’API Géoportail de l’IGN seulement j’ai un problème d’affichage pour les couches provenant du flux WFS.
Je travaille sous QGIS (2.4), J’ai constaté que Qgis n’affiche que 1000 objets par couche, donc les couches apparaissent incomplètes et je ne parviens pas à les afficher entièrement.

Comment enlever la limitation du nombre d’objet à afficher sous QGIS ? J’ai déjà cherché dans les options de paramétrage de Qgis (pour reprendre la solution de la discussion similaire ; décocher l’option «garder en mémoire») mais je n’ai rien trouvé.




De plus, depuis quelques temps plusieurs messages apparaissent lorsque je me connecte au flux WFS de l’IGN.
Lors de la connexion ce message apparait ;



1-Sans cliquer ; ce message apparait



2-Après avoir cliqué sur « Annulé » ce message ;



Puis la connexion au flux se ferme
3-Après avoir cliqué sur « OK »

La couche fini par apparaître mais avec seulement avec 1000 objets par couche, exemple (c’est le même problème pour toutes les couches) ;



Si quelqu’un a une idée ?

-pour l’affichage des couches WFS sans limitation du nombre d’objets ?
-concernant les messages lors de la connexion au flux WFS ?

Merci

[lowzonenose]

pour l’affichage des couches WFS sans limitation du nombre d’objets ?

Le flux WFS est limité à 1000 objets.
(Cette limite est un garde fou, il serait suicidaire de demander à des serveurs d'interroger et d'envoyer de gros volume de données...)

L'idée serait peut être de requêter le service sur des petites emprises (ex. Je n'affiche de la donnée vecteur qu'à partir d'un certain niveau de zoom)

Malheureusement, QGis envoie des requêtes sur le service sans lui passer d'emprise. De plus, il n'interroge le service qu'en version 1.0.0...
Donc, pour pallier ce problème, j'utilise un plugin QGis 'WFS 2.0.0 Client' disponible sur le dépôt officiel.
capture d'écran :



Maintenant, si ton besoin est d'avoir toutes les données sur l'ensemble du territoire,
je t'invite à te tourner vers le service commercial...

concernant les messages lors de la connexion au flux WFS ?

Le mode de connexion sur les flux Geoportail est en HTTPS avec login/password.
J'avoue ne pas rencontrer ce type de problème...
As tu le même problème si tu exécutes cette requête dans un navigateur ?

[acarme]

Bonjour,

concernant les messages d'erreurs, il faut voir avec vos experts réseaux, cela vient probablement de votre proxy.

[mathildeg]

Bonjour,

Merci pour ces réponses.

"Le flux WFS est limité à 1000 objets." Cela signifie que le serveur de l'IGN limite de nombre d'objets à 1000 à chaque connexion au flux WFS?

Qgis n'a pas un nombre limité d'objets à afficher provenant de flux WFS ?


Concernant les messages d'erreurs nous n'avons pas de proxy.
"As tu le même problème si tu exécutes cette requête dans un navigateur ?" Je n'arrive pas a afficher les données dans un navigateur.
Voici le résultat que j'obtiens ; <ExceptionReport><Exception exceptionCode="MissingParameter">Resource parameter not found</Exception></ExceptionReport>

[lowzonenose]

"Le flux WFS est limité à 1000 objets." Cela signifie que le serveur de l'IGN limite de nombre d'objets à 1000 à chaque connexion au flux WFS?

oui

Qgis n'a pas un nombre limité d'objets à afficher provenant de flux WFS ?

non, pas de limite que celle du flux...

Concernant les messages d'erreurs nous n'avons pas de proxy.

pas de proxy d'entreprise ? ok...

Je n'arrive pas a afficher les données dans un navigateur.
Voici le résultat que j'obtiens ; <ExceptionReport><Exception exceptionCode="MissingParameter">Resource parameter not found</Exception></ExceptionReport>

envoie nous ta requête pour que nous la testions...

[little_kevin]

Bonjour,

Je me permet de revenir sur ce ticket, au sujet de la limitations des 1000 premières entités, dans Qgis.

C'est assez problématique lorsque nous avons besoins d'une section cadastrale qui comporte plus de 1000 entités.
Par exemple; une feuille dont j'obtiens des numéros cadastraux qui vont de 1 à 2170. Hors, j'ai aussi besoin de numéros qui sont au-delà de 2170.

Comment avoir les numéros au-delà de 2170? Dans ma requête je ne peux pas mettre

Code :

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

...AND "numeros" > '2170'

, le champ n’étant pas numérique cela ne peux pas fonctionner.

Comment être sûr que les numéros de parcelle manquants ne sont qu’absent de la requête mais existe bien ou sont réellement inexistant?

En fonction, de quel ordre recevons-nous les polygones? Ordre d'un OBJECTID (que l'on ne voit pas dans la table attributaire, mais qui serait requêtable), dans l'ordre des numéros de parcelles ou autre?

Merci

[rwieser]

Bonjour,

Ce type de requête est limité par le nombre d'objets maximum autorisé en une requête par nos serveurs. Pour contourner cette limitation, il possible d'utiliser la pagination du standard WFS 2.0.0. Malheureusement ce standard est encore mal implémenté dans QGIS malgré nos remontés sur le sujet (https://github.com/qgis/QGIS/pull/1665)

Une solution éventuelle consisterait à utiliser le service directement pour télécharger la donnée puis de l'ouvrir sous QGIS.

Exemple pour l’obtention de toute les parcelle de la commune de MINVERSHEIM (67)

Étape 1
Réaliser le filtre wfs pour obtenir tout les objets, dans notre cas (code_dep = 67, code_com = 293) :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
<ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">
 <ogc:And>
  <ogc:PropertyIsEqualTo>
   <ogc:PropertyName>code_dep</ogc:PropertyName>
   <ogc:Literal>67</ogc:Literal>
  </ogc:PropertyIsEqualTo>
  <ogc:PropertyIsEqualTo>
   <ogc:PropertyName>code_com</ogc:PropertyName>
   <ogc:Literal>293</ogc:Literal>
  </ogc:PropertyIsEqualTo>
 </ogc:And>
</ogc:Filter>

Documentation sur les filtres OGC disponible [ici|http://portal.opengeospatial.org/fil...fact_id=39968]

Étape 2
Construction de la requête :

URL de base : https://wxs.ign.fr/key/geoportail/wfs
SERVICE : WFS
VERSION : 2.0.0 (pour la pagination)
REQUEST : GetFeature
TYPENAMES : BDPARCELLAIRE-VECTEUR_WLD_BDD_WGS84G:parcelle (nom de la couche à interroger)
FILTER : Filtre constitué précédemment
SRSNAME : EPSG:2154 (Lambert 93)
RESULTTYPE : hits (uniquement le nombre d'objet répondants aux conditions du filtre)

Ce qui donne :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
https://wxs.ign.fr/key/geoportail/wfs?SERVICE=WFS&VERSION=2.0.0&TYPENAMES=BDPARCELLAIRE-VECTEUR_WLD_BDD_WGS84G:parcelle&FILTER=<ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">
 <ogc:And>
  <ogc:PropertyIsEqualTo>
   <ogc:PropertyName>code_dep</ogc:PropertyName>
   <ogc:Literal>67</ogc:Literal>
  </ogc:PropertyIsEqualTo>
  <ogc:PropertyIsEqualTo>
   <ogc:PropertyName>code_com</ogc:PropertyName>
   <ogc:Literal>293</ogc:Literal>
  </ogc:PropertyIsEqualTo>
 </ogc:And>
</ogc:Filter>&REQUEST=GetFeature&SRSNAME=EPSG:2154&resultType=hits

L’exécution de cette requête nous retourne un fichier XML et dans ce fichier, on peut récupérer l'attribut numberMatched contenant le nombre recherché (ici 2121)

Étape 3
Récupération des objets

3 requêtes sont donc nécessaires :
* 0-999
* 1000-1999
* 2000-2121

Pour cela, nous allons utiliser les paramètres :
* COUNT pour spécifier le nombre d'objets voulut pour chaque requête (1000)
* STARTINDEX pour spécifier l'index de départ de la requête (0, 1000, 2000)
* OUTPUTFORMAT=shape-zip pour récupérer des fichiers Shape

Ce qui nous donne les 3 requêtes suivantes :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
https://wxs.ign.fr/key/geoportail/wfs?SERVICE=WFS&VERSION=2.0.0&TYPENAMES=BDPARCELLAIRE-VECTEUR_WLD_BDD_WGS84G:parcelle&FILTER=<ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">
 <ogc:And>
  <ogc:PropertyIsEqualTo>
   <ogc:PropertyName>code_dep</ogc:PropertyName>
   <ogc:Literal>67</ogc:Literal>
  </ogc:PropertyIsEqualTo>
  <ogc:PropertyIsEqualTo>
   <ogc:PropertyName>code_com</ogc:PropertyName>
   <ogc:Literal>293</ogc:Literal>
  </ogc:PropertyIsEqualTo>
 </ogc:And>
</ogc:Filter>&REQUEST=GetFeature&SRSNAME=EPSG:2154&COUNT=1000&STARTINDEX=0&OUTPUTFORMAT=shape-zip

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
https://wxs.ign.fr/key/geoportail/wfs?SERVICE=WFS&VERSION=2.0.0&TYPENAMES=BDPARCELLAIRE-VECTEUR_WLD_BDD_WGS84G:parcelle&FILTER=<ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">
 <ogc:And>
  <ogc:PropertyIsEqualTo>
   <ogc:PropertyName>code_dep</ogc:PropertyName>
   <ogc:Literal>67</ogc:Literal>
  </ogc:PropertyIsEqualTo>
  <ogc:PropertyIsEqualTo>
   <ogc:PropertyName>code_com</ogc:PropertyName>
   <ogc:Literal>293</ogc:Literal>
  </ogc:PropertyIsEqualTo>
 </ogc:And>
</ogc:Filter>&REQUEST=GetFeature&SRSNAME=EPSG:2154&COUNT=1000&STARTINDEX=1000&OUTPUTFORMAT=shape-zip

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
https://wxs.ign.fr/key/geoportail/wfs?SERVICE=WFS&VERSION=2.0.0&TYPENAMES=BDPARCELLAIRE-VECTEUR_WLD_BDD_WGS84G:parcelle&FILTER=<ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">
 <ogc:And>
  <ogc:PropertyIsEqualTo>
   <ogc:PropertyName>code_dep</ogc:PropertyName>
   <ogc:Literal>67</ogc:Literal>
  </ogc:PropertyIsEqualTo>
  <ogc:PropertyIsEqualTo>
   <ogc:PropertyName>code_com</ogc:PropertyName>
   <ogc:Literal>293</ogc:Literal>
  </ogc:PropertyIsEqualTo>
 </ogc:And>
</ogc:Filter>&REQUEST=GetFeature&SRSNAME=EPSG:2154&COUNT=1000&STARTINDEX=2000&OUTPUTFORMAT=shape-zip

Étape 4

Fusionner les ShapeFile obtenus (avec QGIS ou OGR2OGR par exemple)

Solution artisanale certes, mais il n'y actuellement pas mieux pour QGIS...