资讯详情

GEE-Scholars 数据工具--Sentinel-2遥感植被指数库

沐风-GEE系列专栏

  • Google Earth Engine 系列学习笔记(1)
  • Google Earth Engine 系列学习笔记(2)
  • GEE-Scholars 可视化时间序列农业景观
  • GEE-Scholars 机器学习卫星图像监督分类
  • GEE-Scholars 水体动态监测-河南水灾案例
  • GEE-Scholars 江淮平原水稻面积识别
  • GEE-Scholars 森林火灾的位置和面积
  • GEE-Scholars MODIS地表温度LST时间变化趋势
  • GEE-Scholars 数据工具–Sentinel-2遥感植被指数库

本章目录

  • 沐风-GEE系列专栏
  • 引言
  • 哨兵2号卫星各种遥感植被指数公式


引言

遥感指数是指根据多波段图像计算的图像。 图像强调特定现象,减少其他因素,降低图像效果。 例如,植被指数将索引图像中的健康植被显示为亮色,不健康植被值较低,而裸露区域为暗色。 由于地形变化(山谷)产生的阴影会影响图像强度,指数将通过强调对象的颜色而不是强调对象的强度或亮度来创建。 山谷阴影下健康松树的植被指数值与全光下松树的植被指数值相似。 这些指数通常由波段的加减组合构成,从而产生各种波段比。 它们与电磁光谱特定部分的特定波段有关。 因此,这些指数可能仅对某些传感器或某几类传感器有效,而且在计算中使用正确的波段至关重要。


哨兵2号卫星各种遥感植被指数公式

//RVI function s2_rvi(image){ 
              var rvi = image.expression('NIR / RED', { 
             'NIR' : image.select('B8'),     'RED' : image.select('B4')   }).float();     return image.addBands(rvi.rename("S2_RVI")); }  //EVI function s2_evi(image){ 
              var evi = image.expression('2.5 * ((NIR - RED) / (NIR 6 * RED - 7.5 * BLUE 1))', { 
             'NIR' : image.select('B8'),     'RED' : image.select('B4'),
    'BLUE': image.select('B2')
  }).float();
  return image.addBands(evi.rename('S2_EVI'));
}  
​
//DVI
function s2_dvi(image){ 
        
     var dvi = image.expression('NIR - RED', { 
        
    'NIR' : image.select('B8'),
    'RED' : image.select('B4')
   
  }).float();
  return image.addBands(dvi.rename('S2_DVI'));
}
​
// NDVI——B8,B4
function s2_ndvi_b84(image){ 
        
   return image.addBands(image.normalizedDifference(['B8', 'B4']).rename("S2_NDVI_B84"));
}
​
//LSWI——B8,B11
function s2_lswi_b811(image){ 
        
   return image.addBands(image.normalizedDifference(['B8', 'B11']).rename("S2_LSWI_B811"));
}
​
// LSWI——B8,B12
function s2_lswi_b812(image){ 
        
   return image.addBands(image.normalizedDifference(['B8', 'B12']).rename("S2_LSWI_B812"));
}
​
// LSWI——B8A,B11
function s2_lswi_b8a11(image){ 
        
   return image.addBands(image.normalizedDifference(['B8A', 'B11']).rename("S2_LSWI_B8A11"));
}
​
//LSWI——B8A,B12
function s2_lswi_b8a12(image){ 
        
   return image.addBands(image.normalizedDifference(['B8A', 'B12']).rename("S2_LSWI_B8A12"));
}
​
​
//MNDWI
function s2_mndwi(image){ 
        
   return image.addBands(image.normalizedDifference(['B3', 'B11']).rename("S2_MNDWI"));
}
​
//NDWI
function s2_ndwi(image){ 
        
   return image.addBands(image.normalizedDifference(['B3', 'B8']).rename("S2_NDWI"));
}
​
//GCVI
function s2_gcvi(image){ 
        
     var gcvi = image.expression('NIR / Green - 1', { 
        
    'NIR' : image.select('B8'),
    'Green' : image.select('B3')
  }).float();
  return image.addBands(gcvi.rename('S2_GCVI'));
} 
​
//S2REP——B7,B4,B5,B6
function s2_s2rep(image){ 
        
     var s2rep = image.expression('705 + 35 * ((((NIRn1 + RED) / 2) - RE1) / (RE2 - RE1))', { 
        
    'NIRn1' : image.select('B7'),
    'RE2' : image.select('B6'),
    'RE1' : image.select('B5'),
    'RED' : image.select('B4')
  }).float();
  return image.addBands(s2rep.rename('S2_S2REP'));
}
​
//SAVI
function s2_savi(image) { 
        
  // Add Soil Adjust Vegetation Index (SAVI)
    // using L = 0.5;
    var savi = image.expression('(NIR - RED) * (1 + 0.5)/(NIR + RED + 0.5)', { 
        
    'NIR': image.select('B8'),
    'RED': image.select('B4')
    }).float();
    return image.addBands(savi.rename('S2_SAVI'));
}
 
 //IBI
function s2_ibi(image) { 
        
  // Add Index-Based Built-Up Index (IBI)
  var ibiA = image.expression('2 * SWIR1 / (SWIR1 + NIR)', { 
        
    'SWIR1': image.select('B11'),
    'NIR'  : image.select('B8')
  }).rename(['IBI_A']);
 
  var ibiB = image.expression('(NIR / (NIR + RED)) + (GREEN / (GREEN + SWIR1))', { 
        
    'NIR'  : image.select('B8'),
    'RED'  : image.select('B4'),
    'GREEN': image.select('B3'),
    'SWIR1': image.select('B11')
  }).rename(['IBI_B']);
 
  var ibiAB = ibiA.addBands(ibiB);
  var ibi = ibiAB.normalizedDifference(['IBI_A', 'IBI_B']);
  return image.addBands(ibi.rename(['S2_IBI']));
}
//NDBI
function s2_ndbi(image){ 
        
   return image.addBands(image.normalizedDifference(['B11', 'B8']).rename("S2_NDBI"));
}
​
// NIRv
function s2_nirv(image) { 
        
    
    var nirv = image.expression('((NIR - RED)/(NIR + RED)) * NIR', { 
        
    'NIR': image.select('B8'),
    'RED': image.select('B4')
    }).float();
    return image.addBands(nirv.rename('S2_NIRv'));
}
//CI red edge (CIr) 
​
function s2_cir(image) { 
        
    
    var cir = image.expression('RE3 / RE1 - 1', { 
        
    'RE3': image.select('B7'),
    'RE1': image.select('B5')
    }).float();
    return image.addBands(cir.rename('S2_CIr'));
}
​
//CI green (CIg)
​
function s2_cig(image) { 
        
    
    var cig = image.expression('RE3 / Green - 1', { 
        
    'RE3': image.select('B7'),
    'Green': image.select('B3')
    }).float();
    return image.addBands(cig.rename('S2_CIg'));
}
​
//MTCI
function s2_mtci(image) { 
        
    
    var mtci = image.expression('(RE2 - RE1)/ (RE1 - Red)', { 
        
    'RE2': image.select('B6'),
    'RE1': image.select('B5'),
     'Red': image.select('B4'),
    }).float();
    return image.addBands(mtci.rename('S2_MTCI'));
}
//PSRI
function s2_psri(image) { 
        
      var psri= image.expression(
              '(RED-BLUE)/REDEDRE2',{ 
        
              REDEDRE2:image.select('B6'),
              RED:image.select('B4'),
              BLUE:image.select('B2'),
            }).float().rename('PSRI');  
    return image.addBands(psri);
}

标签: 111系列ndvi传感器组

锐单商城拥有海量元器件数据手册IC替代型号,打造 电子元器件IC百科大全!

锐单商城 - 一站式电子元器件采购平台