时间戳

首页 > 解决方案 > 在 Google Earht Engine 上按时间间隔重新分类 NDVI 栅格

javascript - 在 Google Earht Engine 上按时间间隔重新分类 NDVI 栅格

问题描述

地点：

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

收藏：

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

标签： javascriptgoogle-earth-engine

解决方案

有多种方法可以做到这一点，我更喜欢的方式是使用决策树分类器。从您的问题来看，您想要的范围似乎小于 0.2、02-0.4、0.4-0.6、0.6-0.8（可能错过了问题）并且大于 0.8。我们需要为此构建一个决策树，使用它来创建分类器，然后将其应用于图像。

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);

您可以在此处查看工作示例

或者 您也可以使用逻辑运算符来测试所需范围内的值，然后乘以类别编号以获得每个类别的像素。例如对于第 2 类

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

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

推荐阅读