GEE 统计影像覆盖次数

就是一只白

已于 2023-08-08 17:58:46 修改

阅读量273

点赞数

分类专栏： GEE 文章标签： GEE

于 2023-08-08 11:42:59 首次发布

本文链接：https://blog.csdn.net/xinemachen/article/details/132163653

版权

GEE 专栏收录该内容

2 篇文章 0 订阅

订阅专栏

本文介绍了如何使用GoogleEarthEngine(GEE)对Sentinel-1卫星数据进行月度影像覆盖次数的统计，通过设置ROI(研究区域)和时间过滤器，实现数据筛选并计算HH极化模式的图像数量。

摘要由CSDN通过智能技术生成

以下为利用GEE统计Sentinel-1影像覆盖次数的代码，统计了每月的情况。

注意roi的赋值geometry5是之前就提供的范围变量，可以在代码运行之前就手动编辑。

可以采用

ee.Geometry.BBox(aoi['xmin'], aoi['ymin'], aoi['xmax'], aoi['ymax']),

来赋值，但是我发现如果我用这个下面的Map.addLayer(roi, {'color':'grey'}, 'studyArea');就报错，当鼠标放在上方代码，就可以把它转为一个全局？变量，就是像你自己编辑好的那样，下面的Map之类的我也去掉了，否则还是会报错。

另外我希望只在纬度上控制因此应该是：-180，60，180，90的写法。



// set study area table 就是roi
var roi = geometry5;

// Map.addLayer(roi, {'color':'grey'}, 'studyArea');
// Map.centerObject(roi,2);//可以改变缩放等级

var year_start = 2018; //设定年份
var year_end = 2018;
var month=3;

var imgCol = ee.ImageCollection("COPERNICUS/S1_GRD")
              .filterBounds(roi)
              .filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'HH'))
              .filter(ee.Filter.eq('instrumentMode', 'EW'))
              .filter(ee.Filter.calendarRange(year_start, year_end,'year')).filter(ee.Filter.calendarRange(month,month,'month'));

// compute the image total size
// var imgSize = imgCol.size();
// print("imgSize",imgSize);

/*****************************************************************
compute the image avaibility. There are several ways to achieve that.
1: use the count function
2: use the sum function
******************************************************************/
var imgCount_V1 = imgCol.select("HH").count().clip(roi);


// define image render parameters
var visParam = {
 min: 0,
 max: 10,
 palette: '#00c82d, #15cc2d, #95e11d, #d4d900, #ffae18, #ff6a30, #f53031, #f53031'};

Map.addLayer(imgCount_V1,visParam,"imgCount_V1");

Export.image.toDrive({
    image: imgCount_V1,
    description: "S1_201803_imagecount",
    fileNamePrefix: "S1_201803_imagecount",
    region: imgCount_V1.geometry().getInfo(),
    crs: "EPSG:3413",   //坐标系参数，32650表示WGS_1984_UTM_Zone_50N
    scale: 1000,
    maxPixels: 1e13
});