前言
工具在不断更新,使用完GEE处理MOD10A1数据后,我只想说MRT这次真的可以说“Goodbye”了,这可能是我未来为数不多用这么好的工具了,已经转战新的研究方向了,最近会把之前的一些工作作一些整理放出来了。
那么,GEE的好处就不再详细赘述了,反正用过的人都说“真香”,今天就放出代码大家自己看看吧!
研究区域矢量数据载入
以下均以北美洲为例进行讲解,在Assets下点击new,选择shape files,再点击SELECT,选好要用的shapefile文件,注意时所有的不仅是.shp文件,然后选择一个在GEE中存储的位置。滑到底部选择upload,此时矢量文件就在里面了。最后,点击这个文件,在弹出的窗口中选择import即可在旁边看到这个矢量文件了,到此为止导入矢量文件的工作就结束了。
代码块
接下来就是生成研究区的MOD10A1数据,代码块我放在下面了,这里生成两天的,如果你的需求很多,修改时间RUN就OK了,重要的地方我已经做了注释,如有问题请联系我!
var mod10a1NDSI = ee.ImageCollection("MODIS/006/MOD10A1")
.filterBounds(table7)//适量对应的编号
.filterDate("2019-12-01", "2019-12-03")//time往后写一天
.select("NDSI_Snow_Cover_Class");//需要的产品
print('Collection: ', mod10a1NDSI);
var num = 2;//景数
var list = mod10a1NDSI.toList(num);
print(list);
var count = list.size().getInfo();
print(count)
for(var i=0;i<count;i++)
{
// get the single image from list
var image = ee.Image(list.get(i));
print(i,image);
// get the id property as the part of filenames
var id=image.id();
print(id);
// out name
var name2="NDSI_"+id.getInfo();
print(name2)
Export.image.toDrive({
image: image,
description: 'imageToCOGeoTiffExample',
fileNamePrefix: name2,
scale: 4000,//分辨率
maxPixels: 9999999999,
//crs:'EPSG:32649',
region: table7,
fileFormat: 'GeoTIFF',
formatOptions: {
cloudOptimized: true
}
});
}
结果下载
视窗右侧控制台会显示景数和运行结果,在Task下进行下载,下载结果可存在google drive中,再从google drive中下载到本地即可,因为我尝试直接下载到本地没有成功,你们可以自行尝试,在ENVI中打开看看结果,到此就结束了。
结束语
这不仅仅是虚荣,在我看来,我的面孔是我的存在的意义的全部体现;它是我祖先的记录,世界的一部分;同时,也是我接受这个世界,拼命依附这个世界的方式;也是我向这个世界表白自己的方式。这一切都是只可意会,不可言传的。
——索尔·贝娄 《晃来晃去的人》