PySpark Sedona空间连接报错：LinearRing点未形成闭合线串求解决

问题背景 #

现有两个Spark DataFrame：

• spatial_wi_df（存储POINT类型数据）结构：

root
|-- timestamp: timestamp (nullable = true)
|-- wagonnumber: long (nullable = true)
|-- latitude: double (nullable = true)
|-- longitude: double (nullable = true)
|-- speed: double (nullable = true)
|-- geometry: geometry (nullable = false)
|-- wi_geometry_meter: geometry (nullable = true)

• spatial_station_groups_gdf（存储POLYGON类型数据）结构：

root
|-- geo_name: string (nullable = true)
|-- polygon: geometry (nullable = false)

执行空间连接判断点是否在多边形内时：

spatial_wi_df.createOrReplaceTempView("points")
spatial_station_groups_gdf.createOrReplaceTempView("geofences")
spatial_join_result = spark_sedona.sql("SELECT g.geo_name, p.wagonnumber FROM points AS p, geofences AS g WHERE ST_Contains(g.polygon, p.geometry)")

触发报错：Points of LinearRing do not form a closed line string，且尝试缓冲区、ST_MakeValid、凸包三种修复方法均未解决。

报错原因 #

该错误表明多边形的外环（LinearRing）首尾坐标不重合，不符合OGC几何规范（多边形必须由闭合的LinearRing构成）。之前的修复方法未精准针对“未闭合”问题处理，因此失效：

• ST_Buffer(0)：依赖有效几何结构，对未闭合环无法自动修复
• ST_MakeValid：用法错误（不应搭配LATERAL VIEW），且部分版本对未闭合环修复能力有限
• ST_ConvexHull：生成凸多边形会改变原形状，且无法强制闭合原环

针对性解决方案 #

1. 定位未闭合的多边形 #

先筛选出存在问题的多边形，方便排查：

SELECT geo_name, polygon 
FROM geofences 
WHERE NOT ST_IsClosed(polygon)

2. 强制闭合LinearRing并修复多边形 #

通过提取外环、追加首尾点、重构多边形的方式，强制闭合未合格的几何结构：

CREATE OR REPLACE TEMP VIEW fixed_geofences AS
SELECT 
    geo_name,
    CASE 
        -- 对未闭合的多边形，提取外环后将起点追加到末尾，重新构造多边形
        WHEN NOT ST_IsClosed(polygon) THEN 
            ST_Polygon(
                ST_LineString(
                    ST_Append(ST_ExteriorRing(polygon), ST_StartPoint(ST_ExteriorRing(polygon)))
                )
            )
        -- 已闭合的多边形直接保留
        ELSE polygon
    END AS polygon
FROM geofences;

3. 验证修复效果 #

确认修复后的多边形符合规范：

SELECT 
    geo_name,
    ST_IsClosed(polygon) AS is_closed,
    ST_IsValid(polygon) AS is_valid
FROM fixed_geofences;

4. 重新执行空间连接 #

使用修复后的多边形视图执行原逻辑：

SELECT g.geo_name, p.wagonnumber 
FROM points AS p 
JOIN fixed_geofences AS g 
ON ST_Contains(g.polygon, p.geometry);

额外建议 #

如果修复后仍有异常，可检查：

• 多边形是否存在自相交问题（用ST_IsValid检测），自相交的多边形需要用ST_MakeValid单独处理
• 几何数据的坐标系是否统一，确保点和多边形使用相同的空间参考系

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