Frequência anual de cultivos

// ========================================================
// MapBiomas – Frequência Anual de Cultivos | Coleção 10
// Ano desejado (2017 a 2024)
// ========================================================
var ANO = 2023;
// Classes a analisar (uma ou todas)
// 1 = Um cultivo
// 2 = Dois cultivos
// 3 = Tres cultivos
var CLASSES = [1, 2, 3];

// Região de recorte
// null = Brasil inteiro
var REGIAO = null;

// ========================================================
// Asset
// ========================================================
var ASSET =
'projects/mapbiomas-public/assets/brazil/lulc/collection11/mapbiomas_brazil_collection11_agriculture_number_cycles_v1'

// Nome da banda
var NOME_BANDA = 'classification_' + ANO;

// ========================================================
// Carregar imagem
// ========================================================
var imagem = ee.Image(ASSET).select(NOME_BANDA);

// ========================================================
// Filtrar apenas as classes desejadas
// ========================================================
var mascara = ee.Image(0);
CLASSES.forEach(function(c){
  mascara = mascara.or(imagem.eq(c));
});
var imgFiltrada = imagem.updateMask(mascara);

// ========================================================
// Geometria
// ========================================================
var geometria = REGIAO !== null
  ? REGIAO
  : imagem.geometry();

// ========================================================
// Visualização
// ========================================================
var PALETA = {
  1: '#c7eef2',   // 1 ciclo cultivo
  2: '#fe668c',   // 2 ciclos cultivos
  3: '#8f0620'    // 3 ciclos cultivos
};

var cores = CLASSES.map(function(c){
  return PALETA[c];
});

Map.setCenter(-52, -14, 5);
Map.setOptions('HYBRID');

Map.addLayer(
  imgFiltrada,
  {
    min: Math.min.apply(null, CLASSES),
    max: Math.max.apply(null, CLASSES),
    palette: cores
  },
  'Frequência Anual de Cultivos ' + ANO
);

// ========================================================
// Área por classe (ha)
// ========================================================
var FATOR_HA = 900 / 10000;
var hist = imagem.reduceRegion({
  reducer: ee.Reducer.frequencyHistogram(),
  geometry: geometria,
  scale: 30,
  maxPixels: 1e13,
  bestEffort: true
});

var dictHist = ee.Dictionary(hist.get(NOME_BANDA));
var NOMES = {
  1: 'Um_Cultivo',
  2: 'Dois_Cultivos',
  3: 'Tres_Cultivos'
};

var features = CLASSES.map(function(c){
  var area_ha = ee.Number(
    dictHist.get(String(c), 0)
  ).multiply(FATOR_HA);
  return ee.Feature(null, {
    classe: c,
    frequencia: NOMES[c],
    ano: ANO,
    area_ha: area_ha
  });
});

var tabelaEstat = ee.FeatureCollection(features);
print('Área por frequência de cultivo (ha) - ' + ANO,
      tabelaEstat);