【GEE学习】小范围研究区筛选无云影像

问题描述:

常用的sentinel-2筛选少云影像通常使用Filter

.filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20))

和QA波段去云。

function rmCloudByQA(image) { 
  var qa = image.select('QA60'); 
  var cloudBitMask = 1 << 10; 
  var cirrusBitMask = 1 << 11; 
  var mask = qa.bitwiseAnd(cloudBitMask).eq(0) 
               .and(qa.bitwiseAnd(cirrusBitMask).eq(0)); 
  return image.updateMask(mask); 
} 

但是对于小范围研究区来说,这样的筛选会导致“影像云量多而研究区云量少”、“影像云量少但集中在研究区”等情况,造成时序影像利用效率不高。

思路:

本文拟筛选出研究区内无云影像,但GEE代码水平不高,思路如下:
1、利用qa波段将影像的云掩膜为0值;
2、对研究区内影像进行统计,利用Reducer.allNonZero()检索是否有0值;
3、添加一个属性,有云时值为0,无云时值为1;
4、Filter筛选出该值为1的影像。

代码

function rmCloudByQA(image) {
  var qa = image.select('QA60');

  var cloudBitMask = 1 << 10;
  var cirrusBitMask = 1 << 11;
 
  var mask = qa.bitwiseAnd(cloudBitMask).eq(0)
              .and(qa.bitwiseAnd(cirrusBitMask).eq(0));
  
image = image.updateMask(mask).unmask(0)
  
  var min = image.select("B1")
  .reduceRegion({
  	reducer: ee.Reducer.allNonZero(),
  	geometry:roi,
  	scale: 10,
  	maxPixels: 1e13})
  .get("B1");
  
  return image.set({'cloud_musk':min});
}
var S2Collection = ee.ImageCollection('COPERNICUS/S2')
                  .filterDate(start_date, end_date)
                  .filterBounds(table)
                  .map(rmCloudByQA)
                  .filter(ee.Filter.eq('cloud_musk', 1))

函数说明:

GEE:遥感影像处理中的常见掩膜函数的使用
Google Earth Engine(gee)中的reduce、Reducer
GEE学习:遥感影像设置新的属性并查询

  • 1
    点赞
  • 7
    收藏
    觉得还不错? 一键收藏
  • 4
    评论
geeGoogle Earth Engine(谷歌地球引擎)的简称,是一个用于处理和分析地理空间数据的云端平台。它提供了丰富的遥感影像数据和强大的计算能力,能够帮助我们对遥感影像进行裁剪和研究。 在使用gee裁剪遥感影像时,我们需要首先导入需要处理的影像数据。gee平台拥有大量的遥感影像数据集,包括高分辨率的卫星影像和激光雷达数据等。我们可以选择合适的数据集,并通过代码语言(如JavaScript或Python)将影像数据导入到gee平台中。 接下来,我们可以使用gee提供的裁剪函数对遥感影像进行裁剪。裁剪函数可以根据我们提供的裁剪范围(如矩形范围、多边形范围或栅格掩膜),从原始影像中提取我们感兴趣的域。裁剪后,我们可以将裁剪后的影像数据保存到gee平台中,以便后续的分析和研究。 使用gee裁剪遥感影像可以有多种应用。比如,我们可以裁剪特定地影像,用于监测和分析植被覆盖变化、土地利用变化或城市扩展情况等;我们也可以裁剪海洋影像,用于研究海洋生态系统、海洋污染或海洋气候变化等。 总的来说,gee通过提供丰富的遥感影像数据和强大的计算能力,使得裁剪和研究遥感影像变得更加简便和高效。使用gee对遥感影像进行裁剪,可以方便地提取我们感兴趣的域,为后续的研究和分析提供有价值的数据基础。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值