Le portail francophone de la géomatique

Bon,

visiblement je mélange deux choses,

j'ai pensé que je pourrais faire une boucle car ma requête ne fonctionnait pas en me renvoyant l'erreur précédente.

il se trouve que la requête fonctionne, donc c'est bien un problème de géométrie.
dans l'erreur précédente : Code:

NOTICE:  TopologyException: found non-noded intersection between 
527822 2.56006e+006, 528042 2.55985e+006 
and 527363 2.55994e+006, 527824 2.56006e+006 527823 2.56006e+006

ERROR:  GEOS union() threw an error!

********** Erreur **********

ERROR: GEOS union() threw an error!
État SQL :XX000

527822 2.56006e+006 sont des coordonnées lambert.

je suis allé voir et il y a effectivement une erreur à cet endroit. j'en ai aussi plein d'autre que j'ai découvertes avec la requête :

SELECT pacage, geomunion (the_geom) AS the_geom
FROM rpg76 
WHERE isvalid(the_geom)
GROUP BY pacage;

merci à Xavier-Pierre sur Développez.com.

voilà où ça en est, si ça intéresse quelqu'un, n'hésitez pas à vous manifester.
PS: je ne sais pas s'il est possible de déplacer cette discussion car elle n'a plus rien a faire dans les boucles.

merci à tous,

cordialement,

nicolas.

Bonjour Nicolas,

en effet, j'ai fait la requête isvalid (the_geom) et il m'a renvoyé une quinzaine d'erreur que je suis allé corrigé "ala mano" dans mon logiciel d'information géographique.

mais vous proposez st_isvalid(). quelle est la différence entre st_isvalid et isvalid.

merci pour les ressources, je savais pas que la fonciton st_buffer(the_geom, 0.0) corrigeait les géométrie.

merci pour votre aide, ça m'a été très utile.

Bonjour,

et bien merci pour les explications,

c'est enrichissant.

comment savoir ça sur les fonctions ?

peut-être ai-je survolé de trop haut la doc.

encore merci,

cordialement,

nicolas

et la doc en français sur postgis.fr !!

Y.

Bonjour à tous,

je reprend ce post pour illustrer par image et pourquoi pas relancer la discussion

récemment, j'ai fusionner toutes les parcelles d'une section afin de récupérer le polygone de la section.

 SELECT section.insee,  geomunion(section.the_geom)
   FROM section
  WHERE isvalid(section.the_geom)
  group by section.insee, section.the_geom;

j'avais cependant un polygone invalide (venant pourtant d'une section vectorisée par le cadastre):

j'ai essayé st_buffer(the_geom,0.0) et rien ne fonctionnait
[img]http://www.atelier-cuesta.fr/polygone_a_corriger.gif[/img]

j'ai fait st_snaptogrid(the_geom,1.0) et mon polygone à trou n'en avait plus.
[img]http://www.atelier-cuesta.fr/polygone_corrige.gif[/img]

je me disais que je pouvais peut-être faire les deux en même temps dans une même requête j'ai essayé :

 SELECT section.insee, st_snaptogrid((geomunion(section_yd.the_geom)),1.0)
   FROM section
  WHERE isvalid(section.the_geom)
  group by section.insee, section.the_geom;

mais je suis peut-être un peu optimiste ! ou un peu niais ?!

voilà, pour illustrer par image et vous remercier par la même occasion pour toutes vos interventions.

cordialement,

Nicolas

bonjour Nicolas,

En fait, la requête précédente fonctionne, mais ne me fait pas l'union des géométries,

j'ai essayé de faire le snaptogrid tout seul et le résultat n'est pas convaincant (le 'trous' sont plus grands)

j'ai essayé geomunion tout seul et rien à faire.

Et puis finalement, j'ai fait le geomunion suivi du st_buffer(the_geom,1.0)
cela fonctionne et à comblé les trous alors que st_buffer(the_geom,0.0) ne faisait rien.

j'ai du me trompé la première fois. je n'arrives pas à obtenir le même résultat en faisant d'abord l'union et puis le snaptogrid.

j'ai un peu du mal à m'y retrouver avec tous les tests que j'ai fait. il faut que je reprenne.

par contre la requête suivante :

 SELECT section_yd.commune, st_buffer((geomunion(section_yd.the_geom)),1.0)
   FROM section_yd
  WHERE isvalid(section_yd.the_geom)
  group by section_yd.commune, section_yd.the_geom;

ne fait pas l'union.
pour l'instant j'arrive à avoir ce que je veux en plusieurs étapes.

je ne sais pas si je fais vraiment évoluer le post là, je suis désolé.

en tout cas merci pour les infos,

à bientôt j'espère,

cordialement,