SDC_PBS_Amazon.js

/* Users with access to the repository can add it to the Code Editor using:
https://code.earthengine.google.com/?accept_repo=users/SI_Testing/LULC_EC
*/

/**
 *  * SDC and PBS classification have been modified from Simonetti D, 
 *      https://publications.jrc.ec.europa.eu/repository/handle/JRC99695
 *  * Adding classes such as annual and perennial crops, urban areas.
 *  * Spectral index have been testing by region
 * 
 /**
*/
  //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
//-------------------------------------------------------------------------------------------------------------------------------
//  ----------     Visualization functions                         -----------------------------------------------
//-------------------------------------------------------------------------------------------------------------------------------

function rgb(r,g,b){
          var bin = r << 16 | g << 8 | b;
          return (function(h){
          return new Array(7-h.length).join("0")+h
          })(bin.toString(16).toUpperCase())
}

var CLASS_PBSvis = {
  'bands':'Class',
  'min': 1,
  'max': 10,
  palette: [
            
            rgb(67, 27, 224), //1  Water -----
            rgb(12, 64, 13),    //2  Forest
            rgb(255,153,255),   //3   Annual Crop
            rgb(255,0,127),  //4 Perennial Crop
            rgb(27,229,5),   //5 grassland 
            rgb(255,128,0),  //6 shrub
            rgb(255,213,0),   //7 bare Soil
            rgb(252, 8, 8),   //8 urban  
            rgb(192,192,192), //9 
            rgb(20,20,20)   //     
            
    ] //my palette
};

exports.visPBS = CLASS_PBSvis;

 // ----------------------------------------------------------------------------------------------------------
// ----------------------------------------------------------------------------------------------------------
// ------------------  SINGLE DATE CLASSIFICATION (SDC--------------------------------------------------------
// ----------------------------------------------------------------------------------------------------------
// ----------------------------------------------------------------------------------------------------------

// INPUT 1: Monthly image to be classified 
// INPUT 2: bands conbination (B,G,R,NIR,SWIR1,SWIR2) ['B2','B3','B4','B8','B11','B12']
// OUT :  classified input(same number of layers); Class code is not the same as the PBS  

exports.SDC_dataset =  function (image,BANDS,imageVIIRS){
  
    var th_NDVI_MAX_WATER= 0;
    var BLU=image.select(BANDS[0]);
    var GREEN=image.select(BANDS[1]);
    var RED=image.select(BANDS[2]);
    var NIR=image.select(BANDS[3]);
    var SWIR1=image.select(BANDS[4]);
    var SWIR2=image.select(BANDS[5]);
             
    var OUT=ee.Image(0);
    
    var th_NDVI_SATURATION=0.0040;
    var th_NDVI_MIN_CLOUD_BARE=0.39;//-0.3
    var th_NDVI_MIN_VEGE=0.49;
    
    var th_SHALLOW_WATER=-0.12;
    var th_RANGELAND=0.53;
    var th_GRASS=0.58;
    var th_SHRUB=0.69;
    var th_TREES=0.85 ;
   
    var min123=BLU.min(GREEN).min(RED);
     
    var min1234=min123.min(NIR);
    var min234=GREEN.min(RED).min(NIR);
    
    var max234=GREEN.max(RED).max(NIR);
    var max1234=max234.max(BLU);
    
    var max57=SWIR1.max(SWIR2);
    var max457=max57.max(NIR);
    
    var max123457= max1234.max(max57);
    
    
    var BLUgtGREEN  = BLU.gt(GREEN);
    var BLUgteGREEN = BLU.gte(GREEN);
    var BLUlteNIR   = BLU.lte(NIR);
    
    var GREENgtRED  = GREEN.gt(RED);
    var GREENlteRED = GREEN.lte(RED);
    var GREENgteRED = GREEN.gte(RED);
    var REDlteNIR= RED.lte(NIR);
    
    var REDsubtractGREEN = (RED.subtract(GREEN)).abs();
    var BLUsubtractNIR   = BLU.subtract(NIR)
    
    var BLUgtGREENgtRED=BLUgtGREEN.and(GREENgtRED)
    
    var growing14=(BLU.lte(GREEN)).and(GREENlteRED).and(REDlteNIR);
    var growing15=growing14.and(NIR.lte(SWIR1));
    
    var decreasing2345=(GREENgteRED).and(RED.gte(NIR)).and(NIR.gte(SWIR1));
    
    
    var SATURATION=(max234.subtract(min234)).divide(max234);

    var WETNESS= image.expression('byte(b("'+BANDS[0]+'")*255)*0.2626 + byte(b("'+BANDS[1]+'")*255)*0.21 + byte(b("'+BANDS[2]+'")*255)*0.0926 + byte(b("'+BANDS[3]+'")*255)*0.0656 - byte(b("'+BANDS[4]+'")*255)*0.7629 - byte(b("'+BANDS[5]+'")*255)*0.5388');
    
    var NDVI=(NIR.subtract(RED)).divide(NIR.add(RED));
    var NDSI=(BLU.subtract(SWIR1)).divide(GREEN.add(SWIR1));
    
    var BRIGTSOIL=((BLU.lt(0.168)).and(growing15)).or((BLU.lt(0.168)).and(growing14).and(  ((NIR.subtract(SWIR1)).gt(-0.079)))); 

    var WATERSHAPE= ((BLU.subtract(GREEN)).gt(-0.2)).and(decreasing2345).and(WETNESS.gt(0)); //add other cond
    var OTHERWATERSHAPE= (BLUgteGREEN).and(GREENgteRED).and(NIR.gte(RED)).and(SWIR1.lt(NIR)).and(SWIR2.lte(SWIR1)).and(NIR.lt((RED).multiply(1.3)).and(NIR.lt(0.12)).and(SWIR1.lt(RED)).and(NIR.lte(GREEN)).and(NIR.gt(0.039)).and(WETNESS.gt(0))  ); //add other cond  07/10 (add replaced with and  :) and(NIR.lte(GREEN))
    
    var SNOWSHAPE=(min1234.gt(0.30)).and(NDSI.gt(0.65));
    
    var CLOUDSHAPE = ((SNOWSHAPE.eq(0)).and(BRIGTSOIL.eq(0))).and(
                  ((max123457.gt(0.47)).and(min1234.gt(0.37))).or(
                    
                    ((min123.gt(0.17)).and((SWIR1).gt(min123))).and(
                          ((SATURATION.gte(0.2)).and(SATURATION.lte(0.4)).and(max234.gte(0.35)) ).or ((NDSI.lt(0.65)).and(max1234.gt(0.30)).and( (NIR.divide(RED)).gte(1.3) ).and((NIR.divide(GREEN)).gte(1.3)).and( (NIR.divide(SWIR1)).gte(0.95)  )) 
                                                                   )
                                                                   
                                                                  ) 
                                                              ) 
    
    min123=0
    
    BRIGTSOIL=0
    SATURATION=0
    decreasing2345=0
    var th_NDVI_river = 0.13 //river
    // main groups based on ndvi
    var ndvi_1 = NDVI.lte(th_NDVI_MAX_WATER);
    var ndvi_2 = NDVI.lt(th_NDVI_MIN_VEGE).and(ndvi_1.eq(0)); 
    var ndvi_3 = NDVI.gte(th_NDVI_MIN_VEGE);
    var ndvi_4 = NDVI.lte(th_NDVI_river)
    
    
    //NDBI = (SWIR2 – NIR)/(SWIR2 + NIR)
    var NDBI = (SWIR2.subtract(NIR)).divide(SWIR2.add(NIR));
    var bareThreshold = -0.11;//-0.13
    var waterNdbiThreshold = -0.1;
    var ndbi_th = NDBI.gt(bareThreshold)
    var ndbi_Water_th = NDBI.lte(waterNdbiThreshold);
    var waterNdbiThreshold2 = 0;
    var ndbi_Water_th2 = NDBI.gte(waterNdbiThreshold2);
    var waterNdbiThreshold3 = 0.2;
    var ndbi_Water_th3 = NDBI.lte(waterNdbiThreshold3);
    
    //NDTI
    var NDTI = (SWIR1.subtract(SWIR2)).divide(SWIR1.add(SWIR2));
    var ndtiUrbanThreshold = 0.19;
    var ndti_Urban_th =NDTI.lte(ndtiUrbanThreshold);
    
    
    //UI = (SWIR2 – NIR)/(SWIR2 + NIR)
    var UI = (SWIR2.subtract(NIR)).divide(SWIR2.add(NIR));
    var urbanThreshold = 0; //-0.05;
    var ui_th = UI.lte(urbanThreshold)//UI.gt(urbanThreshold);
    
    ///NDWI
    var NDWI = (GREEN.subtract(NIR)).divide(GREEN.add(NIR));
    var ndwiThreshold = -0.1 ;
    var ndwi_th = NDWI.gte(ndwiThreshold);
    var ndwiThreshold2 = -0.3;//river
    var ndwi_th2 = NDWI.lte(ndwiThreshold2);
    var ndwiSoilsThreshold = -0.16;
    var ndwi_Soils_th = NDWI.lte(ndwiSoilsThreshold);
    var ndwiUrbanThreshold = -0.2
    var ndwi_Urban_th = NDWI.gt(ndwiUrbanThreshold);
    var ndwiUrbanThreshold1 = -0.08
    var ndwi_Urban_th1 = NDWI.lte(ndwiUrbanThreshold1);

      ///---------------------------------------BSI Index------------------------------------------------------//
    var BSI = ((RED.add(SWIR1)).subtract((NIR.add(BLU)))).divide((RED.add(SWIR1)).add((NIR.add(BLU))));
    var bsiUrbanThreshold1 = 0.149;
    var bsi_Urban_th1 = BSI.lte(bsiUrbanThreshold1);
    var bsiUrbanThreshold2 = 0;
    var bsi_Urban_th2 = BSI.gt(bsiUrbanThreshold2);
    var bsiSoilThreshold1 = 0.149;//0.12
    var bsi_Soil_th1 = BSI.gt(bsiSoilThreshold1);
    var bsiWaterThreshold = 0;
    var bsi_Water_Th = BSI.lte(bsiWaterThreshold)
    var bsiWaterThreshold2 = 0.5; //river
    var bsi_Water_Th2 = BSI.lte(bsiWaterThreshold2)
    var bsiWaterThreshold3 = 0.15; //river
    var bsi_Water_Th3 = BSI.gte(bsiWaterThreshold3)

    //GNDVI = (B8-B3/B8+B3)
    var GNDVI = (NIR.subtract(GREEN)).divide(NIR.add(GREEN));
    var greenThreshold = 0.5;
    var gndvi_th = GNDVI.gte(greenThreshold)
    
    //EVI
    var EVI = image.expression(
     '2.5 * ((NIR/10000 - RED/10000) / (NIR/10000 + 6 * RED/10000 - 7.5 * BLUE/10000 + 1))', {
      'RED': image.select('B4'),
      'BLUE': image.select('B2'),
      'NIR': image.select('B8') //calculate NDVI and add an NDVI band
      
    })
    
    
    //SAVI
    var SAVI = image.expression(
        '((NIR-RED) / (RED+NIR+0.5))*(1.5)', {
      'RED': image.select(BANDS[2]),//b4
      'NIR': image.select(BANDS[3]) //calculate SAVI and add an SAVI band b5
        
    })
    
    var SAVI1 = ((NIR.subtract(RED)).divide(RED.add(NIR).add(0.5))).multiply(1.5)
    
    
    
    //-------------------------------------------------------------------------------------------------------------
		//----------------------  SECTION 1 : WATER  ---------------------------------------------------------
		//-------------------------------------------------------------------------------------------------------------
    
    OUT=(ndvi_1.and(SNOWSHAPE)).multiply(3);

    OUT=OUT.where(( (ndvi_1).and(BLU.gt(0.94)).and(GREEN.gt(0.94)).and(RED.gt(0.94)).and(NIR.gt(0.94)) ),1);  // TEST CLOUDS 
                   
		OUT=OUT.where(( (OUT.eq(0)).and(ndvi_4).and(ndbi_Water_th).and(bsi_Water_Th).and(ndwi_th)),2); //water
		
		OUT=OUT.where(( (OUT.eq(0)).and(ndvi_4).and(ndbi_Water_th2).and(ndbi_Water_th3).and(bsi_Water_Th2).and(bsi_Water_Th3).and(ndwi_th2)),2); //river
    
    
    //-------------------------------------------------------------------------------------------------------------
		//---------------------  SECTION 2 : CLOUDS or SOIL  ---------------------------------------------------------
		//------------------------------------------------------------------------------------------------------------
    
     OUT=OUT.where(( (ndvi_2).and(SNOWSHAPE)),3);

     OUT=OUT.where(( (ndvi_2).and(OTHERWATERSHAPE).and(BLU.gt(0.078)).and(max57.lt(0.058))),2 );//4
    
     OUT=OUT.where(( (ndvi_2).and(
                      CLOUDSHAPE.or(
                      (BLUgtGREENgtRED.and(NIR.gt(0.254)).and( BLU.gt(0.165)).and(NDVI.lt(0.40))).or(
                      (BLUgtGREEN.and(BLU.gt(0.27)).and(GREEN.gt(0.21)).and( REDsubtractGREEN.lte(0.1)).and(NIR.gt(0.35)))).or(
                      (BLU.gt(0.94)).and(GREEN.gt(0.94)).and(RED.gt(0.94)).and(NIR.gt(0.94))))
                      
                      )),1);

     CLOUDSHAPE=0
 
     WETNESS=0
   
     OUT=OUT.where(( (OUT.eq(0)).and(NDVI.lte(0.4)).and(NDVI.gte(0)).and(ndti_Urban_th).and(bsi_Urban_th1).and(bsi_Urban_th2)  ),14 );//urban
     
     OUT=OUT.where(( (OUT.eq(0)).and(ndwi_Soils_th).and(bsi_Soil_th1)  ),7 );//soils
     
     OUT=OUT.where(( (ndvi_2).and(BLU.lt(0.13)).and(BLUgtGREENgtRED).and(RED.lt(0.09)).and( BLUsubtractNIR.lt(-0.04))  ),5   );//DARK_VEG
     
     
    
     OUT=OUT.where(( (OUT.eq(0)).and(ndvi_2).and(
                    ((BLU.lt(0.17)).and(BLU.gt(0.133)).and(BLUgtGREENgtRED).and(RED.lt(0.17)).and(NIR.lte(0.27)).and( ((NIR).subtract(BLU)).lte(0.06))).or(
                    ( ((((NIR.subtract(GREEN)).abs().lte(0.06)).add( BLUsubtractNIR.gte(-0.06))).gt(0)).and(BLUgtGREENgtRED).and(NIR.gte(0.22)) )).or(
                    ( (OUT.eq(0)).and(ndvi_2).and(NDVI.lte(0.24)).and(NIR.lte(0.23)).and(GREENlteRED).and(REDlteNIR)) )
                   )),6);//DARK_SOIL
  
     OUT=OUT.where(( (OUT.eq(0)).and(ndvi_2).and(NDVI.lte(0.19)).and(NIR.gt(0.330)).and(growing14)   ),7 );//BRIGHT_SOIL
    
     
     OUT=OUT.where(( (OUT.eq(0)).and(ndvi_2).and(NDVI.gte(0.51)).and(BLUgteGREEN).and(REDsubtractGREEN.lt(0.04))    ),8 );//GRASS
     
     OUT=OUT.where(( (OUT.eq(0)).and(ndvi_2).and(NDVI.gte(0.20)).and( REDsubtractGREEN.lt(0.05))   ),11 );//SPARSE
     
     OUT=OUT.where(( (OUT.eq(0)).and(ndvi_2)),12);//SOIL
     
     OUT=OUT.where(((ui_th.and(ndbi_th).and(ndvi_2)).and(NDVI.lte(0.20)).and(NIR.gt(0.3)).and(growing14) )     ,14);//urban areas
     
     REDsubtractGREEN=0
     BLUgteGREEN=0
     
    //-------------------------------------------------------------------------------------------------------------
		//----------------------  SECTION 3 : VEGETATION  -------------------------------------------------------------
		//-------------------------------------------------------------------------------------------------------------
    
    var MyCOND=(ndvi_3).and(NDVI.lt(th_RANGELAND));
    OUT=OUT.where(( (MyCOND).and(NIR.gte(0.15)) ),8); //GRASS
    OUT=OUT.where(( (MyCOND).and(NIR.lt(0.15))  ),4); //DARK_VEG
    
    MyCOND=(ndvi_3).and(NDVI.lt(th_GRASS));
    OUT=OUT.where(( (MyCOND).and(BLUlteNIR).and(NIR.lt(0.15))  ),4);//DARK_VEG
    OUT=OUT.where(( (OUT.eq(0)).and(
                                    ((MyCOND).and(BLUlteNIR)).or( (NDVI.lt(th_SHRUB) ).and(NIR.gt(0.22)))).and(NDSI.lt(-0.35))  ),16);//SHRUB
  
                                    
    OUT=OUT.where(( (MyCOND).and(BLU.gt(NIR)) ),4);//DARK_VEG
    OUT=OUT.where(( (OUT.eq(0)).and(MyCOND)).and(NDSI.lt(-0.3)),16); //SHRUB                                       
    OUT=OUT.where( (ndvi_3).and(OUT.eq(0)).and(NDSI.gt(-0.25))  ,10);//PURE_FOREST
    
    OUT=OUT.where(((OUT.eq(0)).and(ndvi_3).and(NDVI.gt(th_TREES)) ),9);//BRIGHT_FOREST
    
    OUT=OUT.where(( (OUT.eq(0)).and(NDVI.lt(th_GRASS))),8);//GRASS
    OUT=OUT.where(( (OUT.eq(0)).and(NDSI.lt(-0.25))),13);//DEGRADED_FOREST
    OUT=OUT.where(( (OUT.eq(0))),16);//SHRUB
    
   ///Agriculture
    OUT = OUT.where((SAVI1.gte(0.10).and (SAVI1.lte(0.25)) ).and((GNDVI.gte(0.30)).and ((GNDVI.lte(0.45))) ),15);
   
   

    
    return (OUT.select([0],["Class"]).toByte());
    
   
}   // SINGLE DATE CLASSIFICATION



// ----------------------------------------------------------------------------------------------------------
// ----------------------------------------------------------------------------------------------------------
// ------------------   PHENOLOGY BASED SYNTHESIS  ----------------------------------------------------------
// ----------------------------------------------------------------------------------------------------------
// ----------------------------------------------------------------------------------------------------------

// INPUT: collection of single band image classification (monthly SDC classification)
// OUT :  sinle image with Land Cover labels



exports.PBS_classification =  function (COLL){  

      var tot=COLL.map(function(res){return (res.gt(0)) }).sum().divide(100);
      var relTot=COLL.map(function(res){return (res.gte(2)) }).sum().divide(100);

      print('PBS classification ...');
      
      //----------------------------------------------------------------------------------------------------------
     var CLOUDS=COLL.map(function(res){return res.eq(1)}).sum().divide(1);//Set clouds index
     tot=0;
     var DARK_VEG=COLL.map(function(res){return res.eq(4)}).sum().divide(relTot);
     var LOW_IL = COLL.map(function(res){return res.eq(5)}).sum().divide(relTot);
     
     var DARK_SOIL= COLL.map(function(res){return res.eq(6)}).sum().divide(relTot);
     var BRIGHT_SOIL=COLL.map(function(res){return res.eq(7)}).sum().divide(relTot);
     var SPARSE=COLL.map(function(res){return res.eq(11)}).sum().divide(relTot);
     var SOIL=COLL.map(function(res){return res.eq(12)}).sum().divide(relTot);
     
     var URBAN=COLL.map(function(res){return res.eq(14)}).sum().divide(relTot);
     var BRIGHT_FOREST=COLL.map(function(res){return res.eq(9)}).sum().divide(relTot);
     
     var DEGRADED_FOREST=COLL.map(function(res){return res.eq(13)}).sum().divide(relTot);
     var PURE_FOREST=COLL.map(function(res){return res.eq(10)}).sum().divide(relTot);
     var FOREST=PURE_FOREST.add(DEGRADED_FOREST)
     var SHRUB=COLL.map(function(res){return (res.eq(16))}).sum().divide(relTot);
     
     var GRASS= COLL.map(function(res){return res.eq(8)}).sum().divide(relTot);
     
     var EVERGREEN_FOREST= (BRIGHT_FOREST.add(FOREST).add(DARK_VEG));
     
     var CROP = COLL.map(function(res){return res.eq(15)}).sum().divide(relTot);
     
     var WATER=COLL.map(function(res){return res.eq(2)}).sum().divide(relTot);
     var SNOW=COLL.map(function(res){return res.eq(3)}).sum().divide(relTot);
     
     var DRY = SPARSE.add(SOIL).add(BRIGHT_SOIL).add(DARK_SOIL).add(LOW_IL);
      
     print('Valiables are now set ... ');
      
     var OUT =relTot.multiply(0).select([0],["Class"]).toByte(); 
     
     
     
     OUT=OUT.where( (CLOUDS.gt(9)),10); //no data
     
     OUT=OUT.where( (URBAN.gte(80))  ,8);  // urban
     
     OUT=OUT.where( (WATER.gte(80))  ,1);  // water
     
     OUT=OUT.where( (FOREST.gte(80))  ,2);  // FOREST
     
     
     var GREEN=(EVERGREEN_FOREST.add(SHRUB).add(GRASS));
     var BARE= (BRIGHT_SOIL.add(SOIL).add(DARK_SOIL));
    
     OUT = OUT.where( (GREEN.eq(0)).and(BARE.gt(5)).and(URBAN.gte(50)), 8 ); //urban

     OUT = OUT.where( (BARE.gt(70)).and(URBAN.gt(10)), 7 ); //baresoils
     
     OUT = OUT.where( (BARE.gt(50)).and(URBAN.gt(20)), 7 ); //baresoils
     
     OUT = OUT.where( (BARE.gt(40)).and(URBAN.gt(15)), 7 ); //baresoils
     
     OUT = OUT.where( (BARE.gt(20)).and(URBAN.gt(50)), 8 ); //urban
     
     OUT = OUT.where( (OUT.eq(0)).and(BARE.add(URBAN).gt(65). and (GREEN.lt(25))), 7 ); //baresoil
     
     OUT = OUT.where( (OUT.eq(0)).and(DRY.gt(40). and (GREEN.lt(50))), 7 ); //baresoil
     
     OUT = OUT.where( (DRY.gt(60)).and(GREEN.lt(10)), 8 ); //highway urban
     
     OUT = OUT.where( (DRY.gt(5)).and(URBAN.gt(10)).and(GREEN.lt(25)).and(WATER.lt(5)).and((SPARSE.add(CROP)).lt(50)), 8 ); //highway urban
     
     OUT = OUT.where( (OUT.eq(0)).and(DRY.lt(30). and (WATER.gt(50))), 7 ); //baresoil
    
     OUT = OUT.where( (CLOUDS.gt(9)), 10 ); //nodata
    
     OUT=OUT.where( (OUT.eq(0)).and(GREEN.gt(65)).and((GRASS.add(SHRUB)).gt(20)).and(DRY.lt(35)) ,5);  // DECIDE IF OPEN DRY or HUMID
     
    
     OUT=OUT.where( ((GREEN.gt(90)).and(DRY.lt(5))).and(FOREST.lt(50)),2);//forest
    
     OUT=OUT.where( (GREEN.gt(70)).and(SHRUB.gt(20)).and(FOREST.gt(50)).and((GRASS.add(SHRUB)).gt(30)).and(WATER.lt(5)).and(URBAN.eq(0))   ,2);//forest 
     
    
     OUT=OUT.where( (GREEN.gt(85)).and((SHRUB.add(DEGRADED_FOREST)).gt(70)).and(DRY.lt(5)) ,2);
     
     OUT=OUT.where(  (GREEN.gte(95)).and(FOREST.gte(30)).and(SHRUB.gt(40)).and(EVERGREEN_FOREST.gt(40)).and(WATER.lt(5)).and(URBAN.eq(0))   ,2);  // almost evergreen   
     
     OUT=OUT.where( (FOREST.gt(60)).and(EVERGREEN_FOREST.gt(80)).and(WATER.lt(5)).and(URBAN.eq(0)),2);
     OUT=OUT.where( (FOREST.gt(70)).and(GREEN.eq(FOREST)).and(WATER.lt(5)).and(URBAN.eq(0)),2);
     OUT=OUT.where( (FOREST.gt(50)).and(GREEN.gt(60)).and(((DARK_VEG).add(LOW_IL).add(GREEN)).gt(90)).and(WATER.lt(5)).and(URBAN.eq(0)),2);
     
     OUT=OUT.where( (FOREST.gt(50)).and((EVERGREEN_FOREST.add(SHRUB)).gt(80)).and(DRY.eq(0)).and(WATER.lt(5)).and(URBAN.eq(0)),2);
     
     
     if (EVERGREEN_FOREST.gt(30)){
           var condAC1 = (((EVERGREEN_FOREST.add(SPARSE).add(GRASS)).gt(10)).and((EVERGREEN_FOREST.add(SPARSE).add(GRASS)).lt(60)));
           var condAC2 = ((WATER.gt(10)).and(WATER.lt(50)));
           var condAC3 = ((DRY.add(URBAN).gt(5)).and(DRY.add(URBAN).lt(35)))
           
           OUT=OUT.where(condAC1.and(condAC2).and(condAC3) ,3);//annual crop
     }
     
     
     OUT=OUT.where( (FOREST.gt(90)).and(WATER.lt(5)).and(URBAN.eq(0)),2);
      
     OUT=OUT.where( (WATER.gt(60)),1);
     OUT=OUT.where( (WATER.gt(80)),1);
     
     OUT=OUT.where( (OUT.eq(0)).and((DRY.add(WATER).add(SNOW)).gt(90)).and(WATER.gte(20)) ,1);//water
     
     OUT=OUT.where(((DARK_SOIL.add(LOW_IL)).gte(90)),8);//bare soil
     
     
     var MYCOND=(OUT.eq(0)).and((DRY.add(GRASS)).gt(65)).and(EVERGREEN_FOREST.eq(0));
     OUT=OUT.where( MYCOND.and(GRASS.gte(DRY)).and(SHRUB.lte(5)) ,5); //  bush + grass
     OUT=OUT.where( MYCOND.and(DRY.gt(GRASS)).and(SHRUB.eq(0)) ,5);  // grass 
     MYCOND=0;
     OUT=OUT.where( (OUT.eq(0)).and(GRASS.gte(80)).and((GRASS.add(DRY)).gt(95)) ,5);//grassland
     
   
     
     OUT=OUT.where( (DRY.gte(90)) ,7);//baresoil
     OUT=OUT.where( (SOIL.gte(90)) ,7);
     OUT=OUT.where((BRIGHT_SOIL.gte(90)) ,7);
     OUT=OUT.where( (SPARSE.gte(90)) ,7);
     OUT=OUT.where((URBAN.gte(50)).and(DRY.gt(20)).and(GREEN.eq(0)),8);
     
    if (EVERGREEN_FOREST.gt(30)){
        OUT=OUT.where((CLOUDS.lt(4)).and((BARE.add(SPARSE)).gt(40)).and((CROP.add(EVERGREEN_FOREST)).gt(30)).and(WATER.gte(0)) ,3);//annual crop ***
    }
    
    OUT=OUT.where( (SNOW.gte(80)),9);
     
    OUT=OUT.where(OUT.eq(0).and(FOREST.gte(45)).and(WATER.lt(5)).and(URBAN.eq(0)),2);
    OUT=OUT.where(OUT.eq(0).and(FOREST.gte(25)).and(WATER.lt(5)).and(URBAN.eq(0)),2);
    OUT=OUT.where(OUT.eq(0),5);
    if (EVERGREEN_FOREST.gt(30)){
        OUT=OUT.where((OUT.eq(5)).and((EVERGREEN_FOREST.add(SPARSE)).gt(70)).and(WATER.gte(15)) ,3);//annual crop
        OUT=OUT.where((OUT.eq(5)).and((EVERGREEN_FOREST.add(SPARSE)).gt(30)).and(WATER.gte(15)).and(URBAN.add(DRY).gte(15)) ,3);//annual crop
        OUT=OUT.where((OUT.eq(5)).and((EVERGREEN_FOREST.add(SPARSE)).gt(30)).and(WATER.lte(10)).and(URBAN.add(DRY).gte(15)) ,3);//annual crop
        
        OUT=OUT.where((OUT.eq(5)).and((EVERGREEN_FOREST.add(SPARSE)).gt(15)).and(WATER.gte(5)).and(URBAN.add(DRY).gte(5)).and(CROP.gt(5)) ,3);//annual crop
        OUT=OUT.where((OUT.eq(5)).and((EVERGREEN_FOREST.add(SPARSE)).gt(50)).and(WATER.gte(5)).and(URBAN.add(DRY).gte(5)).and(CROP.gt(0)) ,3);//annual crop
        OUT=OUT.where((OUT.eq(5)).and((EVERGREEN_FOREST.add(SPARSE)).gt(15)).and(WATER.eq(0)).and(URBAN.add(DRY).gte(5)).and(CROP.lt(10)) ,3);//annual crop
     
        /////////manage two image
     
        OUT=OUT.where((CLOUDS.gt(9)).and((EVERGREEN_FOREST.add(SPARSE)).gt(70)).and(WATER.gt(0)) ,3);//annual crop
        OUT=OUT.where((CLOUDS.lte(5)).and(FOREST.gte(DEGRADED_FOREST)).and(DEGRADED_FOREST.gte(20)).and(WATER.eq(0))
                    .and((URBAN.add(DRY)).gte(0)).and(GRASS.gte(10)).and(PURE_FOREST.lt(10)).and(SHRUB.gt(0)) ,3);//annual crop
     }
     
      OUT = OUT.where( (CLOUDS.gt(7)), 10 ); //nodata  
      
     
         
      ///Perennial crop
      
      OUT=OUT.where((CLOUDS.lte(5)).and(FOREST.gte(DEGRADED_FOREST)).and(DEGRADED_FOREST.gte(30)).and(WATER.eq(0))
                    .and(URBAN.eq(0)).and(GRASS.gt(0)).and(PURE_FOREST.lt(22)).and(SHRUB.eq(0)) ,4);//perennial crop crop
      OUT=OUT.where((CLOUDS.lte(5)).and(DEGRADED_FOREST.gte(FOREST)).and(DEGRADED_FOREST.gte(30)).and(WATER.eq(0))
                    .and((DRY.add(URBAN)).lt(20)).and(GRASS.gt(0)).and(SHRUB.gte(0)) ,4);//perennial crop crop
      OUT=OUT.where((CLOUDS.lte(5)).and(FOREST.gte(10)).and(DEGRADED_FOREST.gte(65)).and(WATER.eq(0))
                    .and((DRY.add(URBAN)).lt(20)).and(GRASS.gt(0)).and(SHRUB.gte(0)) ,4);//perennial crop crop              

      //Forest (and dry forest)
      
      var condF1=(EVERGREEN_FOREST.gt(0)).and(EVERGREEN_FOREST.lte(60));
      var condF2=((DEGRADED_FOREST.add(SPARSE)).gt(5)).and((DEGRADED_FOREST.add(SPARSE)).lte(60));
      var condF3=((SOIL.add(DARK_SOIL).add(URBAN)).gt(0)).and((SOIL.add(DARK_SOIL).add(URBAN)).lte(25));
      var condF4=(DEGRADED_FOREST.gt(0)).and(DEGRADED_FOREST.lte(50));
      var condF5=((PURE_FOREST).gt(10)).and((PURE_FOREST).lte(80));
      OUT=OUT.where((CLOUDS.lte(4)).and(PURE_FOREST.gt(50)).and(WATER.eq(0)).and((SOIL.add(DARK_SOIL).add(URBAN)).lte(5)) ,2);//Forest
      
    if (EVERGREEN_FOREST.gt(25)){
       //Shrub
       var condSh1 = (CROP.add(DEGRADED_FOREST.add(DARK_VEG)).gt(0)).and(CROP.add(DEGRADED_FOREST.add(DARK_VEG)).lt(30));
       var condSh2 = ((SOIL.add(DARK_SOIL)).gt(0)).and((SOIL.add(DARK_SOIL)).lte(50));
       var condSh3 = (SHRUB.gt(5)).and(SHRUB.lt(70));
       var condSh4 = ((DEGRADED_FOREST.add(SPARSE)).gt(20)).and((DEGRADED_FOREST.add(SPARSE)).lte(60));
       var condSh5 = (DEGRADED_FOREST.gt(20)).and(DEGRADED_FOREST.lte(50));
       var condSh6 = (URBAN.gt(10)).and(URBAN.lte(56));
       
       
       OUT = OUT.where((CLOUDS.lte(4)).and(condSh1).and(condSh2).and(WATER.eq(0).and(condSh3)),6); //shrubs
       OUT = OUT.where((CLOUDS.lte(4)).and(condSh2).and(condSh4).and(WATER.eq(0)),3);//annual crop ---> Shrubs
       OUT = OUT.where((CLOUDS.lte(4)).and(condSh1).and(condSh2).and(condSh4).and(WATER.eq(0)),6);//annual crop ---> Shrubs
       OUT = OUT.where((CLOUDS.lte(3)).and(OUT.eq(3)).and(condSh2).and(condSh4).and(WATER.eq(0)),6);//annual crop ---> Shrubs
       OUT = OUT.where((CLOUDS.lte(3)).and(OUT.eq(3)).and(condSh2).and(condSh3).and(condSh5).and(WATER.eq(0)),6);//annual crop ---> Shrubs
       OUT = OUT.where((CLOUDS.lte(6)).and(OUT.eq(3)).and(condSh5).and(condSh6).and(WATER.eq(0)),6);//annual crop ---> Shrubs
     
     
    }
     
     return OUT
}  //  END PBS