GEE:变化检测算法,干扰检测,LandTrendr代码与说明(变化年份和扰动幅度示例)

本文主要是介绍GEE:变化检测算法,干扰检测,LandTrendr代码与说明(变化年份和扰动幅度示例),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

LandTrendr是什么?

LandTrendr是一组频谱时间分割算法,可用于中等分辨率卫星图像(主要是Landsat)的时间序列中的变化检测以及生成很大程度上不存在年际信号噪声的基于轨迹的频谱时间序列数据。LT最初是用IDL(交互式数据语言)实现的,但是在Google工程师的帮助下,它已移植到GEE平台。


LandTrendr官方说明网站:

https://emapr.github.io/LT-GEE/index.html
该网站介绍了在Google Earth Engine(GEE)中运行LandTrendr(LT)的基础知识。
它描述了LT概念框架,运行LT-GEE所需的内容,如何运行它,输出是什么以及如何格式化它们。它还提供了三个示例脚本,用于LT操作的一系列用户界面以及一个有目的但易于使用的API。
GEE框架几乎消除了IDL实施中繁琐的数据管理和图像预处理方面。它比IDL实施要快几光年,后者的计算时间以分钟而不是天为单位。(GEE谷歌云平台动用全球的服务器运算)


以下介绍动用API,运行LandTrendr的代码。输出结果为最大的扰动发生的幅度和年份。
完整代码:https://code.earthengine.google.com/65d284842e453de5731810be16bfccf3?noload=true
扰动幅度结果如下图所示,(不同颜色代表扰动强度不同)
在这里插入图片描述
扰动发生的年份,结果如下图所示:(不同颜色代表扰动发生的不同年份)
在这里插入图片描述

// 定义影像集参数
var startYear = 1984;
var endYear = 2019;
var startDay = '06-20';
var endDay = '09-01';
//var aoi = ee.Geometry.Point(-122.8848, 43.7929);
var index = 'NBR';
var maskThese = ['cloud', 'shadow', 'snow', 'water'];// 定义landtrendr参数
var runParams = { maxSegments:            6,spikeThreshold:         0.9,vertexCountOvershoot:   3,preventOneYearRecovery: true,recoveryThreshold:      0.25,pvalThreshold:          0.05,bestModelProportion:    0.75,minObservationsNeeded:  6
};// 定义变化参数
var changeParams = {delta:  'loss',sort:   'greatest',year:   {checked:true, start:1986, end:2019},mag:    {checked:true, value:200,  operator:'>'},dur:    {checked:true, value:4,    operator:'<'},preval: {checked:true, value:300,  operator:'>'},mmu:    {checked:true, value:11},};//各项参数设置好了,就可以开始运行LandTrendr了!// 加载API
var ltgee = require('users/emaprlab/public:Modules/LandTrendr.js'); // add index to changeParams object
changeParams.index = index;// 运行landtrendr
var lt = ltgee.runLT(startYear, endYear, startDay, endDay, aoi, index, [], runParams, maskThese);// 得到 change map layers
var changeImg = ltgee.getChangeMap(lt, changeParams);// 可视化字典
var palette = ['#9400D3', '#4B0082', '#0000FF', '#00FF00', '#FFFF00', '#FF7F00', '#FF0000'];
var yodVizParms = {min: startYear,max: endYear,palette: palette
};// 可视化参数
var magVizParms = {min: 200,max: 800,palette: palette
};// 将扰动图可视化
Map.centerObject(aoi, 11);
Map.addLayer(changeImg.select(['mag']), magVizParms, 'Magnitude of Change');
Map.addLayer(changeImg.select(['yod']), yodVizParms, 'Year of Detection');//设置显示样式:color代表边界颜色;fillcolor代表填充颜色
var styling={color:'red',fillColor:'00000000'}
//JX就是矢量边界
Map.addLayer(aoi.style(styling),{},"JX")// 将结果下载到谷歌云盘var exportImg = changeImg.clip(aoi).unmask(0).short();
Export.image.toDrive({image: exportImg, description: 'lt-gee_disturbance_map', folder: 'lt-gee_disturbance_map_test', fileNamePrefix: 'lt-gee_disturbance_map', region: aoi, scale: 30, crs: 'EPSG:5070', maxPixels: 1e13
});

这篇关于GEE:变化检测算法,干扰检测,LandTrendr代码与说明(变化年份和扰动幅度示例)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



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

相关文章

pandas中位数填充空值的实现示例

《pandas中位数填充空值的实现示例》中位数填充是一种简单而有效的方法,用于填充数据集中缺失的值,本文就来介绍一下pandas中位数填充空值的实现,具有一定的参考价值,感兴趣的可以了解一下... 目录什么是中位数填充?为什么选择中位数填充?示例数据结果分析完整代码总结在数据分析和机器学习过程中,处理缺失数

Pandas统计每行数据中的空值的方法示例

《Pandas统计每行数据中的空值的方法示例》处理缺失数据(NaN值)是一个非常常见的问题,本文主要介绍了Pandas统计每行数据中的空值的方法示例,具有一定的参考价值,感兴趣的可以了解一下... 目录什么是空值?为什么要统计空值?准备工作创建示例数据统计每行空值数量进一步分析www.chinasem.cn处

利用Python调试串口的示例代码

《利用Python调试串口的示例代码》在嵌入式开发、物联网设备调试过程中,串口通信是最基础的调试手段本文将带你用Python+ttkbootstrap打造一款高颜值、多功能的串口调试助手,需要的可以了... 目录概述:为什么需要专业的串口调试工具项目架构设计1.1 技术栈选型1.2 关键类说明1.3 线程模

Python Transformers库(NLP处理库)案例代码讲解

《PythonTransformers库(NLP处理库)案例代码讲解》本文介绍transformers库的全面讲解,包含基础知识、高级用法、案例代码及学习路径,内容经过组织,适合不同阶段的学习者,对... 目录一、基础知识1. Transformers 库简介2. 安装与环境配置3. 快速上手示例二、核心模

Python使用getopt处理命令行参数示例解析(最佳实践)

《Python使用getopt处理命令行参数示例解析(最佳实践)》getopt模块是Python标准库中一个简单但强大的命令行参数处理工具,它特别适合那些需要快速实现基本命令行参数解析的场景,或者需要... 目录为什么需要处理命令行参数?getopt模块基础实际应用示例与其他参数处理方式的比较常见问http

Java的栈与队列实现代码解析

《Java的栈与队列实现代码解析》栈是常见的线性数据结构,栈的特点是以先进后出的形式,后进先出,先进后出,分为栈底和栈顶,栈应用于内存的分配,表达式求值,存储临时的数据和方法的调用等,本文给大家介绍J... 目录栈的概念(Stack)栈的实现代码队列(Queue)模拟实现队列(双链表实现)循环队列(循环数组

Android实现在线预览office文档的示例详解

《Android实现在线预览office文档的示例详解》在移动端展示在线Office文档(如Word、Excel、PPT)是一项常见需求,这篇文章为大家重点介绍了两种方案的实现方法,希望对大家有一定的... 目录一、项目概述二、相关技术知识三、实现思路3.1 方案一:WebView + Office Onl

Mysql用户授权(GRANT)语法及示例解读

《Mysql用户授权(GRANT)语法及示例解读》:本文主要介绍Mysql用户授权(GRANT)语法及示例,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录mysql用户授权(GRANT)语法授予用户权限语法GRANT语句中的<权限类型>的使用WITH GRANT

Java中Scanner的用法示例小结

《Java中Scanner的用法示例小结》有时候我们在编写代码的时候可能会使用输入和输出,那Java也有自己的输入和输出,今天我们来探究一下,对JavaScanner用法相关知识感兴趣的朋友一起看看吧... 目录前言一 输出二 输入Scanner的使用多组输入三 综合练习:猜数字游戏猜数字前言有时候我们在

html5的响应式布局的方法示例详解

《html5的响应式布局的方法示例详解》:本文主要介绍了HTML5中使用媒体查询和Flexbox进行响应式布局的方法,简要介绍了CSSGrid布局的基础知识和如何实现自动换行的网格布局,详细内容请阅读本文,希望能对你有所帮助... 一 使用媒体查询响应式布局        使用的参数@media这是常用的