MULTILINESTRING 附近的 rgeo 点

Question

有一个几何列wkb_geometry，srid 4326 是MULTILINESTRING 我想确定这些记录中的哪些在几何对象POINT 的预定距离内（比如5000m）

虽然下面的方法允许确定一个多边形是否包含一个点

def self.containing_latlon(lat,lon, polygon)
  ewkb = EWKB.generate(RGeo::Geographic.simple_mercator_factory.point(lon, lat).projection)
  where("ST_Intersects(polygon, ST_GeomFromEWKB(E'\\\\x#{ewkb}'))")
end

ST_Intersects 显然不是一个选项，因为它适用于“空间的任何部分然后它们相交”。
我还没有找到文档来确定一条线是否在一个点的 X 距离内。但问题可能反过来了？如果问题不是由 MULTILINESTRING 和缓冲区定义的多边形内的点。

上面的方法需要如何修改才能执行？

Answer 1

请改用ST_DWithin。

对于使用米的距离，将参数转换为geography，例如5公里：

SELECT * FROM t
WHERE ST_DWithin('POINT(7.00 51.82)'::geography,geom::geography,5000);

如果您对 SRS 的测量单位感到满意，请坚持使用“几何”

SELECT * FROM t
WHERE ST_DWithin('POINT(7.00 51.82)'::geometry,geom,42);

WKT 文字后面的 :: 是一种用于转换数据类型的 postgres 语法。但是按照 postgres 的惯例，有很多方法可以做同样的事情。以下示例使用不同的技术将 WKT 文字转换为 geometry：

SELECT 
  CAST('SRID=4326;POINT(1 2)' AS geometry),
  'SRID=4326;POINT(1 2)'::geometry,
  ST_GeomFromText('SRID=4326;POINT(1 2)'),
  ST_SetSRID(ST_MakePoint(1,2),4326);

-[ RECORD 1 ]---+---------------------------------------------------
geometry        | 0101000020E6100000000000000000F03F0000000000000040
geometry        | 0101000020E6100000000000000000F03F0000000000000040
st_geomfromtext | 0101000020E6100000000000000000F03F0000000000000040
st_setsrid      | 0101000020E6100000000000000000F03F0000000000000040

延伸阅读：Getting all Buildings in range of 5 miles from specified coordinates