From d46588f4f15310260f2ac30fccc9add1085e54a2 Mon Sep 17 00:00:00 2001 From: Jeremy Kidwell Date: Fri, 8 Feb 2019 09:44:08 +0000 Subject: [PATCH] adding routines --- testing_sf_plot_times.R | 44 ++++++++++++++++++++++++++++++++--------- 1 file changed, 35 insertions(+), 9 deletions(-) diff --git a/testing_sf_plot_times.R b/testing_sf_plot_times.R index b6aa97a..b0541c3 100644 --- a/testing_sf_plot_times.R +++ b/testing_sf_plot_times.R @@ -3,24 +3,50 @@ require(sf) # new simplefeature data class, supercedes sp in many ways # using GEOS 3.6.1, GDAL 2.1.3, PROJ 4.9.3 require(sp) # needed for proj4string, deprecated by sf() require(rgdal) # version version: 1.3-6 +require(devtools) - - -if (file.exists("data/National_Forest_Inventory_Woodland_Scotland_2017.shp") == FALSE) { - download.file("https://opendata.arcgis.com/datasets/3cb1abc185a247a48b9d53e4c4a8be87_0.zip?outSR=%7B%22wkid%22%3A27700%2C%22latestWkid%22%3A27700%7D", - destfile = "data/National_Forest_Inventory_Woodland_Scotland_2017.zip") - unzip("data/National_Forest_Inventory_Woodland_Scotland_2017.zip", exdir = "data") +# Set up local workspace and load data +if (dir.exists("data") == FALSE) { + dir.create("data") } -forest_inventory <- st_read("data/National_Forest_Inventory_Woodland_Scotland_2017.shp") +transition_wgs <- read.csv(text=getURL("https://zenodo.org/record/165519/files/SCCAN_1.4.csv")) +coordinates(transition_wgs) <- c("X", "Y") +proj4string(transition_wgs) <- CRS(wgs84) +transition_sp <- spTransform(transition_wgs, bng) +transition_sf <- st_as_sf(transition, coords = c("X", "Y"), crs=27700) + + +if (file.exists("data/SSSI_SCOTLAND.shp") == FALSE) { +# TODO: get reliable URL for data download +# http://gateway.snh.gov.uk/natural-spaces/dataset.jsp?dsid=SSSI +# download.file("", destfile = "data/SSSI_SCOTLAND_ESRI.zip") +unzip("data/SSSI_SCOTLAND_ESRI.zip", exdir = "data") +} + +if (file.exists("data/National_Forest_Inventory_Woodland_Scotland_2017.shp") == FALSE) { +download.file("https://opendata.arcgis.com/datasets/3cb1abc185a247a48b9d53e4c4a8be87_0.zip", destfile = "data/National_Forest_Inventory_Woodland_Scotland_2017.zip") + +unzip("data/National_Forest_Inventory_Woodland_Scotland_2017.zip", exdir = "data") +} + +forest_inventory_sf <- st_read("data/National_Forest_Inventory_Woodland_Scotland_2017.shp") forest_inventory_sp <- readOGR("./data", "National_Forest_Inventory_Woodland_Scotland_2017") +st_crs(sssi) <- 27700 +st_crs(ecs_sf) <- 27700 + system.time( ggplot() + - geom_sf(data = forest_inventory) + geom_sf(data = forest_inventory_sf) ) system.time( ggplot() + geom_polygon(data = forest_inventory_sp) - ) \ No newline at end of file + ) + +count_data <- sum(apply(st_within(points_sf, st_buffer(sssi, dist = 50), sparse=FALSE), 1, any)) + + +sessioninfo::session_info() \ No newline at end of file