postgresql 必须用 using gist 创建空间索引,因b-tree无法处理几何类型;mysql spatial索引受限于引擎和字段类型;sqlite需手动维护rtree;st_dwithin失效常因srid或类型不匹配。
postgresql 的 create index 为什么加了 using gist 才生效
空间索引不是普通 b-tree 能处理的——地理数据是二维(甚至三维)结构,point、polygon 这些值没法直接比大小。postgresql 默认的 create index 用的是 b-tree,对 geometry 列直接建索引会报错:data type geometry has no default operator class。
必须显式指定空间索引类型:USING GIST(Generalized Search Tree)是 PostgreSQL 官方支持的空间索引方法,它把几何对象映射成可比较的 bounding box,再组织成树形结构。
-
CREATE INDEX idx_locations_geom ON locations USING GIST (geom);—— 正确写法,geom是GEOMETRY或GEOGRAPHY类型列 - 别写成
USING BTREE,哪怕加了ST_AsText(geom)转字符串也不行——失去空间关系语义,查询ST_Within时照样全表扫 - 如果列是
GEOGRAPHY类型(推荐用于经纬度真实距离计算),GIST同样适用,但内部用球面算法,性能略低于GEOMETRY+ 平面投影
MySQL 8.0+ 的 SPATIAL 索引只支持 MyISAM 和 InnoDB,但有硬限制
MySQL 的 SPATIAL 索引看起来简单:ALTER TABLE places ADD SPATIAL INDEX(geom);,但它背后卡着两个关键前提:存储引擎和字段类型。
不满足任一条件就会静默失败或报错:ERROR 1464 (HY000): Table storage engine for 'places' doesn't support SPATIAL indexes(常见于 Aria、Memory 引擎),或者 ERROR 1170 (42000): BLOB/TEXT column 'geom' used in key specification without a key length(字段类型不对)。
- 必须用
MyISAM或InnoDB(MySQL 5.7+ 支持 InnoDB 的 SPATIAL);MyISAM支持更早,但不支持事务 -
geom字段必须是POINT、LINESTRING、POLYGON等具体地理类型,不能是GEOMETRY通用类型(MySQL 不支持) - 索引只能建在单个地理字段上,不能联合其他普通字段;也不能在表达式上建,比如
ST_Transform(geom, 4326)不行
ST_DWithin 查询没走空间索引?检查 GEOMETRY vs GEOGRAPHY 和 SRID 是否一致
写了 WHERE ST_DWithin(geom, ST_Point(-73.9857, 40.7484), 1000) 却发现执行计划里是 Seq Scan,大概率是类型或坐标系不匹配。
PostgreSQL 中,ST_DWithin 对 GEOMETRY 默认按“投影单位”算距离(比如 Web Mercator 的米),而 GEOGRAPHY 才按真实“米”计算。如果 geom 列是 GEOMETRY 但没设 SRID,或设了 SRID=4326 却没用 ST_SetSRID 显式标注,优化器会放弃使用索引。
- 确认列定义:
geom GEOMETRY(POINT, 4326)或geom GEOGRAPHY(POINT, 4326)—— SRID 必须明确且与查询中一致 - 查询里用
ST_SetSRID(ST_MakePoint(-73.9857, 40.7484), 4326),而不是裸写ST_Point(后者不带 SRID) - 如果列是
GEOMETRY且 SRID=4326,又想按真实距离查,别硬套ST_DWithin,先用ST_Transform投影到米制坐标系(如 EPSG:3857),否则索引失效且结果不准
SQLite / SpatiaLite 的 rtree 索引需要手动维护,且不自动覆盖更新
SpatiaLite 用的是 R*Tree(rtree 表),它不是传统意义上的“索引”,而是一张独立的虚拟表,得自己建、自己同步。
执行 INSERT INTO mytable VALUES (..., GeomFromText('POINT(-73.9857 40.7484)')); 后,rtree_mytable_geom 不会自动更新——查询 SELECT * FROM mytable WHERE id IN (SELECT id FROM rtree_mytable_geom WHERE minx -73.9) 可能漏掉刚插入的数据。
- 建索引要两步:
SELECT CreateSpatialIndex('mytable', 'geom');创建rtree表,再SELECT UpdateLayerStatistics('mytable', 'geom');初始化统计信息 - 每次增删改后,必须显式调用
INSERT INTO rtree_mytable_geom VALUES (...)或用触发器维护(SpatiaLite 提供CreateRTriggers辅助函数) -
rtree只加速范围查询(ST_Contains,ST_Intersects的 bbox 部分),不加速距离函数;想用ST_Distance排序,仍需额外排序开销
空间索引本身不难建,难的是让每一步都对得上类型、SRID、引擎、维护方式——少一个环节,就回到全表扫描。










