GEE 利用Landsat-8影像计算年内最大ndvi

1. 导入内蒙古自治区矢量边界

直接上传边界到assets , import

2. 编译代码

var l8 = ee.ImageCollection("LANDSAT/LC08/C01/T1_RT_TOA");


//计算2022年的NDVI

var year = 2022



var AOI = ee.FeatureCollection(table2)
print(AOI)

var startyear = year+"-01-01"

var endyear = year+"-12-31"

var landsat = l8

    .filterBounds(AOI)

    .filterDate(startyear, endyear);


print (landsat,'landsat');

 

//NDVI计算函数 landsat5 利用 B4 and B3 波段landsat8 利用 B5 and B4波段

var get_NDVI = function(image) {

        var NDVI=image.normalizedDifference(['B5','B4']).rename(['NDVI']);

        image=image.addBands(NDVI)

        return image.select("NDVI")

      };

//针对影像集中的每一张影像裁剪ROI区域、计算NDVI、导出相应波段

var NDVI_Collection = ee.ImageCollection(landsat)

          .map(get_NDVI);

//print (NDVI_Collection,'NDVI_Collection');

 
// 求年内最大值
var greenest = NDVI_Collection.qualityMosaic('NDVI').clip(AOI).multiply(10000).toInt(); 

print(greenest)

 

//查看大体效果

Export.image.toAsset({

    image: greenest,

    description: "ndvi2022max" + year,

    // folder : 'ndvi'+year,

    // fileNamePrefix:  "ndvi2022max",

    scale: 30,

    maxPixels: 1e13,
    
    crs: "EPSG:4326", // 地理坐标系
    // crs:'EPSG:32649',

    region: AOI.geometry(),

    // fileFormat: 'GeoTIFF',

    // formatOptions: {

    // cloudOptimized: true

    // }

});

Map.centerObject(AOI, 5);

Map.addLayer(greenest,{palette: ['black', 'green', 'red']},'greenest')

// index ++;

V2更新----------------------------------------------------

https://code.earthengine.google.com/1666852ef42295870460f52b96a30a4c

// print(l8)

//计算2022年的NDVI

var year = 2022



// var AOI = ee.FeatureCollection(table2)
var AOI = ee.FeatureCollection(table2)
print(AOI)

var startyear = year+"-01-01"

var endyear = year+"-12-31"

var landsat = l8

    .filterBounds(AOI)

    .filterDate(startyear, endyear);
    
    // .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 10));


print (landsat,'landsat');
 


 

//NDVI计算函数 landsat5 利用 B4 and B3 波段landsat8 利用 B5 and B4波段

var get_NDVI = function(image) {

        var NDVI=image.normalizedDifference(['B5','B4']).rename(['NDVI']);

        image=image.addBands(NDVI)

        return image.select("NDVI")

      };

//针对影像集中的每一张影像裁剪ROI区域、计算NDVI、导出相应波段

var NDVI_Collection = ee.ImageCollection(landsat)

          .map(get_NDVI);

//print (NDVI_Collection,'NDVI_Collection');

 
// 求年内最大值
var greenest = NDVI_Collection.qualityMosaic('NDVI').clip(AOI).multiply(10000).toInt(); 

print(greenest)

 

//查看大体效果

/*Export.image.toAsset({

    image: greenest,

    description: "ndvi2022max" + year,

    // folder : 'ndvi'+year,

    // fileNamePrefix:  "ndvi2022max",

    scale: 30,

    maxPixels: 1e13,
    
    crs: "EPSG:4326", // 地理坐标系
    // crs:'EPSG:32649',

    region: AOI.geometry(),

    // fileFormat: 'GeoTIFF',

    // formatOptions: {

    // cloudOptimized: true

    // }

});*/

Export.image.toDrive({

    image: greenest,
    // image: greenset.select("NDVI"),

    description: "l9_ndvi2022max",

    folder : 'l9_ndvi'+year,

    // fileNamePrefix: "ndvi2022max",

    scale: 30,

    maxPixels: 1e13,

    // crs:'EPSG:32649',

    region: AOI.geometry(),

    fileFormat: 'GeoTIFF',

    // formatOptions: {

    //   cloudOptimized: true

    // }

});


Map.centerObject(AOI, 5);

Map.addLayer(greenest,{palette: ['black', 'green', 'red']},'greenest')

Map.addLayer(landsat.select("B4"),{palette: ['black', 'green', 'red']},'landsat')

下载到本地用toDrive()