Pengantar data Forest Monitoring for Action (FORMA)

FORMA adalah sistem peringatan deforestasi dua kali sebulan berbasis MODIS 500 x 500 meter untuk hutan tropis lembap. Dataset FORMA 500 di Earth Engine adalah gambar dengan peringatan yang dimulai pada Januari 2006 dan diperbarui setiap bulan. Setiap pemberitahuan memiliki waktu yang terkait dengannya dalam satu rentang 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 pemberitahuan 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:

Code Editor (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 ditutupi).

Menghitung Peringatan FORMA di Area Minat

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

Code Editor (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 Area yang Diminati

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

Code Editor (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 lain. Perhatikan bahwa setiap region dalam kumpulan kawasan lindung Republik Kongo kini memiliki properti tambahan, sum, yang dinamai sesuai dengan reducer. Nilai properti ini adalah output reducer, atau jumlah pemberitahuan tahun 2012 di kawasan lindung.

Membandingkan Set Data FORMA dan Hansen et al.

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

Code Editor (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)'
);

Kita telah sampai pada akhir ringkasan set data perubahan hutan di Earth Engine. Kami tidak sabar ingin melihat hasil kreasi Anda dengan fitur ini.