基于GEE使用Landsat 8和Landsat 5影像计算RSEI

一、使用2018年landsat 8影像数据计算遥感生态指数RSEI

//could masking，GEE上有现有的去云函数
function maskL8sr(image) {
  // Bit 0 - Fill
  // Bit 1 - Dilated Cloud
  // Bit 2 - Cirrus
  // Bit 3 - Cloud
  // Bit 4 - Cloud Shadow
  var qaMask = image.select('QA_PIXEL').bitwiseAnd(parseInt('11111', 2)).eq(0);
  var saturationMask = image.select('QA_RADSAT').eq(0);

  // Apply the scaling factors to the appropriate bands.
  var opticalBands = image.select('SR_B.').multiply(0.0000275).add(-0.2);
  var thermalBands = image.select('ST_B.*').multiply(0.00341802).add(149.0);

  // Replace the original bands with the scaled ones and apply the masks.
  return image.addBands(opticalBands, null, true)
      .addBands(thermalBands, null, true)
      .updateMask(qaMask)
      .updateMask(saturationMask);
}

var dataset1 = ee.ImageCollection('LANDSAT/LC08/C02/T1_L2')
    .filterDate('2018-x-x', '2018-x-x')
    .filterBounds(roi)
    .filter(ee.Filter.lte('CLOUD_COVER',5))
    .sort('CLOUD_COVER')
    .map(maskL8sr)
    //.mean()
print(dataset1)

var image1=dataset1.mosaic().clip(roi)
print(image1,'maskL8sr')
//Map.addLayer(image1,{bands:["SR_B4","SR_B3","SR_B2"],gamma:1.3,max:108,min:15},'Image1');
//Map.addLayer(roi,{color:'yellow',fillColor: "00000000", width: 1},'roi');

var jrc = jrcwater18.select('b1')
var image = image1.updateMask(jrc.eq(1))
//Map.addLayer(image,{bands:["SR_B4","SR_B3","SR_B2"],gamma:1.3,max:108,min:15},'Image');
print(image,'water and cloud')

//计算SI

function SI_cal(img) {
 var blue = img.select("SR_B2");
 var red = img.select("SR_B4");
 var nir = img.select("SR_B5");
 var swir1 = img.select("SR_B6");
 var swir2 = img.select("SR_B7");
 var SI_temp =((swir1.add(red)).subtract(blue.add(nir)))
              .divide((swir1.add(red)).add(blue.add(nir)))
 //print(SI_temp);
 return SI_temp;
}

//计算IBI
function IBI_cal(img) {
 var green = img.select("SR_B3");
 var red = img.select("SR_B4");
 var nir = img.select("SR_B5");
 var swir1 = img.select("SR_B6");
 var IBI_temp =(((swir1.multiply(2.0)).divide(swir1.add(nir)))
               .subtract((nir.divide(nir.add(red))).add(green.divide(green.add(swir1)))))
               .divide((((swir1.multiply(2.0)).divide(swir1.add(nir)))
               .add((nir.divide(nir.add(red))).add(green.divide(green.add(swir1))))))
 //print(IBI_temp);
 return IBI_temp;
}

//计算湿度指数Wet

function Wet_cal(img) {
 var blue = img.select("SR_B2");
 var green = img.select("SR_B3");
 var red = img.select("SR_B4");
 var nir = img.select("SR_B5");
 var swir1 = img.select("SR_B6");
 var swir2 = img.select("SR_B7");
 var wet_temp =blue.multiply(0.1511)
              .add(green.multiply(0.1973))
              .add(red.multiply(0.3283))
              .add(nir.multiply(0.3407))
              .add(swir1.multiply(-0.7117))
              .add(swir2.multiply(-0.4559))
 //print(wet_temp);
 return wet_temp;
}

//外部导入的NDVI重采样