积雪制图利用可见光与短波红外差异识别雪盖

新雪在可见光波段通常反射较强，在短波红外波段吸收明显，这种光谱差异支持归一化差异积雪指数和多条件分类。卫星积雪产品可描述雪盖范围和季节变化，但云、林冠、山地阴影、薄雪和水体会造成缺测或混淆。

积雪范围、积雪比例、雪深和雪水当量是不同变量。光学影像主要提供表面雪盖证据，不能仅凭白色像元直接估计雪深或水资源量。

NDSI 突出雪与多数云和岩土的差异

NDSI 通常使用绿色和短波红外反射率之差除以两者之和。雪在绿色波段亮而在短波红外较暗，因此指数偏高。阈值只是起点，还需结合近红外、可见光反射率、热红外或水体掩膜，减少水、云和明亮裸地混淆。

不同传感器波段响应、太阳角、积雪粒径、杂质和液态水含量会改变指数。固定阈值不应跨传感器和区域直接套用。森林下积雪被冠层遮挡，像元记录的是树冠与雪地混合信号，开阔地精度不能代表林区。

NASA ARSET 的遥感基础材料提供反射率、波段与空间分辨率的背景：https://appliedsciences.nasa.gov/sites/default/files/Fundamentals_RS_Session2_Land_Final.pdf 。具体积雪识别应以产品算法和质量说明为准。

MODIS/Terra MOD10A1 Collection 6.1 提供每日 500 米全球网格积雪产品及 NDSI、质量和算法标志：https://nsidc.org/data/mod10a1/versions/61 。产品以正弦投影分块，并记录云、极夜及其他无法正常判定的情况。

质量标志不是可选附件。云覆盖是光学积雪制图的主要限制，极夜使高纬冬季无法依靠太阳反射观测，薄雪和短暂积雪也可能漏检。分析时应把云、无数据、夜间和无雪分开，不能统一编码为零。

二值雪盖便于面积统计，NDSI 或分数雪盖保留更多混合像元信息。把指数直接转为面积比例需要区域验证；坡地、林冠和传感器点扩散函数都会影响像元内比例。

坡度和坡向控制太阳入射，阴坡反射率低，积雪可能被误判为非雪。地形校正可以减少照度差异，但数字高程误差、复杂散射和阴影区低信噪比限制效果。陡峭山谷还会出现几何遮挡和像元混合。

应在统一投影和像元网格上计算雪面积，并考虑像元实际地表面积。经纬度网格中不同纬度像元面积不同，直接计数会产生偏差。流域统计要处理边界像元和无效观测比例，避免把云多的日期解释为雪骤减。

多源融合可用雷达、微波或气象模型填补光学缺测，但各传感器测量机制和空间尺度不同。融合结果应标识观测来源和估计部分，不能把插补像元与直接观测混为同等证据。

每日产品中的云会打断积雪序列。时间滤波可利用前后日期填补短空洞，但快速融雪或降雪期间，简单向前填充会延迟变化。云隙填补应给出最大时间窗口，并保留是否填补及来源日期。

积雪初日、终日和持续时间对单日误判敏感，可要求连续多日满足条件，或用概率状态模型结合温度和降水。规则会平滑噪声，也会漏掉短暂积雪，选择时要围绕气候监测、水文模拟或交通应用的需求。

传感器和算法版本变化可能造成时间序列断点。长期分析应使用一致集合，比较重叠版本，并按年记录有效观测数量。数据重新处理后，应重新生成派生指标而不是拼接不同版本结果。

高分辨率影像、地面雪尺、自动站和现场照片可作为参考，但它们的时间和空间支持不同。验证样本应覆盖开阔地、森林、城市、阴影、不同海拔和雪线附近，不能只选择易识别的大片雪区。

分类评价要报告雪类与非雪类的用户精度、生产者精度、漏检和误报，并把云导致的不可判定单独统计。面积产品还应比较区域雪盖比例和边界位置。产品的全球总体精度不能替代具体流域和季节的本地验证。

面向水文应用时，应把雪盖产品与高程带、温度、降水和流域边界结合，但这些辅助变量不能反向掩盖观测缺失。模型同化后的完整雪面属于融合估计，应与卫星直接有效观测分别存储，并通过敏感性试验说明插补对融雪径流结论的影响。

积雪制图可靠性的关键，是区分“没有雪”和“没有有效观测”，并把云、地形、林冠和时间填补带来的不确定性传递到面积和趋势结论。标准产品提供稳定基础，具体应用仍需匹配尺度和独立验证。