在 Google Earth Engine 上按间隔重新分类 NDVI 栅格

Reclassify NDVI raster in intervals on Google Earht Engine

地点:

var roi = /* color: #d63000 */ee.Geometry.Point([-71.97203347683796, -13.529827050320447]);

Collection:

var collection = ee.ImageCollection('COPERNICUS/S2') 
  .filter(ee.Filter.lt("CLOUDY_PIXEL_PERCENTAGE", 10)) 
  .filterDate('2018-01-1' ,'2018-12-31')
  .filterBounds(roi)

计算 NDVI

function addNDVI(image) {
  var a = image.normalizedDifference(['B8', 'B4']);
  return image.addBands(a);
}
var ndvi2 = collection.map(addNDVI)
var ndvi2 = ndvi2.qualityMosaic('nd');

现在我想在内部 [-1-0.2]、[0.2-0.4]、[0.4-0.6]、[0.6-0.8]、[0.8-1] 中重新分类 NDVI 栅格,我试试这个代码

var ndvireclass = ndvi2.select("nd").divide(10).ceil();
Map.addLayer(ndvireclass, {bands:'nd', min: 0, max: 1, gamma: 1.5}, 'NDVI reclass');

但结果图像只有2个类-1和1

有多种方法可以做到这一点,我更喜欢的方法是使用决策树 classifier。从你的问题来看,你想要的范围似乎是 小于 0.2、02-0.4、0.4-0.6、0.6-0.8(可能遗漏)和大于 0.8。 我们需要为此构建一个决策树,用它来创建一个 classifier 然后将其应用于图像。

var DTstring = ['1) root 9999 9999 9999',
'2) nd<=0.2 9999 9999 1 *',
'3) nd>0.2 9999 9999 9999',
'6) nd<=0.4 9999 9999 2 *',
'7) nd>0.4 9999 9999 9999',
'14) nd<=0.6 9999 9999 3 *',
'15) nd>0.6 9999 9999 9999',
'30) nd<=0.8 9999 9999 4 *',
'31) nd>0.8 9999 9999 5 *'].join("\n");

var classifier = ee.Classifier.decisionTree(DTstring);
var reclassifiedImage = ndvi2.select('nd').classify(classifier);

你可以看到工作示例here

您还可以使用逻辑运算符测试所需范围内的值,然后乘以 class 数字以获得每个 class 的像素。 例如 class 2

var nd = ndvi2.select('nd');
var c2 = nd.gt(0.2).and(nd.lte(0.4)).multiply(2);

如果您对 c1、c3、c4、c5 执行类似操作,您应该拥有仅包含值为 0 和 class 数字的像素的栅格。如果你添加所有这些层你应该得到你想要的