Postgis在某些查询期间断开连接

Question

我有两个表，一个在SRID4326中（包含MULTIPOLYGON），另一个在SRID32188中（包含multilestring）。

我试图从t1中选择所有距离t2中的任何行100米以内的多边形。

这是我提出的问题：

SELECT * FROM table1 AS t1, table2 AS t2 
WHERE ST_DWithin(
        ST_Transform(t1.geom, 3857),
        ST_Transform(t2.geom, 3857),
        100);

不幸的是，它不起作用，而且pgAdmin4没有给我一个有用的错误。相反，只要我运行查询，它就会与数据库断开连接（因此这不是超时情况）：

我还尝试了在同一SRID中使用两个表，但仍然存在问题。。。我做错了什么？

编辑：为了尝试缩小问题的范围，我采用了@JGH的建议，并将其调整为针对点测试ST_DWITHIN，而不是表2的多行字符串，它工作正常。

SELECT * FROM table1 t1
WHERE EXISTS (
    SELECT FROM table2 t2
    WHERE ST_DWithin(t1.geom, ST_MakePoint(-73.63,45.52)::geography,100)
);

然而，我使用只包含一条直线的表2进行了测试（尽管仍然是多行线几何体），客户端崩溃。是否与此相关：https://gis.stackexchange.com/questions/31923/why-is-the-database-connection-lost-on-every-query ?

Answer 1

你的客户可能在计时。

实际上，您正在两个表之间进行交叉连接，因此如果t1中的多边形位于t2中10个要素的100m范围内，则将返回10次。您可能希望从t1中获得不同的多边形。要执行此操作，请仅在t1上选择，并在t2上应用的条件满足时保留该行。

第二个问题——不幸的是，这不会引发任何错误——是使用3857计算距离是错误的。在魁北克省，长度畸变约为40%，这是巨大的。对于这两种几何形状，您必须使用合适的局部投影（如MTM或UTM），或者，如果不可能，则转换为地理来进行距离计算。

由于您确实将两个几何转换为新的CRS，因此不能仅使用几何上的空间索引，并且距离计算将花费很长很长的时间，导致您的客户端超时。

SELECT *
FROM   table1 t1
WHERE  EXISTS (
   SELECT FROM table2 t2
   WHERE  ST_DWithin(t1.geom::geography, ST_Transform(t2.geom, 4326)::geography, 100)
   );

不要忘记在地理位置（和/或转换后的几何图形）上创建空间索引

CREATE INDEX t1_geog ON t1 USING gist(geography(t1.geom);

或

CREATE INDEX t1_32188 ON t1 USING gist(st_transform(geom,32188));