像元二分法计算植被覆盖度

通过MODIS的NDVI数据集MOD13Q1和像元二分法计算植被覆盖度

var multi_NDVI = ee.ImageCollection('MODIS/006/MOD13Q1').filterDate('2015-06-01', '2016-09-01').select('NDVI').max().divide(10000).clip(geometry);var ndviVis = {min: 0.0,max: 1,palette: ['FFFFFF', 'CE7E45', 'DF923D', 'F1B555', 'FCD163', '99B718', '74A901','66A000', '529400', '3E8601', '207401', '056201', '004C00', '023B01','012E01', '011D01', '011301'],
};//
Map.centerObject(geometry);//统计影像中对应分位数的值
var num = multi_NDVI.reduceRegion({reducer:ee.Reducer.percentile([5,95]),geometry:geometry,scale:250,maxPixels:1e13
});
var min = ee.Number(num.get("NDVI_p5"));
var max = ee.Number(num.get("NDVI_p95"));
print(min)
print(max)//根据统计值通过像元二分法计算植被覆盖度
var NDVI = ee.ImageCollection('MODIS/006/MOD13Q1').filterDate('2015-06-01', '2016-07-01').select('NDVI').max().divide(10000).clip(geometry);var greaterPart = NDVI.gt(max);
var lessPart    = NDVI.lt(min);
var middlePart  = ee.Image(1).subtract(greaterPart).subtract(lessPart);//计算FVC
var tempf1= NDVI.subtract(min).divide(max.subtract(min));var FVC = ee.Image(1).multiply(greaterPart).add(ee.Image(0).multiply(lessPart)).add(tempf1.multiply(middlePart));Map.addLayer(FVC, ndviVis, 'FVC');