简介
MODIS 表面反射率数据 为 MODIS 数据产品中的 MOD09 系列。包括:Level 2 (MOD09)、Level 2G (MOD09GHK, MOD09GQK, MOD09GST)、Level 2G-Lite daily (MOD09GA, MOD09GQ)、Level 3, 8-day composited products (MOD09A1, MOD09Q1)、daily Level 3 CMG (MOD09CMG, MOD09CMA) 等产品。
产品处理等级(Level 0 ~ 5)介绍见:MODIS 及其 数据、产品_ww ^的CSDN博客
注释:
(1) 2GL 代表 level 2G-lite。
res. 即 resolutions;1表示1km,5 表示 0.05°,H 表示 500m,Q 表示 250m。
(2) collection 是指重新处理过的 MODIS data archive。后面的 collection 会替代前面的。
(3) 用户手册中的 MOD 不是特指 Terra MODIS,而是泛指所有 MODIS 数据(Terra 和 Aqua)。
波段
MODIS 表面反射率产品 是一个 7 波段产品,由 MODIS Level 1B 的陆地波段计算而来。该产品提供了每个波段的地表光谱反射率(surface spectral reflectance)估测结果。该产品在对大气中的气体和气溶胶进行校正后,于地面尺度上测量,即假设没有大气散射和吸收的作用。
波段 | 波长 | 光谱 |
bands 1 | 620-670 nm | 红 |
bands 2 | 841-876 nm | 近红外(NIR) |
bands 3 | 459-479 nm | 青 |
bands 4 | 545-565 nm | 绿 |
bands 5 | 1230-1250 nm | 短波红外(SWIR) |
bands 6 | 1628-1652 nm | 短波红外(SWIR) |
bands 7 | 2105-2155 nm | 短波红外(SWIR) |
产品分类
Product Name | Terra Prod ID | Aqua Prod ID |
Surface Reflectance 8-Day L3 Global 500m | MOD09A1 | MYD09A1 |
Surface Reflectance 8-Day L3 Global 250m | MOD09Q1 | MYD09Q1 |
Surface Reflectance Daily L2G Global 1km and 500m | MOD09GA | MYD09GA |
Surface Reflectance Daily L2G Global 250m | MOD09GQ | MYD09GQ |
Surface Reflectance Daily L3 Global 0.05Deg CMG | MOD09CMG | MYD09CMG |
在 8 天合成产品中,每个像元的值都是在这 8 天中最好的 L2G 观测值,选择的指标包括:高观测覆盖率、低观测角度、没有云和云阴影、没有气溶胶负荷。
产品处理算法
Products of MOD_PR09.exe --- levels 2 and 3
该程序会对所有Level 1B 数据(经过辐射校正)运行,即整个轨道、每个可用的分辨率。对 Daytime 数据进行大气的气体和气溶胶等校正。具体来说,即对 波段 1和2 的250m、波段1 ~ 7 的 500m、波段1 ~ 16 的 1km 等数据进行校正,为每个波段生产出 去除了大气吸收和散射影响的、近似于地表上测量的 表面光谱反射率。同时,为每个波段可用的每个分辨率生产了质量控制信息(辅助数据是否可用、L1B 数据有无谬误 ......)。
Level 2 MOD09 的输出包括:所有经过大气校正的波段及其质量数据、气溶胶检索数据(aerosol retrieval data)、用于评估气溶胶检索算法的质量的数据(热红外波段 20、31、32的亮度温度数据、水蒸气数据、径辐射数据 ......)。这些数据将作为科学数据集(Scientific Data Sets, SDS )写入输出文件。 band quality SDS 包含关于每个像元的大气校正质量的信息。 State QA SDS 包含关于每个像元的状态(state)的信息,即每个像元的与波段或分辨率无关的特性。每个State QA 像元都包含了如 该像元是否被标记为土地/深海/浅海、该像元是否包含云/高气溶胶/低气溶胶/雪/火 之类的信息。State QA 数据反映了像元本身的质量,而非任何表面反射率数据的质量。因此,1km State QA SDS 是另一个重要的输出数据集。
MOD09 的另一个输出是 MOD09 Intermediate Surface Reflectance datasets(MOD09IDN、-IDT、-IDS)。在该数据集中,每个轨道的所有表面反射率数据和波段质量数据都以 5km(0.05°) 分辨率被地理定位到一个线性的经纬度投影上。这些文件中的数据会被平均化。这些文件是中间文件,将被用作 MOD_PR09C.exe 的输入,MOD_PR09C.exe 会把 MOD09IDN、-IDT、-IDS 合成为 daily MOD09CMG、MOD09CMA。
Products of MOD_PRMGR.exe --- level 2G
该程序每天都会对 MODIS sinusoidal grid 中的每个 tile 运行(对当天映射到该 tile 上的所有 MOD09 level 2 granules 运行)。如 figure 1 所示。
该程序生成的数据可以看成一个三维立方体,每个立方体的深度由其对应像元的观测数量决定。每个像元上的观测数量不仅由该地点的环绕圈数决定(赤道处环绕一圈,极地环绕圈数最多,为15圈),也会受到 观测覆盖范围 影响。
L2G 数据会分成两个部分写入输出文件:第一部分由 “first layer” 数据组成,“first layer” 即在立方体深度为 0 处的数据,该部分的形式是 二维 SDSs。第二部分要么是 “full format” 数据,即第一部分之外的所有数据,形式为 三维 SDSs;要么是 “compact format” 数据,删除第一部分之外的所有观测值,形式为 一维 SDSs。本项目生成的是 “compact format”。
L2G 的输出包括:MOD09GHK (500 m surface reflectance and band quality data)、 MOD09GQK (250 m surface reflectance and band quality data)、MOD09GST (1 km State QA data)。除了经过地理校正,此 1km State QA 与 MOD09 中的 State QA一样。这些数据一直存档到 collection 5,后被 L2G-lite 产品取代。
Products of MOD_PR09G.exe --- level 2G-lite
在 L2G 的立方体中,观测值的顺序没有意义,只是按照读取 L2文件的顺序排列的。因此,实际上"first layer" 并不会比 “compact format” 更有用。
但是,随着时间推移,大多数用户更关注 L2G 输出的"first layer"。因此,现在在 L2G-lite 格式中,本项目通过类似于 MOD_PR09A.exe(下一节)的方法对观测值进行排序,从而尽量提高 “first layer” 数据的质量。这种对观测值的排序,也消除了对 level 2G 指针文件中多个 SDS 的需求,从而减少了数据总量。
Products of MOD_PR09A.exe --- level 3
此程序 每 8 天 对 MODIS sinusoidal grid 中的每个 tile 运行。截止(包含) collection 4,该程序在每个tile 的 所有8天period 的 所有MOD09GHK、-GQK、-GST 产品上运行。collection 5,则在每个tile 的 所有8天period 的 所有 MOD09GA、MOD09GQ 上运行。
该程序会选择每个像元的最佳观测值进行合成。在输入文件中、投影到某个特定像元的所有数据,都被称为 观测值(observations)。选择规则为:观测覆盖范围(observational coverage)、观测角度(view angle)、该观测值是否被标记为 cloudy / in cloud shadow / clear / containing high aerosol or low aerosol。
对于每个像元,合成步骤如下:
1)通过分析观测覆盖范围,对相同轨道的观测值进行合成。仅保留观测覆盖范围最广的观测值,其他都丢弃。这个步骤会为每个轨道生成一个观测结果列表。
2)为每个轨道的观测值赋分。
分值 | 描述 | |
0 | 值为 fill value | |
1 | BAD | 该数据对应的 L1B 像元 是 faulty 或 poorly crrected 的 |
2 | HighView | 观测视角 ≥ 60° |
3 | LowSun | 太阳天顶角(solar zenith angle) ≥ 85° |
4 | Cloudy | 该数据被云遮盖 (cloudy / adjacent to cloud) |
5 | Shadow | 该数据包含了云阴影(cloud shadow) |
6 | UnCorrected | 该数据被标记为 uncorrected |
7 | ClimAerosol | 该数据包含了一些气溶胶(default level of aerosols) |
8 | HighAerosol | 该数据包含了很多气溶胶(highest level of aerosols) |
9 | Snow | 该数据被雪覆盖(snow) |
10 | GOOD | 没有上述任何问题的数据 |
该程序会选择分数最高、观测视角最低的观测值作为 MOD09A1 和 MOD09Q1 的输出。MOD09A1 也包含了由所有8天输入数据合成的 1 km State QA 的 500m 版本。
Products of MOD_PR09C.exe --- level 3 CMG
MODIS Surface Reflectance Climate Modeling Grid (CMG) 是 Level 3 数据,且覆盖全球;该数据的投影是一个 线性的纬度和经度(Plate Carre),分辨率为 0.05°,如 figure 2 所示。
产品关系总结
产品细节
MOD09
MODIS Terra/Aqua Surface Reflectance 5-minute L2 Swath
bands 1 - 2 (250m)、bands 1 – 7 (500 m)、bands 1 – 16 (1 km resolution)、multiresolution QA、 1 km observation statistics
Product Granule ID 例子:MOD09.A2011026.0035.005.2011027195038.hdf
MOD09GQ
MODIS Terra/Aqua Surface Reflectance Daily L2G Global 250 m
Product Granule ID 例子:MOD09GQ.A2000339.h10v05.005.2006339053418.hdf.
MOD09GA
MODIS Terra/Aqua Surface Reflectance Daily L2G Global 500 m and 1 km
Product Granule ID 例子:MOD09GA.A2000340.h11v05.005.2006339102700.hdf
GQ 分辨率更高(250m),但是只有波段1-2;GA分辨率低(500m),但波段更全(1-7)。
MOD09Q1
MODIS Terra/Aqua Surface Reflectance 8-Day L3 Global 250 m
Product Granule ID 例子:MOD09Q1.A2000337.h10v05.005.2006342044337.hdf
MOD09A1
MODIS Terra/Aqua Surface Reflectance 8-Day L3 Global 500 m
Product Granule ID 例子:MOD09A1.A2000337.h11v05.005.2006342055602.hdf
MOD09CMG
MODIS Terra/Aqua Surface Reflectance Daily L3 Global 0.05 Deg CMG
band 1-7(0.05-degree resolution)
Product Granule ID 例子:MOD09CMG.A2000341.005.2006347161131.hdf
产品质量
Band QA
Band QA 反映了每个像元的表面反射率数据的质量,主要反映大气校正质量的信息。
State QA
State QA 反映了像元本身的质量,比如:该像元是否被标记为土地/深海/浅海、该像元是否包含云/高气溶胶/低气溶胶/雪/火 ......
Internal CM
Number Mapping
Geolocation flags
Scan value information
Orbit and coverage
参考资料
MODIS Surface Reflectance User’s Guide Collection 6(Version 1.4,2015)
MODIS Algorithm Technical Background Document
MODIS 官网:MODIS Web (nasa.gov)
MODIS 产品验证:MODIS Land Team Home Page (nasa.gov)
DOI:E. Vermote. (2015). MOD09A1 MODIS Surface Reflectance 8-Day L3 Global 500m SIN Grid V006. NASA EOSDIS Land Processes DAAC.
http://doi.org/10.5067/MODIS/MOD09A1.006 (Terra)
http://doi.org/10.5067/MODIS/MYD09A1.006 (Aqua)