分区统计把栅格像元汇总到空间单元

分区统计以一个区域集合为分组边界，对栅格值计算计数、总和、均值、最小最大值或类别面积。行政区内的降水均值、地块内的植被指数、流域内的高程分布都属于这种问题。操作看似只是“区域内取平均”，结果却受网格、边界、无数据和面积单位共同影响。

区域可以来自矢量多边形，也可以由另一张分类栅格中的区号表示。两种输入最终都要明确每个像元属于哪个分区以及贡献多少权重。

矢量区域先与栅格网格建立关系

矢量多边形与栅格像元边界通常不重合。最简单规则按像元中心是否落入多边形决定归属，执行快且结果明确，但窄小区域可能没有任何中心点，边界附近也会产生阶梯误差。

“所有触碰像元”规则提高覆盖，却会把一个边界像元完整计入多个相邻区，造成面积和总量重复。按像元与区域的覆盖比例加权更接近面积积分，但计算成本更高，还依赖几何精度和等面积语义。

选择哪种规则取决于变量含义。温度均值可以按覆盖面积近似，点事件栅格的计数则不能随意分割。分类像元代表主导类别时，用分数面积拆分也未必反映像元内部真实组成。

区域很小时，应报告参与像元数和有效覆盖比例。一个地块只有一两个粗像元时，均值虽然可算，空间代表性却很弱。

坐标参考决定面积和相交是否可靠

矢量与栅格必须位于可比较的坐标参考中。临时重投影区域或栅格都会改变边界与网格，需要记录采用的目标 CRS、分辨率和重采样方法。

经纬度网格的像元面积随纬度变化。若统计总面积、总量或面积权重，不能把像元个数乘一个固定平方度。可以使用等面积投影、逐像元面积模型或测地计算。

重投影连续栅格会插值数值，分类栅格通常使用最近邻或面积保持方法。为了分区统计而多次重投影会积累误差，最好选定统一分析网格，并保存其来源和变换参数。

跨日期比较必须使用一致区域和网格。行政边界修订后，同名区域的统计变化可能来自边界，不是栅格现象。边界版本应与结果一起保存。

无数据不是数值零

云、传感器缺测、覆盖范围外和质量掩膜通常表示为无数据。均值计算应排除无数据，同时报告有效像元数和有效面积；把它们当作零会系统性降低结果。

不同波段或日期的有效像元集合可能不同。比较两个均值时，差异可能来自采样位置变化。严格比较可以使用共同有效区域，或至少报告各期覆盖率。

有些变量中零是合法观测，例如无降水或无植被指数信号。文件元数据、质量层和业务规则要共同确定无数据，不能仅凭数值大小猜测。

区域内全部无数据时，结果应为空并携带原因，而不是返回零。零意味着观测到的量为零，会影响下游汇总和趋势判断。

统计量必须符合变量语义

连续变量可以计算均值、标准差和分位数；类别变量更适合频数、比例和主导类别。对土地覆盖编码求平均通常没有解释意义，即使软件可以计算。

总和是否有意义取决于像元值。每像元总量可以直接求和，单位面积密度需要乘有效面积，点估计或插值场则可能需要积分模型。必须从产品文档确认值代表瞬时、累计、密度还是类别。

极值容易受到孤立噪声影响，分位数或截尾统计可能更稳定，但也会改变问题定义。输出不应只给一个均值，至少保留计数、有效面积和必要的离散程度。

类别面积可通过每类像元面积累加。边界权重、像元面积和分类误差都会进入结果。若分类图有精度评估，可以进一步报告面积不确定性，而不是把地图计数视为无误差真值。

重叠区域需要明确是否允许重复计算

矢量分区可能相互重叠，例如保护区与行政区。每个区域独立统计时，同一像元可以计入多个结果，这是合法的多组汇总；若结果要合计为总体，重复就会造成超额。

互斥分区应先验证拓扑，检查重叠和空洞。按优先级消解重叠会改变区域面积，规则必须可追溯。多部件和带孔多边形也应正确处理。

栅格区号通常每个像元只能有一个值，适合互斥分区。转换矢量到区号栅格时，细小或重叠区域的处理顺序会决定最终归属。

大规模计算要按数据布局优化

逐区域裁剪并读取整张大栅格会产生大量重复 I/O。可以按栅格块读取、建立区域空间索引，并只处理相交候选。云优化格式和内部瓦片有助于范围读取，但前提是访问窗口与存储布局匹配。

区域非常多时，先栅格化区号再做一次分组累计可能更快；需要精确覆盖权重时则要保留矢量相交。性能优化不能悄然改变边界规则。

并行分块时，同一区域可能跨多个块，需要合并可结合的中间统计。均值不能直接平均各块均值，应累计有效权重和加权总和；方差也需要正确的合并公式。

验证从简单几何与守恒关系开始

可以构造与像元完全对齐的矩形，手工计算计数和均值，验证中心点、触碰和权重规则。再测试窄多边形、孔洞、跨越无数据和区域重叠。

面积型结果应检查分区面积之和与有效覆盖面积是否一致。按类别汇总时，各类面积加无数据面积应接近区域总面积，差异要能由边界近似解释。

生产结果需要记录栅格版本、区域版本、分析网格、像元纳入规则、无数据规则和统计单位。分区统计不是一个脱离语义的聚合函数，而是把连续网格映射到管理单元的测量过程。

资料依据

• OSGeo GDAL Project：GDAL Documentation，https://gdal.org/en/stable/
• Open Geospatial Consortium：Simple Feature Access，https://www.ogc.org/standards/sfa/