feat: function to convert data into sf - fix #38

NAMESPACE

@@ -4,6 +4,7 @@ export(compute_abundances)
 export(compute_concentration_bins_CPR_n)
 export(compute_concentrations)
 export(compute_frequencies)
+export(data_to_sf)
 export(filter_by_bbox)
 export(filter_by_month)
 export(filter_by_ocean)

R/data_to_sf.R

@@ -0,0 +1,42 @@
+#' Convert a data frame into an sf object
+#'
+#' @description
+#' This function can be used to convert a data frame into an sf object.
+#' 
+#' @param data a `data.frame`, i.e. a FORCIS dataset or the output of a 
+#'   `filter_*()` function.
+#'
+#' @return An `sf POINTS` object.
+#' 
+#' @export
+#'
+#' @examples
+#' ## ADD EXAMPLE ----
+
+data_to_sf <- function(data) {
+
+  ## Check data object ----
+
+  check_if_not_df(data)
+
+  check_field_in_data(data, "site_lon_start_decimal")
+  check_field_in_data(data, "site_lat_start_decimal")
+
+
+  ## Convert data into sf object -----
+
+  data <- data[!is.na(data$"site_lon_start_decimal"), ]
+  data <- data[!is.na(data$"site_lat_start_decimal"), ]
+
+  data_sf <- sf::st_as_sf(data, 
+                          coords = c("site_lon_start_decimal", 
+                                     "site_lat_start_decimal"),
+                          crs    = sf::st_crs(4326))
+
+
+  ## Project spatial objects into Robinson system ----
+
+  data_sf <- sf::st_transform(data_sf, sf::st_crs(crs_robinson()))
+
+  data_sf
+}

R/map_distribution.R

@@ -19,28 +19,7 @@
 
 map_distribution <- function(data, col = "red", ...) {
 
-  ## Check data object ----
-
-  check_if_not_df(data)
-
-  check_field_in_data(data, "site_lon_start_decimal")
-  check_field_in_data(data, "site_lat_start_decimal")
-
-
-  ## Convert data into sf object -----
-
-  data <- data[!is.na(data$"site_lon_start_decimal"), ]
-  data <- data[!is.na(data$"site_lat_start_decimal"), ]
-
-  data_sf <- sf::st_as_sf(data, 
-                          coords = c("site_lon_start_decimal", 
-                                     "site_lat_start_decimal"),
-                          crs    = sf::st_crs(4326))
-
-
-  ## Project spatial objects into Robinson system ----
-
-  data_sf <- sf::st_transform(data_sf, sf::st_crs(crs_robinson()))
+  data_sf <- data_to_sf(data)
 
 
   ## Map ----