Ujian Akhir Praktikum
Penginderaan Jauh
Land Use Land Cover
(LULC)
Kabupaten Sumedang,
Jawa Barat Tahun 2018 - 2023
Ketik Sentinel 2 dan pilih
yang Level 2A lalu import
dan ganti nama variabel
menjadi s2
Bagi praktikan yang memiliki file shp
dapat diupload kemudian di import
dengan simbol roi
Apabila tidak ada maka bisa
menggunakan data dari FAO (2015)
var roi =
[Link]("FAO/GAUL/2015/level2")
.filter([Link]('ADM2_NAME', ‘Sumedang'));
[Link](roi, , ‘Kabupaten Sumedang',
false);
[Link](roi, 9);
false mengartikan saat
memunculkan layer pada gee
menjadi off
// Cloud masking function
function cloudMask(image){
var scl = [Link]('SCL');
var mask =
[Link](3).or([Link](8).and([Link](10))).eq(
0);
return
[Link](['B.*']).updateMask(mask);
}
// Composite an imagery
var image =
[Link](roi).filterDate('2020-01-01',
'2021-01-01')
.map(cloudMask)
.median().divide(10000).clip(roi);
[Link](image, { min: [0.05, 0.1, 0],
max: [0.3, 0.5, 0.2], bands: ['B11', 'B8',
'B4']}, 'Image');
// Band map for index calculation
var bandMap = {
BLUE: [Link]('B2'),
GREEN: [Link]('B3'),
RED: [Link]('B4'),
NIR: [Link]('B8'),
SWIR1: [Link]('B11'),
SWIR2: [Link]('B12')
};
Modified Normalized Difference Water Index (MNDWI)
(GREEN - SWIR1) / (GREEN + SWIR1)
Enhanced Vegetation Index (EVI)
(2,5 x (NIR - RED)) / (NIR + 6 x RED - 7,5 x BLUE + 1)
Normalized Difference Built-up Index (NDBI)
(SWIR1 - NIR) / (SWIR1 + NIR)
Normalized Difference Bareness Index (NDBaI)
(SWIR1 - SWIR2) / (SWIR1 + SWIR2)
Normalized Burn Ratio
(NIR - SWIR2) / (NIR + SWIR2)
// *******Indeks*******
// MNDWI (Modified Normalized Difference Water
Index)
var mndwi = [Link]('(GREEN - SWIR1) /
(GREEN + SWIR1)', bandMap).rename('MNDWI');
[Link](mndwi, { min: -1, max: 1, palette:
['red', 'white', 'blue'] }, 'MNDWI', false);
// EVI (Enhanced Vegetation Index)
var evi = [Link]('(2.5 * (NIR - RED)) /
(NIR + 6 * RED - 7.5 * BLUE + 1)',
bandMap).rename('EVI');
[Link](evi, { min: -1, max: 1, palette:
['blue', 'white', 'green'] }, 'EVI', false);
// NDBI (Normamlized Difference Built-up Index)
var ndbi = [Link]('(SWIR1 - NIR) /
(SWIR1 + NIR)', bandMap).rename('NDBI');
[Link](ndbi, { min: -1, max: 1, palette:
['blue', 'white', 'red'] }, 'NDBI', false);
// NDBaI (Normamlized Difference Bareness Index)
var ndbai = [Link]('(SWIR1 - SWIR2) /
(SWIR1 + SWIR2)', bandMap).rename('NDBaI');
[Link](ndbai, { min: -1, max: 1, palette:
['blue', 'white', 'red'] }, 'NDBaI', false);
// NBR (Normamlized Burn Ratio)
var nbr = [Link]('(NIR - SWIR2) / (NIR +
SWIR2)', bandMap).rename('NBR');
[Link](nbr, { min: -1, max: 1, palette:
['red', 'white', 'green'] }, 'NBR', false);
// *******Klasifikasi LC******* // Dry vegetation
// Water var dryVegetation = [Link](vegetation);
var water = [Link](0).rename('wt'); [Link]([Link](), { palette: 'limegreen' }, 'Dry
[Link]([Link](), { palette: 'lightskyblue' }, vegetation', false);
'Water', false);
var luasAir = [Link]([Link]()).divide(10000); // Wet vegetation
// Konversi m2 ke hektar var wetVegetation = [Link](vegetation);
var areaAir = [Link]({ [Link]([Link](), { palette: '3CB371' }, 'Wet
reducer: [Link](), vegetation', false);
geometry: roi,
scale: 200, //Jika ada tulisan 'user memory limit' angka scale // Low vegetation
ditambahkan var lowVegetation = [Link]([Link](0.6));
maxPixels: 1e9 [Link]([Link](), { palette: 'lightgreen' },
}); 'Low vegetation', false);
print('Luas Air (hektar):', [Link]('wt'));
// High vegetation
// Land var highVegetation = [Link]([Link](0.6));
var land = [Link](0); [Link]([Link](), { palette: 'darkgreen' },
[Link]([Link](), { palette: 'tan' }, 'Land', false); 'High vegetation', false);
// Wetland (Lahan Basah) // Marsh (Sawah)
var wetland = [Link]([Link](-0.2)); var marsh = [Link](wetland).rename('mr');
[Link]([Link](), { palette: 'olive' }, 'Wetland', [Link]([Link](), { palette: 'yellowgreen' }, 'Marsh',
false); false);
var luasMR = [Link]([Link]()).divide(10000); //
// Dryland (Lahan Kering) Konversi m2 ke hektar
var dryland = [Link]([Link](0)); var areaMR = [Link]({
[Link]([Link](), { palette: 'chocolate' }, reducer: [Link](),
'Dryland', false); geometry: roi,
scale: 200,
// Vegetation maxPixels: 1e9
var vegetation = [Link]([Link](0.4)); });
[Link]([Link](), { palette: 'green' }, print('Luas Sawah (hektar):', [Link]('mr'));
'Vegetation', false);
// Swamp (Rawa) // Forest
var swamp = [Link](wetland).rename('rw'); var forest = [Link](dryland).rename('fr');
[Link]([Link](), { palette: 'olive' }, 'Swamp', [Link]([Link](), { palette: 'darkgreen' }, 'Forest',
false); false);
var luasRW = [Link]([Link]()).divide(10000); // var luasFR = [Link]([Link]()).divide(10000); //
Konversi m2 ke hektar Konversi m2 ke hektar
var areaRW = [Link]({ var areaFR = [Link]({
reducer: [Link](), reducer: [Link](),
geometry: roi, geometry: roi,
scale: 200, scale: 200,
maxPixels: 1e9 maxPixels: 1e9
}); });
print('Luas Rawa (hektar):', [Link]('rw')); print('Luas Hutan (hektar):', [Link]('fr'));
// Shrub (Belukar) // Built-up
var shrub = [Link](dryland).rename('sh'); var builtup = [Link]([Link](-
[Link]([Link](), { palette: 'lightgreen' }, 'Shrub', 0.15)).and([Link](0.15)).rename('bu');
false); [Link]([Link](), { palette: 'lightcoral' }, 'Built-
var luasSH = [Link]([Link]()).divide(10000); // up', false);
Konversi m2 ke hektar var luasbu = [Link]([Link]()).divide(10000); //
var areaSH = [Link]({ Konversi m2 ke hektar
reducer: [Link](), var areabu = [Link]({
geometry: roi, reducer: [Link](),
scale: 200, geometry: roi,
maxPixels: 1e9 scale: 200,
}); maxPixels: 1e9
print('Luas Belukar (hektar):', [Link]('sh')); });
print('Luas Lahan Terbangun (hektar):', [Link]('bu'));
// Bare
var bare = [Link]([Link](0)).and([Link](0));
[Link]([Link](), { palette: 'tan' }, 'Bare', false);
// Dry bareland
var dryBareland = [Link](bare).rename('dl');
[Link]([Link](), { palette: 'tan' }, 'Dry
bareland', false);
var luasDL =
[Link]([Link]()).divide(10000); //
Konversi m2 ke hektar
var areaDL = [Link]({
reducer: [Link](),
geometry: roi,
scale: 200,
maxPixels: 1e9
});
print('Luas Lahan Kering (hektar):', [Link]('dl'));
// Wet bareland
var wetBareland = [Link](bare).rename('wl');
[Link]([Link](), { palette: 'teal' }, 'Wet
bareland', false);
var luasWL =
[Link]([Link]()).divide(10000); //
Konversi m2 ke hektar
var areaWL = [Link]({
reducer: [Link](),
geometry: roi,
scale: 200,
maxPixels: 1e9
});
print('Luas Lahan Basah (hektar):', [Link]('wl'));
// *******LULC Visual*******
// Visual
var vis = {
'lc_class_values': [11, 12, 21, 22, 30, 41, 42, 50],
'lc_class_names': ['Natural/Plantation forest', 'Shrub/Crop', 'Swamp', 'Marsh/ricefield', 'Built-
up', 'Dry bareland', 'Wet bareland', 'Water'],
'lc_class_palette': ['228B22', '90EE90', '808000', '9ACD32', 'F08080', 'D2B48C', '008080', '87CEFA']
};
// Land cover
var landcover = [Link](0).where([Link](1), 11)
.where([Link](1), 12)
.where([Link](1), 21)
.where([Link](1), 22)
.where([Link](1), 30)
.where([Link](1), 41)
.where([Link](1), 42)
.where([Link](1), 50)
.rename('lc')
.set(vis)
.clip(roi);
var luasLC = [Link]([Link]()).divide(10000); // Konversi m2 ke hektar
var areaLC = [Link]({
reducer: [Link](),
geometry: roi,
scale: 200,
maxPixels: 1e9
});
print('LuasLC (hektar):', [Link]('lc'));
// Show the land cover
[Link](landcover, {}, 'Land cover');
// Show legend
var panel = [Link]([[Link]('Legend')], [Link]('vertical'), { position: 'bottom-left' });
[Link](panel);
vis.lc_class_values.map(function(value, index){
[Link](
[Link](
[
[Link]('', { backgroundColor: vis.lc_class_palette[index], height: '20px', width: '30px' }),
[Link](value, { height: '20px' }),
[Link](vis.lc_class_names[index], { height: '20px' }),
],
[Link]('horizontal')
)
);
});
// *******Export Data Sample*******
//Export Hasil Data
[Link]({
image: landcover,
description: "Sumedang_2019_S2", //Nama Task
region: roi,
scale: 30, //Skala Pixel
maxPixels: 600000000,
fileFormat: 'GeoTIFF'
});
2023 2022 2021 2020 2019
470,96
533,87
624,59
562,61
503,60
37769,34
16434,98
26715,53
27422,01
60840,67
8512,97
5993,51
6656,23
6901,54
9616,27
18477,76
24727,51
14462,76
8027,00
5990,86
87354,81
104360,22
103889,43
109651,97
76549,35
0,00
6,22
0,00
2,28
0,00
41,56
417,21
208,43
242,05
54,35
2927,97
3081,84
2998,40
2745,91
2000,27
Luas (Ha)