-
Notifications
You must be signed in to change notification settings - Fork 0
/
SDC_PBS_Coast
525 lines (357 loc) · 24.6 KB
/
SDC_PBS_Coast
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
/* 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;
var th_NDVI_MIN_VEGE=0.4;
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.80 ;
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.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.05;
var ui_th = 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)
var bsiWay1 = -0.1;
var bsiWay2 = 0.1;
var bsiWay_th1 = BSI.gt(bsiWay1);
var bsiWay_th2 = BSI.lt(bsiWay2);
///---------------------------------------GNDVI Index------------------------------------------------------//
//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 Index------------------------------------------------------//
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')
})
///---------------------------------------SAVI Index------------------------------------------------------//
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 b8
})
var SAVI1 = ((NIR.subtract(RED)).divide(RED.add(NIR).add(0.5))).multiply(1.5)
var saviThreshold1 = -0.1;
var savi_th1 = SAVI1.lte(saviThreshold1);//MinForest
var saviThreshold2 = -0.3;
var savi_th2 = SAVI1.gte(saviThreshold2);//MaxForest
//-------------------------------------------------------------------------------------------------------------
//---------------------- 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(bsiWay_th1).and(bsiWay_th2) ),14 );//Highway
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) ) .and(NDVI.gte(0.20)).and( REDsubtractGREEN.lt(0.05))
.and( (ndvi_2).and(BLU.lt(0.13)).and(BLUgtGREENgtRED).and(RED.lt(0.05)).and( BLUsubtractNIR.lt(0.04))) ,14);//urban areas
REDsubtractGREEN=0
BLUgteGREEN=0
//-------------------------------------------------------------------------------------------------------------
//---------------------- SECTION 3 : VEGETATION -------------------------------------------------------------
//-------------------------------------------------------------------------------------------------------------
OUT=OUT.where( (NDVI.gte(0.55)).and(OUT.eq(0)).and(NDSI.gt(-0.25)) ,10);//PURE_FOREST
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(((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
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( (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( (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);
OUT=OUT.where( (GREEN.gt(40)).and(SHRUB.add(GRASS).gt(20)).and(WATER.eq(0)) ,4)//perennial crops
OUT=OUT.where( (SHRUB.add(GRASS).gt(20).and(CROP).gt(10)) ,4)//perennial crops
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);//
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)) ,8);//baresoil
OUT=OUT.where( (SOIL.gte(90)) ,8);
OUT=OUT.where((BRIGHT_SOIL.gte(90)) ,8);
OUT=OUT.where( (SPARSE.gte(90)) ,8);
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);
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
OUT = OUT.where( (CLOUDS.gt(7)), 10 ); //nodata
OUT=OUT.where((CLOUDS.gt(9)).and((EVERGREEN_FOREST.add(SPARSE)).gt(70)).and(WATER.gt(0)) ,3);//annual crop
///Perennial crop
OUT=OUT.where((CLOUDS.lte(3)).and((EVERGREEN_FOREST.add(SPARSE).add(CROP)).gt(80)).and(WATER.eq(0)).and(URBAN.eq(0)) ,4);//perennial crop crop
OUT=OUT.where((CLOUDS.lte(3)).and((GREEN.add(SPARSE).add(CROP)).gt(80)).and(WATER.eq(0)).and(DRY.eq(0)).and(URBAN.eq(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(condF1).and(condF2).and(condF3).and(condF4).and(WATER.eq(0)) ,2);//Forest
OUT=OUT.where((CLOUDS.lte(4)).and(EVERGREEN_FOREST.gt(80)).and(WATER.eq(0)).and((SOIL.add(DARK_SOIL).add(URBAN)).lte(5)) ,2);//Forest
OUT=OUT.where((CLOUDS.lte(4)).and(condF5).and(WATER.eq(0)).and(condF3) ,2);//Forest (could be perennials)
//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)),6);//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
OUT=OUT.addBands(WATER.select([0],["WATER"]));
return OUT
} // END PBS