Couches EFS non ajoutables dans QGIS 3.28

[SimonSAMPERE]

Détecté ici

Avec les service https://carto.isogeo.net/server/rest/services/scan_services_1/EMS_EFS_WMS_WFS/FeatureServer/ on peut seulement ajouter les couches dans QGIS 3.22. Pour info c'est la même version de Python (3.9.5) pour QGIS 3.22 et QGIS 3.28.
J'ai rien trouvé à ce sujet dans les différents Changelogs :

• 3.24
• 3.26
• 3.28

Le code suivant permet effectivement d'ajouter une couche EFS dans QGIS 3.22.* mais pas dans QGIS 3.28.4 ni dans QGIS 3.28.8 :

# -*- coding: utf-8 -*-
from qgis.core import QgsDataSourceUri, QgsProject, QgsVectorLayer
from owslib.wfs import WebFeatureService

wfs_url_base = "crs='EPSG:2154' filter='' url='https://carto.isogeo.net/server/rest/services/scan_services_1/EMS_EFS_WMS_WFS/FeatureServer/2' table='' sql=''"
lyr = QgsVectorLayer(wfs_url_base, baseName="efs_arcgisserver", providerLib="arcgisfeatureserver")
QgsProject.instance().addMapLayer(lyr)

En utilisant l'Outil de débogage/développement, je me suis rendu compte que QGIS requête un URL plus complexe en indiquant la liste des Id des entités à requêter. On note également que les paramètres correspondant au CRS sont différentes. Par exemple pour la couche "ISOGEO_TEST.Un_joli_polygone" qui contient une seule entité :
https://carto.isogeo.net/server/rest/services/scan_services_1/EMS_EFS_WMS_WFS/FeatureServer//11/query?f=json&objectIds=1&inSR=2154&outSR=2154&returnGeometry=true&outFields=*&returnM=false&returnZ=false

Et apparemment le morceau de code suivant fonctionne dans QGIS 3.28 :

# -*- coding: utf-8 -*-
from qgis.core import QgsDataSourceUri, QgsProject, QgsVectorLayer
from owslib.wfs import WebFeatureService

wfs_url_base = "crs='EPSG:2154' url='https://carto.isogeo.net/server/rest/services/scan_services_1/EMS_EFS_WMS_WFS/FeatureServer//11/'"
lyr = QgsVectorLayer(wfs_url_base, baseName="efs_arcgisserver", providerLib="arcgisfeatureserver")
QgsProject.instance().addMapLayer(lyr)

Ressources :

• https://gis.stackexchange.com/questions/422207/load-layers-with-pyqgis-from-arcgis-feature-service
• https://gis.stackexchange.com/questions/456167/pyqgis-add-arcgis-feature-service-layer-to-qgis-including-a-query
• REST server not loading in 3.28 Firenze but work fine in previous versions qgis/QGIS#51409

J'ai créé ce ticket dans le dépôt GitHub de QGIS, en espérant avoir plus d'infos : qgis/QGIS#53736
J'ai également posté une question sur GIS StackExchange : https://gis.stackexchange.com/questions/463079/creating-qgsvectorlayer-from-arcgis-feature-service-layer-using-python-doesnt-w

[SimonSAMPERE]

Finalement il suffisait d'enlever filter='', table='' et sql='' de la variable wfs_url_base, j'ai donc répondu à ma question sur GIS StackExchange et commenté mon ticket dans le dépôt GitHub de QGIS.

[SimonSAMPERE]

Recette

Vérifier qu'on peut bien ajouter les couches issues d'EFS en utilisant le plugin Isogeo dans QGIS 3.22 et QGIS 3.28.

[SimonSAMPERE]

https://github.com/isogeo/isogeo-plugin-qgis/releases/tag/3.5.0