-
Notifications
You must be signed in to change notification settings - Fork 77
/
api_source_hls.R
84 lines (84 loc) · 3.6 KB
/
api_source_hls.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
#' @title Create an items object in an HLS cube
#' @keywords internal
#' @noRd
#' @description \code{.source_items_new()} this function is called to create
#' an items object. In case of Web services, this function is responsible for
#' making the Web requests to the server.
#' @param source Name of the STAC provider.
#' @param ... Other parameters to be passed for specific types.
#' @param collection Collection to be searched in the data source.
#' @param stac_query Query that follows the STAC protocol
#' @param tiles Selected tiles (optional)
#' @return An object referring the images of a sits cube.
#' @export
.source_items_new.hls_cube <- function(source, ...,
collection,
stac_query,
tiles = NULL) {
# set caller to show in errors
.check_set_caller(".source_items_new_hls_cube")
# check netrc file
suppressWarnings(.check_netrc_gdal(attributes = .conf("HLS_ACCESS_URL")))
# convert tiles to a valid STAC query
if (!is.null(tiles)) {
roi <- .s2_mgrs_to_roi(tiles)
stac_query[["params"]][["intersects"]] <- NULL
stac_query[["params"]][["bbox"]] <- c(roi[["lon_min"]],
roi[["lat_min"]],
roi[["lon_max"]],
roi[["lat_max"]]
)
} else {
# Convert roi to bbox
lon <- stac_query[["params"]][["intersects"]][["coordinates"]][, , 1]
lat <- stac_query[["params"]][["intersects"]][["coordinates"]][, , 2]
stac_query[["params"]][["intersects"]] <- NULL
stac_query[["params"]][["bbox"]] <- c(min(lon),
min(lat),
max(lon),
max(lat)
)
}
# making the request
items_info <- rstac::post_request(q = stac_query, ...)
.check_stac_items(items_info)
# if more than 2 times items pagination are found the progress bar
# is displayed
progress <- rstac::items_matched(items_info) > 2 *
.conf("rstac_pagination_limit")
# check documentation mode
progress <- .check_documentation(progress)
# fetching all the metadata and updating to upper case instruments
items_info <- rstac::items_fetch(items = items_info, progress = progress)
# checks if the items returned any items
.check_stac_items(items_info)
return(items_info)
}
#' @title Organizes items by tiles for HLS collections
#' @param source Name of the STAC provider.
#' @param ... Other parameters to be passed for specific types.
#' @param items \code{STACItemcollection} object from rstac package.
#' @param collection Collection to be searched in the data source.
#' @return A list of items.
#' @keywords internal
#' @noRd
#' @export
.source_items_tile.hls_cube <- function(source, ...,
items,
collection = NULL) {
tiles <- strsplit(rstac::items_reap(items, field = "id"), "\\.")
tiles <- purrr::map_chr(tiles, function(x) x[[3]])
substr(tiles, 2, 6)
}
#' @noRd
#' @title Configure access.
#' @param source Data source
#' @param collection Image collection
#' @return Called for side effects
.source_configure_access.hls_cube <- function(source, collection = NULL) {
.check_set_caller(".source_configure_access_hls_cube")
# check netrc file
.check_netrc_gdal(attributes = .conf("HLS_ACCESS_URL"))
# done!
return(invisible(source))
}