4つの式の計算（式）を学習するGoogle Earth Engine（GEE）のコード例

タイトル式

この記事では、1つの画像の表現を計算する方法について説明します。

コード例は、MOVI表面反射率画像EVIを計算することです。
正規化された差分インデックスは、上記のGEEによって提供される関数を使用して簡単に計算されます。

//加载MODIS地表反射率数据其单位为*10000，所以这里需要*0.0001
var img = ee.Image('MODIS/006/MOD09GA/2012_03_09').multiply(0.0001);

// 定义自己的表达式
//是在单张影像上进行操作
// variable name to band name in the input image.
var evi = img.expression(
    '2.5 * (nir - red) / (nir + 6 * red - 7.5 * blue + 1)',
    {
        red: img.select('sur_refl_b01'),    // 620-670nm, RED
        nir: img.select('sur_refl_b02'),    // 841-876nm, NIR
        blue: img.select('sur_refl_b03')    // 459-479nm, BLUE
    });
// 地图中心
Map.setCenter(-94.84497, 39.01918, 8);
// 加载影像，真彩色合成
Map.addLayer(img.select(['sur_refl_b01', 'sur_refl_b04', 'sur_refl_b03']),
         {min: 0, max: 0.2}, 'MODIS bands 1/4/3');
//加载计算的EVI，会覆盖在上幅影像上
Map.addLayer(evi, {min: 0, max: 1}, 'EVI');

コード操作の結果

次に、コードを変更して、Landsat8燃焼指数BAIを計算します。

var img = ee.Image('LANDSAT/LC08/C01/T1_SR/LC08_123032_20160113').multiply(0.0001);

var MyExpressionBAI='1/(sqrt(0.1-RED)+sqrt(0.06-NIR))'
var BAI = img.expression(MyExpressionBAI,
    {
        RED: img.select('B4'),    // 620-670nm, RED
        NIR: img.select('B5'),    // 841-876nm, NIR
    });

// Center the map.
Map.setCenter(116, 39.9, 8);

Map.addLayer(BAI, {min: 1.79, max: 3.97}, 'BAI');

運用実績