STIntersects必须写在JOIN的ON子句中才能利用空间索引,放WHERE会导致笛卡尔积;需确保类型一致、SRID相同、索引配置合理,并可用STEnvelope预筛选优化性能。

STIntersects JOIN 必须写在 ON 子句里,不能放 WHERE
SQL Server 优化器只在 ON 条件中直接使用 STIntersects 时才考虑空间索引;一旦挪到 WHERE,就会先做笛卡尔积再过滤,百万级数据直接卡死。
常见错误写法:SELECT * FROM points p, polygons z WHERE p.geom.STIntersects(z.geom) = 1
这是全量配对,哪怕建了空间索引也无效。
正确写法必须显式 JOIN:SELECT p.id, z.name FROM points p JOIN polygons z ON p.geom.STIntersects(z.geom) = 1
- 函数返回值必须是布尔表达式,写成
= 1是 SQL Server 的强制要求(不像 PostGIS 返回 true/false) - 不能包裹:如
COALESCE(p.geom.STIntersects(z.geom), 0) = 1会绕过索引 - 不能混用类型:
geometry列和geography列不能直接比较,会隐式转换并丢弃索引
两个 geometry 列必须 SRID 一致且类型匹配
SQL Server 不会自动对齐坐标系。若 p.geom 是 SRID 4326,而 z.geom 是 SRID 0 或 3857,STIntersects 内部会尝试转换,但此时空间索引完全失效。
检查方法:SELECT TOP 1 geom.STSrid FROM pointsSELECT TOP 1 geom.STSrid FROM polygons
- 不一致时,统一转换(推荐在 ETL 阶段完成):
ALTER TABLE points ALTER COLUMN geom TYPE geometry USING geom.STTransform(4326) - 注意:
STTransform要求源列有有效 SRID;若为 0,需先用STSetSRID显式指定:UPDATE points SET geom = geom.STSetSRID(4326) WHERE geom.STSrid = 0 - 确保两列都是
geometry类型——geography上的STIntersects索引支持有限,尤其对多边形
必须建空间索引,且网格配置要合理
SQL Server 空间索引不是“建了就快”,它依赖四级网格划分。默认 MEDIUM(8×8)适合中小数据,但对高密度点或复杂面,常需调高密度。
创建示例:CREATE SPATIAL INDEX SIndx_points_geom ON points(geom) USING GEOMETRY_GRID WITH (GRIDS = (LEVEL_1 = HIGH, LEVEL_2 = HIGH, LEVEL_3 = MEDIUM, LEVEL_4 = LOW), CELLS_PER_OBJECT = 16)
-
HIGH(16×16)提升细粒度覆盖,适合小目标(如建筑轮廓);LOW(4×4)适合大范围粗筛(如省界) -
CELLS_PER_OBJECT控制每个几何对象最多占多少网格单元,太大导致索引膨胀,太小漏匹配 - 务必在建索引前运行
UPDATE STATISTICS points,否则优化器可能误判选择性而跳过索引 - 用
sys.dm_db_index_usage_stats查看索引是否被真正使用:SELECT user_seeks, user_scans FROM sys.dm_db_index_usage_stats WHERE object_id = OBJECT_ID('points')
JOIN 前加边界框预筛选(&& 等效写法)
SQL Server 没有 PostGIS 那样的 && 操作符,但可用 STEnvelope().STIntersects() 模拟边界框快速剪枝——它走空间索引极快,能筛掉 80%+ 无关记录。
写法示例:SELECT p.id, z.name FROM points p JOIN polygons z ON p.geom.STIntersects(z.geom) = 1 WHERE p.geom.STEnvelope().STIntersects(z.geom.STEnvelope()) = 1
- 这个条件本身不保证精确相交,但能大幅缩小参与
STIntersects计算的组合数 - 如果驱动表是面(比如你要查哪些面覆盖某批点),应在点表上加范围约束:
AND p.geom.STWithin(geometry::STGeomFromText('POLYGON((...))', 4326)) = 1 - 避免在
ON中用STBuffer或STUnion——这些高开销函数会让整个 JOIN 退化为嵌套循环
最易被忽略的是:SQL Server 空间索引对 STIntersects 的加速效果高度依赖数据分布。一个面包含十万点,和十万点各自落在不同小面里,执行计划可能完全不同——别只看单条语句的 EXPLAIN,得用真实数据集跑 STATISTICS IO 和实际耗时验证。










