为什么PostGIS中ST_DWithin查询无法使用已创建的GIST索引？

• 坐标系单位不匹配是核心问题：你使用的SRID=4326（WGS84地理坐标系）的Geometry类型默认单位是度，而ST_DWithin传入的第三个参数62815.14是米为单位的距离，相当于要求查询六万多度范围内的点，覆盖范围接近半个地球，PostgreSQL优化器判断走索引过滤的行数占比极高，反而不如直接全表扫描效率高，因此选择了并行顺序扫描。
• 之前的多边形查询正常走索引，是因为多边形坐标用的是和坐标系匹配的度单位，范围合理，过滤后返回行数少，符合索引生效的条件。

方法1：转换为Geography类型适配米单位（无需修改存量数据） #

直接将查询中的Geometry对象强转为Geography类型，ST_DWithin对Geography类型默认使用米作为距离单位，修改后的查询条件如下：

SELECT distinct alias
FROM position
WHERE
date > '2021-11-28T19:26:18.574Z' AND date < '2021-11-28T20:26:18.574Z'
AND
ST_DWithin (
  point::geography,
  ST_GeomFromText('POINT (14.32983409613371 49.91815471231952)',4326)::geography,
  62815.14152820495
);

修改后查询范围会被修正为62公里左右，过滤后行数大幅减少，优化器会自动选择GIST索引扫描。

方法2：使用投影坐标系（性能更优，适合固定业务区域） #

如果你的业务数据集中在固定区域，可以将point字段转换为对应区域的UTM投影坐标系（单位为米），之后直接使用原生Geometry类型的ST_DWithin查询即可，性能比Geography类型高15%~30%。

可选优化：创建时空复合索引 #

如果你的查询几乎都会同时携带时间范围过滤，可以借助btree_gist扩展创建时空复合索引，进一步提升查询效率：

-- 启用btree_gist扩展
CREATE EXTENSION IF NOT EXISTS btree_gist;
-- 创建时空复合索引
CREATE INDEX position_date_point_idx ON position USING gist (date, point);

内容的提问来源于stack exchange，提问作者Jurij Hladuk