inst/extdata/processing.R

In jacob-ogre/ESAListings: An R Data package for U.S. Endangered Species Act Listed Species

library(ecosscraper)
library(geojsonio)
library(maptools)
library(rgdal)
library(leaflet)

#format and join lat and lon from county_atrib csv file
#county_attrib <- read.csv("inst/extdata/US_counties_attrib.csv")
#county_attrib$GEOID <- as.character(county_attrib$GEOID)
#county_attrib$GEOID <- ifelse(nchar(county_attrib$GEOID) == 4,
#                             paste0("0", county_attrib$GEOID),
#                             county_attrib$GEOID)

source("C:/Users/mevans/repos/ecosscraper/R/options.R")
TECP_date <- get_TECP_table()

county_map <- readShapeSpatial("C:/Users/mevans/GIS/cb_2015_us_county_500k.shp")

#join species counts by county table to dataframe of shapefile
county_map@data <- left_join(county_map@data, counties, by = "GEOID")


county_map@data$style <- pal(county_map@data$count)

test_geo <- geojson_json(county_map)
test_geo <- geojson_read("county_geo.json")

#read topoJSON in list format
test_topo <- jsonlite::fromJSON(readLines("county_topo3.json") %>%
                                  paste(collapse="\n"))

#read topoJSON in character vector format
test_topo <- readLines("county_topo_rb76.json") %>% paste(collapse="\n")

test_geo$features <- lapply(test_geo$features, function(x) {
  x$properties$style <- list(fillColor = palfx(x$properties$count))
  x})

geojson_out <- jsonlite::toJSON(test_geo, pretty = TRUE)
write(geojson_out, file="county_geo2.json")


cur_map <- leaflet() %>%
  addTopoJSON(topo_prac, weight = 0.75, color = "black", fillOpacity = 0.5) %>%
  addProviderTiles("Stamen.TonerLite") %>%
  setView(lng=-95, lat=38, zoom = 4)  %>%
  mapOptions(zoomToLimits = "never") %>%
  addLegend(title = "Number of<br>Listed Species",
            position = "bottomleft",
            pal = palfx,
            values = seq(0, 80, 10))