PostGIS -ST_SHORTESTLINE与LAT LON不正确

Question

我正在使用PostGis st_shortestline来计算一条线和点之间的最短线（一些米）：

SELECT ST_AsText(
        ST_ShortestLine(ST_GeomFromText('POINT(2.33123610021 48.87902639841)', 4326),
                ST_GeomFromText('LINESTRING ( 2.33122725689 48.87902421718, 2.33123229444 48.87901190847)', 4326))
        ) As sline;

我得到的结果似乎不连贯，给定的行不是最短的线：

LINESTRING(2.33123610021 48.87902639841,2.331227760998549 48.87902298544515)

这是一幅图纸结果，使用Mercator投影（JOSM）。

可以解释这一点？

Answer 1

如果您依靠眼睛来确定绘制线是否是最短的线条，那么您可能会误导这一结论。 ST_SHORTESTLINE将返回一条与st_distance完全相同的行，这是两个几何的最小2D笛卡尔距离。这正是正在发生的事情：

WITH j (line,point) AS (
  VALUES ('SRID=4326;POINT(2.33123610021 48.87902639841)',
          'SRID=4326;LINESTRING(2.33122725689 48.87902421718, 2.33123229444 48.87901190847)')
)
SELECT   
  ST_Length(ST_ShortestLine(point,line)), -- length of the shortest line
  ST_Distance(line,point),                -- distance between 'point' and 'line'
  ST_AsEWKT(ST_ShortestLine(point,line))  -- the shortest line as EWKT 
FROM j;

       st_length       |      st_distance      |                                       st_asewkt                                        
-----------------------+-----------------------+----------------------------------------------------------------------------------------
 9.010592472335791e-06 | 9.010592472335791e-06 | SRID=4326;LINESTRING(2.331227760998549 48.87902298544515,2.33123610021 48.87902639841)
(1 row)

也许您分享了您期望的结果，我们可以从那里去。