Le portail francophone de la géomatique

Bonjour Ted,

je ne réponds pas à la question "Que fais-je de mal" mais une piste de solution ici (st_node  + st_polygonize) :

https://georezo.net/forum/viewtopic.php … 78#p250778

Mathieu

Par contre si je veux récupérer mon id_triangle, je fais comment.
J'ai beau ajouter l'entrée, il me dit qu'il faut l'ajouter dans la clause group by??

Ted

Merci Mathieu pour ta réponse, mais le résultat n'est pas positif à 100%. De nombreux polygones ne récupèrent pas l'id du bon triangle ;-(.

C'est quand même étrange, qu'avec la requête avec le With, on ne puisse pas introduire cette information.

Je vais poursuivre mes recherches.

a+

ST_WITHIN n'est peut-être pas l’opérateur spatial adéquat ?

Bonjour,

J'ai réussi en "trichant" un peu. J'ai utilisé des buffers négatifs (vers l'intérieur) des objets et ça fonctionne. Peut-être pas propre, mais ça me convient.
Voici la requête :

create table test_interpo."Trpolygo_Update" as
select distinct a.id_polygone,b.id_triangle,round(st_zmin(a.geom)::decimal, 1) as Zmin_Polygone, round(st_zmin(b.geometryexplose)::decimal, 1) as Zmin_triangle, a.geom
from test_interpo.Trpolygo a, test_interpo."7123_TRIANGLE_DISTINCT_test_interpo" b where st_intersects(st_buffer(b.geometryexplose,-0.0000001),a.geom)

Merci pour votre aide

Ça manque de schéma :-) mais je me suis loupé dans ma réponse : la table bnd est une table de lignes, pas de polygones... il faudrait faire l'update des triangles avec la table de polygones initiale.

je serai preneur d'un jeu de données test pour mieux comprendre le truc.

Merci je veux bien tester moi aussi. Ça me remettra le fonctionnement de st_node() et st_polygonize() en tête.

Bonsoir,

En pièce jointe, des shp pour tester (linéaire et surfacique).

Avant la requête de Nicolas, voici la requête pour construire mes polygones:

create table test_interpo.Trpolygo as 
with bnd as (
    select  id_triangle, geom from test_interpo."7123_ligne_segmentize_Tr2" union all
    select id_triangle,geom from test_interpo.line_intertrianglestep2
), col as (
    select st_collect(geom) as geom from bnd
) select (st_dump(st_polygonize(st_node(geom)))).geom as geom
    from col;

Puis la requête pour récupérer les id_triangle :

alter table  test_interpo.Trpolygo add column id_polygone bigserial;

--Récupération de l'id du triangle DELAUNAY dans chaque polygone et du minZ du Triangle et MinZ du Polygone
create table test_interpo."Trpolygo_Update" as
select distinct a.id_polygone,b.id_triangle,round(st_zmin(a.geom)::decimal, 1) as Zmin_Polygone,
 round(st_zmin(b.geometryexplose)::decimal, 1) as Zmin_triangle, a.geom
from test_interpo.Trpolygo a, test_interpo."7123_TRIANGLE_DISTINCT_test_interpo" b
 where st_intersects(st_buffer(b.geometryexplose,-0.0000001),a.geom);

J'obtiens 77 polygones.
@ Nicolas, j'ai executé ça requête en la mettant à jour, mais le résultat n'est pas bon pour les polygones (seulement 52, il est peut-être nécessaire de faire un st_snaptogrid?), par contre on récupère bien pour chaque polygone, l'id_triangle:

create table test_interpo.tmp_Ribot as
with tmp as (
  SELECT st_exteriorRing(geometryexplose) AS geom
  FROM test_interpo."7123_TRIANGLE_DISTINCT_test_interpo"
  UNION ALL
  SELECT geom
  FROM test_interpo.line_intertrianglestep2
), tmp1 as (
  select st_collect(geom) as geom
  from tmp
), tmp2 as (
  select st_node(st_union(geom)) as geom
  from tmp1
), tmp3 as (
    SELECT st_dump(st_polygonize(geom)) AS d
    FROM tmp2
), tmp4 as (
  select (d).path[1] as id, (d).geom
  from tmp3
), tmp5 as (
  select t.id, t.geom, tr.id_triangle
  from tmp4 t join test_interpo."7123_TRIANGLE_DISTINCT_test_interpo" tr on st_contains(tr.geometryexplose, st_pointOnSurface(t.geom))
) SELECT *
  FROM tmp5;

les noms de mes fichiers sont des noms temporaires, car je teste sur une petite portion de mon jeu de données.

Merci

Dernière modification par Ted (Wed 29 March 2017 18:43)

Re Bonjour,

J'ai refait toute une batterie de tests et je ne comprends pas le problème.
J'ai fait des tests en limitant le nombre d'objets :
Avec 10 polygones  tout fonctionne
Avec 20 polygones, certains nouveaux polygones sont bien construits d'autres ne le sont pas.
Sur la totalité, j'ai des objets construits, mais qui ne ressemblent à rien (polygones "aléatoires").
J'ajoute des shape de mes données (échantillon) avec les problèmes.
J'utilise cette requête (fourni par Nicolas Ribot) pour créer mes polygones.

 create table "CTRL_SONDE_7123"."7123_Trg_polygo_non_int_Delau_sel" as 
with bnd as (
    select  id_triangle, geom from "CTRL_SONDE_7123"."7123_ligne_segmentize_Delau_Tr_sel" union all
    select id_triangle,geom from "CTRL_SONDE_7123"."7123_line_intertriangle_Delau_step1_sel"
), col as (
    select st_collect(geom) as geom from bnd
) select (st_dump(st_polygonize(st_node(geom)))).geom as geom
    from col;

Si quelqu'un pouvait m'aider ça serait cool ;-)

Je n'ai pas pu insérer mes copies d'écran.

D'avance merci

Dernière modification par Ted (Thu 11 May 2017 18:56)

Bonjour et merci de ton aide tumasgiu,

J'ai testé ta commande en l'aménageant comme ceci :

WITH t AS (
select p.id_triangle, array_agg(l.geom)||st_exteriorring(p.geom) as t
from "CTRL_SONDE_7123"."7123_TRIANGLE_non_intersect_Delau" p, "CTRL_SONDE_7123"."7123_line_intertriangle_Delau_step1_sel" l
WHERE st_intersects(p.geom, l.geom) group by  p.id_triangle,p.geom),
u AS (SELECT id_triangle, st_dump(st_node(st_collect(t))) d FROM t)

SELECT (st_dump(st_polygonize((d).geom ORDER BY (d).path))).geom FROM u
GROUP BY id_triangle;

Point positif, le calcul est beaucoup plus rapide qu'avec la précédente.
Point négatif, il me manque des résultats (voir zip avec captures d'écran).

D'où cela vient-il? Pourquoi certains polygones sont construits et pas d'autres?
Ne faudrait-il pas utiliser un st_snaptogrid pour une tolérance?

Merci de votre aide

Dernière modification par Ted (Fri 12 May 2017 13:35)

Bon je sais je suis un peu obnubilé par la "chose", même un dimanche matin ;-).

Je viens de réaliser un test en isolant les segments interieurs d'un polygone (ceux qui servent à découper).

J'ai relancé la requête fourni par tumasgiu avec mes aménagements (voir post plus haut) et je n'obtiens pas le même résultat que sur la globalité des données. Il doit y avoir un problème d'intersection, mais le st_snaptogrid ne le résout pas!

Des pistes? pour m'aider!

Bonne fin de week end

Ted