Guest User

Untitled

a guest
Jan 16th, 2019
62
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.09 KB | None | 0 0
  1. var geometry = /* color: #d63000 */ee.Geometry.Point([5.7788, 52.7005]);
  2.  
  3.  
  4. var country = ee.FeatureCollection('ft:1tdSwUL7MVpOauSgRzqVTOwdfy17KDbw-1d9omPw')
  5. .filter(ee.Filter.eq('Country', 'Netherlands'));
  6.  
  7. // Get the VV collection.
  8. var collectionVV = ee.ImageCollection('COPERNICUS/S1_GRD')
  9. .filter(ee.Filter.eq('instrumentMode', 'IW'))
  10. .filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VV'))
  11. .filter(ee.Filter.eq('orbitProperties_pass', 'DESCENDING'))
  12. .select(['VV']);
  13.  
  14.  
  15.  
  16. // Create a 3 band stack by selecting from different periods (months)
  17. var im1 = ee.Image(collectionVV.filterDate('2016-04-01', '2016-05-30').mean());
  18. var im2 = ee.Image(active_collection.filterDate('2016-06-01', '2016-08-31').mean());
  19.  
  20. var roi1 = im1.clip(country)
  21. var roi2 = im2.clip(country)
  22.  
  23.  
  24. var change = im2.subtract(im1).rename('CHANGE');
  25.  
  26.  
  27. // Display the result.
  28. Map.centerObject(geometry, 13);
  29.  
  30. // Change
  31. var nbrParams = {min: -25, max: 0, palette: ['blue', 'red', 'green']};
  32. Map.addLayer(change, nbrParams, 'Change image');
  33.  
  34. // Stack
  35. Map.addLayer(im1.addBands(im2), {min: -25, max: 0}, 'VV stack');
Add Comment
Please, Sign In to add comment