为什么 PostGIS 将所有点堆叠在一起？ （使用 ST_DWithin 查找半径 1000m 内的所有结果）

Question

PostGIS/PostgreSQL 的新手...任何帮助将不胜感激！

我在 postgres 数据库中有两个表，别名分别为 gas 和 ev。 我正在尝试选择一个特定的加油站 (gas.site_id=11949) 并在 1000m 半径内找到所有 EV/替代燃料充电站。 但是，当我运行以下命令时，PostGIS 会返回许多在地图中堆叠在一起的 ev 站（请参见屏幕截图）。

任何人都知道为什么会这样？ 如何让 PostGIS 可视化指定加油站 1000m 半径内的点？

with myplace as (
  SELECT gas.geom 
  from nj_gas gas 
  where gas.site_id = 11949 limit 1)
select myplace.*, ev.* 
from alt_fuel ev, myplace 
where ST_DWithin(ev.geom1, myplace.geom, 1000)

Answer 1

函数ST_DWithin不使用geometry类型参数计算以米为单位的距离。

从文档：

对于几何：距离以几何空间参考系统定义的单位指定。 为了使此函数有意义，源几何图形必须具有相同的坐标投影，具有相同的 SRID。

因此，如果您想以米为单位计算距离，则必须使用数据类型geography ：

对于地理单位以米为单位，测量默认为 use_spheroid=true，为了更快检查，use_spheroid=false 沿球体测量。

话虽如此，您必须转换几何的数据类型。 除此之外，您的查询看起来还不错-考虑到您的数据是正确的:-)

WITH myplace as (
  SELECT gas.geom 
  FROM nj_gas gas 
  WHERE gas.site_id = 11949 LIMIT 1)
SELECT myplace.*, ev.* 
FROM alt_fuel ev, myplace 
WHERE ST_DWithin(ev.geom1::GEOGRAPHY, myplace.geom::GEOGRAPHY, 1000)

样本数据：

CREATE TABLE t1 (id INT, geom GEOGRAPHY);
INSERT INTO t1 VALUES (1,'POINT(-4.47 54.22)');

CREATE TABLE t2 (geom GEOGRAPHY);
INSERT INTO t2 VALUES ('POINT(-4.48 54.22)'),('POINT(-4.41 54.18)');

询问

WITH j AS (
  SELECT geom FROM t1 WHERE id = 1 LIMIT 1)
SELECT ST_AsText(t2.geom) 
FROM j,t2 WHERE ST_DWithin(t2.geom, j.geom, 1000);

     st_astext      
--------------------
 POINT(-4.48 54.22)
(1 Zeile)

Answer 2

您正在交叉加入这些表，并让 PostgreSQL 在选择myplace.* & ev.*时返回两者的笛卡尔积。

因此，虽然myplace只有一行， myplace它的geom将与alt_fuel每一行合并（即结果集将在两者的每种可能组合中包含两个表的所有列）； 由于结果集因此有两个几何列，您的客户端应用程序可能会选择第一个或名为geom （而不是alt_fuel.geom1 ）的列来显示！

我看不出你对结果集中的myplace.geom感兴趣，所以我建议运行

WITH
  myplace as (
    SELECT gas.geom 
    FROM   nj_gas gas 
    WHERE  gas.site_id = 11949
    LIMIT  1
  )

SELECT ev.* 
FROM   alt_fuel AS ev
JOIN   myplace AS mp
  ON   ST_DWithin(ev.geom1, mp.geom, 1000)  -- ST_DWithin(ev.geom1::GEOGRAPHY, mp.geom::GEOGRAPHY, 1000)
;

如果出于某种原因，您还想将myplace.geom与站一起显示，则必须在myplace上使用SELECT *对上述内容进行UNION[ ALL] ； 请注意，您还必须在该SELECT与alt_fuel.*相同的列列表和结构（相同的数据类型！）（或者更好， UNION[ ALL]的另一侧）！

---

请注意@JimJones关于单位的建议； 如果您的数据不是在基于米的 CRS 中投影的（而是在地理参考系统中；“LonLat”），请使用GEOGRAPHY ST_DWithin转换让ST_DWithin将输入视为米（并使用球体代数而不是平面（欧几里得）进行计算） ！

Answer 3

通过使用解决：

WITH
  myplace as (
    SELECT geom as g
    FROM   nj_gas  
    WHERE  site_id = 11949 OR site_id = 11099 OR site_id = 11679 or site_id = 480522
  ), myresults AS (
SELECT ev.* 
FROM   alt_fuel AS ev
JOIN   myplace AS mp
  ON   ST_DWithin(ev.geom, mp.g, 0.1))
select * from myresults```

Thanks so much for your help @ThingumaBob and @JimJones ! Greatly appreciate it.