python:布伊山德U检验(Buishand U test,BUT)突变点检测(以NDVI时间序列为例)

作者:CSDN @ _养乐多_

本文将介绍布伊山德U检验(Buishand U test,BUT)突变点检测代码。以 NDVI 时间序列为例。输入数据可以是csv,一列NDVI值,一列时间。代码可以扩展到遥感时间序列突变检测(突变年份、突变幅度等)中。

结果如下图所示,
在这里插入图片描述


文章目录

      • 一、准备数据
      • 二、BUT介绍和代码
          • 2.1 原理和步骤
          • 2.2 核心函数
      • 三、读取csv格式时序数据的示例


一、准备数据

测试数据(0积分下载):https://download.csdn.net/download/qq_35591253/88895803

该数据是GEE上提取的,参考博客《GEE:基于Landsat5/7/8/9数据提取一个点的NDVI时间序列(1986-2024)》

二、BUT介绍和代码

Buishand U test突变点检测中文名为布伊山德U检验,而其原理是基于正态分布变量的单一变点模型。

2.1 原理和步骤

Buishand U test是用于检测时间序列数据中是否存在一个突变点,即某个时刻数据的特性发生了显著的变化。这种检测方法在气象、水文、生态和经济学等领域具有广泛的应用。

Buishand U test的原理是考虑一个正态随机变量X,并假设存在一个单变点将数据集分为两部分。在变点之前,数据遵循均值为μ和方差σ²的正态分布;在变点之后,数据的均值变为μ+δ。测试的零假设H₀是δ=0,即没有变点;备择假设H₁是δ≠0,即存在变点。为了进行检测,计算调整后的局部和统计量Sₖ,使用样本标准偏差D(x),然后构造U统计量:

U = 1 n ∗ ( n + 1 ) ∗ ∑ k = 1 n − 1 ( S [ k ] − D x ) 2 U = \frac{1}{n * (n + 1)} * \sum_{k=1}^{n-1} (S[k] - Dx)^2 U=n(n+1)1

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/734867.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

基于Token的身份验证:安全与效率的结合

🤍 前端开发工程师、技术日更博主、已过CET6 🍨 阿珊和她的猫_CSDN博客专家、23年度博客之星前端领域TOP1 🕠 牛客高级专题作者、打造专栏《前端面试必备》 、《2024面试高频手撕题》 🍚 蓝桥云课签约作者、上架课程《Vue.js 和 E…

前端文件上传

文件上传方式 前端文件上传有两种方式,第一种通过二进制blob传输(formData传输),第二种是通过base64传输 文件相关的对象 file对象其实是blob的子类 blob对象的第一个参数必须是一个数组,你可以把一个file对象放进去…

C语言连接【MySQL】

稍等更新图片。。。。 文章目录 安装 MySQL 库连接 MySQLMYSQL 类创建 MySQL 对象连接数据库关闭数据库连接示例 发送命令设置编码格式插入、删除或修改记录查询记录示例 参考资料 安装 MySQL 库 在 CentOS7 下,使用命令安装 MySQL: yum install mysq…

【蓝桥 2021】扫雷

扫雷 题目描述 在一个 n 行 m 列的方格图上有一些位置有地雷,另外一些位置为空。 请为每个空位置标一个整数,表示周围八个相邻的方格中有多少个地雷。 输入描述 输入的第一行包含两个整数 n,m。 第 2 行到第 n1 行每行包含 m 个整数,相…

关于查看 CentOS7虚拟机的 ip地址

1. 启动网卡 1.1 打开网卡配置文件。 vi /etc/sysconfig/network-scripts/ifcfg-eth01.2 启动网卡 修改为下图中的ONBOOTyes 2. 重启网络服务 sudo service network restart3. 查看ip地址 ip addr

【C/C++ 学习笔记】数组

【C/C 学习笔记】数组 视频地址: Bilibili 一维数组 数据类型 数组名[数组长度];数据类型 数组名[数组长度] { 值1, 值2, … }数据类型 数组名[] { 值1, 值2, … } 特点: 放在一块连续的内存空间数组中每个元素都是相同数据类型 数组名: 可以统计整个数组在内存中的长…

MySQL8.0数据库开窗函数

简介 数据库开窗函数是一种在SQL中使用的函数,它可以用来对结果集中的数据进行分组和排序,以便更好地分析和处理数据。开窗函数与聚合函数不同,它不会将多行数据聚合成一行,而是保留每一行数据,并对其进行分组和排序。…

HellaSwag数据集分享

来源: AINLPer公众号(每日干货分享!!) 编辑: ShuYini 校稿: ShuYini 时间: 2024-3-10 该数据集是由斯坦福大学研究人员提出的,用于评估NLP模型在常识自然语言推理(NLI)任务上的性能,…

凌鲨微应用架构

微应用是静态网页加上凌鲨提供的扩展能力而形成的一种应用,主要特点是开发便捷,安全。 微应用架构 组件说明 名称 说明 微应用 webview窗口,显示web服务器上的页面 接口过滤器 根据权限配置,屏蔽非授权接口访问 接口提供者 tauri注入…

一个开源免费的实时AI绘画软件Krita

Krita是由Krita插件ComfyUILCM插件结合的一个实时可预览的AI绘画工具,是一款自由开源、免费的专业级绘画软件,Krita采用的是GNU GPL许可证确保它将一直保持自由开源的本色。 一:搭建Krita创作平台 1、Krita客户端软件,该软件下载…

乡村治理深度解析:策略、挑战与解决方案

毋庸置疑,在今天这个崭新的时代,乡村治理的过程已然向我们发出了挑战。为了迎难而上,我们必须摒弃陈旧观念,勇敢迎接并大胆尝试探索与实践新的思路!为了达到这一宏伟目标,我们需要首先廓清如下关键概念&…

云计算 3月7号 (crontab-计划任务 日志轮转 免密登录)

一、计划任务 计划任务概念解析 在Linux操作系统中,除了用户即时执行的命令操作以外,还可以配置在指定的时间、指定的日期执行预先计划好的系统管理任务(如定期备份、定期采集监测数据)。RHEL6系统中默认已安装了at、crontab软件…

使用Rust开发小型搜索引擎

一、概述 用Rust创建搜索引擎是探索该语言在性能和安全性方面具有优势的绝佳方式。 这个项目将索引和搜索概念转移到Rust的生态系统中,由于Rust独特的语法和范式,这是一个挑战,但也是有益的。 二、构建搜索引擎 步骤1,创建项目…

新一代电话机器人正式上线

演示账户 代理商演示账户 地址:http://119.23.229.15:8080 用户名:c0508 密码:123456 功能:包含AI外呼管理,话术管理,CRM管理,坐席管理等功能。 管理员演示账户: 地址:h…

LeetCode 0299.猜数字游戏:计数

【LetMeFly】299.猜数字游戏:计数 力扣题目链接:https://leetcode.cn/problems/bulls-and-cows/ 你在和朋友一起玩 猜数字(Bulls and Cows)游戏,该游戏规则如下: 写出一个秘密数字,并请朋友猜…

随机森林原理sklearn实现

原理 定义 随机森林就是通过集成学习的思想将多棵树集成的一种算法,它的基本单元是决策树, 而它的本质属于机器学习的一大分支——集成学习(Ensemble Learning)方法。 随机森林的名称中有两个关键词,一个是“随机”&a…

拥塞控制的作用和方法

目录 1.概念 2.慢启动 3.快速恢复 4.快速重传 1.概念 TCP数据流的控制. 为了保证可靠传输和提高网络通信质量,内核需要对TCP数据流进行控制 TCP数据流控制的两个方面:超时重传和拥塞控制, 在计算机网络中的链路容量(即带宽)、交换结点中的缓存和处理机等,都是…

复合查询【MySQL】

文章目录 复合查询测试表 单表查询多表查询子查询单行子查询多行子查询IN 关键字ALL 关键字ANY 关键字 多列子查询 合并查询 复合查询 测试表 雇员信息表中包含三张表,分别是员工表(emp)、部门表(dept)和工资等级表&…

11 数据库优化

文章目录 数据库优化数据库设计范式MySQL存储引擎字段数据类型和键的选择explain语句SQL优化表的拆分 数据库优化 数据库设计范式 设计关系数据库时,遵从不同的规范要求,设计出合理的关系型数据库,这些不同的规范要求被称为不同的范式。 目…

命令行hive的基本操作

使用SQL语句在命令行创建Hive库: CREATE DATABASE database_name;其中,database_name是要创建的数据库的名称。例如: CREATE DATABASE mydatabase;这将在Hive中创建一个名为"mydatabase"的数据库。 使用SQL语句在命令行创建Hive…