GEE学习记录（三）：Google Earth Engine处理MODIS产品——以MOD10A1积雪产品为例

前言

工具在不断更新，使用完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中打开看看结果，到此就结束了。

结束语

这不仅仅是虚荣，在我看来，我的面孔是我的存在的意义的全部体现；它是我祖先的记录，世界的一部分；同时，也是我接受这个世界，拼命依附这个世界的方式；也是我向这个世界表白自己的方式。这一切都是只可意会，不可言传的。
——索尔·贝娄 《晃来晃去的人》