GEE批量下载Sentinel2数据

浩瀚地学

已于 2023-11-03 00:08:47 修改

阅读量3.9k

点赞数 56

分类专栏： Google earth engine 数据下载文章标签： 1024程序员节经验分享学习笔记 arcgis js sentinel

于 2023-10-24 21:57:04 首次发布

本文链接：https://blog.csdn.net/zbh13859825167/article/details/134022781

版权

数据下载同时被 2 个专栏收录

6 篇文章 4 订阅

订阅专栏

Google earth engine

1 篇文章 0 订阅

订阅专栏

GEE批量下载数据

Google Earth Engine https://earthengine.google.com网站批量下载数据，需要有谷歌账号和科学上网

界面功能

导入研究区边界

1.学会导入roi的矢量文件，在【Assets】下选择【New】，导入shp文件

2.在弹出的界面可以选择自己的文件进行上传，除了在gee上的新建的数据集文件名需要修改之外，其他保持默认即可，最后点击【Upload】上传

3.创建好之后就可以去copy别人的代码进行批量下载了，有两个数据可以需要自己手动导入，网上的教程有时没讲这部分，通过【Import】导入roi文件和数据集。其中，roi就是刚刚的矢量文件，在旁边的目录中有，点击之后在弹出的界面找到【Import】就可以导入了；数据集在上面的搜索栏中可以找到，比如下载Sentinel-2 L2A数据，就可以直接搜索“Sentinel”，点击搜索，之后找到对应的数据集，点击之后找到【Import】就可以导入

4.在导入这两个之后，就会再代码部分最上面出现如下三行（这三行不算代码，代码行数的计数从这三行之后开始），点击可以将”table”和“imageCollection”直接改名，因此在网上的一些代码中出现了未定义变量的报错，可能就是这两个错误；也可以通过在之后通过定义新变量来指代这两个，如下所示，这样名称就不用改了

var roi = table;
var s2_collection = imageCollection;

数据下载

以哨兵数据批量下载为例（无云掩膜）

var roi = table;

function filterImageCollection(collection) {
  // 日期 
  var startDate = '2021-01-01';
  var endDate = '2021-06-30';
  var filteredCollection = collection.filterDate(startDate, endDate);
  filteredCollection = filteredCollection.filterBounds(roi);
  // 云量
  filteredCollection = filteredCollection.filter(ee.Filter.rangeContains("CLOUDY_PIXEL_PERCENTAGE", 0, 30));
  return filteredCollection;
}
  var s2_collection = ee.ImageCollection('COPERNICUS/S2_SR');
  var filteredCollection = filterImageCollection(s2_collection);
  // 获取影像列表
  var imageList = filteredCollection.toList(filteredCollection.size());
  // 获取影像数量
  var count = imageList.size().getInfo();
  // 使用 for 循环进行迭代
for (var i = 0; i < count; i++) {
  // 通过索引获取影像
  var image = ee.Image(imageList.get(i));
  // 获取原始影像的文件名
  var originalFilename = ee.String(image.get('system:index')).getInfo();
  // 对影像进行研究区域裁剪和去云处理
  var maskedImage = image.clip(roi);
  // 选择需要的波段
  var composite = maskedImage.select(['B2', 'B3', 'B4', 'B8', 'B11']);
  // 导出波段合成影像
  Export.image.toDrive({
    image: composite,
    description: 'S2L2A_' + originalFilename,
    folder: 'S2L2A',
    region: roi,
    scale: 10,
    crs: "EPSG:4326",
    maxPixels: 1e13
  });
}
print('OK');