GEE(Google Earth Engine)计算全球的NDVI和SAVI!

2024-02-01 19:55:33 浏览数 (2)

 2002年11月NDVI均值

此次案例主要是通过Landsat7来进行2002年11月-12月的NDVI和SAVI指数的计算,而且是通过表达式,这中表达是在于,以后你也想通过波段计算的话,可以改改公式即可使用,话不多说!上代码:

代码语言:javascript复制
// Filter the L7 collection to a single month.
var collection = ee.ImageCollection('LANDSAT/LE07/C01/T1_TOA')
    .filterDate('2002-11-01', '2002-12-01');

// A function to compute NDVI.
var NDVI = function(image) {
  return image.expression('float(b("B4") - b("B3")) / (b("B4")   b("B3"))');
};

// A function to compute Soil Adjusted Vegetation Index.
var SAVI = function(image) {
  return image.expression(
      '(1   L) * float(nir - red)/ (nir   red   L)',
      {
        'nir': image.select('B4'),
        'red': image.select('B3'),
        'L': 0.2
      });
};

// Shared visualization parameters.这个条带颜色可以保存起来,以后根据你的影像值修改min,max就行了,要不然以后自己写的时候会比较麻烦
var vis = {
  min: 0,
  max: 1,
  palette: [
      'FFFFFF', 'CE7E45', 'DF923D', 'F1B555', 'FCD163', '99B718',
      '74A901', '66A000', '529400', '3E8601', '207401', '056201',
      '004C00', '023B01', '012E01', '011D01', '011301'
  ]
};

Map.setCenter(-93.7848, 30.3252, 11);

// Map the functions over the collection, reduce to mean and display.
//这里只是展示图像,如若想要求出具体的平均值或者打印出来,需要用到reducer
Map.addLayer(collection.map(NDVI).mean(), vis, 'Mean NDVI');
Map.addLayer(collection.map(SAVI).mean(), vis, 'Mean SAVI');

 2002年11月SAVI均值

0 人点赞