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

GEE入门教程:详细步骤教你下载影像数据

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

GEE入门教程:详细步骤教你下载影像数据

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

一、基本流程

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

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

C. 数据信息查看:查看所选数据集的详细描述、波段信息、影像属性等,并了解如何通过代码访问这些数据。

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


图1:GEE平台界面示例


图2:数据检索界面示例

二、下载流程

1. 导入研究区域

首先需要在ArcMap中导出研究区域的矢量数据(如县级行政区划数据),然后将其上传到GEE的个人资产中。使用ee.FeatureCollection代码导入研究区域。

var roi = ee.FeatureCollection("users/exampleuser/researcharea");

2. 数据调用和预处理

通过ee.ImageCollection代码访问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);

3. 查看影像

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

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");

4. 合成影像和下载

由于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: "ResearchAreaESA",  
      scale: 10,  
      region: roi,  
      maxPixels: 1e13,  
      crs: "EPSG:4326"  
  });
// 其他分辨率的影像下载代码省略,结构类似

5. 验证下载结果

将下载的影像数据在ENVI等软件中加载,检查影像质量。

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