if (FALSE) {
library(rgee)
library(sf)
ee_Initialize()
# Case 1: Geometry*
geom1 <- ee$Geometry$Point(list(-73.53, -15.75))
Map$centerObject(geom1, zoom = 8)
m1 <- Map$addLayer(
eeObject = geom1,
visParams = list(
pointRadius = 10,
color = "FF0000"
),
name = "Geometry-Arequipa"
)
# Case 2: Feature
feature_arq <- ee$Feature(ee$Geometry$Point(list(-72.53, -15.75)))
m2 <- Map$addLayer(
eeObject = feature_arq,
name = "Feature-Arequipa"
)
m2 + m1
# Case 4: Image
image <- ee$Image("LANDSAT/LC08/C01/T1/LC08_044034_20140318")
Map$centerObject(image)
m4 <- Map$addLayer(
eeObject = image,
visParams = list(
bands = c("B4", "B3", "B2"),
max = 10000
),
name = "SF"
)
# Case 5: ImageCollection
nc <- st_read(system.file("shape/nc.shp", package = "sf")) %>%
st_transform(4326) %>%
sf_as_ee()
ee_s2 <- ee$ImageCollection("COPERNICUS/S2")$
filterDate("2016-01-01", "2016-01-31")$
filterBounds(nc)
ee_s2 <- ee$ImageCollection(ee_s2$toList(2))
Map$centerObject(nc$geometry())
m5 <- Map$addLayers(ee_s2)
m5
# Case 6: Map comparison
image <- ee$Image("LANDSAT/LC08/C01/T1/LC08_044034_20140318")
Map$centerObject(image)
m_ndvi <- Map$addLayer(
eeObject = image$normalizedDifference(list("B5", "B4")),
visParams = list(min = 0, max = 0.7),
name = "SF_NDVI"
) + Map$addLegend(list(min = 0, max = 0.7), name = "NDVI", position = "bottomright", bins = 4)
m6 <- m4 | m_ndvi
m6
# Case 7: digging up the metadata
m6$rgee$tokens
m5$rgee$tokens
# Case 8: COG support
# See parameters here: https://api.cogeo.xyz/docs
server <- "https://storage.googleapis.com/pdd-stac/disasters/"
file <- "hurricane-harvey/0831/20170831_172754_101c_3B_AnalyticMS.tif"
resource <- paste0(server, file)
visParams <- list(bands = c("B3", "B2", "B1"), min = 3000, max = 13500, nodata = 0)
Map$centerObject(resource)
Map$addLayer(resource, visParams = visParams, shown = TRUE)
}
Run the code above in your browser using DataLab