VNP43IA3：VIIRS BRDF/反照率日值全球 500 米数据集（V002）
VNP43IA3：VIIRS BRDF/反照率日值全球 500 米数据集（V002）
引言
地表反照率（Albedo）是描述地表对太阳辐射反射能力的关键物理参数，它直接影响地表能量收支平衡，是气候模型、陆面过程模拟和遥感反演中不可或缺的输入变量。从极地冰盖到热带雨林，从荒漠到农田，反照率的时空变化深刻影响着区域乃至全球的气候格局。
VNP43IA3是 NASA/NOAA 联合发布的 VIIRS（可见光红外成像辐射仪套件）BRDF/反照率日值全球 500 米数据集（第 2 版）。它继承并延续了 NASA MODIS BRDF/反照率数据产品的技术路线，利用 Suomi NPP 卫星的 VIIRS 传感器观测数据，以每日为时间分辨率、500 米为空间分辨率，提供全球陆表反照率估计。
该数据集的核心价值在于：通过 BRDF（双向反射率分布函数）模型对地表各向异性反射特性进行精确描述，进而反演得到黑空反照率（BSA）和白空反照率（WSA）两个关键物理量，分别代表直射光和漫射光条件下的地表反照率。
数据集概述
| 属性 | 说明 |
|---|---|
| 数据集ID | NASA/VIIRS/002/VNP43IA3 |
| 时间范围 | 2012-01-17 至 2026-07-25（持续更新） |
| 空间分辨率 | 500 米（所有波段） |
| 时间分辨率 | 1 天 |
| 空间范围 | 全球 |
| 数据生产者 | NASA 陆地过程分布式主动存档中心（LP DAAC）/ USGS 地球资源观测与科学中心 |
| 数据版本 | Version 002 |
数据生产方法
VNP43IA3 产品的生产融合了物理模型与多天观测的协同反演策略：
1. 多天合成策略
产品并非由单日观测直接反演，而是每天使用 16 天的 VIIRS 数据进行生产，并在时间维度上加权至第 9 天——这一时间重心信息体现在文件名中。这种多天合成策略有效降低了单日观测中云层、气溶胶等因素带来的噪声，提高了反照率估计的稳定性。
2. BRDF 模型反演
算法采用RossThick/Li-Sparse-Reciprocal（RTLSR）半经验核驱动 BRDF 模型。该模型通过三个核权重（来源于 VNP43IA1 产品）来重建地表各向异性反射效应，将不同观测几何条件下的方向反射率校正为统一的观测几何（VNP43IA4），进而计算出地方太阳正午时分的综合黑空反照率（BSA）和白空反照率（WSA）。
3. 波段覆盖
产品提供三个 VIIRS 影像波段（I1、I2 和 I3）的反照率估计，分别对应可见光、近红外和短波红外谱段，满足不同应用场景的需求。
扩展应用：研究人员可以将 BRDF 模型参数与简单的多项式结合，以获得任意太阳光照角度下的黑空反照率。同样，BSA 和 WSA 波段也可与多项式结合，手动估计瞬时实际反照率（晴空反照率）。
波段信息
所有波段的空间分辨率均为500 米。
| 波段名称 | 说明 |
|---|---|
| BRDF_Albedo_Band_Mandatory_Quality_I1 | 波段 I1 的 BRDF/反照率强制质量标记 |
| BRDF_Albedo_Band_Mandatory_Quality_I2 | 波段 I2 的 BRDF/反照率强制质量标记 |
| BRDF_Albedo_Band_Mandatory_Quality_I3 | 波段 I3 的 BRDF/反照率强制质量标记 |
| Albedo_BSA_I1 | 波段 I1 的黑空反照率（Black-Sky Albedo） |
| Albedo_BSA_I2 | 波段 I2 的黑空反照率 |
| Albedo_BSA_I3 | 波段 I3 的黑空反照率 |
| Albedo_WSA_I1 | 波段 I1 的白空反照率（White-Sky Albedo） |
| Albedo_WSA_I2 | 波段 I2 的白空反照率 |
| Albedo_WSA_I3 | 波段 I3 的白空反照率 |
物理含义：
- 黑空反照率（BSA）：仅考虑直射太阳辐射条件下的反照率，随太阳天顶角变化
- 白空反照率（WSA）：仅考虑漫射辐射条件下的反照率，不随太阳天顶角变化
- 两者结合可估算任意天空条件下的实际反照率
Earth Engine 代码示例
JavaScript（Code Editor）
// 加载 VNP43IA3 数据集vardataset=ee.ImageCollection('NASA/VIIRS/002/VNP43IA3');// 查看数据集信息print('Dataset:',dataset);print('Number of images:',dataset.size());// 按时间过滤（例如：2024 年 7 月）varfiltered=dataset.filterDate('2024-07-01','2024-07-31');// 选择白空反照率（WSA）波段 I1，计算月均值varwsa=filtered.select('Albedo_WSA_I1').mean();// 可视化参数（反照率值范围 0-1）varvisParams={min:0,max:0.5,palette:['white','yellow','orange','red','darkred']};// 添加到地图Map.centerObject(wsa,4);Map.addLayer(wsa,visParams,'WSA I1 (Monthly Mean)');多波段分析与合成（JavaScript）
// 加载最新一期数据varlatest=dataset.sort('system:time_start',false).first();// 提取三个波段的白空反照率varwsa_i1=latest.select('Albedo_WSA_I1');varwsa_i2=latest.select('Albedo_WSA_I2');varwsa_i3=latest.select('Albedo_WSA_I3');// 计算可见光-近红外组合反照率（简单平均）varcompositeAlbedo=wsa_i1.add(wsa_i2).add(wsa_i3).divide(3);// 可视化Map.addLayer(compositeAlbedo,{min:0,max:0.5,palette:['white','yellow','orange','red']},'Composite Albedo');提取时间序列（JavaScript）
// 定义感兴趣点varpoint=ee.Geometry.Point([-100,40]);// 提取白空反照率时间序列varts=dataset.select('Albedo_WSA_I1').map(function(image){varval=image.reduceRegion({reducer:ee.Reducer.mean(),geometry:point,scale:500,maxPixels:1e9});returnee.Feature(null,{'date':image.date().format('YYYY-MM-dd'),'WSA_I1':ee.Number(val.get('Albedo_WSA_I1'))});});print('WSA I1 Time Series:',ts);全部代码
vardataset=ee.ImageCollection('NASA/VIIRS/002/VNP43IA3').filter(ee.Filter.date('2021-06-01','2021-06-03'));varblackSkyAlbedo=dataset.select('Albedo_BSA_I1').first();varpalette=['000080','0000d9','4000ff','8000ff','0080ff','00ffff','00ff80','80ff00','daff00','ffff00','fff500','ffda00','ffb000','ffa400','ff4f00','ff2500','ff0a00','ff00ff',];varvisParams={min:0,max:400,palette:palette,};// cadetbluevarbackground=ee.Image.rgb(95,158,160).visualize({min:0,max:255});varimage=blackSkyAlbedo.visualize(visParams);varlon=-8;varlat=60;Map.addLayer(background,{},'background');Map.addLayer(image,{},'Black-Sky Albedo I1');Map.setCenter(lon,lat,3);可视化结果
Python（geemap / ee）
importeeimportgeemap# 初始化 Earth Engineee.Initialize()# 加载数据集dataset=ee.ImageCollection('NASA/VIIRS/002/VNP43IA3')# 按时间过滤filtered=dataset.filterDate('2024-07-01','2024-07-31')wsa=filtered.select('Albedo_WSA_I1').mean()# 创建交互式地图Map=geemap.Map()Map.centerObject(wsa,4)# 可视化参数vis_params={'min':0,'max':0.5,'palette':['white','yellow','orange','red','darkred']}Map.addLayer(wsa,vis_params,'WSA I1 (Monthly Mean)')Map.addLayerControl()Map应用场景
VNP43IA3 数据集在以下领域具有重要的应用价值：
- 气候与能量平衡研究：反照率是地表能量收支的核心参数，直接影响地表温度、感热通量和潜热通量的模拟
- 陆面过程与水文建模：为 Noah-MP、CLM 等陆面模型提供高分辨率反照率输入
- 城市热岛效应研究：监测城市下垫面反照率变化，评估热岛效应的缓解策略（如白色屋顶）
- 农业与干旱监测：反照率变化可作为植被健康状态和土壤湿度的指示因子
- 极地与冰雪监测：追踪冰盖和积雪反照率变化，评估冰川消融速率
使用条款与引用
数据使用政策：
LP DAAC NASA 数据可免费访问；不过，如果作者发布这些数据或基于这些数据创作作品，则需要在出版物的正文中引用相应数据集，并在引用列表中添加对这些数据集的引用。
建议引用格式：
Schaaf, C., Z. Wang, A. Erb, I. Paynter. VIIRS/NPP BRDF/Albedo Albedo Daily L3 Global 500m SIN Grid V002. 2024, NASA EOSDIS Land Processes Distributed Active Archive Center. https://doi.org/10.5067/VIIRS/VNP43IA3.002
注意事项
数据版本：当前为Version 002，V001 版本已于 2025 年 7 月 31 日停用。用户应使用 V002 版本。
时间合成：产品虽为“日值”，但并非单日观测，而是基于 16 天观测窗口加权合成至第 9 天。在分析快速变化的地表过程时需注意这一时间平滑效应。
质量标记：使用前建议检查
BRDF_Albedo_Band_Mandatory_Quality质量波段，筛选高质量像元。波段选择：I1、I2、I3 三个波段分别对应不同的光谱区间，用户应根据具体应用选择合适的波段或组合使用。
MODIS 延续性：该产品旨在促进 NASA MODIS BRDF/反照率数据产品套件的连续性，可与 MODIS 产品（MCD43A3）进行时间序列衔接分析。
总结
VNP43IA3（V002）是 VIIRS 时代全球陆表反照率监测的旗舰级数据集。它以每日 500 米分辨率，提供黑空反照率（BSA）和白空反照率（WSA）两种关键物理量，覆盖三个重要光谱波段。通过继承 MODIS 成熟的 BRDF 模型反演算法，该数据集实现了从 MODIS 到 VIIRS 的无缝延续，为气候研究、陆面过程模拟和生态环境监测提供了高质量、长时间序列的反照率数据支撑。借助 Google Earth Engine 平台，用户可以便捷地访问、处理和分析这一宝贵数据资源。