使用rgee快速获取河南洪涝地区Sentinel-2 10m分辨率影像

Sentinel-2是目前公开的最高分辨率的数据，而且有着最短5天的重访周期，对于灾后数据应急有着非常重要的作用。

前面的推文中讲过Sentinel-2数据下载和处理的方法，但是数据需要从哨兵分发网站下载到本地，速度慢，拼接繁琐，在这里提供一种在Google Earth Engine支持下，使用rgee快速获取Sentinel-2影像的方法。

rgee快速获取Sentinel-2 10米分辨率数据

• Sentinel-2十米的波段有四个，B8/B4/B3/B2
• 为了快速获取灾区影像，使用了Sentinel-2 L1C数据
• 由于洪涝后仍有大量云覆盖，云筛选没有生效

library(rgee)
library(sf)

ee_Initialize(drive = TRUE)

# Define a region of interest with sf
ee_roi <-  read_sf("D:/Henan/SHP/XX.shp")%>%  #选择研究区SHP
  sf_as_ee()

s2_sr <- ee$ImageCollection('COPERNICUS/S2')$filterBounds(ee_roi)$
  filterDate('2021-07-25', '2021-07-27')#$   时间筛选
 # filter(ee$Filter$lt('CLOUD_COVERAGE_ASSESSMENT', 10))   #云量筛选

mosaic = s2_sr$select(c('B8','B4', 'B3', 'B2'))$mosaic()$clip(ee_roi)   #选取四个10米分辨率的波段

#预览数据
imageVisParam <- list(bands = c('B8', 'B4', 'B3'), min = 300, max = 2000)
Map$addLayer(mosaic, imageVisParam, 'mosaic')+
  Map$addLayer(ee_roi, visParams = list(palette = "red"), "roi")
Map$centerObject(ee_roi)

#数据输出
s2r <- ee_as_raster(
  image = mosaic,
  region = ee_roi$geometry(),
  scale = 10,
  dsn = './S2r/s2b.tif',
  via = "drive"
)

郑州Sentinel-2影像20210726，标准假彩色合成

新乡Sentinel-2影像20210726，标准假彩色合成

• 可以查看影像数据集时间

ee_get_date_ic(s2_sr)

查看影像时间

郑州和新乡洪涝灾后数据下载

分辨率：10米

波段：B8/B4/B3/B2，包含近红外和红绿蓝四个波段

成像时间：郑州2021年7月24日和26日

                    新乡2021年7月26日

数据下载链接：https://pan.baidu.com/s/1rDZU47fj477tsr_PLpN_2A 提取码：56ua

参考文献