Advertisement
Guest User

Untitled

a guest
Feb 22nd, 2019
76
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.99 KB | None | 0 0
  1. var image = ee.Image(ee.ImageCollection('LANDSAT/LC8_L1T_TOA')
  2. .filterBounds(roi)
  3. .filterDate('2015-01-01', '2015-12-31')
  4. .sort('CLOUD_COVER')
  5. .first());
  6. print(image);
  7.  
  8.  
  9. Map.addLayer(image, {bands: ['B4', 'B3', 'B2'], max: 0.3}, 'image_L8');
  10.  
  11.  
  12. var clip = image.clip(table); //cliping image to shapefile
  13.  
  14.  
  15. var newclassified =
  16. vegetation.merge(water).merge(agri_bare).merge(builtarea);
  17. print(newclassified);
  18. var bands = ['B2', 'B3', 'B4', 'B5', 'B6', 'B7', 'B10', 'B11'];
  19.  
  20.  
  21. var training = clip.select(bands).sampleRegions({
  22. collection: newclassified,
  23. properties: ['landcover'],
  24. scale: 30
  25. });
  26.  
  27.  
  28. var classifier = ee.Classifier.cart().train({
  29. features: training,
  30. classProperty: 'landcover',
  31. inputProperties: bands
  32. });
  33.  
  34. var classified = clip.select(bands).classify(classifier);
  35.  
  36. Map.centerObject(newclassified, 11);
  37. Map.addLayer(clip,
  38. {bands: ['B4', 'B3', 'B2'], max: 0.4},
  39. 'Landsat image');
  40.  
  41. Map.addLayer(classified,
  42. {min: 1, max: 4, palette: ['green', 'blue', 'yellow', 'red']},
  43. 'classification');
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement