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

该博客介绍了如何使用Sentinel-1的C波段双极化合成孔径雷达数据进行处理。内容包括Sentinel-1任务的数据特性,如分辨率、极化组合和仪器模式。通过代码展示了如何过滤并选择具有VV和VH极化模式的干涉宽幅模式图像,以及如何根据上升和下降轨道获取不同视角的影像。最后,利用2019年春季的数据创建了极化和后向散射特征的复合图像。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

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. 结果请添加图片描述

GEE(Google Earth Engine)是一个强大的云平台,用于处理地理空间数据。在这个平台上,用户可以访问大量的遥感数据集并进行分析。哨兵(Sentinel)影像数据就是其中之一,这些数据由欧洲航天局(ESA)的哨兵卫星任务提供,包括哨兵-1(S1)合成孔径雷达(SAR)数据。 要从GEE下载VV+VH偏振的哨兵影像,可以按照以下步骤操作: 1. 访问Google Earth Engine的代码编辑器(https://code.earthengine.google.com/)。 2. 在代码编辑器中,选择Sentinel-1影像集,通常使用以下代码进行筛选: ```javascript var startDate = 'YYYY-MM-DD'; // 设置开始日期 var endDate = 'YYYY-MM-DD'; // 设置结束日期 var sentinel1 = ee.ImageCollection('COPERNICUS/S1_GRD') .filterBounds(geometry) // 使用自己的几何形状或坐标边界来筛选影像 .filter(ee.Filter.date(startDate, endDate)) // 设置日期范围 .filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VV')) // 筛选VV偏振影像 .filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VH')) // 筛选VH偏振影像 .filter(ee.Filter.eq('instrumentMode', 'IW')); // 使用条带式干涉测量模式 ``` 3. 在筛选结果中选择所需的影像,然后使用`.getInfo()`方法或导出工具下载数据。例如: ```javascript var image = sentinel1.first(); // 选择第一个影像 varVV = image.select('VV'); // 选择VV偏振数据 varVH = image.select('VH'); // 选择VH偏振数据 Export.image.toDrive({ // 导出影像到Google Drive image:VV, description: 'VV_image', scale: 10, // 设置影像的分辨率 region: geometry // 设置导出区域 }); Export.image.toDrive({ // 导出VH偏振数据 image:VH, description: 'VH_image', scale: 10, region: geometry }); ``` 4. 点击代码编辑器界面上的“运行”按钮,并在运行后点击“导出”按钮来启动数据导出流程。 请注意,根据影像的大小和GEE使用政策,可能需要一些时间来完成数据的下载。务必确保遵守GEE的使用条款和数据导出限制。
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值