Oracle SDO_INTERSECTION在SRID4326下有交集却返回NULL问题

核心原因 #

当多边形包含镂空（内环）时，Oracle的SDO_GEOM.SDO_INTERSECTION在地理坐标系（SRID4326）下，可能因拓扑计算的精度限制或内环处理逻辑问题返回NULL，即便SDO_GEOM.RELATE能检测到二者存在交互。

可行解决方案 #

1. 临时转换投影坐标系计算，再转回地理坐标系 #

投影坐标系（如对应区域的UTM分带坐标系）的拓扑运算精度更高，适合处理带内环的多边形：

-- 先将几何图形转换为投影坐标系（示例用UTM 33N，SRID=32633）
WITH projected_geoms AS (
    SELECT
        SDO_CS.TRANSFORM(geom1, 32633) AS proj_geom1,
        SDO_CS.TRANSFORM(geom2, 32633) AS proj_geom2
    FROM your_table
    WHERE id = your_id
)
-- 在投影坐标系下计算交集，再转回SRID4326
SELECT
    SDO_CS.TRANSFORM(
        SDO_GEOM.SDO_INTERSECTION(proj_geom1, proj_geom2, 0.005),
        4326
    ) AS intersection_geom
FROM projected_geoms;

注意：需根据数据所在区域选择对应投影SRID，确保图形完全落在该投影有效范围内；0.005为容差参数，可依据数据精度调整。

2. 修复镂空多边形的拓扑结构 #

若镂空多边形存在微小拓扑错误（如内环与外环过近、顶点重合），可通过SDO_UTIL.REPAIR_GEOMETRY修复后再计算：

-- 修复几何后计算交集
SELECT
    SDO_GEOM.SDO_INTERSECTION(
        geom1,
        SDO_UTIL.REPAIR_GEOMETRY(geom2, 0.0000001),
        0.0000001
    ) AS intersection_geom
FROM your_table
WHERE id = your_id;

SRID4326下单位为度，容差需用极小数值（如0.0000001度，约对应1厘米），匹配数据精度。

3. 拆分镂空多边形为简单多边形再计算 #

将带内环的多边形拆分为多个无内环的简单多边形，分别计算交集后合并结果：

-- 拆分带内环的多边形为简单多边形集合
WITH split_geoms AS (
    SELECT SDO_UTIL.POLYGONTOSURFACE(geom2) AS simple_geom
    FROM your_table
    WHERE id = your_id
)
-- 计算每个简单多边形与geom1的交集，再合并
SELECT
    SDO_AGGR_UNION(
        SDOAGGRTYPE(
            SDO_GEOM.SDO_INTERSECTION(geom1, simple_geom, 0.0000001),
            0.0000001
        )
    ) AS intersection_geom
FROM your_table, split_geoms
WHERE id = your_id;

注意事项 #

• 所有方案中的容差参数需与数据精度匹配，避免过度简化或计算错误。
• 转换投影坐标系时，必须确认数据所在区域对应的正确投影SRID，防止图形变形或计算偏差。

内容的提问来源于stack exchange，提问作者Manny K SoSo