Advertisement
Dragon1573

GEE Code

Mar 4th, 2025
134
0
6 days
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
JavaScript 1.91 KB | Source Code | 0 0
  1. var district = ee.FeatureCollection("projects/gee-cindy-250304/assets/ws");
  2.  
  3. var dsize = district.size();
  4. print(dsize);
  5.  
  6. var district_geometry = district.geometry();
  7.  
  8. Map.centerObject(district_geometry,12);
  9. Map.addLayer(district);
  10.  
  11. function maskS2clouds(image) {
  12.   var qa = image.select('QA60');
  13.  
  14.   // Bits 10 and 11 are clouds and cirrus, respectively.
  15.   var cloudBitMask = 1 << 10;
  16.   var cirrusBitMask = 1 << 11;
  17.  
  18.   // Both flags should be set to zero, indicating clear conditions.
  19.   var mask = qa.bitwiseAnd(cloudBitMask).eq(0)
  20.       .and(qa.bitwiseAnd(cirrusBitMask).eq(0));
  21.  
  22.   return image.updateMask(mask).divide(10000);
  23. }
  24.  
  25.  
  26. // Map the function over one year of data.
  27. // Load Sentinel-2 TOA reflectance data.
  28. var dataset = ee.ImageCollection("COPERNICUS/S2_HARMONIZED")
  29.                   .filterBounds(district_geometry)
  30.                   .filterDate('2019-01-01', '2020-01-01')
  31.                   // Pre-filter to get less cloudy granules.
  32.                   .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 10))
  33.                   .map(maskS2clouds);
  34.  
  35. var rgbVis = {
  36.   min: 0.0,
  37.   max: 0.3,
  38.   bands: ['B4', 'B3', 'B2'],
  39. };
  40.  
  41.  
  42. Map.addLayer(dataset.median(), rgbVis, 'RGB');
  43.  
  44. //export data
  45. var exportdataset =  ee.ImageCollection('COPERNICUS/S2_HARMONIZED')
  46.                   .filterBounds(district_geometry)
  47.                   .filterDate('2019-01-01', '2020-01-01')
  48.                   // Pre-filter to get less cloudy granules.
  49.                   .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 10))
  50.                   .map(maskS2clouds)
  51.                   .select(['B4', 'B3', 'B2']);
  52.                  
  53.   var mosaic = exportdataset.mosaic();
  54.  
  55.   Export.image.toDrive({
  56.       image: mosaic,
  57.       description:'Yunshan-Yunnan',
  58.       scale: 10,
  59.       maxPixels: 1e13,
  60.       region:district_geometry,
  61.       fileFormat: 'GeoTIFF',
  62.       formatOptions: {
  63.         cloudOptimized: true
  64.       }
  65.     });
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement