বার্ষিক বনভূমি ধ্বংসের চার্ট

বার্ষিক বনভূমি ধ্বংসের হিসেব করা

আগের বিভাগে আপনি শিখেছেন কীভাবে reduceRegion পদ্ধতি ব্যবহার করে প্রদত্ত অঞ্চলে মোট হারানো বনাঞ্চল গণনা করতে হয়। মোট ক্ষতির পরিমাণ হিসেব করার পরিবর্তে, প্রতি বছরের ক্ষতির পরিমাণ হিসেব করলে তা সহায়ক হবে। Earth Engine-এ এটি করার উপায় হল গ্রুপ করা রিডিউসার ব্যবহার করা।

reduceRegion()-এর আউটপুট গ্রুপ করতে, আপনি একটি গ্রুপিং ব্যান্ড নির্দিষ্ট করতে পারেন যা পূর্ণসংখ্যা পিক্সেল মান দ্বারা গ্রুপ সংজ্ঞায়িত করে। নিম্নলিখিত উদাহরণে, আমরা আগের কোডটি সামান্য পরিবর্তন করেছি এবং আসল ছবিতে lossYear ব্যান্ড যোগ করেছি। lossYear ব্যান্ডের প্রতিটি পিক্সেলে ০ থেকে ১৪ পর্যন্ত ভ্যালু থাকে - যা নির্দেশ করে যে কোন বছরে ক্ষতি হয়েছে। এছাড়াও, আমরা রিডিউসারকে গ্রুপ করা রিডিউসারে পরিবর্তন করি, গ্রুপিং ব্যান্ডের (১) ব্যান্ড ইনডেক্স নির্দিষ্ট করে যাতে পিক্সেল এরিয়া যোগ করা হবে এবং lossYear ব্যান্ডের ভ্যালু অনুযায়ী গ্রুপ করা হবে।

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

// Load country boundaries from LSIB.
var countries = ee.FeatureCollection('USDOS/LSIB_SIMPLE/2017');
// Get a feature collection with just the Congo feature.
var congo = countries.filter(ee.Filter.eq('country_co', 'CF'));

// Get the loss image.
// This dataset is updated yearly, so we get the latest version.
var gfc2017 = ee.Image('UMD/hansen/global_forest_change_2017_v1_5');
var lossImage = gfc2017.select(['loss']);
var lossAreaImage = lossImage.multiply(ee.Image.pixelArea());

var lossYear = gfc2017.select(['lossyear']);
var lossByYear = lossAreaImage.addBands(lossYear).reduceRegion({
  reducer: ee.Reducer.sum().group({
    groupField: 1
    }),
  geometry: congo,
  scale: 30,
  maxPixels: 1e9
});
print(lossByYear);

উপরের কোডটি রান করার পরে, আপনি groups নামক নেস্টেড তালিকায় বার্ষিক বনভূমি হ্রাসের ক্ষেত্রফল দেখতে পাবেন। আমরা আউটপুটটি একটু ফর্ম্যাট করতে পারি যাতে ফলাফলটি একটি অভিধান হয়, যেখানে বছরটি কী এবং ক্ষতির ক্ষেত্রটি মান হিসাবে থাকে। লক্ষ্য করো যে আমরা format() পদ্ধতি ব্যবহার করে বছরের মান ০-১৪ থেকে ২০০০-২০১৪-এ পরিবর্তন করছি।

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

var statsFormatted = ee.List(lossByYear.get('groups'))
  .map(function(el) {
    var d = ee.Dictionary(el);
    return [ee.Number(d.get('group')).format("20%02d"), d.get('sum')];
  });
var statsDictionary = ee.Dictionary(statsFormatted.flatten());
print(statsDictionary);

চার্ট তৈরি করা

এখন আমাদের কাছে বার্ষিক ক্ষতির সংখ্যা আছে, আমরা একটি চার্ট তৈরি করতে প্রস্তুত। আমরা ui.Chart.array.values() পদ্ধতি ব্যবহার করব। এই পদ্ধতিতে ইনপুট ভ্যালুর একটি অ্যারে (বা তালিকা) এবং X-অ্যাক্সিসের জন্য লেবেলের একটি অ্যারে (বা তালিকা) নেওয়া হয়।

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

var chart = ui.Chart.array.values({
  array: statsDictionary.values(),
  axis: 0,
  xLabels: statsDictionary.keys()
}).setChartType('ColumnChart')
  .setOptions({
    title: 'Yearly Forest Loss',
    hAxis: {title: 'Year', format: '####'},
    vAxis: {title: 'Area (square meters)'},
    legend: { position: "none" },
    lineWidth: 1,
    pointSize: 3
  });
print(chart);

ফলাফলটি নিচের চার্টের মতো দেখতে হতে হবে।

ছবি ১.

-এর মাধ্যমে বছর অনুযায়ী বনভূমি ধ্বংসের চার্ট

পরবর্তী বিভাগে, আপনি বন উজাড় মনিটর করার আরেকটি ডেটাসেট, FORMA সম্পর্কে জানবেন এবং Hansen et al. ডেটার সাথে এটির তুলনা করবেন।