问小白 wenxiaobai
资讯
历史
科技
环境与自然
成长
游戏
财经
文学与艺术
美食
健康
家居
文化
情感
汽车
三农
军事
旅行
运动
教育
生活
星座命理

GEE入门教程:遥感影像数据下载完整指南

创作时间:
作者:
@小白创作中心

GEE入门教程:遥感影像数据下载完整指南

引用
CSDN
1.
https://blog.csdn.net/Promising_GEO/article/details/137526315

本文将详细介绍如何使用Google Earth Engine (GEE)平台下载遥感影像数据。通过本文,你将学会从数据检索、代码编写到最终下载的完整流程。

一、基本流程

A. 平台进入:访问GEE官网(https://developers.google.com/earth-engine),使用邮箱注册并登录。

B. 数据检索:在平台上点击“Dataset”,可以查看所有可用的数据集。选择“Landsat”数据,进一步选择“Landsat9OLI/TIRS”数据集,以及大气层顶反射率数据集“Top of Atmosphere”。

C. 数据信息查看:查看所选数据集的详细描述、波段信息、影像属性等,并了解平台上的代码调用方式。

D. 数据下载:GEE支持通过JavaScript和Python两种接口进行代码编辑。下载流程通常包括:导入研究区域→访问影像数据集→空间、时间和云量筛选→查看影像结果→代码调整→影像下载。


图1:GEE平台界面示例


图2:数据检索界面示例

二、下载流程

① 导入研究区域

在ArcMap中导出县级行政区划数据,并在个人资产上上传研究区域的矢量shp文件。使用以下代码导入研究区域:

var roi = ee.FeatureCollection("users/hesuixinya511/rongchang");

② 数据调用和预处理

通过GEE平台访问Sentinel-2数据,根据研究区域边界和时间范围筛选数据,并进行去云处理:

function maskS2clouds(image) {
  var qa = image.select('QA60');
  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);
}

var dataset = ee.ImageCollection('COPERNICUS/S2_HARMONIZED')
                  .filterDate('2020-01-01', '2020-12-31')
                  .filterBounds(roi)
                  .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 0.05))
                  .map(maskS2clouds);
print(dataset);

③ 查看影像

生成影像列表并可视化显示,选择合适的影像:

var list = dataset.toList(dataset.size());
var downImage = ee.Image(list.get(3));
var upImage = ee.Image(list.get(5));
var visualization = {
  min: 0.0,
  max: 0.3,
  bands: ['B4', 'B3', 'B2'],
};
Map.centerObject(roi,6);
Map.addLayer(upImage,visualization,"upImage");
Map.addLayer(downImage, visualization, 'downImage');
Map.addLayer(roi,{"color":"red"},"ROI");

④ 合成影像和下载

由于Sentinel-2数据的波段具有不同分辨率,需要先进行影像合成,然后设置下载参数:

var downImage10 = downImage.select("B2","B3","B4","B8");
var downImage20 = downImage.select("B5","B6","B7","B8A","B11","B12");
var downImage60 = downImage.select("B1","B9","B10");
var upImage10 = upImage.select("B2","B3","B4","B8");
var upImage20 = upImage.select("B5","B6","B7","B8A","B11","B12");
var upImage60 = upImage.select("B1","B9","B10");

Export.image.toDrive({  
      image: downImage10,  
      description: "downImage10",  
      fileNamePrefix: "Sentinel-2down10",  
      folder: "RongchangESA",  
      scale: 10,  
      region: roi,  
      maxPixels: 1e13,  
      crs: "EPSG:4326"  
  });

// 其他分辨率影像的下载代码省略,格式类似

⑤ 验证下载结果

将下载的Sentinel-2影像数据和Landsat8数据在ENVI中加载,检查影像质量。

© 2023 北京元石科技有限公司 ◎ 京公网安备 11010802042949号