Forest Monitoring for Action (FORMA) ডেটার ভূমিকা

FORMA হল MODIS ভিত্তিক ৫০০ x ৫০০ মিটার আর্দ্র গ্রীষ্মমণ্ডলীয় বনের জন্য মাসে দুবার বন উজাড়ের সতর্কতা সিস্টেম। Earth Engine-এ FORMA 500 ডেটাসেট হল জানুয়ারি ২০০৬ থেকে শুরু হওয়া অ্যালার্ট সহ একটি ছবি এবং এটি প্রতি মাসে আপডেট করা হয়। প্রতিটি সতর্কতার সাথে alert_date নামের একটি ব্যান্ডে ইপক সেকেন্ড ইউনিটে সময় যুক্ত থাকে। FORMA ডেটাসেট ব্যবহার করে আপনি যেসব গুরুত্বপূর্ণ কাজ করতে পারেন তার মধ্যে দুটি হল তারিখ অনুযায়ী FORMA ফিল্টার করা এবং আগ্রহের এলাকার মধ্যে অ্যালার্ট গণনা করা।

তারিখ অনুযায়ী FORMA ফিল্টার করা হচ্ছে

শুধুমাত্র ২০১২ সালে ঘটা সতর্কতাগুলি দেখাতে, ১ জানুয়ারি, ১৯৭০ মধ্যরাত থেকে সেকেন্ডে প্রকাশ করা, ২০১২ সালের প্রথম দিন এবং ২০১৩ সালের প্রথম দিনের মধ্যে সময় আছে এমন পিক্সেলগুলি খুঁজুন:

কোড এডিটর (জাভাস্ক্রিপ্ট)

// 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'
);

এই উদাহরণে, forma2012 হল একটি বাইনারি ছবি যাতে শুধুমাত্র সেইসব পিক্সেল আছে যেগুলি ২০১২ সালে ঘটেছে (অর্থাৎ, অন্য সব পিক্সেল মাস্ক করা আছে)।

আগ্রহের অঞ্চলে FORMA সতর্কতা গণনা করা

আগের বিভাগে আমরা যেভাবে Hansen et al. ডেটা ব্যবহার করেছি, সেইভাবে আমরা আগ্রহের কোনও এলাকায় FORMA অ্যালার্টের (পিক্সেল) সংখ্যা গণনা করে শুরু করতে পারি। যেমন, ২০১২ সালে কঙ্গো প্রজাতন্ত্রের সুরক্ষিত এলাকায় অ্যালার্টের সংখ্যা গণনা করতে, আগের উদাহরণের উপর ভিত্তি করে নিম্নলিখিত কোডটি লিখুন:

কোড এডিটর (জাভাস্ক্রিপ্ট)

// 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'));

আগ্রহের বিভিন্ন অঞ্চলে FORMA অ্যালার্ট গণনা করা হচ্ছে

এখনও পর্যন্ত, আমরা একবারে একটি অঞ্চলের পরিসংখ্যান গণনা করছি। একসাথে একাধিক অঞ্চলে পরিসংখ্যান গণনা করার জন্য, আপনি reduceRegions() ব্যবহার করতে পারেন। আবার আগের উদাহরণের ভিত্তিতে:

কোড এডিটর (জাভাস্ক্রিপ্ট)

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

কনসোলে প্রিন্ট করা অবজেক্টটি পরীক্ষা করুন এবং লক্ষ্য করুন যে reduceRegions()-এর আউটপুট হল আরেকটি FeatureCollection। মনে রাখবেন যে কঙ্গো প্রজাতন্ত্রের সুরক্ষিত এলাকার সংগ্রহের প্রতিটি অঞ্চলে এখন একটি অতিরিক্ত প্রপার্টি আছে, sum, রিডিউসারের নামানুসারে নামকরণ করা হয়েছে। এই প্রপার্টির ভ্যালু হল রিডিউসারের আউটপুট অথবা সুরক্ষিত এলাকায় ২০১২ সালের সতর্কতার সংখ্যা।

FORMA এবং Hansen et al. ডেটাসেটের তুলনা করা

FORMA এবং Hansen et al. ডেটাসেট তুলনা করতে, আপনি লজিক্যাল অপারেটর ব্যবহার করতে পারেন। (লজিক্যাল অপারেশন সম্পর্কে আরও জানুন)। বিশেষ করে, আমরা এমন একটি ছবি তৈরি করতে চাই যেখানে FORMA এবং Hansen et al. ডেটা, উভয় দ্বারা বন উজাড় হিসেবে চিহ্নিত পিক্সেলগুলি ১ এবং বাকিগুলি শূন্য। এই কোড ২০১২ সালের জন্য এই ধরনের ছবি তৈরি করে এবং অন্যান্য পূর্বাভাসিত বন উজাড়ের লেয়ারের সাথে এটি দেখায়:

কোড এডিটর (জাভাস্ক্রিপ্ট)

// 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)'
);

Earth Engine-এ বনের পরিবর্তন সংক্রান্ত ডেটাসেটের ওভারভিউ এখানে শেষ হচ্ছে। আপনি এগুলি দিয়ে কী করতে পারেন তা দেখার জন্য আমরা অপেক্ষা করছি!