SQL Server Spatial多面多点匹配求助：判断点是否在多边形内

嘿，作为刚接触SQL Server Spatial的新手，碰到这种百万级点和多边形匹配的效率坑太正常了！你的核心问题是原来的查询用了交叉连接（笛卡尔积），导致100万点×800多边形直接生成8亿条中间记录，不仅结果不符合需求（你只需要每个点是否在任意多边形内），还把数据库拖垮了。下面给你几个实用的解决方案：

一、用EXISTS子查询实现布尔判断（最直接） #

EXISTS是半连接逻辑，它会为每个点去查找是否存在包含它的多边形，一旦找到第一个匹配项就停止查询，完全避免了全量笛卡尔积的灾难，效率提升非常明显。

1. 直接查询结果 #

SELECT 
    points_id,
    CASE 
        WHEN EXISTS (
            SELECT 1 
            FROM [polydb] p2 
            WHERE p1.GEOM.STWithin(p2.GEOM) = 1
        ) THEN 'yes' 
        ELSE 'no' 
    END AS results
FROM [pointsdb] p1;

2. 给points库新增布尔字段并更新 #

如果需要把结果持久化到表中，可以先加字段再批量更新：

-- 新增布尔字段（BIT类型比字符串更高效）
ALTER TABLE [pointsdb] ADD IsInPolygon BIT;

-- 批量更新字段值
UPDATE p1
SET IsInPolygon = CASE 
    WHEN EXISTS (
        SELECT 1 
        FROM [polydb] p2 
        WHERE p1.GEOM.STWithin(p2.GEOM) = 1
    ) THEN 1 
    ELSE 0 
END
FROM [pointsdb] p1;

二、必须加空间索引！（性能提升的关键） #

百万级数据下，没有空间索引的空间查询基本都是灾难。SQL Server需要空间索引来快速筛选出可能包含目标点的多边形，而不是逐个遍历所有800个多边形。

给多边形库的GEOM字段创建空间索引：

CREATE SPATIAL INDEX SIX_polydb_GEOM 
ON [polydb] (GEOM)
USING GEOMETRY_AUTO_GRID;

注：GEOMETRY_AUTO_GRID是SQL Server推荐的自动网格索引，适合大多数场景。如果你的多边形分布有特定规律，也可以手动指定网格层级，但自动网格足够用了。

另外要检查两个库的空间参考系（SRID）是否一致！如果pointsdb和polydb的GEOM字段SRID不同，STWithin会返回NULL，导致判断错误。可以用下面的语句检查：

-- 检查points库的SRID
SELECT TOP 1 GEOM.STSrid FROM [pointsdb];
-- 检查多边形库的SRID
SELECT TOP 1 GEOM.STSrid FROM [polydb];

如果不一致，需要用STTransform转换其中一方的SRID，比如把点转换为多边形的SRID：

WHERE p1.GEOM.STTransform(4326).STWithin(p2.GEOM) = 1 -- 假设多边形SRID是4326

三、如果需要知道具体落入哪个多边形（可选） #

如果之后你不仅需要布尔判断，还想知道点属于哪个多边形，可以用OUTER APPLY只返回第一个匹配的多边形（避免重复记录）：

SELECT 
    p1.points_id,
    CASE WHEN p2.poly_id IS NOT NULL THEN 'yes' ELSE 'no' END AS results,
    p2.poly_id -- 可选，显示匹配的多边形ID
FROM [pointsdb] p1
OUTER APPLY (
    SELECT TOP 1 poly_id 
    FROM [polydb] p2 
    WHERE p1.GEOM.STWithin(p2.GEOM) = 1
) p2;

额外小技巧 #

• 测试时可以先用TOP 1000点验证逻辑，比如SELECT TOP 1000 points_id...，避免一次性跑百万数据卡壳。
• 查看执行计划，确认空间索引是否被使用（执行计划中会显示Spatial Index Seek），如果没用到，可能需要检查索引是否正确创建，或者查询语句是否有问题。

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