GEE详细教程之:将Landsat8与Landsat9影像合成一个影像

2024-02-09 13:28

本文主要是介绍GEE详细教程之:将Landsat8与Landsat9影像合成一个影像,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

1.前言

        因项目需求,需要获取一个研究区的Landsat8影像,但Landsat8重复周期长,加之天气的影响,很难获取影像质量较好的影像。Landsat4/5/7的波段顺序与landsat8不同,除此之外,landsat7影像还需要工具进行条带修复,因此只考虑最近发射的landsat9进行融合。

        Landsat9,就号称可以和Landsat8协同,降低重访周期。考虑到使用同一系列卫星的场景比较多,所以这篇博文就使用Landsat8与Landsat9卫星的影像合成。

2.思路

        (1)首先,要获取一个包含需要的时间段/影像质量/区域范围的Landsat8系列影像集合。

        (2)然后,再获取一个包含需要的时间段/影像质量/区域范围的Landsat9系列影像集合。

        (3)最后,将两个集合的影像通过merge工具进行拼接。

        效果图如下所示:

(1)只使用landsat8

(2)只使用landsat9

(3)融合landsat8与landsat9

        可以看出,融合后的数据相当于是两者数据的并集,可以将缺失的数据填充上,形成完整的一张图。

3.代码

var table = ee.FeatureCollection("注意!替换成你自己的矢量边界");// 导入 Landsat 8 和 9 的影像集合
var landsat8 = ee.ImageCollection('LANDSAT/LC08/C02/T1_L2');
var landsat9 = ee.ImageCollection('LANDSAT/LC09/C02/T1_L2');// 设置时间范围
var startDate = '2023-11-14';
var endDate = '2023-11-28';// 选择时间范围内的影像
landsat8 = landsat8.filterDate(startDate, endDate).filter(ee.Filter.lte('CLOUD_COVER',30)).filterBounds(table).map(applyScaleFactors);
landsat9 = landsat9.filterDate(startDate, endDate).filter(ee.Filter.lte('CLOUD_COVER',30)).filterBounds(table).map(applyScaleFactors);function applyScaleFactors(image) {var opticalBands = image.select('SR_B.').multiply(0.0000275).add(-0.2);var thermalBands = image.select('ST_B.*').multiply(0.00341802).add(149.0);return image.addBands(opticalBands, null, true).addBands(thermalBands, null, true);
}
print(landsat8);
print(landsat9);//landsat8=landsat8.select(['SR_B4', 'SR_B3', 'SR_B2']);
//print(landsat8)// 转换为 ImageCollection
var mergedImageCollection = landsat8.merge(landsat9);// 打印合并后的影像集合
print(mergedImageCollection);var mergedImage = mergedImageCollection.median().clip(table);
print("mergedImage",mergedImage);// 在地图上显示合并后的影像集合
Map.addLayer(mergedImage, {bands: ['SR_B4', 'SR_B3', 'SR_B2'], min: 0, max: 0.3}, 'Merged RGB');// 设置地图视图
Map.centerObject(table, 6);

4.结尾

        此时,如果遇到云雨比较频繁的地区云、贵、川等地区,可能会出现Landsat8与Landsat9融合后仍出现缺失的问题,这就需要把Landsat与哨兵二号进行融合。

        本文后续会将landsat与sentinel-2 融合的资源代码整理出来,提供给有需要的小伙伴。代码为JS格式,直接复制到你的GEE code Editor,并替换成自己的矢量区域即可使用~

这篇关于GEE详细教程之:将Landsat8与Landsat9影像合成一个影像的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



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

相关文章

全网最全Tomcat完全卸载重装教程小结

《全网最全Tomcat完全卸载重装教程小结》windows系统卸载Tomcat重新通过ZIP方式安装Tomcat,优点是灵活可控,适合开发者自定义配置,手动配置环境变量后,可通过命令行快速启动和管理... 目录一、完全卸载Tomcat1. 停止Tomcat服务2. 通过控制面板卸载3. 手动删除残留文件4.

Python中isinstance()函数原理解释及详细用法示例

《Python中isinstance()函数原理解释及详细用法示例》isinstance()是Python内置的一个非常有用的函数,用于检查一个对象是否属于指定的类型或类型元组中的某一个类型,它是Py... 目录python中isinstance()函数原理解释及详细用法指南一、isinstance()函数

Python的pandas库基础知识超详细教程

《Python的pandas库基础知识超详细教程》Pandas是Python数据处理核心库,提供Series和DataFrame结构,支持CSV/Excel/SQL等数据源导入及清洗、合并、统计等功能... 目录一、配置环境二、序列和数据表2.1 初始化2.2  获取数值2.3 获取索引2.4 索引取内容2

uni-app小程序项目中实现前端图片压缩实现方式(附详细代码)

《uni-app小程序项目中实现前端图片压缩实现方式(附详细代码)》在uni-app开发中,文件上传和图片处理是很常见的需求,但也经常会遇到各种问题,下面:本文主要介绍uni-app小程序项目中实... 目录方式一:使用<canvas>实现图片压缩(推荐,兼容性好)示例代码(小程序平台):方式二:使用uni

python依赖管理工具UV的安装和使用教程

《python依赖管理工具UV的安装和使用教程》UV是一个用Rust编写的Python包安装和依赖管理工具,比传统工具(如pip)有着更快、更高效的体验,:本文主要介绍python依赖管理工具UV... 目录前言一、命令安装uv二、手动编译安装2.1在archlinux安装uv的依赖工具2.2从github

Python屏幕抓取和录制的详细代码示例

《Python屏幕抓取和录制的详细代码示例》随着现代计算机性能的提高和网络速度的加快,越来越多的用户需要对他们的屏幕进行录制,:本文主要介绍Python屏幕抓取和录制的相关资料,需要的朋友可以参考... 目录一、常用 python 屏幕抓取库二、pyautogui 截屏示例三、mss 高性能截图四、Pill

C#实现SHP文件读取与地图显示的完整教程

《C#实现SHP文件读取与地图显示的完整教程》在地理信息系统(GIS)开发中,SHP文件是一种常见的矢量数据格式,本文将详细介绍如何使用C#读取SHP文件并实现地图显示功能,包括坐标转换、图形渲染、平... 目录概述功能特点核心代码解析1. 文件读取与初始化2. 坐标转换3. 图形绘制4. 地图交互功能缩放

java时区时间转为UTC的代码示例和详细解释

《java时区时间转为UTC的代码示例和详细解释》作为一名经验丰富的开发者,我经常被问到如何将Java中的时间转换为UTC时间,:本文主要介绍java时区时间转为UTC的代码示例和详细解释,文中通... 目录前言步骤一:导入必要的Java包步骤二:获取指定时区的时间步骤三:将指定时区的时间转换为UTC时间步

MySQL批量替换数据库字符集的实用方法(附详细代码)

《MySQL批量替换数据库字符集的实用方法(附详细代码)》当需要修改数据库编码和字符集时,通常需要对其下属的所有表及表中所有字段进行修改,下面:本文主要介绍MySQL批量替换数据库字符集的实用方法... 目录前言为什么要批量修改字符集?整体脚本脚本逻辑解析1. 设置目标参数2. 生成修改表默认字符集的语句3

MySQL使用EXISTS检查记录是否存在的详细过程

《MySQL使用EXISTS检查记录是否存在的详细过程》EXISTS是SQL中用于检查子查询是否返回至少一条记录的运算符,它通常用于测试是否存在满足特定条件的记录,从而在主查询中进行相应操作,本文给大... 目录基本语法示例数据库和表结构1. 使用 EXISTS 在 SELECT 语句中2. 使用 EXIS