当前位置：首页 > news >正文

GEE计算 RSEI（遥感生态指数）

news 来源：原创 2025/6/20 1:27:43

🛰️ 什么是 RSEI？为什么要用它评估生态环境？

RSEI（遥感生态指数，Remote Sensing Ecological Index） 是一种通过遥感数据计算得到的、综合反映区域生态环境质量的指标体系。

它的设计初衷是用最少的变量，概括最多的生态信息，便于大尺度、自动化地评估城市或自然区域的生态健康状态。

🔍 RSEI 的四大核心因子：

因子	代表意义	常用遥感指数/数据源
NDVI	植被覆盖度，代表“绿意”	NIR与红光构建
Wet	土壤/植被含水量，代表“湿润”	Tasseled Cap 湿度因子
NDBSI	裸土+建成区比例，代表“干扰”	反映裸地和城市建设强度
LST	地表温度，代表“热压力”	Landsat 热红外波段反演地表温度

这些因子经过标准化处理后，通过**主成分分析（PCA）**提取出主成分，再归一化成 0–1 范围的指数值，就是最终的 RSEI。

值越大，代表生态质量越好；越接近 0，说明该区域可能存在生态退化、城市扩张、土地干扰等问题。

🛠️ 如何用 GEE 实现 RSEI 计算？（完整流程）

我们以 Landsat 5 2000 年夏季影像为例，结合四川西部地区，进行了如下步骤：

1️⃣ 数据准备

使用 GEE 调用 Landsat 5 L2 SR 产品（含大气校正）
设置研究时间为 6~9月（植被生长旺季），增强 NDVI 信噪比
设置研究区域（AOI）为上传的 chuanxi 区域文件

2️⃣ 遥感因子计算

NDVI：NIR 与 RED 构建标准植被指数
Wet：用蓝、绿、红、NIR、SWIR 构建湿度因子（Tasseled Cap 线性组合）
NDBSI：结合 SI、IBI 等多个裸地指数构建综合干扰指标
LST：基于热红外波段，计算辐射亮温并反演地表温度，考虑地表发射率影响

所有因子都经过 标准差标准化（非简单归一化），以增强 PCA 分析的鲁棒性

3️⃣ 主成分分析（PCA）

将 4 个标准化指标合并为一个影像堆栈
使用 GEE 的 centeredCovariance 计算协方差矩阵
求解特征值与特征向量，获取 PC1（第一主成分）作为生态质量主因子
反转 PC1，生成原始的生态指数影像（RSEI0）

4️⃣ 归一化处理

提取 PC1 的 1%-99% 分位值范围
将 RSEI0 线性归一化到 0~1 区间，生成最终 RSEI

5️⃣ 可视化与导出

使用配色方案（绿好、红差）展示生态分布差异
导出 GeoTIFF 文件，支持 ArcGIS/QGIS 进一步分析

仅部分代码，需完整版代码关注全域智图后台私信

仅部分代码，需完整版代码关注全域智图工作室后台私信
var year = '2000';
var startDate = ee.Date(year + '-06-01');
var endDate = ee.Date(year + '-09-30');
var aoi = ee.FeatureCollection("projects/ee-xu135150/assets/chuanxi");
Map.centerObject(aoi, 6);// 1. Landsat5 预处理 
function maskL5sr(image) {var qa = image.select('QA_PIXEL');var mask = qa.bitwiseAnd(1 << 4).eq(0).and(qa.bitwiseAnd(1 << 3).eq(0));return image.updateMask(mask).select(['SR_B1','SR_B2','SR_B3','SR_B4','SR_B5','SR_B7','ST_B6']).copyProperties(image, ['system:time_start']);
}var l5 = ee.ImageCollection("LANDSAT/LT05/C02/T1_L2").filterDate(startDate, endDate).filterBounds(aoi).map(maskL5sr).median().clip(aoi);// 检查影像像素覆盖
print("Landsat 5 影像像素计数", l5.reduceRegion({reducer: ee.Reducer.count(),geometry: aoi,scale: 200,maxPixels: 1e13
}));// 波段转换
var B = l5.select('SR_B1').multiply(0.0000275).add(-0.2);
var G = l5.select('SR_B2').multiply(0.0000275).add(-0.2);
var R = l5.select('SR_B3').multiply(0.0000275).add(-0.2);
var NIR = l5.select('SR_B4').multiply(0.0000275).add(-0.2);
var SWIR1 = l5.select('SR_B5').multiply(0.0000275).add(-0.2);-----------------------------------部分代码----------------------------------------// 打印 RSEI 统计
print('RSEI 值域:', RSEI.reduceRegion({reducer: ee.Reducer.minMax(), geometry: aoi, scale: 200, maxPixels: 1e13}));
print('RSEI 均值:', RSEI.reduceRegion({reducer: ee.Reducer.mean(), geometry: aoi, scale: 200, maxPixels: 1e13}));// 5. 可视化 
Map.addLayer(NDVI, {min: -2, max: 2, // 标准化后的值域palette: ['#d9f0d3', '#addd8e', '#31a354']
}, 'NDVI (植被指数)');Map.addLayer(Wet, {min: -2, max: 2,palette: ['#f7fbff', '#6baed6', '#08306b']
}, 'Wet (湿度指数)');Map.addLayer(NDBSI, {min: -2, max: 2,palette: ['#ffffcc', '#fd8d3c', '#bd0026']
}, 'NDBSI (裸土指数)');Map.addLayer(LST, {min: -2, max: 2,palette: ['#ffffcc', '#fc8d59', '#d73027']
}, 'LST (温度)');Map.addLayer(RSEI.updateMask(RSEI), {min: 0, max: 1,palette: ['d7191c', 'fdae61', 'ffffbf', 'a6d96a', '1a9641']
}, 'RSEI (生态质量指数)');// 6. 导出 RSEI 
Export.image.toDrive({image: RSEI,description: 'RSEI_' + year,folder: 'RSEI_Export',scale: 200,region: aoi,maxPixels: 1e13,fileFormat: 'GeoTIFF'
});

GEE计算 RSEI（遥感生态指数）

🛰️ 什么是 RSEI？为什么要用它评估生态环境？ RSEI（遥感生态指数，Remote Sensing Ecological Index） 是一种通过遥感数据计算得到的、综合反映区域生态环境质量的指标体系。它的设计初衷是用最少的变量&…...

编程日记 2025/6/20 1:27:43

k8s监控方案实践补充（二）：使用kube-state-metrics获取资源状态指标

k8s监控方案实践补充（二）：使用kube-state-metrics获取资源状态指标文章目录 k8s监控方案实践补充（二）：使用kube-state-metrics获取资源状态指标一、Metrics Server简介二、kube-state-metrics实战部署1. 创…...

编程日记 2025/6/20 0:52:00

JavaScript：PC端特效--元素可视区client系列

一、client系列 client翻译过来就是客户端，我们使用client系列的相关属性来获取元素可视区的相关信息。通过client系列的相关属性可以动态的得到该元素的边框大小、元素大小等。 client系列属性作用element.clientTop返回元素上边框的大小element.clientLeft返回元…...

编程日记 2025/6/20 1:11:34

Centos7 中 Docker运行配置Apache

1、拉去httpd镜像（不加版本号，默认拉最新版本） docker pull httpd 2、运行httpd docker run -di --name httpd-test -p 8080:80 httpd 3、创建文件夹后边做映射 mkdir -p /Docker/apache/www /Docker/apache/logs /Docker/apache/conf 4、…...

编程日记 2025/6/20 1:13:43

PostgreSQL中的全页写

一、概述在PGSQL数据库中，默认的页面大小为8KB，但是磁盘buffer的大小为4KB，扇区大小为512B。这就导致在操作系统的角度看数据库的写操作，其实并不是一种原子操作。如果操作系统发生了系统级别的故障，此时正好操作系统…...

编程日记 2025/5/27 7:38:01

对称二叉树的判定：双端队列的精妙应用

一、题目解析题目描述给定一个二叉树，检查它是否是镜像对称的。例如，二叉树 [1,2,2,3,4,4,3] 是对称的： 1/ \2 2/ \ / \ 3 4 4 3而 [1,2,2,null,3,null,3] 则不是镜像对称的： 1/ \2 2\ \3 3问题本质判断一棵二叉…...

编程日记 2025/5/26 14:54:06

Redis + ABP vNext 构建分布式高可用缓存架构

🚀 Redis ABP vNext 构建分布式高可用缓存架构 🔧 环境准备开发环境 .NET 8.0 SDKVisual Studio 2022 / VS CodeDocker & Docker Compose NuGet 包 Volo.Abp.Caching.StackExchangeRedis v8.1.5Volo.Abp.DistributedLocking.StackExchangeRedis v…...

编程日记 2025/5/27 7:38:49

jvm第一篇《内存与垃圾回收》学习笔记第一章jvm初始

jvm是虚拟机的通称。 java实际默认的应用是hotspot（基于栈的指令集架构） 注：注意区分寄存器的指令集和栈指令集的架构。（大概理解java移植性好就是因为是栈指令集） jvm虚拟机，具有跨语言功能&#xff0…...

编程日记 2025/6/1 9:24:35

MySQL——3、数据类型

数据类型 1、数据类型分类2、数值类型2.1、tinyint类型2.2、bit类型2.3、小数类型2.3.1、float2.3.2、decimal 3、字符串类型3.1、char3.2、varchar3.3、char和varchar比较3.4、日期和时间类型3.5、enum和set 1、数据类型分类 2、数值类型 2.1、tinyint类型首先创建t1表&…...

编程日记 2025/6/4 9:44:51

Flutter - 集成三方库:日志(logger)

日志使用print方法时,会提示添加依赖 $ flutter pub add logger下载依赖 $ flutter pub get使用打印 import package:logger/logger.dart;var logger Logger(); logger.d("debug"); logger.e("error"); logger.i("info"); logger.f(&qu…...

编程日记 2025/5/29 17:14:24

第五部分：第五节 - Express 路由与中间件进阶：厨房的分工与异常处理

随着你的 Express 应用变得越来越大，所有的路由和中间件都写在一个文件里会变得难以管理。这时候就需要将代码进行拆分和组织。此外，一个健壮的后端应用必须能够优雅地处理错误和一些常见的 Web 开发问题，比如跨域。路由模块化 (express.Ro…...

编程日记 2025/5/27 6:16:46

国标GB/T 12536-90滑行试验全解析：纯电动轻卡行驶阻力模型参数精准标定

摘要本文以国标GB/T 12536-90为核心框架，深度解析纯电动轻卡滑行试验的完整流程与数据建模方法，提供： 法规级试验规范：从环境要求到数据采集全流程详解行驶阻力模型精准标定：最小二乘法求解 ( FAv^2BvC ) 的MATLAB实…...

编程日记 2025/5/26 10:41:45

组件导航 (Navigation)+flutter项目搭建-混合开发+分栏

组件导航 (Navigation)flutter项目搭建接上一章flutter项目的环境变量配置并运行flutter 上一章面熟了搭建flutter并用编辑器运行了ohos项目，这章主要是对项目的工程化改造先创建flutter项目，再配置Navigation 1.在开发视图的resources/base/profi…...

编程日记 2025/5/29 18:35:56

物联网中的WiFi模式解析：AP、STA与混合模式

物联网现在还是比较火的，各种设备都要联网，那么WiFi已成为设备联网的“标配”。但你是否想过，为什么有的设备能自己创建WiFi热点，有的只能连接路由器？为什么有些网关既能收数据又能传数据？ 主要还是因为Wi…...

编程日记 2025/5/28 14:35:02

spring cloud gateway 源码解析

参考：Spring Cloud Gateway SpringCloud gateway源码走读（顺带聊聊响应式） - 掘金 1，原理图还是从starter 开始看要实现网关的核心概念，肯定是需要接受请求的server ,从上面的截图看 starter-gateway 只负责了包的依赖，并没有定义自动配置，他依赖了starter-webf…...

编程日记 2025/5/27 20:24:01

游戏引擎学习第286天:开始解耦实体行为

回顾并为今天的内容定下基调我们目前正在进入实体系统的一个新阶段，之前我们已经让实体的移动系统变得更加灵活，现在我们想把这个思路继续延伸到实体系统的更深层次。今天的重点，是重新审视我们处理实体类型（entity type&#x…...

编程日记 2025/5/23 16:17:46

【论文阅读】KIMI K1.5: SCALING REINFORCEMENT LEARNING WITH LLMS

KIMI K1.5: SCALING REINFORCEMENT LEARNING WITH LLMS Scaling的解释： 通过系统性的方法扩展强化学习算法的能力，使其能够处理更复杂的问题、更大的状态/动作空间、更长的训练周期或更高效的资源利用原文摘要： 研究背景与问题定位传统预训…...

编程日记 2025/5/27 5:03:57

大语言模型与多模态模型比较

一、核心差异：输入数据类型与模态融合输入数据类型 LLM：仅处理文本数据，例如文本分类、机器翻译、问答等任务，通过大规模语料库学习语言规律。 LMM：支持文本、图像、音频、视频等多种模态输入，例如根据图…...

编程日记 2025/5/22 19:14:52

vscode debug node + 前端

方法 2：调试全栈（Node 前端） 如果需同时调试后端和前端： 分别启动两个调试会话一个配置调试 Node.js 后端（server.js）。另一个配置调试浏览器前端（如上）。 {// Use IntelliS…...

编程日记 2025/5/25 9:17:41

RK3568-鸿蒙5.1与原生固件-扇区对比分析

编译生成的固件目录地址 ../openharmony/out/rk3568/packages/phone/images鸿蒙OS RK3568固件分析通过查看提供的信息，分析RK3568开发板固件的各个组件及其用途： 主要固件组件根据终端输出的文件列表，RK3568固件包含以下关键组件&#x…...

编程日记 2025/5/25 13:18:13

Java线程池（Thread Pool）性能优化解析

在高性能、高并发的Java应用开发中，线程池（Thread Pool）是不可或缺的组件。它通过复用线程，避免了线程频繁创建和销毁带来的资源开销，提高了系统的响应速度和稳定性。然而，不合理的线程池配置和使用方式也可能成为系统性能瓶颈的根源。本文旨在深入解析Java线程池的性能…...

编程日记 2025/5/25 14:01:12

AI重塑未来学者：研究生教育的“进化论”与“数字化生存指南

目录：一、引言：AI浪潮下的“象牙塔”新挑战与新机遇二、AI的“双刃剑”：深度剖析对研究生教育的颠覆性影响 1. 研究范式的革新：从“人工”到“智能” 2. 知识获取与传授方式的重塑 3. 创新能力与批判性思维的再定义 4. 伦理困境与学术诚信的新考验三、他山之石：发达国家…...

编程日记 2025/5/25 11:33:44

IHttpHandler和Tcp Listener的web服务器接收上传文件有什么区别

IHttpHandler和Tcp Listener的web服务器接收上传文件有什么区别 IHttpHandler 与 TCP Listener 处理文件上传的核心区别 IHttpHandler 和 TcpListener 是ASP.NET中处理 HTTP 请求的两种不同抽象层级，它们在文件上传处理上存在以下关键区别： 1. 抽象层…...

编程日记 2025/5/31 5:57:31

C++ --- new与delete

new与delete 一、回顾1.malloc2.calloc3.realloc4.free 二、new与delete的特殊之处（1）（2） 三、new与delete的底层原理四、总结一、回顾在C语言阶段我们学习了动态内存管理：malloc,calloc,realloc,free。 1.malloc …...

编程日记 2025/5/26 5:50:31

Guided Filtering相关记录

一、背景介绍以前折腾保边滤波时候，刷了一些Guided Filtering相关资料。这里主要是对它们做个算法效果复现和资料简单整理。二、Guided Filtering 1、基本原理原版Guided Filtering的提出，主要是为了改善双边滤波做保边平滑滤波器时候的梯度翻转伪影…...

编程日记 2025/5/25 12:38:17

Rust 学习笔记：关于 String 的练习题

Rust 学习笔记：关于 String 的练习题 Rust 学习笔记：关于 String 的练习题选出描述正确的那一个。该程序最多可能发生多少次堆的内存分配？哪种说法最能解释为什么 Rust 不允许字符串索引？哪种说法最能描述字符串切片 &str 和字…...

编程日记 2025/5/27 4:51:50

AG-UI 协议：重构多模态交互，开启智能应用新纪元

一、协议诞生的时代背景：填补 AI 生态最后一块拼图在人工智能技术飞速发展的今天，AI 代理（Agent）作为能够主动执行复杂任务的智能实体，正从实验室走向生产环境，重塑各个行业的工作流程。然而，…...

编程日记 2025/5/30 15:42:41

网络世界的“百变身份“：动态IP让连接更自由

深夜的程序调试凌晨两点，我盯着电脑屏幕上的报错信息："Connection timed out"。这是本周第三次测试服务器响应时被拒绝访问了——只因为之前同一个IP地址尝试登录太过频繁。正在改代码的朋友小王凑过来看了眼："老兄&…...

编程日记 2025/5/31 14:50:22

【学习笔记】因果推理导论第1课

因果推理导论第1课为何要做因果推理一、辛普森悖论一个例子二、相关不代表因果性三、什么揭示因果四、观测研究五、结论本节课通过一、辛普森悖论一个例子书中举了一个疫情两种治疗方法A,B,分析哪一个治疗方法更好的例子. 首先已知B治疗方法更稀缺,因此观测数据样本上…...

编程日记 2025/5/28 15:34:33

Android 中权限分类及申请方式

在 Android 中，权限被分为几个不同的类别，每个类别有不同的申请和管理方式。一、普通权限（Normal Permissions）普通权限通常不会对用户隐私或设备安全造成太大风险。这些权限在应用安装时自动授予，无需用户在运行时手动授权。 android.permission.INTERNETandroid.pe…...

编程日记 2025/5/29 9:17:27

深度学习算法介绍

深度学习算法是一种基于人工神经网络结构的机器学习方法，其核心理念是通过多层次的神经元组成的神经网络来模拟人类大脑的工作原理。以下是几种常见的深度学习算法及其简要介绍： 多层感知器（Multilayer Perceptron, MLP）&#xff…...

编程日记 2025/5/19 22:56:12

Java【13_1】final、初始化块、继承(测试题)

测试题 1、简述final修饰符的功能 ① 修饰类该类不能被继承 ② 修饰变量该变量就是常量(一旦被初始化，就不可以修改) ③ 修饰方法该方法不能被重写 2、写出程序结果 (仔细认真) public class MyClass { static int x,y; static{ …...

编程日记 2025/5/29 17:58:59

小结：JavaScript 模块化工具链

JavaScript 模块化工具链是现代前端开发中用于组织、管理和优化模块化代码的核心工具集合。以下是关于 JS 模块化工具链的概述，包括关键工具、作用和常见工作流程：** **1. **模块化的背景 JavaScript 模块化是为了解决代码组织、依赖管理和作用域隔离…...

编程日记 2025/6/4 15:39:21

RabbitMQ ④-持久化 || 死信队列 || 延迟队列 || 事务

消息确认机制简单介绍 RabbitMQ Broker 发送消息给消费者后，消费者处理该消息时可能会发生异常，导致消费失败。如果 Broker 在发送消息后就直接删了，就会导致消息的丢失。为了保证消息可靠到达消费者并且成功处理了该消息，…...

编程日记 2025/5/31 7:14:04

十一、Hive JOIN 连接查询

作者：IvanCodes 日期：2025年5月16日专栏：Hive教程在数据分析的江湖中，数据往往分散在不同的“门派”（表）之中。要洞察数据间的深层联系，就需要JOIN这把利器，将相关联的数据串联起来…...

编程日记 2025/5/28 13:51:47

启用rvzi可视化自己的机器人发现joint state publisher gui没有滑块

启用rvzi可视化自己的机器人发现joint state publisher gui没有滑块？ 解决方法： 一：查看urdf中joint type定义是不是revolute。二：查看urdf关节限制是不是正确，如果是0到0则不正确。 <joint name"joint_…...

编程日记 2025/5/29 10:35:16

Gitee DevOps：中国企业数字化转型的“本土化加速器“

在数字化浪潮席卷全球的当下，DevOps已经从技术热词转变为企业的核心生产力工具。根据IDC最新报告，到2025年中国DevOps市场规模将达到15亿美元，年复合增长率高达25%。在这一快速增长的市场中，一个显著趋势正在显现：越来…...

编程日记 2025/5/29 11:12:34

RKNN开发环境搭建（ubuntu22.04）

以下情况在RV1106G3的平台上验证正常。 1、conda安装 1）conda --version//确认是否安装 2）创建一个安装目录，进行下一步 3）wget https://mirrors.tuna.tsinghua.edu.cn/anaconda/miniconda/Miniconda3-4.6.14-Linux-x…...

编程日记 2025/5/29 17:10:57

2025年上软考考试时间+准考证打印全攻略

各位同学好呀！我是你们的老朋友大老李~之前咱们聊过软考报名和机考模拟系统的使用技巧。今天要和大家分享一下，最近好多同学私信的问题："老李！准考证到底哪天能打印？""考试时间会不会有变动&#xff1f…...

编程日记 2025/5/31 11:31:10

vllm量化05—FP8 W8A8

本系列基于Qwen2.5-7B，学习如何使用vllm量化，并使用benchmark_serving.py、lm_eval 测试模型性能和评估模型准确度。测试环境为： OS: centos 7 GPU: nvidia l40 driver: 550.54.15 CUDA: 12.3本文是该系列第5篇——FP8 W8A8 一、量化 fro…...

编程日记 2025/5/27 6:13:35

OpenCV 图像透视变换详解

在计算机视觉领域，图像的视角问题常常会影响后续的分析与处理。例如，从倾斜角度拍摄的文档、带有畸变的场景图像等，都需要通过特定的方法进行矫正。OpenCV 作为计算机视觉领域的重要库，提供了强大的图像透视变换功能，能…...

编程日记 2025/5/26 8:04:53

Vue主题色切换实现方案(CSS 变量 + 类名切换)

定义主题变量 // styles/themes.scss :root {--primary-color: #314099;--secondary-color: #1E3968;--text-color: #2c44ce; }[data-theme"红系主题"] {--primary-color: #d74146;--secondary-color: #c20707;--text-color: #db3b3b; }在组件中使用变量 <templ…...

编程日记 2025/5/29 11:53:09

MySQL 高可用

目录一、概述MySQL高可用 1. 什么是高可用性（High Availability, HA）？ 2. MySQL 高可用常见方案 3. 高可用核心机制 4. 选型建议注意事项二、案例环境组成 MySQL 主主复制 keepalived HAProxy 的高可用负载均衡架构三、案例架构…...

编程日记 2025/5/29 9:34:28

C#学习教程（附电子书资料）

概述 C#（读作"C Sharp"）是一种由微软开发的现代编程语言，结合了C的高效性和Java的简洁性，专为.NET框架设计。以下是其核心特性和应用领域的详细介绍电子书资料：https://pan.quark.cn/s/6fe772420f95 一、语…...

编程日记 2025/5/25 10:35:55

MySQL性能优化

目录一、索引优化 1、慢查询日志分析 2、EXPLAIN 执行计划分析 3、索引类型选择 4、索引使用原则 5、常见索引失效场景二、SQL语句优化 1、避免低效操作符 2、减少数据扫描量 3、子查询优化 4、其他高频优化技巧三、表设计优化 1、数据类型优化四、架构设计优…...

编程日记 2025/5/27 7:49:42

Mendix 中的XPath 令牌（XPath Tokens）详解

在 Mendix 中，XPath 令牌（XPath Tokens） 是一种特殊的动态参数化查询技术，允许你在 XPath 表达式中使用变量或上下文相关的值，从而实现更灵活的查询逻辑。 1. 什么是 XPath 令牌？ XPath 令牌是 Mendix 提…...

编程日记 2025/5/26 10:51:32

Feign异步模式丢失上下文问题

Feign异步模式丢失上下文问题问题描述当我们使用异步对我们代码进行操作优化时，代码中使用了RequestContextHolder去获取上下文的数据，当我们执行原来可以执行的业务时发现报了空指针异常或数据为空，这是为什么呢？ 原理解释 …...

编程日记 2025/5/26 9:12:45

保姆教程-----安装MySQL全过程

1.电脑从未安装过mysql的，先找到mysql官网：MySQL :: Download MySQL Community Server 然后下载完成后，找到文件，然后双击打开 2. 选择安装的产品和功能依次点开“MySQL Servers”、“MySQL Servers”、“MySQL Servers 5.7”、…...

编程日记 2025/5/26 6:16:21