中国30m逐年NDVI最大值数据集是基于Google Earth Engine（GEE）遥感云计算平台，利用美国陆地卫星Landsat 5/8遥感影像计算的1986年以来的逐年度NDVI最大值数据集。计算过程中对每年度全年所有的Landsat5/8遥感影像计算NDVI，然后获取每个像元所在位置的一年中的NDVI最大值，从而生成逐年度NDVI最大值数据。其中Landsat5 NDVI计算利用 B4和B3波段，Landsat8 NDVI计算利用 B5和B4波段。该数据集受Landsat影像数量影响明显，有的年份局部区域(主要是新疆、西藏)数据缺失，不同影像图幅间可能会有整体性条纹状差异。
数据空间分辨率为30m，时间分辨率为1986-2021年逐年度。数据格式为GeoTIFF和arcgis grid格式。数据数值范围为-2000-10000（由-0.2-1放大10000倍）。
GEE 遥感云计算平台计算及数据下载参考代码程序如下：
//2013年及以后NDVI计算调用Landsat8影像库，之前的计算则调用Landsat5影像库"LANDSAT/LT05/C01/T1_TOA"
var l8 =
ee.ImageCollection("LANDSAT/LC08/C01/T1_RT_TOA");
//计算2021年的NDVI
var year = 2021
//采取分省计算策略，"anhui"为上传到云平台的安徽省边界shape文件，
//也可以增加多个省界同时计算
var provinces = ["anhui"]
print(provinces);
var index = 0;
print("provinces.length is " +
provinces.length);
while (index < provinces.length)
{
print("while loop is " +
provinces[index]);
var qu = provinces[index]
var AOI =
ee.FeatureCollection("users/xuxl/"+qu);
var startyear = year+"-01-01"
var endyear = year+"-12-31"
var landsat = l8
.filterBounds(AOI)
.filterDate(startyear, endyear);
print (qu,'qu');
print (landsat,'landsat');
//NDVI计算函数
landsat5 利用 B4 and B3 波段landsat8 利用 B5 and B4波段
var get_NDVI = function(image) {
var NDVI=image.normalizedDifference(['B5','B4']).rename(['NDVI']);
image=image.addBands(NDVI)
return image.select("NDVI")
};
//针对影像集中的每一张影像裁剪ROI区域、计算NDVI、导出相应波段
var NDVI_Collection =
ee.ImageCollection(landsat)
.map(get_NDVI);
//print
(NDVI_Collection,'NDVI_Collection');
var greenest = NDVI_Collection.qualityMosaic('NDVI').clip(AOI).multiply(10000).toInt();
print(greenest)
//查看大体效果
Export.image.toDrive({
image: greenest,
description: qu,
folder : 'ndvi'+year,
fileNamePrefix: qu,
scale: 30,
maxPixels: 1e13,
// crs:'EPSG:32649',
region: AOI.geometry(),
fileFormat: 'GeoTIFF',
// formatOptions: {
// cloudOptimized: true
// }
});
Map.centerObject(AOI, 5);
Map.addLayer(greenest,{palette: ['black',
'green', 'red']},'greenest')
index ++;
}