5

How to find the closest line vertex from a linestring? ST_ClosestPoint will give a point which is not from the linestring.

RyanKDalton
  • 23,068
  • 17
  • 110
  • 178
user32876
  • 81
  • 3
  • Could you draw a sketch about your case with the linestrings and mark which point you would like to find. Attach also the SQL query you have tried. – user30184 Aug 26 '14 at 08:17

1 Answers1

3

ST_ClosestPoint will find the closest point on a linestring. For example, locating a point on a road that is closest to a point-of-interest off to one side:

SELECT ST_AsText(ST_ClosestPoint(road, poi))
FROM (SELECT
  'LINESTRING (10 10, 20 20)'::geometry AS road,
  'POINT (14 10)'::geometry AS poi
) AS f;
  st_astext
--------------
 POINT(12 12)
(1 row)

There are theoretically an infinite number of points that are on a line, and this location is exactly on the road linestring:

SELECT ST_Intersects('LINESTRING (10 10, 20 20)'::geometry, 'POINT (12 12)'::geometry);
 st_intersects
---------------
 t
(1 row)

If the objective is to get the closest vertex or coordinate from the linestring, the road needs to be converted to a MultiPoint geometry. Here's a function that does that:

CREATE OR REPLACE FUNCTION ST_AsMultiPoint(geometry) RETURNS geometry AS
'SELECT ST_Union((d).geom) FROM ST_DumpPoints($1) AS d;'
LANGUAGE sql IMMUTABLE STRICT COST 10;

Then use it similar to before:

SELECT ST_AsText(ST_ClosestPoint(ST_AsMultiPoint(road), poi))
FROM (SELECT
  'LINESTRING (10 10, 20 20)'::geometry AS road,
  'POINT (14 10)'::geometry AS poi
) AS f;
  st_astext
--------------
 POINT(10 10)
(1 row)
Mike T
  • 42,095
  • 10
  • 126
  • 187