// Definir región de interés var countries = ee.Filter.expression('ADM2_CODE == 31823') var paises = ee.FeatureCollection("FAO/GAUL/2015/level2").filter(countries) var roi = paises.geometry() //Seleccionar Colección de imágenes utilizada como base para el remuestreo var dataset = ee.ImageCollection("MODIS/061/MOD13Q1").filter(ee.Filter.date('2002-05-01','2023-05-25')).filterBounds(roi); //Definir máscara agrícola var v200 = ee.ImageCollection("ESA/WorldCover/v200").filterBounds(roi).median() var AgMask = v200.eq(40).clip(roi) // Cambiar proyección de la máscara, reducir su resolución espacial y remuestrear var crsDS = dataset.first().projection() var AgMask_250 = AgMask.setDefaultProjection({crs :crsDS, scale: 10}).reduceResolution({ reducer: ee.Reducer.mean(), maxPixels: 625 }) .reproject({ crs: crsDS }).gte(0.5); // Definir parametros visualizacion var empty = ee.Image().byte(); var outline = empty.paint({ featureCollection: paises, color: 1, width: 1 }); //Agregar capas al mapa Map.setOptions("SATELLITE") Map.centerObject(roi, 12); Map.addLayer(dataset.first().select('NDVI').clip(roi), {min: 5000, max: 8000}, "NDVI", false); Map.addLayer(AgMask, {}, "AgMask_10m", false); Map.addLayer(AgMask_250, {}, "AgMask_250m", false); Map.addLayer(outline, {palette: 'ff0000'}, 'Deptos'); //Map.addLayer(paises, {}, 'paises') //Exportar la máscara crsDS.evaluate(function(proj_obj) { Export.image.toAsset({ image: AgMask_250, description: 'AgMask_250', assetId: 'AgMask_250', // modificar ID region: roi, //scale: 250, // No definir la escala porque reproyecta la imagen y la modifica crs: proj_obj.crs, crsTransform: proj_obj.transform }); });
Preview:
downloadDownload PNG
downloadDownload JPEG
downloadDownload SVG
Tip: You can change the style, width & colours of the snippet with the inspect tool before clicking Download!
Click to optimize width for Twitter