Гребля Каховської гідроелектростанції знищено армією Росії 6 червня 2023 року. Як виглядає Каховське водосховище рік потому? Використаємо для візуалізації та аналізу платформу Google Earth Engine, супутникові знімки Sentinel-2 та індекс NDWI. Нижче наведено приклад коду для отримання наведених вище карт.
var geometry =
ee.Geometry.Polygon(
[[[33.0019701001224, 48.02863226215282],
[33.0019701001224, 46.58353286204405],
[35.5288255688724, 46.58353286204405],
[35.5288255688724, 48.02863226215282]]], null, false);
var addNDWI = function (image) {
var nd = image.normalizedDifference(['B3', 'B8']).rename('ndwi');
return image.addBands(nd);
}
var s2 = ee.ImageCollection("COPERNICUS/S2_HARMONIZED");
var csPlus = ee.ImageCollection('GOOGLE/CLOUD_SCORE_PLUS/V1/S2_HARMONIZED');
var QA_BAND = 'cs';
var CLEAR_THRESHOLD = 0.5;
var composite1 = s2
.filterBounds(geometry)
.filterDate('2022-06-01', '2022-08-01')
.linkCollection(csPlus, [QA_BAND])
.map(function(img) {
return img.updateMask(img.select(QA_BAND).gte(CLEAR_THRESHOLD));
})
.map(addNDWI)
.mean().clip(geometry);
var composite2 = s2
.filterBounds(geometry)
.filterDate('2024-06-01', '2024-08-01')
.linkCollection(csPlus, [QA_BAND])
.map(function(img) {
return img.updateMask(img.select(QA_BAND).gte(CLEAR_THRESHOLD));
})
.map(addNDWI)
.mean().clip(geometry);
var s2Viz = {bands: ['B4', 'B3', 'B2'], min: 0, max: 2000};
Map.addLayer(composite1, s2Viz, 'S2 composite1');
Map.addLayer(composite2, s2Viz, 'S2 composite2');
Map.centerObject(geometry);
var ndwiParams = {min: -0.7, max: 0.4, palette: ['red', 'yellow', 'green', 'cyan', 'blue']};
var S2_ndwi1 = composite1.select('ndwi');
var S2_ndwi2 = composite2.select('ndwi');
Map.addLayer(S2_ndwi1, ndwiParams, 'NDWI 1');
Map.addLayer(S2_ndwi2, ndwiParams, 'NDWI 2');
var ndwiThreshold1 = S2_ndwi1.gte(0.0);
var ndwiThreshold2 = S2_ndwi2.gte(0.0);
var ndwiMask1 = ndwiThreshold1.updateMask(ndwiThreshold1);
var ndwiMask2 = ndwiThreshold2.updateMask(ndwiThreshold2);
Map.addLayer(ndwiMask1, {palette:['blue']}, 'NDWI 1 Mask');
Map.addLayer(ndwiMask2, {palette:['blue']}, 'NDWI 2 Mask');
Для співставлення обрано два періоди: 2022 та 2024 роки. При усередненні для кожного року обрано період з 1 червня до 1 серпня, аби виключити хмарні частини знімків.
Спочатку наведені усереднені зображення знімків Sentinel-2 у видимому спектрі (з каналами ‘B4’, ‘B3’, ‘B2’ для RGB) для 2022 р. та 2024 р. відповідно:
На цих знімках видно масштаби змін на територіях Каховського водосховища та зміна русла р. Дніпро. По-друге, для аналізу за цей самий період та ці колекції знімків Sentinel-2 наведені візуалізації водного індексу NDWI для 2022 р. та 2024 р. відповідно:
Також у використаному коді є можливість навести маски для водних площ, для чого встановлено поріг водного індексу понад значення 0: var ndwiThreshold1 = S2_ndwi1.gte(0.0). Використовуючи цей поріг можна розрахувати площі води та її зміни. Детальніше з цим можна ознайомитись в статті: Francisco J. Peña. Water Detection Using NDWI on Google Earth Engine. https://medium.com/@melqkiades/water-detection-using-ndwi-on-google-earth-engine-2919a9bf1951
Залишити відповідь