kidwellj/mapping_environmental_action
1
0
Fork 0
mirror of https://github.com/kidwellj/mapping_environmental_action.git synced 2024-10-31 23:42:20 +00:00
Code Issues Packages Projects Releases Wiki Activity

adding routines

Browse source
This commit is contained in:
Jeremy Kidwell 2019-02-08 09:44:08 +00:00
parent 47605c147f
commit d46588f4f1
1 changed files with 35 additions and 9 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

44
testing_sf_plot_times.R
View file

@ -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)
)
)
count_data <- sum(apply(st_within(points_sf, st_buffer(sssi, dist = 50), sparse=FALSE), 1, any))
sessioninfo::session_info()