nationalparkservice · RobLBaker · Feb 28, 2023 · Feb 8, 2023 · Feb 8, 2023 · Feb 8, 2023
diff --git a/DESCRIPTION b/DESCRIPTION
@@ -1,7 +1,7 @@
 Package: NPSutils
 Type: Package
 Title: Collection of Functions to read/write information from the NPS Data Store
-Version: 0.2.0.0
+Version: 0.2.0.3
 Authors@R: c(
     person(given = "Robert", family = "Baker", email = "robert_baker@nps.gov",
            role = c("aut", "cre"),
@@ -29,8 +29,9 @@ Imports:
     readr,
     magrittr,
     crayon,
-    stringr
-RoxygenNote: 7.2.1
+    stringr,
+    leaflet
+RoxygenNote: 7.2.3
 Suggests: 
     knitr,
     rmarkdown

diff --git a/NAMESPACE b/NAMESPACE
@@ -13,5 +13,6 @@ export(loadDataPackages)
 export(load_data_package)
 export(load_domains)
 export(load_metadata)
+export(map_wkt)
 export(validate_data_package)
 importFrom(magrittr,"%>%")
diff --git a/NEWS.md b/NEWS.md
@@ -1,3 +1,6 @@
+# NPSutils 0.2.0.3
+  * added `map_wkt()` function to map points, polygons, or both from Well Known Text coordinates (WKT). WKT is used in place to GPS coordinates when sensitive species locations have been "fuzzed". In this case, providing a polygon rather than the an exact (albeit fuzzed) is preferable as it is clear that the location is not exact. WKT is an efficient way to store geographic shapes such as polygons in flat files such as .csv.
+
 # NPSutils 0.2.0.2
 
   * deprecated `get_data_package()` in favor of `getDataPackages()`. This is just a name change and does not change any functionality.

diff --git a/R/map_wkt.R b/R/map_wkt.R
@@ -0,0 +1,68 @@
+#' Map WKT geometry (points and polygons)
+#'
+#' @description `map_wkt()` takes a well-known text (WKT) geometry column and maps points and polygons onto a gray leaflet map. All NA geometry is dropped before mapping.
+#'
+#' @details Define your dataframe, the column that contains WKT, and an option to map specific geometry types.
+#'
+#' @param df - The name of the data frame that contains WKT geometry.
+#' @param wellknowntext - The name of the specific column within the data frame that contains the WKT geometry. This parameter is currently not fully implemented and defaults to the Darwin Core 'footprintsWKT'.
+#' @param type -  Pick one from "points", "polygons", or "all" to map specific geometry types.
+#'
+#' @return The function returns a dynamic, zoomable leaflet map with the specific geometry plotted.
+#' @importFrom magrittr %>%
+#' @export
+#'
+#' @examples
+#' \dontrun{
+#' #map species observations
+#' map_wkt(my_NPS_species_obs)
+#' 
+#' #map species observations - points only
+#' map_wkt(my_NPS_species_obs, wellknowntext = "footprintWKT", type = "points")
+#' }
+map_wkt <- function(df, wellknowntext = "footprintWKT", type = "all") {
+  #filter to just wellknowntext column:
+  wkt_grepl <- paste0('\\b', wellknowntext, '\\b')
+  df <- df[grepl(wkt_grepl, colnames(df))]
+  #omit NAs - not important for plotting anyway:
+  df <- na.omit(df)
+  #convert to geographic object:
+  df <- sf::st_as_sf(df, wkt = wellknowntext)
+  #new column in data frame for the geometry type
+  df$geometry_types <- sf::st_geometry_type(df)
+  #use the geometry_type column to filter only for POINT or POLYGON
+  df_pts <- df[df$geometry_types == "POINT",]
+  df_polys <- df[df$geometry_types == "POLYGON",]
+
+  #only map what is requested
+  if(type == "points") {
+    map <- leaflet::leaflet(df, 
+              options = leafletOptions(preferCanvas = TRUE)) %>%
+      #addTiles(group = "OSM (default)"); prevent unwanted map updates:
+              leaflet::addProviderTiles(providers$Esri.WorldGrayCanvas, 
+                                        options = providerTileOptions(
+                                          updateWhenZooming = FALSE,
+                                          updateWhenIdle = TRUE)) %>%
+              leaflet::addCircles(data = df_pts, color = "blue",) #odd stray ,
+  } else if(type == "polygons") {
+    map <- leaflet::leaflet(df, 
+              options = leafletOptions(preferCanvas = TRUE)) %>%
+      #addTiles(group = "OSM (default)"); prevent unwanted map updates:
+              leaflet::addProviderTiles(providers$Esri.WorldGrayCanvas, 
+                                        options = providerTileOptions(
+                                          updateWhenZooming = FALSE,
+                                          updateWhenIdle = TRUE)) %>%
+              leaflet::addPolygons(data = df_polys, color = "red",) #odd stray ,
+  } else if(type == "all") {
+    map <- leaflet::leaflet(df, 
+              options = leafletOptions(preferCanvas = TRUE)) %>%
+      #addTiles(group = "OSM (default)") %>%; prevent unwatned map updates:
+              leaflet::addProviderTiles(providers$Esri.WorldGrayCanvas,
+                                        options = providerTileOptions(
+                                          updateWhenZooming = FALSE,
+                                          updateWhenIdle = TRUE)) %>%
+              leaflet::addCircles(data = df_pts, color = "blue",) %>%
+              leaflet::addPolygons(data = df, color = "red",)  #odd stray ,
+  }
+  return(map)
+}
diff --git a/docs/404.html b/docs/404.html
diff --git a/docs/LICENSE-text.html b/docs/LICENSE-text.html
diff --git a/docs/LICENSE.html b/docs/LICENSE.html
diff --git a/docs/articles/NPSutils.html b/docs/articles/NPSutils.html
diff --git a/docs/articles/index.html b/docs/articles/index.html
diff --git a/docs/authors.html b/docs/authors.html
diff --git a/docs/index.html b/docs/index.html
diff --git a/docs/news/index.html b/docs/news/index.html
diff --git a/docs/pkgdown.yml b/docs/pkgdown.yml
@@ -1,7 +1,6 @@
-pandoc: 2.17.1.1
+pandoc: 2.19.2
 pkgdown: 2.0.6.9000
 pkgdown_sha: 5bb671b1a8bac8e926479868681936f7c9ee41b3
 articles:
   NPSutils: NPSutils.html
-last_built: 2023-02-08T18:49Z
-
+last_built: 2023-02-28T17:14Z
diff --git a/docs/reference/getDataPackages.html b/docs/reference/getDataPackages.html
diff --git a/docs/reference/get_park_code.html b/docs/reference/get_park_code.html
diff --git a/docs/reference/get_park_taxon_citations.html b/docs/reference/get_park_taxon_citations.html
diff --git a/docs/reference/get_park_taxon_refs.html b/docs/reference/get_park_taxon_refs.html
diff --git a/docs/reference/get_park_taxon_url.html b/docs/reference/get_park_taxon_url.html
diff --git a/docs/reference/get_ref_info.html b/docs/reference/get_ref_info.html
diff --git a/docs/reference/get_unit_code.html b/docs/reference/get_unit_code.html
diff --git a/docs/reference/get_unit_code_info.html b/docs/reference/get_unit_code_info.html
diff --git a/docs/reference/get_unit_info.html b/docs/reference/get_unit_info.html
diff --git a/docs/reference/index.html b/docs/reference/index.html
diff --git a/docs/reference/loadDataPackages.html b/docs/reference/loadDataPackages.html
diff --git a/docs/reference/load_data_package.html b/docs/reference/load_data_package.html
diff --git a/docs/reference/load_domains.html b/docs/reference/load_domains.html
diff --git a/docs/reference/load_metadata.html b/docs/reference/load_metadata.html