4. 矿山评价区边界确定与遥感影像预处理

NDVI时序评价的第一步并不是下载影像，而是把“评价谁”这件事在空间上钉死。矿山生态修复的评价对象通常不是一个完整的行政县，也不是一整景卫星影像，而是一个具有明确权属、责任主体和修复任务的采矿权范围或修复规划单元。边界一旦含糊，后续所有NDVI均值、趋势斜率、恢复等级统计都会失去可比性——同一片复垦草地被算进两个评价单元，或者把周边未扰动的天然林地误纳入评价区，都会让“植被恢复趋势”这一结论失真。因此本节要解决两个前后衔接的问题：评价区边界如何确定，以及影像进入计算前必须完成哪些预处理。

一、评价区边界的三种确定路径

实际工作中，评价区边界来源主要有三类，各有适用条件，也常组合使用。

采矿权边界是最具法定效力的依据。可从采矿许可证拐点坐标或自然资源部门提供的采矿权范围矢量数据（常见为Shapefile或GeoJSON）直接获得。它的优点是权属清晰、责任明确，适合以“矿山企业修复责任履行”为评价目标的场景。但要注意：采矿权范围往往包含工业场地、道路、选矿厂等非植被恢复目标区域，直接整块统计NDVI会稀释植被信号，通常需要从中剔除建设用地或保留为单独子单元。

修复规划图适合已编制矿山地质环境恢复治理与土地复垦方案的矿山。规划图上的复垦区、边坡治理区、排土场平台、沉陷区等分区，本身就是按修复措施划分的，直接作为评价单元比采矿权边界更贴合“修复效果”这一评价目的。使用时应将规划图矢量化，并核对坐标系与影像是否一致。

高分辨率影像勾绘用于缺少矢量资料或边界已发生实际变化的情形。以亚米级或米级影像（如高分二号、资源三号或无人机正射影像）为底图，人工或半自动勾绘评价单元边界。勾绘时应遵循“同质单元”原则：同一单元内土地利用类型、修复措施和地形条件尽量一致，例如把复垦耕地、复垦林草地、裸露边坡分别成图。勾绘精度建议控制在影像2～3个像元以内，并保留勾绘底图与操作记录，便于复核。

二、评价单元矢量化的关键控制点

边界确定后要落成可计算的矢量面。几个容易出错的细节值得强调：其一，统一坐标系。遥感影像常为UTM或Albers等投影坐标系，而采矿权拐点可能是地理坐标（经纬度），必须统一到同一投影坐标系后再叠加，否则会出现数十米乃至上百米的偏移。其二，检查拓扑。相邻评价单元之间不应重叠，也不应留缝隙，可用GIS的拓扑检查或“融合—再分割”方式处理。其三，为每个单元赋予唯一编号和属性字段，如unit_id、area_ha、repair_type（复垦耕地/复垦林草地/边坡/建设用地），这些字段将贯穿整个时序统计流程。

三、遥感影像预处理的四道工序

边界是“空间框架”，预处理则是让影像像素值具备物理可比性。NDVI由近红外与红光波段反射率计算，若不做辐射与大气处理，不同时相影像的NDVI差异可能主要来自大气状况和传感器状态，而非植被真实变化。四道工序如下。

辐射定标把影像的DN值转换为辐亮度或表观反射率。以Landsat 8 OLI为例，使用元数据中的RADIANCE_MULT_BAND_x与RADIANCE_ADD_BAND_x：

L = DN × RADIANCE_MULT_BAND_x + RADIANCE_ADD_BAND_x
ρ_TOA = (π × L × d²) / (ESUN_x × cosθ)

其中d为日地距离（天文单位），θ为太阳天顶角，ESUN为波段大气外太阳辐照度。多数处理软件（如ENVI、Google Earth Engine、Sen2Cor配套工具）可自动读取元数据完成，但需确认输出是辐亮度还是表观反射率，二者不可混用。

大气校正把表观反射率还原为地表反射率，消除气溶胶、水汽散射影响。常用方法有FLAASH、6S模型和暗像元法。对NDVI时序评价而言，大气校正的意义在于让不同日期影像的NDVI处于同一物理基准上。若研究区跨越多个季节且气溶胶变化大，建议统一采用基于辐射传输模型的大气校正，而不是简单暗像元减法。

几何精校正保证多时相影像与评价单元矢量严格套合。Landsat、Sentinel-2等产品通常已做过正射校正，但仍有数米到十余米的几何残差。可用同一区域的高分辨率正射影像或控制点做二次精校正，控制点应均匀分布、数量不少于15～20个，校正后中误差建议控制在1个像元以内。几何误差会直接表现为评价单元边缘的“混入像元”，对面积较小的复垦地块影响尤其明显。

裁剪按评价单元矢量对影像做掩膜裁剪。裁剪时建议向外缓冲1～2个像元再裁，避免边缘像元因重采样被截断；统计时再按单元内部有效像元计算，并记录每单元的有效像元数与云掩膜比例，作为后续数据质量判断依据。

四、结果如何判断是否可用

预处理完成后不能直接进入趋势计算，应先做三项检查。第一，几何套合检查：将裁剪后的影像与评价单元矢量叠加，观察单元边界是否与影像地物吻合，重点看复垦地块边缘是否出现明显错位。第二，数值合理性检查：抽取复垦林草地单元，其地表反射率在近红外波段应明显高于红光波段，对应NDVI一般在0.4～0.8之间（示例参考值，实际随季节和植被类型变化）；若某时相NDVI普遍偏低且红光反射率异常升高，往往提示大气校正不足或云污染。第三，有效像元检查：统计每个单元在各时相的有效像元比例，若某单元长期低于70%，说明该单元面积过小或云掩膜过严，应考虑合并单元或改用更高分辨率数据。

把边界与预处理做扎实，后续的NDVI时序统计才有意义。一个常见的工程教训是：评价单元边界与影像存在系统性偏移时，NDVI趋势曲线看起来“平滑合理”，但实际统计的可能是道路或裸地边缘，结论完全不可信。因此建议在正式计算前，保留一份边界叠加底图与预处理参数记录表，作为整个评价过程可追溯的依据。