docs: more guidance for big layers

DESCRIPTION

@@ -52,8 +52,8 @@ Config/testthat/edition: 3
 Encoding: UTF-8
 LazyData: true
 Roxygen: list(markdown = TRUE)
-RoxygenNote: 7.3.1
+RoxygenNote: 7.3.2
 SystemRequirements: C++11, GDAL (>= 2.0.1), GEOS (>= 3.4.0), PROJ (>=
     4.8.0)
-Remotes:  
+Remotes: 
     eblondel/ows4R

R/layers.R

@@ -28,7 +28,9 @@
 #' @param reduce_layers whether to reduce output layers to a single `sf` object.
 #' @param ... additional vendor parameter arguments passed to
 #' [`ows4R::GetFeature()`](https://docs.geoserver.org/stable/en/user/services/wfs/reference.html#getfeature).# nolint
-#' For example, including `count=1` returns the first available feature.
+#' For example, including `count = 1` returns the first available feature.
+#' Or `outputFormat = "CSV"` (or `outputFormat = "JSON"`) might help downloading
+#' bigger datasets.
 #' @return If `reduce_layers = FALSE` (default), a list of `sf`
 #' objects, one element for each layer. Any layers for which download was
 #' unsuccessful will be NULL. If `reduce_layers = TRUE`, all layers are
@@ -37,8 +39,7 @@
 #' an `sf` out of a single layer request instead of a list of length 1.
 #' @export
 #'
-#' @examples
-#' \dontrun{
+#' @examplesIf interactive()
 #' # Layers as character vector
 #' emodnet_get_layers(
 #'   service = "biology",
@@ -58,7 +59,27 @@
 #'   layers = "mediseh_zostera_m_pnt",
 #'   count = 1
 #' )
-#' }
+#'
+#' # Usage of csv output
+#' data <- emodnet_get_layers(
+#'     service = "biology_occurrence_data",
+#'     layers = "abiotic_observations",
+#'     outputFormat = "CSV"
+#' )
+#' str(data[["abiotic_observations"]])
+#'
+#' @section Big downloads:
+#'
+# nolint start: line_length_linter
+#' If a layer is really big (like `"abiotic_observations"` of the
+#' `"biology_occurrence_data"` service),
+#' you might consider a combination of these ideas:
+#' - using [`outputFormat = "CSV"`](https://docs.geoserver.org/stable/en/user/services/wfs/reference.html#getfeature);
+#' - filtering using [`cql_filters`](https://emodnet.github.io/EMODnetWFS/articles/ecql_filtering.html) or
+#' [bounding boxes](https://emodnet.github.io/EMODnetWFS/articles/request-params.html#limit-spatial-extent-using-a-boundary-box)
+#' (possibly splitting the area of interests into several requests);
+#' - Using [EMODnet's download toolbox](https://emodnet.ec.europa.eu/geoviewer/).
+# nolint end
 emodnet_get_layers <- function(wfs = NULL,
                                service = NULL,
                                service_version = NULL,
@@ -137,7 +158,7 @@
   # if reduce_layers = T, reduce to single sf --------------------------------
   if (reduce_layers) {
     tryCatch(
       out <- purrr::reduce(out, rbind),
       error = function(e) {
         cli::cli_abort(
           c(