GEE案例:利用谷歌地球引擎中的哨兵-1 号合成孔径雷达数据,以无监督方式跟踪烧毁区域的进展情况

目录

简介

代码:数据预处理

数据预处理 

代码2:

代码3

 


简介

野火在全球范围内日益频繁,在过去二十年中造成了三分之一的森林损失。 使用传统的光学遥感技术跟踪烧毁面积的变化会受到云层和烟雾覆盖的影响。 因此,本研究采用了不受大气影响的多时合成孔径雷达(SAR)卫星数据。 该研究以 2021 年希腊的四场野火为重点,对位时和多时合成孔径雷达极坐标特征进行了无监督 k-means 聚类。 评估了不同移动核窗口大小的输入特征平滑化对提高准确性的影响。 使用这些平滑特征可大幅提高所有四个检测区域的准确性,而选择 19x19 的窗口大小则是为了在保留精细细节和减少斑点之间取得适当的平衡。 此外,在聚类后添加滤波器以去除小于 2 公顷的区域,还能进一步提高准确度,尤其是佣金误差。 使用所定义设置的结果显示,F1 分数为 0.75-0.88,总体准确率为 81-94%,遗漏/误差分别为 33-16% 和 14-3%。 在农业区占很大比例的地区观察到了挑战,而地形效应对结果没有实质性影响。 利用高分辨率图像进行的视觉验证证明了合成孔径雷达主要对较大的结构变化敏感的假设。 此外,还利用所介绍的方法开发了一个谷歌地球引擎(GEE)工具箱 "哨兵-1 烧毁区域进展"(Sentinel-1 Burned Area Progression,S1-BAP),该工具箱可在 GitHub 上免费供科学界使用。

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
GEE,可以使用Sentinel-2影像计算每月的NDVI,并且可以使用哨兵影像的数据来实现这个目的。下面是一个计算2018-2020年每月NDVI的代码实现: ```javascript // 定义计算NDVI函数 var addNDVI = function(image){ var ndvi = image.normalizedDifference(['B8', 'B4']).rename('NDVI'); return image.addBands(ndvi); }; // 定义日期范围 var startDate = '2018-01-01'; var endDate = '2020-12-31'; var startYear = 2018; var endYear = 2020; // 加载哨兵影像 var sentinel2 = ee.ImageCollection('COPERNICUS/S2_SR') .filterBounds(geometry) .filterDate(startDate, endDate) .map(addNDVI); // 将每个月的影像合并成一个影像 var months = ee.List.sequence(1, 12); var years = ee.List.sequence(startYear, endYear); var byMonthYear = ee.ImageCollection.fromImages( years.map(function (y) { return months.map(function (m) { var monthImage = sentinel2.filter(ee.Filter.calendarRange(y, y, 'year')) .filter(ee.Filter.calendarRange(m, m, 'month')) .reduce(ee.Reducer.mean()); return monthImage.set('year', y).set('month', m).set('system:time_start', ee.Date.fromYMD(y, m, 1)); }); }).flatten()); // 将每个月的NDVI作为一个图层添加到地图上 var ndviPalette = ['FFFFFF', 'CE7E45', 'DF923D', 'F1B555', 'FCD163', '99B718', '74A901', '66A000', '529400', '3E8601', '207401', '056201', '004C00', '023B01', '012E01', '011D01', '011301']; var ndviParams = {min: -0.2, max: 1, palette: ndviPalette}; byMonthYear = byMonthYear.map(function(image){ return image.clip(geometry); }); Map.addLayer(byMonthYear.select('NDVI'), ndviParams, 'NDVI'); ``` 这个代码实现首先定义了一个计算NDVI的函数`addNDVI`,然后加载了哨兵影像,将每个月的影像合并成一个影像,并将每个月的NDVI作为一个图层添加到地图上。其,`ndviPalette`是一个NDVI的颜色调色板,用于将NDVI值转换为对应的颜色。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

此星光明

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值