Le portail francophone de la géomatique

Bonjour,

Je souhaite retourner le point d'une polyligne en entrée ("PathLines"/Multilinestring) le plus proche d'un point en entrée ("s17_pt") et obtenir la mesure en kilomètre de ce point en sortie par rapport à l'origine de la polyligne.

Pour cela :
J'utilise ST_Line_Locate_Point avec comme entrée les alias foo.the line et foo2.points des requêtes

(SELECT "the_geom" AS points FROM "s17_pt" ) as foo_2,


(SELECT ST_LineMerge(

    ST_Union( "Geometry" )
    )
    
    as the_line FROM "PathLines" ) as foo

ST_Union est utilisé pour que l'ensemble de la géométrie de la polyligne soit dans un enregistrement. Ceci afin que le calcul de ST_Line_Locate_Point s'opère non pas sur chaque segment de la polyligne mais sur l'ensemble.

ST_LineMerge quant à elle convertit la géométrie de multi à simple linestring.
En effet, bien qu'ayant la version 1.5.3 de Postgis et L'info "Availability: 1.5.1 - support for MULTILINESTRINGS was introduced." dans ST_Line_Locate_Point, il me faut une linestring pour que cela marche...

(SELECT sum(ST_Length(
  ST_GeographyFromText(AsText("Geometry"))))/1000 As total_length FROM "PathLines") as foo3

Retourne la longueur totale de la polyligne ou de l'ensemble des segments.

SELECT (ST_Line_Locate_Point(foo.the_line, foo_2.points))* foo3.total_length as kp

Multipliée à la fraction entre 0 et 1 retournée par ST_Line_Locate_Point, j'obtiens la mesure kilométrique de chaque point projeté sur la polyligne.

Mais les valeurs sont trop importantes et l'écart entre la valeur souhaitée et celle obtenue s'accroit plus on s'éloigne de l'origine de la polyligne.

geographyFromText semble fonctionner car la valeur de sum(ST_Length(...)) est correcte.

Je pense qu'il y a peut être une erreur de projection.

Les points s'affichent bien dans QGIS après avoir été importés dans pgsql. Pour la route (polyligne) je ne sais pas car c'est du multilinestring et QGIS ne les lit pas.

La projection est Mercator srid 4326

Merci

Dernière modification par Squid (Fri 20 January 2012 11:33)

Merci Nicolas,

Mes données sont mondiales, avez-vous une suggestion pour la reprojection?

Dois-je utiliser ST_Transform()?

Merci

_ST_BestSRID(geography)

Elle n'est cependant pas vraiment faite pour être utilisée en dehors des "bidouilles" interne de traitement des types geography. Mais ça marche.

Vincent

Bonjour,
concernant le rendu en multilinestring de st_linemerge, j'ai récemment trouvé la solution. Il faut ordonner les objets selon le code de route et le pk en entrée, puis faire un ST_collect et pas un st_union avant de lancer un linemerge..
Si les objets sont bien connexes, on obtient alors bien des LINESTRING, puis des LINESTRING M avec le Add_mesure..

Voili voilou.. ça mériterait une fonction toute faite dans postgis du type make_route('monschema', 'matable', 'monid_route', 'monchamppk', 'multilineOutputYesNo').. et qui envoie des sets d'erreurs pour les objets renvoyant des multilinestring  , ou un bon blog pour documenter ça..

Régis

Finalement ça marche si je sélectionne le tronçon proche des points (3.5% de la route totale) et que j'utilise la zone UTM qui va bien.

Merci à vous.

J'ai plus qu' à développer une fonction qui découpe ma route en tronçons suivant les zones UTM en pensant aux extrémités...

Je mettrai mon code lorsque j'aurai terminé

Squid

Une version provisoire de calculer des KP. Elle n'est valable que si l'on reste dans la zone couverte par le SRID.

CREATE OR REPLACE FUNCTION calib2(table_name_line text,table_name_point text,column_line_geometry text,column_point_geometry text, srid integer) RETURNS SETOF double precision AS $$
DECLARE
    rec record;
    ma_var double precision;
BEGIN

FOR rec IN EXECUTE ('SELECT (ST_Line_Locate_Point(foo.the_line, foo_2.points))* foo3.total_length as kp
FROM
(SELECT ST_Transform(' ||column_point_geometry|| ', '||srid||') AS points FROM ' || table_name_point || ' ) as foo_2,
(SELECT 
ST_LineMerge(
    ST_collect( ST_Transform( foo4.geom , '||srid||' ) ) )as the_line FROM (SELECT ' ||column_line_geometry|| ' as geom FROM ' || table_name_line || '  )as foo4)as foo  ,
(SELECT sum(ST_Length(
  ST_GeographyFromText(AsText(' ||column_line_geometry|| '))))/1000 As total_length FROM ' || table_name_line || ' ) as foo3;') LOOP
    ma_var:=rec.kp;
    RETURN NEXT ma_var;
  END LOOP;

END;
$$ LANGUAGE plpgsql;

BEGIN;
SELECT calib2(quote_ident('votre_nom_de_relation_route') ,quote_ident('votre_nom_de_relation_point'),quote_ident('votre_nom_de_colonne_geometrie_route'),quote_ident('votre_nom_de_colonne_geometrie_point'), votre_SRID_integer );

Merci de m'indiquer si vous voyez des améliorations possibles.

D'autre part je réfléchis sur la possibilité d'utiliser cette base de script avec plusieurs SRID suivant la zone UTM où l'on se trouve. (_ST_BestSRID(geography))

Merci.

Dernière modification par Squid (Tue 24 January 2012 11:55)