从GEE上提取波段信息到CSV文件

最新推荐文章于 2024-07-05 10:32:37 发布

一个金牛座的GISer

最新推荐文章于 2024-07-05 10:32:37 发布

阅读量195

点赞数 2

分类专栏： GEE 文章标签：云计算

本文链接：https://blog.csdn.net/weixin_44586490/article/details/137041734

版权

GEE 专栏收录该内容

1 篇文章 0 订阅

订阅专栏

GEE上可以提取遥感影像的特征到CSV文件，方便后面数据分析和数据清洗，然后使用机器学习算法处理


```javascript
// 定义时间范围和区域（使用与Sentinel-2相同的时间范围和区域）
var startDate = '2021-10-01';
var endDate = '2021-12-31';

// 加载Sentinel-1A数据
var sentinel1 = ee.ImageCollection('COPERNICUS/S1_GRD')
  .filterDate(startDate, endDate)
  .filter(ee.Filter.eq('instrumentMode', 'IW'))  // Interferometric Wide swath mode
  .filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VH')) // 确保数据包含VH极化
  .select('VH')  // 选择VH波段
  .median();  // 创建一个中值合成影像

// Make a cloud-free Landsat 8 surface reflectance composite.
var image = ee.ImageCollection('COPERNICUS/S2')
  .filterDate('2021-10-01', '2021-12-31')
  .median();
  
// 将Sentinel-1的VH波段添加到Sentinel-2的影像中
var combinedImage = image.addBands(sentinel1);  

// 创建NDVI
var ndvi = combinedImage.normalizedDifference(['B8', 'B4']).rename('NDVI');
var ndwi = combinedImage.normalizedDifference(['B3', 'B8']).rename('NDWI');
  
// 将NDVI和NDWI波段添加到组合影像中
combinedImage = combinedImage.addBands([ndvi, ndwi]);

// 选择需要的波段
var bands = ['NDVI', 'NDWI', 'B11', 'B12', 'VH'];

combinedImage = combinedImage.select(bands);

var pointsDataset = ee.FeatureCollection('projects/ee-zhangan28886666/assets/jiangsu');
var pointsDateset2 = ee.FeatureCollection('projects/ee-zhangan28886666/assets/jiangsu-n')
var combineDataset = pointsDataset.merge(pointsDateset2);
//var pointsDataset = ee.FeatureCollection('projects/ee-zhangan28886666/assets/train4');
//var pointsDataset = ee.FeatureCollection('projects/ee-zhangan28886666/assets/training_2020_autumn');

/*var reprojectedDataset = pointsDataset.map(function(feature) {
  return feature.setGeometry(feature.geometry().transform('EPSG:32650'));
});*/
// 选择需要的波段
var selectedImage = combinedImage.select(bands);

// 对每个点应用reduceRegion并提取值
var pointValues = combineDataset.map(function(feature) {
  return feature.set(selectedImage.reduceRegion({
    reducer: ee.Reducer.mean(), // 使用first获取该位置的值
    geometry: feature.geometry(),
    scale: 10 // Sentinel-2的分辨率通常是10米至20米不等，根据需要选择
  }));
});

// 打印输出结果
//print(pointValues);
// 导出结果

Export.table.toDrive({
  collection: pointValues,
  description: 'exported_point_values',
  fileFormat: 'CSV'
});