Guest User

Untitled

a guest
Jul 21st, 2018
106
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.93 KB | None | 0 0
  1. var geometry = /* color: #98ff00 */ee.Geometry.Polygon(
  2. [[[73.85016678395664, 14.304113379426035],
  3. [73.85016583974448, 14.048667693939866],
  4. [74.0477686190103, 14.155139129529527]]]),
  5. point1 = /* color: #d63000 */ee.Geometry.Point([73.91532897949219, 14.20980182205916]),
  6. point2 = /* color: #98ff00 */ee.Geometry.MultiPoint(
  7. [[73.9105224609375, 14.169860112219006],
  8. [73.87619018554688, 14.186503346854852],
  9. [73.87550354003906, 14.233097898500935],
  10. [73.93730163574219, 14.199817054773131],
  11. [73.95378112792969, 14.181177644664091],
  12. [73.88717651367188, 14.155213055929545],
  13. [73.92906188964844, 14.141230893912757],
  14. [73.94691467285156, 14.155878851710582],
  15. [73.97300720214844, 14.161205147702967],
  16. [73.98399353027344, 14.165865554216062],
  17. [73.9764404296875, 14.197820048470204]]);
  18.  
  19. var image = ee.ImageCollection("LANDSAT/LC8_L1T_TOA").filterBounds(geometry)
  20. .filterDate('2013-01-01', '2013-04-14');
  21.  
  22. var addNDWI = function(image) {
  23. var ndwi = image.normalizedDifference(['B3', 'B6']).rename('NDWI');
  24. return image.addBands(ndwi);
  25. };
  26. var image1 = image.map(addNDWI);
  27.  
  28. var fg_points =
  29. [ee.Feature(geometry),ee.Feature(point1),ee.Feature(point2)];
  30. Map.addLayer(ee.Feature(geometry), {}, "geometry");
  31. Map.addLayer(ee.Feature(point1), {}, "point1");
  32. Map.addLayer(ee.Feature(point2), {}, "point2");
  33.  
  34.  
  35. var mapfunc = function(feat) {
  36. var geom = feat.geometry()
  37. var addProp = function(img, f) {
  38. var newf = ee.Feature(f)
  39. var date = img.date().format()
  40. var value = img.reduceRegion(ee.Reducer.first(), geom, 30).get('NDWI')
  41. return ee.Feature(ee.Algorithms.If(value,
  42. newf.set(date, ee.String(value)),
  43. newf.set(date, ee.String('No data'))))
  44. }
  45. var newfeat = ee.Feature(image1.iterate(addProp, feat))
  46. return newfeat
  47. };
  48.  
  49.  
  50. var newft = fg_points.map(mapfunc);
  51.  
  52. Export.table.toDrive(newft,
  53. "export_Points",
  54. "export_Points",
  55. "export_Points");
Add Comment
Please, Sign In to add comment