Guest User

Untitled

a guest
Mar 19th, 2018
91
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.76 KB | None | 0 0
  1. var image=ee.Image('LANDSAT/LC08/C01/T1_SR/LC08_169062_20180205');
  2.  
  3. //define parameters
  4.  
  5. var green = image.select('B3');
  6. var swi = image.select('B6');
  7. var mndwi = green.subtract(swi).divide(green.add(swi)).rename('MNDWI');
  8.  
  9. var mndwipara = {min: 0, max: 0.7, palette: ['white', 'blue']};
  10.  
  11.  
  12. Map.addLayer(mndwi, mndwipara, 'MNDWI');
  13.  
  14.  
  15. //apply threshold to select only positive pixels
  16. var lake_mask= mndwi.gt(0);
  17. Map.addLayer(lake_mask);
  18.  
  19.  
  20. //classification attempt - scum
  21. var scum=image.expression(
  22. '(nir)*0 + (nir < 0.018)*1', {
  23. 'nir': image.select('B5'),
  24. 'red': image.select ('B4'),
  25. 'green': image.select('B3'),
  26. });
  27.  
  28. //classification attempt - sediment
  29. var sediment=image.expression(
  30. '(nir)*0 + (nir > 0.018 && red > green && green/nir < 1.3 && red/green>
  31. nir/red)*1', {
  32. 'nir': image.select('B5'),
  33. 'red': image.select ('B4'),
  34. 'green': image.select('B3'),
  35. });
  36.  
  37. //classification attempt - lowbiomass
  38. var lowbiomass=image.expression(
  39. '(nir)*0 + (nir > 0.018 && red < green && green > 0.065 && nir < red)*1', {
  40. 'nir': image.select('B5'),
  41. 'red': image.select ('B4'),
  42. 'green': image.select('B3'),
  43. });
  44.  
  45.  
  46. //classification attempt-scum
  47. var scum=image.expression(
  48. '(nir)*0 + (nir > 0.018 && red < green && green > 0.065 && nir > red && nir/red > 4.07)*1', {
  49. 'nir': image.select('B5'),
  50. 'red': image.select ('B4'),
  51. 'green': image.select('B3'),
  52. });
  53.  
  54. //classifcation attempt bleached scum
  55. var bleached_scum=image.expression(
  56. '(nir)*0 + (nir > 0.018 && red > green && green/nir > 1.3 && green > 0.065)*1', {
  57. 'nir': image.select('B5'),
  58. 'red': image.select ('B4'),
  59. 'green': image.select('B3'),
  60. });
  61.  
  62. //classification attempt - microphytobenthos
  63. var microphytobenthos=image.expression(
  64. '(nir)*0 + (nir > 0.018 && red > green && green/nir < 1.3 && red/green< nir/red)*1', {
  65. 'nir': image.select('B5'),
  66. 'red': image.select ('B4'),
  67. 'green': image.select('B3'),
  68. });
  69.  
  70. //classification attempt - highbiomass
  71. var highbiomass=image.expression(
  72. '(nir)*0 + (nir > 0.018 && red < green && green > 0.065 && nir > red && nir/red < 4.07)*1', {
  73. 'nir': image.select('B5'),
  74. 'red': image.select ('B4'),
  75. 'green': image.select('B3'),
  76. });
  77.  
  78. /*
  79. lake = 1
  80. scum = 2
  81. bleached_scum = 3
  82. microphytobenthos = 4
  83. highbiomass = 5
  84. sediment
  85.  
  86. */
  87.  
  88.  
  89. var lake_mask_ = ee.Image(1).mask(lake_mask).toInt();
  90. var scum_ = ee.Image(2).mask(lake_mask.mask(scum)).toInt();
  91. var bleached_scum_ =
  92. ee.Image(3).mask(lake_mask.mask(bleached_scum)).toInt();
  93. var microphytobenthos_ =
  94. ee.Image(4).mask(lake_mask.mask(microphytobenthos)).toInt();
  95. var biomass_ = ee.Image(5).mask(lake_mask.mask(highbiomass)).toInt();
  96.  
  97. var test = ee.ImageCollection([lake_mask_,scum_,bleached_scum_,microphytobenthos_,biomass_,]);
  98.  
  99. var test = test.reduce(ee.Reducer.max());
  100.  
  101. print(test);
  102. Map.addLayer(test,{min: 1, max: 5, palette: ['0000FF',
  103. '1BCFFF','66ff33','F4FF0B','FA0007']},'Classification');
Add Comment
Please, Sign In to add comment