// 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
});
});
Comments