GeoSpark项目中RS_ZonalStats功能增强：支持全接触像素统计

在空间数据分析领域，栅格数据与矢量数据的交互分析是一个常见需求。GeoSpark作为开源空间计算框架，其RS_ZonalStats函数提供了基于多边形的栅格统计功能，但在处理小面积要素时存在统计盲区。

现有问题分析

当前RS_ZonalStats实现采用严格的空间包含关系判断，只有当栅格像素中心点完全位于多边形内部时才会被纳入统计。这种机制在处理建筑物等小面积要素时会出现统计遗漏，特别是当要素尺寸小于栅格分辨率时，极易导致返回NaN值。

技术解决方案

通过引入all_touched参数可以完善该功能：

1. 当all_touched=false时（默认值），保持现有中心点判断逻辑
2. 当all_touched=true时，采用边界接触判断，任何与多边形边界相交的像素都将被纳入统计

实现原理

该增强功能需要修改栅格-矢量相交判断逻辑：

• 对于每个候选像素，构建其边界多边形
• 使用空间关系谓词判断像素多边形与目标多边形的相交关系
• 根据all_touched参数选择严格包含或边界接触判断模式

应用价值

该改进将显著提升小要素场景下的分析可靠性：

• 建筑物轮廓分析：确保小尺寸建筑的完整统计
• 精细农业应用：准确获取田块边界特征值
• 生态保护评估：完整统计保护区内边缘生境

技术展望

未来可进一步扩展为可配置的接触判断策略：

• 中心点模式（当前默认）
• 全接触模式（本改进）
• 面积占比模式（按覆盖面积比例加权统计）

该改进将使得GeoSpark在精细化空间分析场景中具备更强的竞争力，为使用者提供更灵活、更准确的空间统计能力。