# NOT RUN {
library(mapview)
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"
)
m4
# Case 5: mapview + EarthEnginemap
library(sf)
nc <- st_read(system.file("shp/arequipa.shp", package="rgee"))
mapview(nc) + m2
m2 + mapview(nc)
# Case 6: mapedit
library(mapedit)
# my_geometry <- m2 %>% ee_as_mapview() %>% editMap()
# Case 7: 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_get(0:4)
Map$centerObject(nc$geometry())
Map$addLayers(ee_s2)
# }
Run the code above in your browser using DataLab