在GEE平台提取Sentinel-1 SAR GRD的VV+VH波段

1. 数据描述

The Sentinel-1 mission provides data from a dual-polarization C-band Synthetic Aperture Radar (SAR) instrument at 5.405GHz (C band). This collection includes the S1 Ground Range Detected (GRD) scenes, processed using the Sentinel-1 Toolbox to generate a calibrated, ortho-corrected product. The collection is updated daily. New assets are ingested within two days after they become available.

This collection contains all of the GRD scenes. Each scene has one of 3 resolutions (10, 25 or 40 meters), 4 band combinations (corresponding to scene polarization) and 3 instrument modes. Use of the collection in a mosaic context will likely require filtering down to a homogeneous set of bands and parameters. See this article for details of collection use and preprocessing. Each scene contains either 1 or 2 out of 4 possible polarization bands, depending on the instrument’s polarization settings. The possible combinations are single band VV or HH, and dual band VV+VH and HH+HV:

VV: single co-polarization, vertical transmit/vertical receive
HH: single co-polarization, horizontal transmit/horizontal receive
VV + VH: dual-band cross-polarization, vertical transmit/horizontal receive
HH + HV: dual-band cross-polarization, horizontal transmit/vertical receive

2. 代码实现

// Load the Sentinel-1 ImageCollection.
var sentinel1 = ee.ImageCollection('COPERNICUS/S1_GRD');

// Filter by metadata properties.
var vh = sentinel1
  // Filter to get images with VV and VH dual polarization.
  .filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VV'))
  .filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VH'))
  // Filter to get images collected in interferometric wide swath mode.
  .filter(ee.Filter.eq('instrumentMode', 'IW'))
 .map(function(image) {
          var edge = image.lt(-30.0);
          var maskedImage = image.mask().and(edge.not());
          return image.updateMask(maskedImage).clip(jianghan);
        });

// Filter to get images from different look angles.
var vhAscending = vh.filter(ee.Filter.eq('orbitProperties_pass', 'ASCENDING'));
var vhDescending = vh.filter(ee.Filter.eq('orbitProperties_pass', 'DESCENDING'));
var spring = ee.Filter.date('2019-02-23', '2019-12-31');
// Create a composite from means at different polarizations and look angles.
var composite = ee.Image.cat([
  // vhAscending.select('VH').mean(),
  ee.ImageCollection(vhAscending.select('VH').merge(vhDescending.select('VH'))).filter(spring).mean(),
  ee.ImageCollection(vhAscending.select('VV').merge(vhDescending.select('VV'))).filter(spring).mean(),
  // vhDescending.select('VH').mean()
]).focal_median();
print('composite',composite);
// Display as a composite of polarization and backscattering characteristics.
Map.addLayer(composite, {min: -20, max: 0}, 'composite');

3. 结果