使用GEE或PIE下载Landsat8影像

最新推荐文章于 2024-11-04 16:13:10 发布

陨星落云

最新推荐文章于 2024-11-04 16:13:10 发布

阅读量1.1w

点赞数 20

分类专栏： GEE与PIE 文章标签：深度学习人工智能

本文链接：https://blog.csdn.net/qq_28368377/article/details/121929529

版权

GEE与PIE 专栏收录该内容

7 篇文章

订阅专栏

本文对比了使用Google Earth Engine (GEE)与Pie-Engine (PIE)下载Landsat8影像的过程，包括数据选择、区域裁剪、影像显示及导出步骤，帮助读者理解两种平台的异同和操作方法。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

使用遥感云平台下载Landsat8影像

一、GEE下载Landsat8影像

GEE官网：

https://code.earthengine.google.com/

上传下载区域的矢量边界
GEE代码：

//首先，导入自己的研究区矢量数据
var ROI = ee.FeatureCollection("users/yunxingluoyun/Heze");
var ROI = ROI.geometry();

//然后，选择数据集和合成数据的时间范围
var dataset = ee.ImageCollection('LANDSAT/LC08/C02/T1')
    .filterDate('2021-11-01', '2021-12-01').filter(ee.Filter.lt("CLOUD_COVER",5));


//用clip（）函数里面的行政边界切割
var ROI_collection = ee.Algorithms.Landsat.simpleComposite(dataset).clip(ROI)

print(ROI_collection)
Map.centerObject(ROI, 7);
//显示菏泽的shp边界
Map.addLayer(ROI,{color:'yellow',fillColor: "00000000", width: 1},'Heze Boundary');

//显示ROI影像并且以321合成，后面命名
Map.addLayer(ROI_collection,{bands:["B4","B3","B2"],gamma:1.3,max:108,min:15},'Heze Image');

//下载处理好的影像
Export.image.toDrive({
  image:ROI_collection,
  description: 'Landsat8_2021_11',
  fileNamePrefix:'Landsat8_heze',
  region:ROI,
  scale:30,
  crs: "EPSG:4326",
  maxPixels:1e13
})

结果：

二、PIE下载Landsat8影像

PIE-Engine官网：

https://engine.piesat.cn/workbench/sign-in

例：下载济南地区Landsat8影像

PIE代码如下：

//确定区域
//首先，导入自己的研究区矢量数据
var roi = pie.FeatureCollection("NGCC/CHINA_CITY_BOUNDARY")
    .filter(pie.Filter.eq("name","济南市"))
    .first()
    .geometry();


//pie然后，选择数据集和合成数据的时间范围
//加载Landsat8 Collection2 SR数据集并日期、云量筛选
var img = pie.ImageCollection('LC08/02/SR')
            .filterBounds(roi)
            .filterDate("2021-1-1", "2021-12-5")
            .filter(pie.Filter.lt("cloud_cover", 2))
            .sort()
            .map(function(image){
                var qa = image.select('QA_PIXEL');
                var cloudMask = qa.bitwiseAnd(1<<3).eq(0);
                return image.updateMask(cloudMask);})
            .select(["B2", "B3", "B4"])
            .mosaic()
            .clip(roi);
print(img);

//按影像定位地图中心
Map.centerObject(roi, 7);
//显示HY的shp边界，后面‘’里面命名
Map.addLayer(roi,{color:'yellow', fillColor: "00000000", width: 1},'济南边界');
//加载显示影像
Map.addLayer(img, {min: 6000, max: 20000, bands: ["B4", "B3", "B2"]}, "img");

//导出影像到云盘
Export.imageToCloud({
  image:img,
  description: 'Landsat8',
  assetId:'Landsat8_Jinan1',
  region:roi,
  scale:30,
  crs: "EPSG:4326",
  maxPixels:1e13
})

/**
 * 导出影像，参数分别为：
 *  image：导出影像
 *  description：导出任务描述也就是任务名称
 *  assetId：导出的《资源》中数据集合路径
 *  region：导出范围
 *  scale：分辨率30米
 * */
/*
//导出到资源下
Export.image({
    image: img,
    description: "landsat8",
    assetId: "Landsat8_Jinan",
    region: roi,
    scale: 30
});
*/

结果：