Le portail francophone de la géomatique

Bonjour,
nous avons été confronté au même problème et nous n'avons pas de solution directe avec les produits Esri.
Suite à nos tests et recherches voici un script python nécessitant ogr et arcpy qui semble règler le problème des géométries incompatibles Esri/Ogr

import sys
import arcpy
import ogr
import os

#Chemin du repertoire contenant le lot des fichiers shape a analyser et traiter
#pathDirShape= r"d:/......."
#ici le script est parametre pour fonctionner avec un script dans ARCGIS PRO ou Arcgis Desktop si on installe ogr.

pathDirShape= arcpy.GetParameterAsText(0)
## dev
#recuperation de tous les fichiers shape du repertoire
arcpy.env.workspace = pathDirShape
listShp = arcpy.ListFeatureClasses('*.shp')

#boucle sur les fichiers shape decouverts
for nameShape in listShp:
    arcpy.AddMessage(nameShape)
    #reconstitution du chemin du fichier shapre en cours d'ananlyse
    pathShape = os.path.join(pathDirShape,nameShape)
    # creation d'un curseur de mise a jour
    cursShp = arcpy.UpdateCursor(pathShape)
    #boucle sur chaque les entites du shape en cours d'analyse
    for cur in cursShp:
        # recuperation de la geometrie
        shapeV =cur.getValue("shape")
        # recuperation de le geometreie en wkt pour d'interoperabilite avec ogr
        wktStr = shapeV.WKT
        # creation d'un objet geometrique au sens ogr
        objGeom = ogr.CreateGeometryFromWkt(wktStr)
        # test de la validite au norme ogr
        if not objGeom.IsValid():
            # si invalide 
            arcpy.AddWarning("""
            ERREUR DE GEOMETRIE
                COUCHE      : %s
                OBJECTID    : %i """%(nameShape,cur.getValue("FID")))
            # comme indique dans la doc de geoportail on applique un buffer à 0
            objGeom =objGeom.Buffer(0) 
            # si l'objet est alors valide suit au passage du buffer on met a jour
            # la geometrie du shape
            if objGeom.IsValid():
                featwkt = arcpy.FromWKT(objGeom.ExportToWkt())
                cur.setValue("shape",featwkt)
                cursShp.updateRow(cur)
                arcpy.AddWarning(""">>> Reparation effectue <<<""")
            else:
                arcpy.AddWarning(""">>><<< Reparation NON effectuee >>><<<""")
            
        del shapeV
        del objGeom
        del wktStr
    
    if "cursShp" in globals() :del cursShp
    if "cur" in globals() : del cur

pour l'instant nous sommes en phase de test mais cela semble passer au detecteur de mensonge du geoportail

bon courage

Bonjour,

Je reviens à la charge sur ce post.

Je passe un temps fou à corriger mes saisies de zonage de PLU(i) pour que les fichiers soient propre et valide lors du dépôt sur le GPU.
J'utilise Qgis (3.34) et l'un des 2 formats tolérés par le GPU : shp
J'utilise bien évidemment la version la plus récente des cadastres que j'importe via l'extension "cadastre".

Je me suis rendu compte de plusieurs choses problématiques dans la saisie de polygones :
- le fichier parcelle du cadastre est loin d'être "propre" (des noeuds en doublons, des polygones avec des artefacts qui ne seraient pas validés par le GPU), aussi, lorsque j'utilise l'outil d'accrochage automatique aux sommets pour saisir par exemple mon zonage sur le tracé des parcelles... mon tracé se retrouve avec ces mêmes défauts;
-  sur ma propre couche "zonage" : qgis ajoute (sans bien sur que je ne lui demande !) des sommets qui n'existe pas sur la couche parcelle, de sorte que mes polygones limitrophes ne sont plus jointifs
- quand j'ai fini de saisir et que je corrige la topologie de mon fichier, il arrive quasi tout le temps que lorsque je supprime un sommet en doublon, qgis en recréé plein d'autres (ouinnnn !); pareil quand je coupe un polygone en 2, bien souvent des nouveaux sommets en doublon sont créés :-(

Je me doute qu'il faut faire attention aux projections :
- je télécharge le cadastre dans sa version "système RGF 93, projection Lambert 93"
- mes couches zonage and co sont également créées dans ce système de projection
- les orthophotos que je mets dessous aussi (me semble-t-il !)

Mes questions :
- pensez-vous que de saisir mes données dans le format GPKG puis d'exporter ensuite mes couches en shp limiterait ces problèmes topologiques ?
   --> cf discussion shp / gpkg ici : https://georezo.net/forum/viewtopic.php … sie+edigeo
- y'aurait-il quelque part dans qgis une option pour limiter une "fausse précision" des couches géographique : faire par exemple que les sommets de mes polygones jointifs se positionnent se superposent exactement quand bien même le cadastre n'aurait pas ses polygones jointifs au 0,001 cm ! Si les coordonnées de mes polygones sont ajustées même à 10 cm prés, si cela permet une superposition parfaite, cela me suffit largement.
- vaut-il mieux travailler sur le cadastre dans sa version "base de données", ou comme je le fais (c'est peut-être un des problèmes), en important les couches au format shp ou gpkg ?

Bien cordialement.
Elise