Dans cet article nous verrons plus sp\u00e9cifiquement comment d\u00e9tecter et corriger les probl\u00e8mes g\u00e9om\u00e9triques dans une table PostgreSQL\/Postgis en utilisant des requ\u00eates SQL.<\/p>\n
V\u00e9rification des erreurs avec QGis.<\/h2>\n
Vous pouvez d\u00e9tecter rapidement les erreurs g\u00e9om\u00e9triques en utilisant QGis. Il suffit de charger la couche et de lancer le menu Vecteur -> Outils de g\u00e9om\u00e9trie -> V\u00e9rifier la validit\u00e9<\/strong><\/em> La m\u00e9thode GEOS est plus rapide mais n’indique que la premi\u00e8re erreur rencontr\u00e9e pour chaque objet. Elle renvoie un fichier \u2018Sortie invalide<\/em>\u2019 qui est une couche de polygones compl\u00e9t\u00e9e par une colonne\u00a0 _errors<\/em> \u00a0:<\/p>\n Par contre si des erreurs existent je passe directement sur pgAdmin 4 (ou \u00e0 d\u00e9faut dans le gestionnaire de base de donn\u00e9es de QGis) pour travailler en SQL.<\/p>\n La premi\u00e8re \u00e9tape est l’identification des erreurs. Notons que PostGIS<\/a> est conforme au standard OGC OpenGIS Specifications, et donc la validit\u00e9 est v\u00e9rifi\u00e9e par rapport \u00e0 GEOS.<\/p>\n Nous allons donc r\u00e9aliser une requ\u00eate qui recherchera les trois types d’erreurs principaux:<\/p>\n Un dernier type d’erreur, en g\u00e9n\u00e9ral sans cons\u00e9quence, est la pr\u00e9sence de n\u0153uds doubles (deux points de l’entit\u00e9 situ\u00e9s exactement au m\u00eame endroit). Ils ne perturbent pas les op\u00e9rateurs spatiaux, mais ce n’est pas une raison pour ne pas les corriger.<\/p>\n Voici une requ\u00eate qui vous permet de d\u00e9tecter les trois premiers types d’erreurs en un seul passage:<\/p>\n
\n![\"\"](\"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2018\/01\/g1.png?resize=640%2C428\")
<\/a>
\nVous avez le choix entre trois m\u00e9thodes :<\/p>\n\n
![\"\"](\"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2018\/01\/g3.png?resize=640%2C170\")
<\/a>
\nLa m\u00e9thode Qgis peut renvoyer plusieurs erreurs par objet. Cette m\u00e9thode renvoie dans le fichier \u2018Sortie invalide<\/em>\u2019 une couche de polygones compl\u00e9t\u00e9e par une colonne _errors<\/em> avec un message moins normalis\u00e9 que par la m\u00e9thode GEOS :<\/p>\n![\"\"](\"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2018\/01\/g4.png?resize=640%2C238\")
<\/a>Personnellement, j’utilise cet outil pour v\u00e9rifier s’il y a des erreurs ou pas. J’ex\u00e9cute les deux m\u00e9thodes et, si elles ne renvoient aucune erreur je consid\u00e8re que la couche est valide.<\/p>\nD\u00e9tection des erreurs g\u00e9om\u00e9triques avec Postgis<\/h2>\n
\n
<\/div>
\nUNION
\nSELECT ‘geom nulle’ AS nb, count(*) FROM table WHERE the_geom is null
\nUNION
\nSELECT ‘collection’ AS nb, count(*) FROM table WHERE not ST_IsValid(the_geom)
\nAND ST_GeometryType(ST_MakeValid(the_geom))=’ST_GeometryCollection’;<\/div><\/div>\n![\"\"](\"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2018\/01\/g5.png?resize=640%2C352\")
<\/a>Vous pouvez avoir le d\u00e9tail des erreurs GEOS avec la requ\u00eate suivante:<\/p>\n
![\"\"](\"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2018\/01\/g6.png?resize=640%2C455\")
<\/a><\/p>\n