Advertisement
Guest User

Untitled

a guest
Jul 18th, 2019
130
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.13 KB | None | 0 0
  1. var l5 = ee.ImageCollection("LANDSAT/LT05/C01/T1_SR"),
  2. roi = /* color: #d63000 */ee.Geometry.Polygon(
  3. [[[77.68018052033392, 31.0091863423649],
  4. [80.36084458283392, 29.5863837026937],
  5. [80.98706528595892, 30.233944725529557],
  6. [78.85571762970892, 31.57247908795104]]]),
  7. elevation = ee.Image("USGS/SRTMGL1_003"),
  8. l7 = ee.ImageCollection("LANDSAT/LE07/C01/T1_RT_TOA");
  9.  
  10. // Defining a filter that gives results for months September to November
  11. var mf = ee.Filter.calendarRange(9, 11, 'month');
  12.  
  13. Map.centerObject(roi, 9);
  14.  
  15. // Definig the image collections
  16. var img1 = ee.ImageCollection(l5
  17. .filterDate('1984-01-01','1999-12-31')
  18. .filterBounds(roi)
  19. .filter(mf));
  20. var img2 = ee.ImageCollection(l7
  21. .filterDate('2000-01-01','2018-12-31')
  22. .filterBounds(roi)
  23. .filter(mf));
  24.  
  25. // add NDWI band function
  26. var addNDWI = function(image) {
  27. return image
  28. // NDWI
  29. .addBands(image.normalizedDifference(['B2', 'B4']).rename('NDWI'))
  30. };
  31.  
  32. // map add ndwi function to both image collections
  33. var image1 = img1.map(addNDWI);
  34. var image2 = img2.map(addNDWI);
  35.  
  36. var test2 = image2.median();
  37. Map.addLayer(test2.clip(roi), {}, 'image2');
  38.  
  39. // masking function for elevation
  40. var mask = function(image) {
  41. var elevMask = elevation.updateMask(elevation.gt(1500));
  42. return image.updateMask(elevMask);
  43. };
  44.  
  45. // masking function for NDWI
  46. var mask1 = function(image){
  47. var ndwiMask = image.select('NDWI').gte(0.4);
  48. return image.updateMask(ndwiMask);
  49. };
  50.  
  51. // map masking functions
  52. var maskedImg = image1.map(mask).map(mask1); //For 1984-1999
  53. var maskedImg2 = image2.map(mask).map(mask1); //For 2000-2018
  54.  
  55. // getting rid of bad pixels
  56. function trimImage(image){
  57. return image.clip(image.geometry().buffer(-5000));
  58. }
  59.  
  60. maskedImg2 = maskedImg2.map(trimImage);
  61. print(maskedImg);
  62. print(maskedImg2);
  63.  
  64. //set up params to view NDWI band
  65. var ndwiViz = {bands: 'NDWI',min: 0.3, palette: ['00FFFF', '0000FF']};
  66.  
  67. var test = maskedImg.median();
  68. var test1 = maskedImg2.median();
  69.  
  70. Map.addLayer(test.clip(roi), ndwiViz, 'NDWI 1984-1999');
  71. Map.addLayer(test1.clip(roi), ndwiViz, 'NDWI 2000-2018');
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement