GEE题目：筛选指定地区的Sentinel-1影像缺乏VH等波段

目录
题目简介
原始代码
原始代码
题目解析
修改后的代码

---

题目简介

亲爱的同事们，我正在尝试使用 SAR 图像 - Sentinel-1 来改进使用呆板学习的地上生物量建模。我想处置惩罚 Sentinel 图像并将它们作为波段插入以增强模子。通过阅读文档，可用的极化（已在 GRD 中处置惩罚）是 HH、HV、VV 和 VH。但是，出现一条消息，指出只有 VV 波段可用。有人知道其他极化是否不可用吗？
原始代码

https://code.earthengine.google.com/00d394a0be6a840eb4a260a78cabca7e
原始代码

1. var table = ee.FeatureCollection("projects/ee-selper/assets/aoi_sample");
2. // Define period of analysis
3. var start = '2021-01-01';
4. var end = '2021-12-31';
5. var season = ee.Filter.date(start,end);
6. print(season);
9. // Import Sentinel-1 collection
10. var sentinel1 =  ee.ImageCollection('COPERNICUS/S1_GRD');
11. // Filter Sentinel-1 collection for study area, date ranges and polarization components
12. var sCollection =  sentinel1
13.                     //filter by aoi and time
14.                     .filterBounds(table)
15.                     .filter(season)
16.                     // Filter to get images with VV and VH dual polarization
17.                     .filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VV'))
18.                     .filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VH'))
19.                     // Filter to get images collected in interferometric wide swath mode.
20.                     .filter(ee.Filter.eq('instrumentMode', 'IW'));
22. var desc = sCollection.filter(ee.Filter.eq('orbitProperties_pass', 'DESCENDING'));
23. var asc = sCollection.filter(ee.Filter.eq('orbitProperties_pass', 'ASCENDING'));
25. // Contar o número de imagens em cada coleção
26. print('Número de Imagens Ascendente (VH)', asc.filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VH')).size());
27. print('Número de Imagens Ascendente (VV)', asc.filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VV')).size());
28. print('Número de Imagens Descendente (VH)', desc.filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VH')).size());
29. print('Número de Imagens Descendente (VH)', desc.filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VV')).size());
30. // Criar um compósito a partir da média das polarizações e modos de órbita
31. var vh_asc = asc.select('VH').mean().rename('VH_asc');
32. var vv_asc = asc.select('VV').mean().rename('VV_asc');
33. var vh_desc = desc.select('VH').mean().rename('VH_desc');
34. var vv_desc = desc.select('VV').mean().rename('VV_desc');
37. // Combinar as bandas em uma única imagem
38. var composite = ee.Image.cat([vh_desc]).focal_median();
40. // Definir parâmetros de visualização para compósito (usar valores individuais para min e max)
41. var visParams = {
42.   bands: ['VH_desc'],  // bandas para os canais R, G, B
43.   min: -25,  // valor mínimo para todas as bandas
44.   max: 0     // valor máximo para todas as bandas
45. };
47. // Exibir no mapa
48. Map.centerObject(table, 10);  // Centralizar o mapa na área de interesse
49. Map.addLayer(composite.clip(table), visParams, 'composite');
51. // Exportar a imagem para o asset
52. Export.image.toAsset({
53.   image: composite.clip(table), // Recortar para a área de interesse
54.   description: 'Sentinel1', // Nome da tarefa de exportação
55.   assetId: 'projects/ee-selper/assets/Sentinel1', // Substitua pelo caminho correto do seu asset
56.   scale: 30, // Resolução da imagem em metros
57.   region: table, // Área de interesse para recorte
58.   maxPixels: 1e13 // Número máximo de pixels permitidos (ajuste conforme necessário)
59. });

复制代码

题目解析

Sentinel-1卫星是由欧洲空间局（ESA）运行的合成孔径雷达（SAR）卫星系统。它们以天基的方式提供大量的雷达图像数据，用于监测地球表面的变化和活动。
在Sentinel-1数据中，地区只有Descending（下行）或只有Ascending（上行）数据是由于卫星系统的运行模式和地球自转的缘故原由。
在多数情况下，Sentinel-1卫星会同时网络Descending和Ascending数据覆盖一个地区。如许做可以提供更全面的地面覆盖，减少数据缺失和提高数据质量。同时，Descending和Ascending数据还可以用于获取更准确的地表变化信息，比方地表沉降等。
然而，由于技能和资源限制，有些地区大概只有一种数据可用。这大概是因为在某些时间地区内，只有Descending或只有Ascending数据网络，这取决于卫星任务安排和优先级。此外，在一些地区，由于地形、地貌等因素，只有一种数据方向可以有效地传输信号，而另一种方向大概碰到干扰或信号丧失。
总结起来，Sentinel-1数据只有Descending或只有Ascending覆盖某个地区是由于卫星运行模式、地球自转和地理因素等多种缘故原由造成的。这些数据的不同方向可以提供更全面和准确的地表观测，并支持各种应用，如环境监测、农业、都会规划等。
修改后的代码

1. var table = ee.FeatureCollection("projects/ee-selper/assets/aoi_sample");
2. // Define period of analysis
3. var start = '2021-01-01';
4. var end = '2021-12-31';
5. var season = ee.Filter.date(start,end);
6. print(season);
9. // Import Sentinel-1 collection
10. var sentinel1 =  ee.ImageCollection('COPERNICUS/S1_GRD');
11. // Filter Sentinel-1 collection for study area, date ranges and polarization components
12. var sCollection =  sentinel1
13.                     //filter by aoi and time
14.                     .filterBounds(table)
15.                     .filter(season)
16.                     // Filter to get images with VV and VH dual polarization
17.                     .filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VH'))
18.                     // Filter to get images collected in interferometric wide swath mode.
19.                     .filter(ee.Filter.eq('instrumentMode', 'IW'));
21. var desc = sCollection.filter(ee.Filter.eq('orbitProperties_pass', 'DESCENDING'));
22. var asc = sCollection.filter(ee.Filter.eq('orbitProperties_pass', 'ASCENDING'));
24. // Contar o número de imagens em cada coleção
25. print('Número de Imagens Ascendente (VH)', asc.filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VH')).size());
26. print('Número de Imagens Ascendente (VV)', asc.filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VV')).size());
27. print('Número de Imagens Descendente (VH)', desc.filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VH')).size());
28. print('Número de Imagens Descendente (VH)', desc.filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VV')).size());
29. // Criar um compósito a partir da média das polarizações e modos de órbita
30. var vh_asc = asc.select('VH').mean().rename('VH_asc');
31. var vv_asc = asc.select('VV').mean().rename('VV_asc');
32. var vh_desc = desc.select('VH').mean().rename('VH_desc');
33. var vv_desc = desc.select('VV').mean().rename('VV_desc');
36. // Combinar as bandas em uma única imagem
37. var composite = ee.Image.cat([vh_desc]).focal_median();
39. // Definir parâmetros de visualização para compósito (usar valores individuais para min e max)
40. var visParams = {
41.   bands: ['VH_desc'],  // bandas para os canais R, G, B
42.   min: -25,  // valor mínimo para todas as bandas
43.   max: 0     // valor máximo para todas as bandas
44. };
46. // Exibir no mapa
47. Map.centerObject(table, 10);  // Centralizar o mapa na área de interesse
48. Map.addLayer(composite.clip(table), visParams, 'composite');
50. // Exportar a imagem para o asset
51. Export.image.toAsset({
52.   image: composite.clip(table), // Recortar para a área de interesse
53.   description: 'Sentinel1', // Nome da tarefa de exportação
54.   assetId: 'projects/ee-selper/assets/Sentinel1', // Substitua pelo caminho correto do seu asset
55.   scale: 30, // Resolução da imagem em metros
56.   region: table, // Área de interesse para recorte
57.   maxPixels: 1e13 // Número máximo de pixels permitidos (ajuste conforme necessário)
58. });

复制代码

免责声明：如果侵犯了您的权益，请联系站长，我们会及时删除侵权内容，谢谢合作！更多信息从访问主页：qidao123.com:ToB企服之家，中国第一个企服评测及商务社交产业平台。