这是我用来估算NDVI的所有代码。

初始值

var roi =/* color: #98ff00 *//* displayProperties: [{"type": "rectangle"}] */
    ee.Geometry.Polygon(
        [[[-72.3734130976651, -13.430548306629259],
          [-72.3734130976651, -14.326448420293916],
          [-70.7254638789151, -14.326448420293916],
          [-70.7254638789151, -13.430548306629259]]], null, false);
var startDate = ee.Date('2018-08-23');
var endDate = ee.Date('2018-12-21');
var cloud = 20


选择卫星
0:(2015年6月23日至2019年7月17日)前哨2级1C

var selector = 0;
var Sentinel = "";


NDVI

function getNDVI(image){
    var isLandsat75 = (selector !== 0);
    var NDVI = image.expression('(nir - red)/(nir + red)',
    {'nir':image.select(isLandsat75 ? 'B4' : 'B8'),
    'red':image.select(isLandsat75 ? 'B8' : 'B4')});
    //print(NDVI);
    return NDVI;
}


选择卫星

switch(selector){ case 0: Sentinel =  'COPERNICUS/S2'}


采集

var collection = ee.ImageCollection(Sentinel)
  .filter(ee.Filter.lt("CLOUDY_PIXEL_PERCENTAGE", cloud))
  .filterDate(startDate, endDate)
  .filterBounds(roi);
print(collection)


计算NDVI

var ndviFiltered = collection.map(getNDVI).qualityMosaic('B8').clip(roi);


重新分类栅格

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

var classifier = ee.Classifier.decisionTree(DTstring);
var reclass = ndviFiltered.select('B8').classify(classifier);
print(reclass)

Map.addLayer(reclass,{min:1,max:5}, 'Reclass');


转换为向量

var vectors = ndviFiltered.reduceToVectors({
  geometry: roi,
  scale: 1000,
  geometryType: 'polygon',
  reducer: ee.Reducer.countEvery()
});
print(vectors)


在这里,我得到这个错误“ FeatureCollection(错误)
Image.reduceToVectors:图像的第一个波段('B8')必须是整数”
频段B8已经是必不可少的

最佳答案

我正面临类似的问题。
NDVI值为连续值,因此无法以这种格式将其转换为向量。
您必须定义离散区域并用空白图像掩盖。
参见此link

09-25 20:02