ToB企服应用市场:ToB评测及商务社交产业平台

标题: GEE题目：筛选指定地区的Sentinel-1影像缺乏VH等波段 [打印本页]
作者: 诗林    时间: 2024-10-11 22:40
标题: 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企服之家，中国第一个企服评测及商务社交产业平台。



欢迎光临 ToB企服应用市场:ToB评测及商务社交产业平台 (https://dis.qidao123.com/) Powered by Discuz! X3.4