Pengantar data Forest Monitoring for Action (FORMA)

FORMA adalah sistem peringatan deforestasi dua bulanan berukuran 500 x 500 meter berbasis MODIS untuk hutan tropis lembap. Set data FORMA 500 di Earth Engine adalah gambar dengan pemberitahuan yang dimulai pada Januari 2006 dan diperbarui setiap bulan. Setiap pemberitahuan memiliki waktu yang terkait dengannya dalam satu band bernama alert_date dalam satuan detik epoch. Memfilter FORMA menurut tanggal dan menghitung pemberitahuan dalam area minat adalah dua hal terpenting yang dapat Anda lakukan dengan set data FORMA.

Memfilter FORMA menurut Tanggal

Untuk menampilkan hanya peringatan yang terjadi pada tahun 2012, temukan piksel yang memiliki waktu antara hari pertama tahun 2012 dan hari pertama tahun 2013, yang dinyatakan dalam detik sejak tengah malam, 1 Januari 1970:

Editor Kode (JavaScript)

// Convert dates from milliseconds to seconds.
var start = ee.Date('2012-01-01').millis().divide(1000);
var end = ee.Date('2013-01-01').millis().divide(1000);

// Load the FORMA 500 dataset.
var forma = ee.Image('FORMA/FORMA_500m');

// Create a binary layer from the dates of interest.
var forma2012 = forma.gte(start).and(forma.lte(end));

Map.setCenter(15.87, -0.391, 7);
Map.addLayer(
  forma2012.mask(forma2012),
  {palette: ['FF0000']},
  'FORMA alerts in 2012'
);

Dalam contoh ini, forma2012 adalah gambar biner yang hanya berisi piksel yang memiliki waktu yang terjadi pada tahun 2012 (yaitu, semua piksel lainnya disamarkan).

Menghitung Peringatan FORMA di Region of Interest

Seperti yang kita lakukan di bagian sebelumnya dengan data Hansen et al., kita dapat memulai dengan menghitung jumlah pemberitahuan FORMA (piksel) di area minat. Misalnya, untuk menghitung jumlah pemberitahuan di kawasan lindung Republik Kongo pada tahun 2012, kembangkan contoh sebelumnya sebagai berikut:

Editor Kode (JavaScript)

// Load country features from Large Scale International Boundary (LSIB) dataset.
var countries = ee.FeatureCollection('USDOS/LSIB_SIMPLE/2017');

// Subset the Congo Republic feature from countries.
var congo = ee.Feature(
  countries
    .filter(ee.Filter.eq('country_na', 'Rep of the Congo'))
    .first()
);

// Subset protected areas to the bounds of the congo feature
// and other criteria. Clip to the intersection with congo.
var protectedAreas = ee.FeatureCollection('WCMC/WDPA/current/polygons')
  .filter(ee.Filter.and(
    ee.Filter.bounds(congo.geometry()),
    ee.Filter.neq('IUCN_CAT', 'VI'),
    ee.Filter.neq('STATUS', 'proposed'),
    ee.Filter.lt('STATUS_YR', 2010)
  ))
  .map(function(feat){
    return congo.intersection(feat);
  });

// Display protected areas on the map.
Map.addLayer(
  protectedAreas,
  {color: '000000'},
  'Congo Republic protected areas'
);

// Calculate the number of FORMA pixels in protected
// areas of the Congo Republic, 2012.
var stats = forma2012.reduceRegion({
  reducer: ee.Reducer.sum(),
  geometry: protectedAreas.geometry(),
  scale: 500
});
print('Number of FORMA pixels, 2012: ', stats.get('constant'));

Menghitung Notifikasi FORMA di Beberapa Region yang Diminati

Sejauh ini, kita telah menghitung statistik di satu region dalam satu waktu. Untuk menghitung statistik di beberapa region sekaligus, Anda dapat menggunakan reduceRegions(). Sekali lagi, berdasarkan contoh sebelumnya:

Editor Kode (JavaScript)

var regionsStats = forma2012.reduceRegions({
  collection: protectedAreas,
  reducer: ee.Reducer.sum(),
  scale: forma2012.projection().nominalScale()
});
print(regionsStats);

Periksa objek yang dicetak ke konsol dan amati bahwa output reduceRegions() adalah FeatureCollection lainnya. Perhatikan bahwa setiap region dalam koleksi kawasan lindung Republik Kongo kini memiliki properti tambahan, sum, yang dinamai sesuai dengan reducer. Nilai properti ini adalah output dari reducer, atau jumlah pemberitahuan tahun 2012 di kawasan lindung.

Membandingkan Set Data FORMA dan Hansen et al.

Untuk membandingkan set data FORMA dan Hansen et al., Anda dapat menggunakan operator logika. (Pelajari lebih lanjut operasi logika). Secara khusus, kita ingin membuat gambar yang pikselnya ditandai oleh data FORMA dan Hansen dkk. sebagai deforestasi adalah 1 dan sisanya adalah 0. Kode ini membuat gambar tersebut untuk tahun 2012 dan menampilkannya bersama dengan lapisan prediksi deforestasi lainnya:

Editor Kode (JavaScript)

// Convert dates from milliseconds to seconds.
var start = ee.Date('2012-01-01').millis().divide(1000);
var end = ee.Date('2013-01-01').millis().divide(1000);
var region = ee.Geometry.Rectangle([-59.81163, -9.43348, -59.27561, -9.22818]);

// Load the FORMA 500 dataset.
var forma = ee.Image('FORMA/FORMA_500m');

// Create a binary layer from the dates of interest.
var forma2012 = forma.gte(start).and(forma.lte(end));

// Load Hansen et al. data and get change in 2012.
var gfc = ee.Image('UMD/hansen/global_forest_change_2015');
var gfc12 = gfc.select(['lossyear']).eq(12);

// Create an image which is one where the datasets
// both show deforestation and zero elsewhere.
var gfc_forma = gfc12.eq(1).and(forma2012.eq(1));

// Display data on the map.
Map.setCenter(-59.58813, -9.36439, 11);
Map.addLayer(forma.updateMask(forma), {palette: '00FF00'}, 'Forma (green)');
Map.addLayer(gfc12.updateMask(gfc12), {palette: 'FF0000'}, 'Hansen (red)');
Map.addLayer(
  gfc_forma.updateMask(gfc_forma),
  {palette: 'FFFF00'},
  'Hansen & FORMA (yellow)'
);

Sekian ringkasan set data perubahan hutan di Earth Engine. Kami menantikan apa yang dapat Anda lakukan dengan fitur-fitur ini.