以下为利用GEE统计Sentinel-1影像覆盖次数的代码,统计了每月的情况。
注意roi的赋值geometry5是之前就提供的范围变量,可以在代码运行之前就手动编辑。
可以采用
ee.Geometry.BBox(aoi['xmin'], aoi['ymin'], aoi['xmax'], aoi['ymax']),
来赋值,但是我发现如果我用这个下面的Map.addLayer(roi, {'color':'grey'}, 'studyArea');就报错,当鼠标放在上方代码,就可以把它转为一个全局?变量,就是像你自己编辑好的那样,下面的Map之类的我也去掉了,否则还是会报错。文章来源:https://www.toymoban.com/news/detail-632596.html
另外我希望只在纬度上控制因此应该是:-180,60,180,90的写法。文章来源地址https://www.toymoban.com/news/detail-632596.html
// 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
});
到了这里,关于GEE 统计影像覆盖次数的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!