之前我以为这个问题是因为研究区很大,会导致影像的自动分片分区,当下载的单个文件超过一定的大小时候,就会产生如下图这样的分区文件并且自动编号,但是这如果研究区小,出现如此多的文件,主要的问题就是没有将影像在一起用mosaic在一起。
原有的下载数据用到的是max(),所以数据就是提取像元的最大值,所以下载的数据再arcgis中就是最大的值。
最后修改后:
代码:
var table = ee.FeatureCollection("users/bqt2000204051/beijing"); var roi = table; //landsat 8 function NDVI_landsat_8(image) { var ndvi = image.normalizedDifference(['B5', 'B4']); return image.addBands(ndvi.rename('NDVI')); } var imageCollection = ee.ImageCollection("LANDSAT/LC08/C01/T1_SR").select(['B5', 'B4']) for(var i=2019;i<=2020;i++){ for(var j=1;j<=12;j++){ var data_collection = imageCollection.filter(ee.Filter.calendarRange(j, j,'month')) .filter(ee.Filter.calendarRange(i, i, 'year')) .map(function(img) { return img.set('year', img.date().get('year')); }).map(NDVI_landsat_8).filterBounds(roi); var HY_collection = data_collection.qualityMosaic('NDVI').clip(roi); print(data_collection,i+"_year_"+j+"_month"); Map.addLayer(HY_collection,{},i+'-'+j); Export.image.toDrive({ image: HY_collection.int8(), description: i+'-'+j, crs: "EPSG:32649", scale: 30, region: roi, maxPixels: 1e13, folder: 'ShanXi' }) } } Map.addLayer(roi,{},'hengyang Boundary');