GEE 依照范围裁剪 下载Sentinel-2数据


0. GEE介绍

Google Earth Engine(GEE) 是由Google开发的一种云端平台,旨在提供强大的地理空间数据处理和分析工具。GEE集成了大量的遥感影像数据和地理空间数据集,以及高性能的计算资源,使用户能够在云端高效地进行大规模地理空间数据处理、分析和可视化,为环境监测、资源管理、自然灾害监测等领域提供了强大的支持。


使用Google Earth Engine(GEE)进行遥感影像分析是一种强大而灵活的方法,它允许用户在云端高效处理大规模的地理空间数据。在Google Earth Engine 中,可以使用 ee.ImageCollection.clip() 等方法来依照范围裁剪 Sentinel-2 数据。


步骤 1:访问Google Earth Engine

  • 打开浏览器,转到 Google Earth Engine 网站。
    • https://code.earthengine.google.com/

  • 如果你尚未登录,使用Google账号登录。


步骤 2:启动Google Earth Engine Code Editor

  • 进入GEE的代码编辑器,按照要求创建项目。

  • 完成Earth引擎代码编辑器


步骤 3:编写和运行代码

在代码编辑器中,可以编写JavaScript代码来处理遥感影像数据。GEE提供了大量的API和函数来操作和分析遥感数据。

  • 通过点击 Assets -> NEW -> Shape files 创建boundary.shp文件,作为文件裁剪范围。
    • 通过点击 SELECT 筛选文件类型,上传边界数据。(上传文件类型见截图)
    • 注意:上传可能有延迟,试着多刷新,稍等一段时间


步骤 4:查看边界数据

  • 在代码编辑器中核验 Cuba 国家行政区划边界数据集boundary.shp
// 从用户的数据集中加载行政区划边界数据集
var district = ee.FeatureCollection("users/xxxx/boundary");// 获取行政区划边界数据集的大小(即要素数量)并打印输出
var dsize = district.size();
print(dsize);// 提取行政区划边界数据集的几何信息
var district_geometry = district.geometry();// 将地图视图中心设置为行政区划边界的几何中心,并缩放到指定级别
Map.centerObject(district_geometry, 7);// 在地图上添加行政区划边界图层
Map.addLayer(district);


步骤 5:下载Sentinel-2影像并导出数据

  • 在代码编辑器中的代码执行过程中,可以在右侧的 Console 面板中查看代码运行的输出结果,包括图像、表格等。
  • 完成了分析,你可以将结果保存在你的 GEE 账户中,以便日后访问和使用。
  • 还可以通过导出功能将结果数据导出至 Google Drive :GeoTIFF、CSV等常见格式,以便在其他GIS软件中使用。
// 定义函数:用于遮蔽Sentinel-2影像中的云
function maskS2clouds(image) {// 选择Sentinel-2影像中的QA60波段,该波段用于云掩蔽var qa = image.select('QA60');// 云和卷云的标志位分别在第10位和第11位var cloudBitMask = 1 << 10;var cirrusBitMask = 1 << 11;// 如果两个标志位都为零,则表示天气晴朗var mask = qa.bitwiseAnd(cloudBitMask).eq(0).and(qa.bitwiseAnd(cirrusBitMask).eq(0));// 更新影像的遮罩并进行归一化return image.updateMask(mask).divide(10000);
}// 加载boundary.shp文件
var boundary = ee.FeatureCollection('users/xxxx/boundary');// 映射函数至一年的数据。
// 加载Sentinel-2 TOA反射数据。
var dataset = ee.ImageCollection('COPERNICUS/S2').filterBounds(boundary)  // 根据区域几何范围筛选数据.filterDate('2020-01-01', '2021-01-01')  // 筛选日期范围// 预先过滤以获取云量较少的数据.filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20)).map(maskS2clouds);  // 应用云掩蔽函数// 定义RGB可视化参数
var rgbVis = {min: 0.0,max: 0.3,bands: ['B4', 'B3', 'B2'],
};// 将处理后的影像数据按照boundary.shp文件裁剪,并添加到地图中
var clippedImage = dataset.median().clip(boundary);
Map.addLayer(clippedImage, rgbVis, 'Clipped Image');// 导出裁剪后的影像数据到Google Drive
Export.image.toDrive({image: clippedImage.select(['B4', 'B3', 'B2']),  // 选择RGB波段description: 'Sentinel2_Cropped',  // 导出影像的描述folder: 'GEE_Images',  // 导出到Google Drive中的文件夹名称region: boundary.geometry(),  // 导出影像的区域scale: 10,  // 分辨率maxPixels: 1e13  // 最大像素数量
});


补充材料

  • GEE的学习曲线可能相对陡峭,请参考官方操作文档和示例代码:
    • https://developers.google.com/earth-engine/guides/getstarted
    • https://earthengine.google.com/platform/
  • GEE提供了强大的云端计算资源,但某些操作可能需要一定的时间才能完成,尤其是对大规模数据的处理。
    • 可以考虑基于小行政区划shp文件,进行分区下载。

如果这对您有所帮助,希望点赞支持一下作者! 😊

点击查看原文
https://mp.weixin.qq.com/s?__biz=Mzk0MTU1MjU5Mw==&mid=2247484716&idx=1&sn=1968d0c1bfacb28e11c73e655f7f137b&chksm=c2d1e392f5a66a849a615c94c3a9f63c9ed9d32e06bef1b3d74df2923f9ff10926bb985df2b9#rd

file

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

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

相关文章

GDB调试入门笔记

文章目录 What&#xff1f;WhyHow安装GDB安装命令查看是否安装成功调试简单的程序预备一个程序调试 使用breakinfolistnextprintstep一些小技巧在gdb前shell日志功能watch point| catch point 调试core调试一个运行的程序 What&#xff1f; GDB是什么&#xff1f; 全称GNU sym…

1、MQ_介绍、优缺点、类型等

MQ介绍 1. MQ概述 MQ&#xff08;Message Queue&#xff09;&#xff1a;消息队列&#xff0c;是基础数据结构中FIFO&#xff08;first in first out&#xff09;的一种数据结构。一般用来解决流量削峰、应用解耦、异步处理等问题&#xff0c;实现高性能&#xff0c;高可用&a…

鸿蒙实战开发:数据交互【RPC连接】

概述 本示例展示了同一设备中前后台的数据交互&#xff0c;用户前台选择相应的商品与数目&#xff0c;后台计算出结果&#xff0c;回传给前台展示。 样例展示 基础信息 RPC连接 介绍 本示例使用[ohos.rpc]相关接口&#xff0c;实现了一个前台选择商品和数目&#xff0c;后台…

【大数据】-- 创建 Paimon 外部表

如今&#xff0c;在数据湖三剑客&#xff08;delta lake、hudi、iceberg&#xff09;之上&#xff0c;又新出一派&#xff1a; apache paimon。我们恰好在工作中遇到&#xff0c;以下介绍在 dataworks 上&#xff0c;使用 maxcompute odps sql 创建 apache paimon 外部表的一些…

Claude3深夜震撼发布!模型特点分析,附使用教程

Claude3深夜震撼发布&#xff01;模型特点分析&#xff0c;附使用教程 引言 最新发布的Claude3引起了广泛关注&#xff0c;这次发布一举推出了三个不同类型的模型&#xff0c;分别是Claude 3 Haiku、Claude 3 Sonnet和Claude 3 Opus。每个模型都具有独特的特点和能力&#xff…

深色系可视化界面看腻了,来点浅色系?安排,20页来了。

只要不放在大屏上展示&#xff0c;贝格前端工场还是非常推崇浅色系的可视化界面&#xff0c;把它作为配色的首选 。浅色系可视化界面具有以下几个优势&#xff1a; 清晰明了 浅色系界面通常使用明亮的颜色&#xff0c;如白色、浅灰色等&#xff0c;使界面元素更加清晰可见。这…

Python 开发图形界面程序

用 Python 语言开发图形界面的程序&#xff0c;有2种选择&#xff1a; Tkinter 基于Tk的Python库&#xff0c;这是Python官方采用的标准库&#xff0c;优点是作为Python标准库、稳定、发布程序较小&#xff0c;缺点是控件相对较少。 PySide2/PySide6 基于Qt 的Python库&#x…

机器学习-面经(part7、无监督学习)

机器学习面经系列的其他部分如下所示&#xff1a; 机器学习-面经&#xff08;part1&#xff09; 机器学习-面经(part2)-交叉验证、超参数优化、评价指标等内容 机器学习-面经(part3)-正则化、特征工程面试问题与解答合集机器学习-面经(part4)-决策树共5000字的面试问题与解答…

【ArcGIS超级工具】基于ArcPy的矢量数据批量自动化入库工具

最近&#xff0c;有很多做规划的朋友私信我&#xff0c;想让我帮忙开发一款ArcGIS自动化脚本工具&#xff0c;实现点、线、面的自动化入库操作&#xff0c;帮他们在平时的内业数据处理工作中减少机械式重复性的工作&#xff0c;提高工作效率。为此&#xff0c;我详细了解了下目…

这本书太好了!150页就能让你上手大模型应用开发

如果问个问题&#xff1a;有哪些产品曾经创造了伟大的奇迹&#xff1f;ChatGPT 应该会当之无愧入选。仅仅发布 5 天&#xff0c;ChatGPT 就吸引了 100 万用户——当然&#xff0c;数据不是关键&#xff0c;关键是其背后的技术开启了新的 AI 狂潮&#xff0c;成为技术变革的点火…

数据结构与算法学习【算法思想之二分法基础】

文章目录 数据结构与算法学习【算法思想之二分查找基础】本文学习目标或巩固的知识点 最基础的二分查找&#x1f7e2;通过题目可知题解结果验证 数据结构与算法学习【算法思想之二分查找基础】 本文学习目标或巩固的知识点 学习二分法类题目 巩固基础的二分法 提前说明&#…

Jmeter之Ramp-up Period(in seconds)

1、Ramp-up Period概念 &#xff08;in seconds&#xff09;–并发用户启动周期&#xff0c;告知JMeter 要在多长时间内启动全部Vuser用户。 2、为什么需要有“ramp-up period”&#xff0c;立即启动所有的并发用户数不是更好&#xff1f; 对于绝大多数的网址或应用&#xf…

【数据结构】堆的TopK问题

大家好&#xff0c;我是苏貝&#xff0c;本篇博客带大家了解堆的TopK问题&#xff0c;如果你觉得我写的还不错的话&#xff0c;可以给我一个赞&#x1f44d;吗&#xff0c;感谢❤️ 目录 一. 前言二. TopK三. 代码 一. 前言 TOP-K问题&#xff1a;即求数据结合中前K个最大的元…

【Docker】技术架构演变

【Docker】技术架构演变 目录 【Docker】技术架构演变架构中的概念架构演进单机架构相关软件 应用数据分离架构应用服务集群架构相关软件 读写分离/主从分离架构相关软件 引入缓存——冷热分离架构相关软件 垂直分库&#xff08;分布式数据库架构&#xff09;相关软件 业务拆分…

Day14:单元测试、Junit单元测试框架、反射、注解

单元测试 针对最小的功能单元&#xff08;方法&#xff09;进行正确性测试编写正规的单元测试框架传统的无法执行自动化测试&#xff0c;且无法得到测试报告 Junit单元测试框架 Junit的作用&#xff1a; 具体步骤 测试类取名&#xff1a;原类名Test&#xff08;大驼峰&#…

蓝桥杯备战刷题four(自用)

1.砝码称重 #include <iostream> #include <vector> using namespace std; const int N110; const int M100010; int w[N]; int n; int f[N][M]; int m; int ans; //f[i][j]表示到第i个砝码进行放置时的称得的重量为j的方案数 int main() {cin>>n;for(int i1…

DbSchema导出HTML/PDF版表结构

一、连接数据库 登录成功默认显示当前用户的所有资源&#xff08;表、视图、序列、方法、触发器等&#xff09;&#xff0c;如果不操作将导出此用户的全部信息。 至此连接数据库完成 二、表结构导出 本次不想给用户全部导出&#xff0c;只给导出几张&#xff0c;选择需要…

分享7款前端动画特效(附效果图及在线演示)

分享7款好玩的前端动画特效 其中有CSS动画、SVG动画、js小游戏等等 下方效果图可能不是特别的生动 那么你可以点击在线预览进行查看相应的动画特效 同时也是可以下载该资源的 canvas彩色画树特效 基于canvas实现的画树特效 同时还可选择树枝的初始数目进行彩色树生成 以下效果…

企业对接Walmart平台API流程 On-request Reports API(二)

对接On-request Reports API 1、对接指南1.1 报告生成时间1.2 报告保留期1.3 请求限制1.4 报告请求工作流如何申请报告第 1 步&#xff1a;申请取消报告第 2 步&#xff1a;获取报表可用性状态第 3 步&#xff1a;下载报告 URL 2、代码实现2.1、获取访问API的token2.2、构建公共…

【教育部白名单赛事】C语言编程题解析--软件编程邀请赛(决赛)

文章目录 1、保留12位小数的浮点数2、气温统计3.大写字母的判断4、【递归】母鸡的故事5、小白免再排队 1、保留12位小数的浮点数 输入一个双精度浮点数&#xff0c;保留12位小数&#xff0c;输出这个浮点数。 时间限制&#xff1a;1000 内存限制&#xff1a;65536 【输入】 只…