MOD13Q1.061质量控制（去云等）之GEE

一、MOD13Q1.061数据说明（来自GEE）

1. DATA

MOD13Q1.061 Terra Vegetation Indices 16-Day Global 250m

2.DESCRIPTION

The MOD13Q1 V6.1 product provides a Vegetation Index (VI) value at a per pixel basis. There are two primary vegetation layers.
The first is the Normalized Difference Vegetation Index (NDVI) which is referred to as the continuity index to the existing National Oceanic and Atmospheric Administration-Advanced Very High Resolution Radiometer (NOAA-AVHRR) derived NDVI. The second vegetation layer is the Enhanced Vegetation Index (EVI) that minimizes canopy background variations and maintains sensitivity over dense vegetation conditions.
The MODIS NDVI and EVI products are computed from atmospherically corrected bi-directional surface reflectances that have been masked for water, clouds, heavy aerosols, and cloud shadows.

3.BANDS

该数据集共有12个波段：

4.部分波段详细信息（主要涉及QA质量评估波段）

所有波段详细信息可直接前往GEE平台搜索查看。

二、质量控制说明

QA波段中存储的值显示为十进制，该十进制数值需转换成二进制数值使用。
一般采用bitwiseAnd函数进行质量控制，是对二进制进行处理，按照以下方式计算：
bitwiseAnd : 1&1=1，1&0=0，0&1=0，0&0=0
二进制位数说明：
位数 —— 98 7654 3210（从右往左，0~15，16位），如：
1 —— 0000 0000 0000 0001
0 —— 0000 0000 0000 0000
左移符号说明：
例如：
1<<10 —— 0000 0100 0000 0000
需要提前将十进制转为二进制。

实在不懂可以参见以下文章：
遥感影像中的QA波段（质量评估波段）存储方式介绍
GEE质量评估（QA）波段——去云、云量统计

云不仅仅影响被云遮蔽的像元，也会影响邻近像元，一般去云的时候，还会考虑去除冰雪以及气溶胶等，按需求而定。

三、注意事项

MOD13Q1有两个质量控制波段：DetailedQA和SummaryQA，适用于不同波段的质量控制。

DetailedQA波段：涉及的是每个反射率波段的单独QA以及影像的信息，不用来对VI（NDVI，EVI）进行质量控制。用于对反射率波段的质量控制。

SummaryQA波段：适用于对VI进行质量控制。

四、质量控制示例——对NDVI波段进行质量控制

对于MOD13Q1数据集中的NDVI数据，使用SummaryQA波段进行质量控制。质量指标选择Bits:0-1 = 0，即高质量。
Bits 0-1是指第1位和第2位（从右往左）的值为0，则图像质量就是0：Good data，use with confidence.（看数据说明中的表格）
GEE代码示例：

//quality control 
function maskclouds(image) {
  var qa = image.select('SummaryQA');
  
  var Mask1 = 1 << 0;  
  var mask = qa.bitwiseAnd(Mask1).eq(0);
  
  var Mask2 = 1 << 1;
  var Mask = qa.bitwiseAnd(Mask2).eq(0);
  mask = mask.and(Mask);
  
  return image.updateMask(mask);
}

注：在选择数据时，要先使用质量控制函数进行掩膜，再选择自己所需的NDVI波段，才不会报错。否则会出以下错误：