GEE:基于光谱距离度量方法的巴以冲突造成的地表覆盖变化检测

本文主要是介绍GEE:基于光谱距离度量方法的巴以冲突造成的地表覆盖变化检测,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

作者:CSDN @ _养乐多_

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

结果如下图所示,

在这里插入图片描述


文章目录

      • 一、核心函数
          • 1.1 spectralDistance函数
          • 1.2 光谱距离度量方法
      • 二、代码链接
      • 三、完整代码


一、核心函数

1.1 spectralDistance函数

spectralDistance 函数用于计算两个图像之间每个像素的光谱距离。

参数类型默认值描述
image1Image第一个图像
image2Image第二个图像
metricString“sam”光谱距离度量标准。可以是 ‘sam’ (光谱角映射器)、‘sid’ (光谱信息差异)、‘sed’ (平方欧氏距离) 或 ‘emd’ (地球移动者距离) 中的一个。

注意,如果图像是基于数组的,则只使用每个图像的第一个波段;否则,所有波段都参与距离计算。因此,需要这两个图像包含相同数量的波段或具有相同的一维数组长度。

1.2 光谱距离度量方法
方法描述
SAM光谱角制图(Spectral Angle Mapper,SAM)用来计算两个数组之间的相似性,其计算结果可看作两数组之间余弦角。如果两个光谱非常相似,则夹角接近零度;如果它们之间的差异很大,则夹角接近 90 或 180 度。
SID光谱信息散度(spectral information divergence,SID)是一种衡量两个光谱之间差异的指标。它基于信息论的概念,通过计算两个光谱分布之间的 KL 散度来量化它们之间的差异。SID 越高,表示两个光谱之间的差异越大。
SED平方欧几里得距离(Squared Euclidean distance,SED)是欧氏距离的平方。欧氏距离是两个点之间的直线距离,在光谱分析中用于衡量两个光谱向量之间的差异。SED 越大,表示两个光谱之间的差异越大。
EMD陆地移动距离(Earth Mover’s Distance,EMD)是一种衡量两个分布之间差异的度量。它考虑了两个分布之间的距离以及在将一个分布转换成另一个分布时所需的最小成本。在图像处理中,EMD 可以用于衡量两个图像的相似性。

二、代码链接

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'] }, '变化较大的地方');

这篇关于GEE:基于光谱距离度量方法的巴以冲突造成的地表覆盖变化检测的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/928431

相关文章

Window Server2016加入AD域的方法步骤

《WindowServer2016加入AD域的方法步骤》:本文主要介绍WindowServer2016加入AD域的方法步骤,包括配置DNS、检测ping通、更改计算机域、输入账号密码、重启服务... 目录一、 准备条件二、配置ServerB加入ServerA的AD域(test.ly)三、查看加入AD域后的变

Window Server2016 AD域的创建的方法步骤

《WindowServer2016AD域的创建的方法步骤》本文主要介绍了WindowServer2016AD域的创建的方法步骤,文中通过图文介绍的非常详细,对大家的学习或者工作具有一定的参考学习价... 目录一、准备条件二、在ServerA服务器中常见AD域管理器:三、创建AD域,域地址为“test.ly”

NFS实现多服务器文件的共享的方法步骤

《NFS实现多服务器文件的共享的方法步骤》NFS允许网络中的计算机之间共享资源,客户端可以透明地读写远端NFS服务器上的文件,本文就来介绍一下NFS实现多服务器文件的共享的方法步骤,感兴趣的可以了解一... 目录一、简介二、部署1、准备1、服务端和客户端:安装nfs-utils2、服务端:创建共享目录3、服

Java 字符数组转字符串的常用方法

《Java字符数组转字符串的常用方法》文章总结了在Java中将字符数组转换为字符串的几种常用方法,包括使用String构造函数、String.valueOf()方法、StringBuilder以及A... 目录1. 使用String构造函数1.1 基本转换方法1.2 注意事项2. 使用String.valu

Python中使用defaultdict和Counter的方法

《Python中使用defaultdict和Counter的方法》本文深入探讨了Python中的两个强大工具——defaultdict和Counter,并详细介绍了它们的工作原理、应用场景以及在实际编... 目录引言defaultdict的深入应用什么是defaultdictdefaultdict的工作原理

使用Python进行文件读写操作的基本方法

《使用Python进行文件读写操作的基本方法》今天的内容来介绍Python中进行文件读写操作的方法,这在学习Python时是必不可少的技术点,希望可以帮助到正在学习python的小伙伴,以下是Pyth... 目录一、文件读取:二、文件写入:三、文件追加:四、文件读写的二进制模式:五、使用 json 模块读写

Oracle数据库使用 listagg去重删除重复数据的方法汇总

《Oracle数据库使用listagg去重删除重复数据的方法汇总》文章介绍了在Oracle数据库中使用LISTAGG和XMLAGG函数进行字符串聚合并去重的方法,包括去重聚合、使用XML解析和CLO... 目录案例表第一种:使用wm_concat() + distinct去重聚合第二种:使用listagg,

Java后端接口中提取请求头中的Cookie和Token的方法

《Java后端接口中提取请求头中的Cookie和Token的方法》在现代Web开发中,HTTP请求头(Header)是客户端与服务器之间传递信息的重要方式之一,本文将详细介绍如何在Java后端(以Sp... 目录引言1. 背景1.1 什么是 HTTP 请求头?1.2 为什么需要提取请求头?2. 使用 Spr

Java如何通过反射机制获取数据类对象的属性及方法

《Java如何通过反射机制获取数据类对象的属性及方法》文章介绍了如何使用Java反射机制获取类对象的所有属性及其对应的get、set方法,以及如何通过反射机制实现类对象的实例化,感兴趣的朋友跟随小编一... 目录一、通过反射机制获取类对象的所有属性以及相应的get、set方法1.遍历类对象的所有属性2.获取

Java中的Opencv简介与开发环境部署方法

《Java中的Opencv简介与开发环境部署方法》OpenCV是一个开源的计算机视觉和图像处理库,提供了丰富的图像处理算法和工具,它支持多种图像处理和计算机视觉算法,可以用于物体识别与跟踪、图像分割与... 目录1.Opencv简介Opencv的应用2.Java使用OpenCV进行图像操作opencv安装j