草庐IT

GEE学习记录(一)基于GEE利用LANDSAT 8数据计算遥感生态指数(RSEI)

gxwwww 2023-04-19 原文

最近老师让看一下关于GEE的东西,实现大面积的反演、计算地表温度等,也算熟悉一下。参考网上很多大佬的文章,按照自己的思路和想法算出了RSEI,参考的文章都有列出来。

目录

所用数据集

影像数据

影像数据是:LANDSAT/LC08/C02/T1_L2(就是LANDSAT8下面的Collection 2 Level2 T1数据)。总的来说就是Collection 2比Collection 1好,T1比T2好,具体的不同可以见如下链接。数据介绍

以下是此数据集的波段信息,SR表示是地面反射率,不同于TOA顶层反射率,也就是TOA没有经过大气校正,而SR数据是经过大气校正的,更加适合进行地面反演。有些地方谈到不用进行大气校正也可以进行反演,结果也不会差太多,但是用校正过的肯定好一些。

图像包含 5 个可见光和近红外 (VNIR) 波段和 2 个短波红外 (SWIR) 波段( 没有8、9波段是因为USGS没有处理),处理为正射校正表面反射率,以及一个热红外 (TIR) 波段处理为正射校正表面温度。

注意第十波段的介绍:“地表温度”。因此可以直接用此波段进行地表温度的计算,在计算中有一些小点后面会谈到。

矢量数据

内蒙古锡林郭勒盟整个盟。

地表气温数据

MODIS/006/MOD11A2

RSEI计算流程图

并不是完全按照此流程图进行的,实际中就只做了“去云-裁剪-计算指标-掩膜水体-归一化-主成分分析-结果”。

四个指数的计算

Wet、NDVI、NDBSI的计算

以上部分参考如下文章:张元培等.基于遥感生态指数的钟祥市生态变化分析[J].资源环境与工程,2019,33(04):474-480.DOI:10.16536/j.cnki.issn.1671-1211.2019.04.006.

LST的计算

地表温度反演算法主要有以下三种:大气校正法(也称为辐射传输方程:Radiative Transfer Equation——RTE)、单通道算法和分裂窗算法。参考

在上面的影像数据介绍中谈到数据的第十波段,可以直接用来进行反演。以下是官方介绍中对数据进行修正的参数,修正完之后可以直接拿来用。


也比较好理解,对SR_B.也就是5 个可见光和近红外 (VNIR) 波段和 2 个短波红外 (SWIR) 波段进行第一行的处理,ST_B.*就是对热红外波段进行处理。

处理完之后(" ST_B10-273.15 ")即可得到地表温度。参考
得到结果之后发现,计算出来的相比于实际温度是有一点偏高的,而且出现了一些异常值,50℃、60℃肯定是不太正常的。


不想麻烦再用其他算法或者去除异常值,用GEE不就是因为数据集多、方便吗?所以就直接拿MODIS的数据进行计算了,这个肯定没有问题。

对水体进行掩膜

为什么要对水体进行掩膜呢?因为遥感生态指数主要应用于陆地为主的地区,不适宜大面积的水域地区。如果不进行掩膜,会出现异常值(血的教训)。参考:张元培等.基于遥感生态指数的钟祥市生态变化分析[J].资源环境与工程,2019,33(04):474-480.DOI:10.16536/j.cnki.issn.1671-1211.2019.04.006.
本来计算的大部分都是在正常范围之内的,但是突然出现了一个400多,后来一检查就是出现在水体的部分。

归一化

得到四个指数之后,除了温度,其他的指数范围都应该在[-1,1]之间,如果不是,可能出现了问题。之后进行归一化处理。将所有的值都归一化到[0,1]之间。下图是归一化前后的对比。

主成分分析


参考:徐涵秋.城市遥感生态指数的创建及其应用.生态学报,2013,33(24):7853—7862.
主成分分析的结果如上。

出图


原图:

代码

部分代码截图

完整代码
https://code.earthengine.google.com/7f41d4897f872c2bf7ae485db391bf92

参考链接

基本上列出来了,如果哪位大佬的有被参考却没被列出来,请及时私信我!
https://blog.csdn.net/qq_41664159/article/details/120522718
https://zhuanlan.zhihu.com/p/364850569
https://mp.weixin.qq.com/s/spi10OeHDRo6VpIrrpvYHA

有关GEE学习记录(一)基于GEE利用LANDSAT 8数据计算遥感生态指数(RSEI)的更多相关文章

  1. ruby - 解析 RDFa、微数据等的最佳方式是什么,使用统一的模式/词汇(例如 schema.org)存储和显示信息 - 2

    我主要使用Ruby来执行此操作,但到目前为止我的攻击计划如下:使用gemsrdf、rdf-rdfa和rdf-microdata或mida来解析给定任何URI的数据。我认为最好映射到像schema.org这样的统一模式,例如使用这个yaml文件,它试图描述数据词汇表和opengraph到schema.org之间的转换:#SchemaXtoschema.orgconversion#data-vocabularyDV:name:namestreet-address:streetAddressregion:addressRegionlocality:addressLocalityphoto:i

  2. ruby - Sinatra:运行 rspec 测试时记录噪音 - 2

    Sinatra新手;我正在运行一些rspec测试,但在日志中收到了一堆不需要的噪音。如何消除日志中过多的噪音?我仔细检查了环境是否设置为:test,这意味着记录器级别应设置为WARN而不是DEBUG。spec_helper:require"./app"require"sinatra"require"rspec"require"rack/test"require"database_cleaner"require"factory_girl"set:environment,:testFactoryGirl.definition_file_paths=%w{./factories./test/

  3. ruby - Ruby 有 `Pair` 数据类型吗? - 2

    有时我需要处理键/值数据。我不喜欢使用数组,因为它们在大小上没有限制(很容易不小心添加超过2个项目,而且您最终需要稍后验证大小)。此外,0和1的索引变成了魔数(MagicNumber),并且在传达含义方面做得很差(“当我说0时,我的意思是head...”)。散列也不合适,因为可能会不小心添加额外的条目。我写了下面的类来解决这个问题:classPairattr_accessor:head,:taildefinitialize(h,t)@head,@tail=h,tendend它工作得很好并且解决了问题,但我很想知道:Ruby标准库是否已经带有这样一个类? 最佳

  4. ruby-on-rails - Rails 5 Active Record 记录无效错误 - 2

    我有两个Rails模型,即Invoice和Invoice_details。一个Invoice_details属于Invoice,一个Invoice有多个Invoice_details。我无法使用accepts_nested_attributes_forinInvoice通过Invoice模型保存Invoice_details。我收到以下错误:(0.2ms)BEGIN(0.2ms)ROLLBACKCompleted422UnprocessableEntityin25ms(ActiveRecord:4.0ms)ActiveRecord::RecordInvalid(Validationfa

  5. ruby - 我如何添加二进制数据来遏制 POST - 2

    我正在尝试使用Curbgem执行以下POST以解析云curl-XPOST\-H"X-Parse-Application-Id:PARSE_APP_ID"\-H"X-Parse-REST-API-Key:PARSE_API_KEY"\-H"Content-Type:image/jpeg"\--data-binary'@myPicture.jpg'\https://api.parse.com/1/files/pic.jpg用这个:curl=Curl::Easy.new("https://api.parse.com/1/files/lion.jpg")curl.multipart_form_

  6. 世界前沿3D开发引擎HOOPS全面讲解——集3D数据读取、3D图形渲染、3D数据发布于一体的全新3D应用开发工具 - 2

    无论您是想搭建桌面端、WEB端或者移动端APP应用,HOOPSPlatform组件都可以为您提供弹性的3D集成架构,同时,由工业领域3D技术专家组成的HOOPS技术团队也能为您提供技术支持服务。如果您的客户期望有一种在多个平台(桌面/WEB/APP,而且某些客户端是“瘦”客户端)快速、方便地将数据接入到3D应用系统的解决方案,并且当访问数据时,在各个平台上的性能和用户体验保持一致,HOOPSPlatform将帮助您完成。利用HOOPSPlatform,您可以开发在任何环境下的3D基础应用架构。HOOPSPlatform可以帮您打造3D创新型产品,HOOPSSDK包含的技术有:快速且准确的CAD

  7. 叮咚买菜基于 Apache Doris 统一 OLAP 引擎的应用实践 - 2

    导读:随着叮咚买菜业务的发展,不同的业务场景对数据分析提出了不同的需求,他们希望引入一款实时OLAP数据库,构建一个灵活的多维实时查询和分析的平台,统一数据的接入和查询方案,解决各业务线对数据高效实时查询和精细化运营的需求。经过调研选型,最终引入ApacheDoris作为最终的OLAP分析引擎,Doris作为核心的OLAP引擎支持复杂地分析操作、提供多维的数据视图,在叮咚买菜数十个业务场景中广泛应用。作者|叮咚买菜资深数据工程师韩青叮咚买菜创立于2017年5月,是一家专注美好食物的创业公司。叮咚买菜专注吃的事业,为满足更多人“想吃什么”而努力,通过美好食材的供应、美好滋味的开发以及美食品牌的孵

  8. FOHEART H1数据手套驱动Optitrack光学动捕双手运动(Unity3D) - 2

    本教程将在Unity3D中混合Optitrack与数据手套的数据流,在人体运动的基础上,添加双手手指部分的运动。双手手背的角度仍由Optitrack提供,数据手套提供双手手指的角度。 01  客户端软件分别安装MotiveBody与MotionVenus并校准人体与数据手套。MotiveBodyMotionVenus数据手套使用、校准流程参照:https://gitee.com/foheart_1/foheart-h1-data-summary.git02  数据转发打开MotiveBody软件的Streaming,开始向Unity3D广播数据;MotionVenus中设置->选项选择Unit

  9. 使用canal同步MySQL数据到ES - 2

    文章目录一、概述简介原理模块二、配置Mysql使用版本环境要求1.操作系统2.mysql要求三、配置canal-server离线下载在线下载上传解压修改配置单机配置集群配置分库分表配置1.修改全局配置2.实例配置垂直分库水平分库3.修改group-instance.xml4.启动监听四、配置canal-adapter1修改启动配置2配置映射文件3启动ES数据同步查询所有订阅同步数据同步开关启动4.验证五、配置canal-admin一、概述简介canal是Alibaba旗下的一款开源项目,Java开发。基于数据库增量日志解析,提供增量数据订阅&消费。Git地址:https://github.co

  10. 基于C#实现简易绘图工具【100010177】 - 2

    C#实现简易绘图工具一.引言实验目的:通过制作窗体应用程序(C#画图软件),熟悉基本的窗体设计过程以及控件设计,事件处理等,熟悉使用C#的winform窗体进行绘图的基本步骤,对于面向对象编程有更加深刻的体会.Tutorial任务设计一个具有基本功能的画图软件**·包括简单的新建文件,保存,重新绘图等功能**·实现一些基本图形的绘制,包括铅笔和基本形状等,学习橡皮工具的创建**·设计一个合理舒适的UI界面**注明:你可能需要先了解一些关于winform窗体应用程序绘图的基本知识,以及关于GDI+类和结构的知识二.实验环境Windows系统下的visualstudio2017C#窗体应用程序三.

随机推荐