Lampiran G: Contoh kode sumber terbuka dan perhitungan
Demo Google Earth Engine dan akses ke kode ISBM di Savimbo GitHub
Apakah ini membantu?
Apakah ini membantu?
// Muat poligon plot18
var plot18_data = require("users/drea/map:plot18_data");
var plot18 = plot18_data.polygon;
// Hitung luas plot dalam hektare
var plotArea = plot18.area();
var plotAreaHectares = plotArea.divide(10000);
// Muat titik jaguar
var points_jaguar_data = require("users/drea/map:points_jaguar_data");
var puntos = points_jaguar_data.points;
var radios = [];
var sumMultipliedArea = ee.Number(0);
// Tentukan koleksi fitur untuk menyimpan poligon irisan
var intersectionPolygons = ee.FeatureCollection([]);
// Tentukan fungsi assignedArea
var assignedAreaFunction = function(offset) {
var day = startDate.advance(offset, 'day');
var dayString = day.format('YYYY-MM-dd');
var feature = ee.Feature(null, { date: dayString, intersectionArea: intersectionArea });
return feature.set('date_area', ee.String(dayString).cat(' - ').cat(intersectionArea));
};
// Hitung radius dan buffer untuk setiap titik
for (var i = 0; i < puntos.length; i++) {
var point = puntos[i].geometry;
var date = puntos[i].date;
// Hitung tanggal mulai dengan mengurangi 30 hari
var startDate = ee.Date(date).advance(-30, 'day');
// Hitung tanggal akhir dengan menambah 30 hari
var endDate = ee.Date(date).advance(30, 'day');
// Buat fitur dengan geometri titik dan tanggal sebagai properti
var feature = ee.Feature(point, { date: date });
var pointBuffer = feature.buffer(800);
radios.push(pointBuffer);
// Hitung irisan dengan plot18
var intersection = pointBuffer.intersection(plot18);
// Hitung luas dalam hektare
var area = intersection.area().divide(10000);
// Ambil bulan dan tahun dari tanggal
var month = ee.Date(date).get('month');
var year = ee.Date(date).get('year');
// Buat label untuk bulan dan tahun
var monthYearLabel = ee.String(month).cat('-').cat(year).cat(' Hektare N°');
// Hitung irisan dengan plot18 untuk titik saat ini
var intersectionPlot18 = pointBuffer.intersection(plot18);
// Hitung luas dalam hektare untuk irisan dengan plot18
var intersectionArea = intersectionPlot18.area().divide(10000);
var multipliedArea = intersectionArea.multiply(60);
// Tambahkan multipliedArea ke jumlah
sumMultipliedArea = sumMultipliedArea.add(multipliedArea);
// Tentukan rentang tanggal jaguar
var jaguarRange = endDate.difference(startDate, 'day');
// Tetapkan nilai intersectionArea untuk setiap hari dalam rentang jaguar
var assignedArea = ee.FeatureCollection(ee.List.sequence(0, jaguarRange.subtract(1)).map(assignedAreaFunction));
// Tampilkan hasil untuk setiap titik
print('Jaguar:', i + 1);
print('Tanggal pengambilan citra:', date);
print('Tanggal mulai:', startDate.format('YYYY-MM-dd'));
print('Tanggal akhir:', endDate.format('YYYY-MM-dd'));
print('Total luas irisan per hari dalam hektare:', intersectionArea);
print('Luas irisan dikali 60:', multipliedArea);
print('Luas yang ditetapkan per hari:', assignedArea);
print('----------------------');
// Tambahkan buffer titik saat ini ke peta dengan warna biru
Map.addLayer(pointBuffer, { color: 'blue' }, 'Radius ' + (i + 1));
// Tambahkan geometri irisan ke intersectionPolygons
intersectionPolygons = intersectionPolygons.merge(intersection);
}
// Lakukan penggabungan poligon untuk menghindari duplikat
var unionPolygons = intersectionPolygons.union();
// Hitung total luas irisan tanpa duplikat
var totalIntersectionArea = unionPolygons.geometry().area().divide(10000);
// Tampilkan hasil
print('Total luas irisan dalam hektare:', totalIntersectionArea);
print('Luas plot18 dalam hektare:', plotAreaHectares);
print('Jumlah luas irisan dikali 60 hari dalam hektare:', sumMultipliedArea);
// Tambahkan layer plot18 ke peta
Map.addLayer(plot18, { color: 'gold' }, "plot18");
Map.centerObject(plot18);