library(rgee)
library(sf)

ee_Initialize(drive = T)  #启动GEE和Google Drive

ee_roi <- read_sf("./SHP/ChinaSizhi.shp") %>%
  sf_as_ee()

copDEM <- ee$ImageCollection("projects/sat-io/open-datasets/GLO-30")

#数据预览
dem_composite = copDEM$
  mosaic()$
  clip(ee_roi)

demParams <- list(
  min = 0,
  max = 6000,
  palette = c(
    '#006600', '#002200', '#fff700', '#ab7634', '#c4d0ff', '#ffffff'
  ))

Map$centerObject(ee_roi)
Map$addLayer(dem_composite, demParams, "dem")+
  Map$addLayer(ee_roi)+
  Map$addLegend(demParams)

demtif <- ee_as_raster(
  image = dem_composite,
  region = ee_roi$geometry(),
  scale = 30,
  dsn = "./copDEM/dem.tif",
  via = "drive",
  maxPixels = 1e13
)

Logo

瓜分20万奖金 获得内推名额 丰厚实物奖励 易参与易上手

更多推荐