任何一个遥感人不用GEE下载遥感影像我都会伤心的ok?!! 众所周知哨兵二号影像可以在usgs上或者欧空局网站下载,usgs感觉还相对更方便些。但是都还是避免不了下载影像后先用SNAP预处理大气校正重采样再波段组合什么的才能得到L2A级的tif文件。 现在有一个很方便的方法摆在我们面前,只需要十几行代码分分钟下载到预处理后的tif影像,还可以一步解决裁剪、投影、重采样的操作。
GEE界面: 总之就是先导入你研究区的shapefile和哨兵二号L2A,运行下面代码就好了(注释写得很详细了!!)。
var batch = require('users/fitoprincipe/geetools:batch')
var imageCollection = s2_col.filterDate('2020-2-10', '2020-6-01').filterBounds(roi).filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 10));
var imageList = imageCollection.toList(imageCollection.size())
var imageListSize = imageList.size().getInfo()
print(imageList)
print(imageListSize)
var visParams = {bands: ['B4', 'B3', 'B2'],min: 0,max: 3000,gamma: 1.4,};
for(var i=0;i<imageListSize;i++){
var image = ee.Image(imageList.get(i))
Map.addLayer(image, visParams, image.get("system:index").getInfo());
}
Map.centerObject(roi,11);
Map.addLayer(roi,{color:'FFFF00'});
batch.Download.ImageCollection.toDrive(imageCollection,"result", {
scale: 10,
crs: "EPSG:4526",
region: roi,
maxPixels:1e13,
type:"int16" })
打印的结果是这样的: 下载到google drive里是这样的: 把下载的影像在arcmap里显示是这样的: 在arcmap里怎么设置符号系统展示成如上图所示的样子就不说了,之前在一篇笔记里写过方法。
再次感叹!!GEE真是21世纪对遥感人(也可能仅对本人。。)来说最伟大的发明!!!!!!
|