【GEE】基于光谱距离度量方法的加沙地区地表覆盖变化检测

Twcat_tree

发布于 2024-04-25 14:47:30

770

发布于 2024-04-25 14:47:30

文章被收录于专栏：二猫の家二猫の家

以下文章来源于GEEer成长日记，作者_养乐多_

本文将介绍如何在 Google Earth Engine (GEE) 平台中使用光谱距离度量方法进行地表覆盖变化检测，并以加沙地区为例，使用Sentinel2数据展示2023年3月和2024年3月的地表覆盖变化区域。

结果如下图所示，

2. 代码链接

https://code.earthengine.google.com/5c2e0d8d82c1144a9907c58756228be0?noload=true

完整代码

var roi = ee.FeatureCollection('users/949384116/GAZA/1roiGAZA');
var bands = ['B2', 'B3', 'B4', 'B5', 'B6', 'B7', 'B8', 'B8A', 'B11', 'B12'];
Map.centerObject(roi, 11);
var roi1 = ee.Image().toByte().paint({ featureCollection: roi, color: 0, width: 3 });
Map.addLayer(roi1, {}, '边界')
//
var s2 = ee.ImageCollection('COPERNICUS/S2');
var rgbVis = {
    min: 0.0,
    max: 3000,
    bands: ['B4', 'B3', 'B2'],
};
var filtered = s2
    .filter(ee.Filter.bounds(roi))
    .select('B.*');
var csPlus = ee.ImageCollection('GOOGLE/CLOUD_SCORE_PLUS/V1/S2_HARMONIZED');
var csPlusBands = csPlus.first().bandNames();
var filteredS2WithCs = filtered.linkCollection(csPlus, csPlusBands);
function maskLowQA(image) {
    var qaBand = 'cs';
    var clearThreshold = 0.65;
    var mask = image.select(qaBand).gte(clearThreshold);
    return image.updateMask(mask);
}
var filteredMasked = filteredS2WithCs
    .map(maskLowQA);
var before = filteredMasked
    .filterDate('2023-03-01', '2023-03-31')
    .median();
var after = filteredMasked
    .filterDate('2024-03-01', '2024-03-31')
    .median();
Map.addLayer(before.clip(roi), rgbVis, '冲突前-2023年3月');
Map.addLayer(after.clip(roi), rgbVis, '冲突后-2024年3月');
var sed = after.spectralDistance(before, 'sed');
var distance = sed.sqrt();
var mask = distance.clip(roi).gt(1400);
var maskedImage = distance.clip(roi).updateMask(mask);
Map.addLayer(maskedImage, { min: 0, max: 1500, palette: ['white', 'red'] }, '变化较大的地方');