网站首页   |   专业术语  |  0 1级数据标准   |   产品标准   |   共享统计   |   技术论坛
欢迎访问MODIS技术标准网    
MODIS 植被指数数据产品参考规范(讨论稿)
内容索引:

1 主题内容与适用范围

2 植被指数类型

3 术语

4 植被指数数据产品的生产标准
5 MODIS植被指数质量保证与验证 6 参考资料


1 主题内容与适用范围

  • 主题内容
    本标准规定了国家对地观测系统MODIS共享平台植被指数数据产品术语、类型、制作、和验证过程,用以规范我国MODIS植被指数数据产品在产生、保藏、交换和应用中的一致性。

  • 适用范围
    本规范适用于国家科技基础条件平台对地观测系统MODIS共享平台植被指数数据产品及与之相关的数据产品在数据源、数据合成、数据质量检验和数据交换过程中的活动规范。

2 植被指数类型

MODIS植被指数分为归一化植被指数(NDVI)和增强型植被指数(Enhanced Vegetation Index, EVI)二种类型。

3 术语

植被指数

通过地表覆盖物在可见光波谱段的吸收和在近红外波谱段的反射特性,建立的用于描述植被数量和质量的参数。植被指数没有量纲。

地表反射数据

指经过大气校正的 MODIS 1-7 波段数据,即MOD09 产品。

植被指数合成

在多日植被指数中,按照一定标准和规则,选择其中一个植被指数的过程。

植被指数合成期

用于实施合成的时间段。以天、旬、月度计算。

BRDF 合成

双向反射分布函数。指把传感器视角的观测值,统一为星下点观测值;同时把不同太阳高度角统一为有代表性的一个角度。

4 植被指数数据产品的生产标准

  • 4.1 单日植被指数计算
    4.1.1 输入数据:输入去云并且经过大气校正的地面反射数据。MODIS 1-7 波段。其中:1-2波段空间分辨率 250m,3-7波段空间分辨率500m.
    4.1.2 植被指数定义及计算公式
    NDVI = (B2-B1)/(B2+B1)
    EVI =2.5*(B2-B1)/(B2+6*B1–7.5*B3 + 1)
    其中:
    NDVI:归一化植被指数
    EVI:增强型植被指数
    B1:MODIS第1波段
    B2:MODIS第2波段
    B3:MODIS第3波段
    4.1.3 输出数据:NDVI和 EVI,日数据,空间分辨率250米。
  • 4.2 植被指数合成规范
    4.2.1 合成期:旬(10天)合成期和月合成期。
    4.2.2 合成算法使用原则:
    MODIS 植被指数的合成是在像元水平上进行的。 只能采用以下4种合成方法中的一种,根据输入数据的质量,按照以下优先序列决定使用的合成方法。
    (1) BRDF 合成: 如果合成期内,晴天观测值超过30%,采用BRDF 原理合成。
    (2) 有限视角内的最大值合成(CV-MVC):如果合成期内,晴天观测值少于30%,且多于2天,选其中视角最小的2个计算植被指数,取二者中最大值。
    (3) 直接计算植被指数:如果只有1天无云,则直接使用这天数据计算的植被指数。
    (4) 最大值合成(MVC):如果所有观测全部有云,则取合成期所有植被指数的最大值。

5 MODIS植被指数质量保证与验证

  • 5.1 质量保证的定义和范围
    业务化运行过程中的标记性数据产品,表达与预期精度不相符的数据,同时存储对后期处理有用的相关数据和信息。

  • 5.2 质量标记方法
    (1) 像元层次 (Pixel):根据4.2.4的合成原则,将旬、月度合成期有关数据质量的信息存在一个大小为16bit的波段。合成后的数据是由被选中像元所对应的所有原始数据组成。其中,1km 分辨率植被指数数据包括 11 个波段(表9-1)。
表9-1 1km 分辨率植被指数数据的波段组成

 

数据名称

数据类型

尺度转换

有效值域

填充值

1

1 km 10 天 NDVI

INT 16

10000

-2000 ~ 10000

-3000

2

1 km 10 天 EVI

INT 16

10000

-2000 ~ 10000

-3000

3

1 km 10 天 NDVI 质量标记

UINT 16

不适用

0-65534

65535

4

1 km 10 天 EVI 质量标记

UINT 16

不适用

0-65534

65535

5

1 km 10 天 red 反射值

INT 16

10000

0-10000

-1000

6

1 km 10 天 NIR 反射值

INT 16

10000

0-10000

-1000

7

1 km 10 天 blue 反射值

INT 16

10000

0-10000

-1000

8

1 km 10 天 MIR 反射值

INT 16

10000

0-10000

-1000

9

1 km 10 天 平均视角

INT 16

100

-9000 ~ 9000

-10000

10

1 km 10 天平均太阳高度角

INT 16

100

-9000 ~ 9000

-10000

11

1 km 10 天平均方位角

INT 16

10

-3600 ~ 3600

-4000

其中,3和4是根据合成前的各种质量信息的总结,由16 bit 组成。各bit 表示的质量信息如表9-2所示。

表 9-2 各 bit 表示的质量信息

16 bit 序号

各 bit 表达的质量信息

0-1

NDVI 质量信息: 4
(00) 已经制作 NDVI ,质量可靠;
(01) 已经制作 NDVI ,但使用时需要检查质量标记;
(10 ) 没有制作 NDVI, 原因时云覆盖;
(11) 没有制作 NDVI ,其他原因。

2-5

植被指数可用性分级: 16
(0) 质量最高; ………. (7) 质量中等 ……
(13) 没有经过大气校正;
(14) 质量太差,无法使用;
(15) 由于其他原因,数据无法使用。

6-7

大气气溶胶影响: (00) 当地气候基准值水平 (01) 气溶胶水平低 (10) 平均 (11) 高

8

是否进行邻边校正: (1) 是 , (0) 否

9

在大气校正时是否使用了双向反射分布函数: (1) 是 (0) 否

10

是否存在混合云: (1) 是 (0) 否

11-12

海陆标记: (00) 海洋 , (01) 海岸 t, (10) 湿地 , (11) 陆地

13

冰雪可能性: (1) 是 (0) 否

14

云影可能性: (1) 是 (0) 否

15

NDVI 合成方法: (0) BRDF 合成,把 VI 换算为星下点; (1) CV-MVC.

(2)文件层次 (Tile):通过元数据文件描述,可以是定量的(如云覆盖率),或定性的(如,根据检验结果,描述数据整体质量)。

  • 5.3 植被指数质量验证方法
    使用较高空间分辨率的遥感数据,或使用实地观测数据,通过对比研究,验证MODIS 植被指数的可靠性。
    (1) 样地代表性原则:选择样地应该包括多种生态类型。既要有植被覆盖很少的样地,也要有植被覆盖很高的样地,以便检查植被指数在不同生态类型稳定性和灵敏性。
    (2) 地面最小面积:样地所在的均一区域,应该为所检验像元的2-3倍以上。如,如果检验空间分辨率250m 的像元,样地所代表的区域应该在500x500m 以上。
  • 5.4 验证指标
    (1) 定性分类能力:是否能够反映生态系统时间和空间变化。时间指标主要表现为对土地覆盖季节性的表达;空间指标指对土地覆盖空间分异性的表达。
    (2) 定量能力:与实地观测生物量的关联程度。
    (3) MODIS-NDVI 与AVHRR-NDVI 比较验证:检验MODIS-NDVI是否能够延续AVHRR-NDVI的时间序列。

6 参考文献

[1] Alfredo Huete,Chris Justice,Wim van Leeuwen. 1999.MODIS VEGETATION INDEX (MOD 13) ALGORITHM THEORETICAL BASIS DOCUMENT,Version 3. (April 30, 1999).
[2] A. Huete, K. Didan, T. Miura, E.P. Rodriguez, X. Gao, L.G. Ferreira. 2002. Overview of the radiometric and biophysical performance of the MODIS vegetation indices. Remote Sensing of Environment. 83 (2002) 195–213.

 

(执笔:王正兴、刘闯、赵冰茹、朱晓华)

 

 

 
联系信息: 电子信箱:wangzx@igsnrr.ac.cn  电话:010-64889073
 
    © 中国科学院地理科学与资源研究所全球变化信息研究中心 2004
    
 www.nfiieos.cn