Le portail francophone de la géomatique

Bon j'ai pas mal fait de test sans succès.
J'ai réglé le problème du doublon dans les tables, certainement une table temporaire qui n'est pas effacée.
J'ai essayé de travailler directement avec du wkt, mais la couche n'est toujours pas parfaite.
Voici l'état d'avancement de la commande :

CREATE TABLE pt_move AS
SELECT DISTINCT ON (pt_id)
    pt_id,
    ln_id, 
    ST_AsEWKT(ST_LineInterpolatePoint(
            ln_geom, 
            ST_LineLocatePoint(ln_geom, pt_geom)
        )
        )
   
FROM
(
    SELECT 
        ln.geom AS ln_geom,
        pt.geom AS pt_geom, 
        ln.gid AS ln_id, 
        pt.gid AS pt_id
    FROM
        pt, 
        ln 
    WHERE 
        ST_DWithin(pt.geom, ln.geom, 10.00) 
    ORDER BY
pt_id
) AS subquery;

ALTER TABLE pt_move
ADD geom geometry(point,32719);
UPDATE pt_move SET geom=st_asewkt;
ALTER TABLE pt_move
DROP st_asewkt;

Merci pour votre aide,
Franck

Salut,

les calculs d'interpolation étant effectuées avec des nombres
à virgule flottante, il se peut que des erreurs d'arrondi surviennent,
car elles sont inhérente à l'arithmétique les concernant.

Selon la précision que vous utilisez pour valider la topologie de vos données,
il y aura probablement des erreurs, même si elles ne sont pas visibles "à l'oeil nu".

Les coordonnées générées par les fonctions d'interpolation de postgis ne sont elles-même non
"consistentes" avec ses propres prédicats topologiques, par exemple :

    WITH foo as (SELECT 'LINESTRING(0 0 , 10 99)'::geometry g)
    SELECT st_intersects(ST_LineInterpolatePoint(g,0.15) ,g) 
    FROM foo

renvoie faux.

Pour la partie postgis, il est apparemment prévu d'ajouter dans les prochaines versions
la possibilité de choisir une tolérance pour les prédicats.
Pour le moment, on peut s'en sortir en snappant au préalable les géomètres à une grille .

Pour QGIS, si c'est le plugin de vérification topologique que vous utilisez, vous pouvez
dans la fenêtre principale régler la tolérance à appliquer .

Dernière modification par tumasgiu (Mon 12 February 2018 17:08)

La façon classique de travailler avec ces données est d'utiliser st_dwithin(geom1, geom2, dist), plutot que les prédicats OGC st_touches, st_intersects, etc.
Pas de souci particuliers avec du récursif.
St_dwithin permet de simuler cette précision à laquelle on veut travailler (précision des données géo).

Sinon, vous pouvez construire une topologie Postgis, où là les relations entre objets seront explicites.

Comme le disait tumasgiu, "vraiment parfaite" ca n'existe pas en informatique lorsque les réels sont représentés en double précision: il y a tjs des arrondis. (Par ex un point se trouvant a 1/3 de la longueur d'une ligne).
La question est plutot: "vraiment parfaite à la précision près".

Nicolas

Dernière modification par Nicolas Ribot (Mon 12 February 2018 16:21)