inst/doc/openeo-02-sample_data.R

In openeo: Client Interface for 'openEO' Servers

## ----setup,echo=FALSE, eval=FALSE---------------------------------------------
#  library(openeo)

## ---- eval=FALSE--------------------------------------------------------------
#  con = connect("https://openeo.cloud")
#  login()

## ---- eval=FALSE--------------------------------------------------------------
#  p = processes()
#  coll = list_collections()
#  f = list_file_formats()
#  
#  bounding_box = list(west=6.75,south=51.85,east=7.25,north=52.15)
#  
#  data = p$load_collection(id = coll$SENTINEL2_L2A,
#                           bands=c("B04","B08"),
#                           spatial_extent = bounding_box,
#                           temporal_extent = list("2021-03-01","2021-07-15"))
#  
#  ndvi = p$reduce_dimension(data=data, dimension = "bands",reducer = function(x, context) {
#    b04 = x[1]
#    b08 = x[2]
#    (b08-b04)/(b08+b04)
#  })

## ---- eval=FALSE--------------------------------------------------------------
#  res = p$save_result(ndvi, format = f$output$netCDF)

## ---- eval=FALSE--------------------------------------------------------------
#  as(res,"Process")

## ---- eval=FALSE--------------------------------------------------------------
#  center = list(lon=mean(c(bounding_box$west, bounding_box$east)), lat=mean(c(bounding_box$south, bounding_box$north)))
#  diff = 0.0003
#  sample_bbox = list(west=center$lon-diff,south=center$lat-diff,east=center$lon+diff,north=center$lat+diff)

## ---- eval=FALSE--------------------------------------------------------------
#  data$parameters$spatial_extent = sample_bbox

## ---- eval=FALSE--------------------------------------------------------------
#  data

## ---- eval=FALSE--------------------------------------------------------------
#  manual_file = "test_manual.nc"
#  file = compute_result(res,output_file = manual_file)

## ---- eval=FALSE,results='asis', warning=FALSE--------------------------------
#  library(stars)
#  obj = read_stars(manual_file,driver = detect.driver(manual_file))
#  obj

## ---- eval=FALSE, warning=FALSE-----------------------------------------------
#  library(lubridate)
#  library(stars)
#  library(ggplot2)
#  
#  dates = as_date(st_get_dimension_values(obj,"t"))
#  obj = st_set_dimensions(obj,which="t",values = dates)
#  obj

## ---- eval=FALSE,out.height='400px',out.width='600px'-------------------------
#  ggplot() + geom_stars(data=obj) + facet_wrap(~t) + coord_equal() + theme_void()

## ---- eval=FALSE--------------------------------------------------------------
#  data$parameters$spatial_extent = bounding_box
#  data

## ---- eval=FALSE, warning=FALSE-----------------------------------------------
#  library(mapview)
#  library(sf)
#  
#  initial_bbox = as.numeric(bounding_box)
#  names(initial_bbox) = c("xmin","ymin","xmax","ymax")
#  sample_bbox = as.numeric(sample_bbox)
#  names(sample_bbox) = c("xmin","ymin","xmax","ymax")
#  
#  initial_bbox = st_as_sfc(st_bbox(initial_bbox,crs=4326))
#  sample_bbox = st_as_sfc(st_bbox(sample_bbox,crs=4326))
#  
#  received_data_bbox = st_as_sfc(st_bbox(obj))
#  
#  mapview(list(initial_bbox,sample_bbox,received_data_bbox))

## ---- eval=FALSE--------------------------------------------------------------
#  ?get_sample

## ---- eval=FALSE--------------------------------------------------------------
#  ?compute_result

## ---- eval=FALSE--------------------------------------------------------------
#  f = list_file_formats()

## ---- eval=FALSE, warning=FALSE-----------------------------------------------
#  sample_file = "test_get_sample.nc"
#  obj = get_sample(ndvi,
#                   as_stars=TRUE,
#                   output_file=sample_file,
#                   format = f$output$netCDF)

## ---- eval=FALSE--------------------------------------------------------------
#  obj

## ---- eval=FALSE--------------------------------------------------------------
#  library(lubridate)
#  library(stars)
#  library(ggplot2)
#  
#  dates = as_date(st_get_dimension_values(obj,"t"))
#  obj = st_set_dimensions(obj,which="t",values = dates)

## ---- eval=FALSE, out.height="400px",out.width="600px"------------------------
#  ggplot() + geom_stars(data=obj) + facet_wrap(~t) + coord_equal() + theme_void()

## ---- eval=FALSE--------------------------------------------------------------
#  library(mapview)
#  library(sf)
#  
#  received_data_bbox = st_as_sfc(st_bbox(obj))
#  
#  mapview(list(initial_bbox,received_data_bbox))