Commit 8ca8ef89 authored by LDalby's avatar LDalby
Browse files

Don't read the map from o

parent ca0c35d6
......@@ -30,9 +30,9 @@ library(fs)
library(colorblindr)
library(rebus)
library(patchwork)
library(here)
resultpath <- path_home("almass", "run-directory")
o_drive <- "/Volumes/dfs/ST_GooseProject/"
git_dir <- path_home("git")
runinfo = readLines(path(resultpath, 'BatchALMaSS.ini'))
simyears = as.numeric(runinfo[length(runinfo) - 1])
......@@ -54,8 +54,7 @@ Date stamp on ALMaSS exe file: `r file.info(paste0(resultpath, 'ALMaSS_CmdLine.e
## Notes to run
`r Sys.Date()`.
+ See the config table below for the settings used for the goose model.
Standard run
## Available forage
```{r forage, fig.width=14, fig.height=16}
......@@ -333,7 +332,9 @@ forage %>%
utm32 <- "+proj=utm +zone=32 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"
# Define AOR grid size
grid_size = 1000
ext <- read_sf(path(o_drive, "/ALMaSS/GIS/vejlerne-outline_2018-07-04.shp")) %>%
ext <- read_sf(here::here("data", "vejlerne-outline.gpkg")) %>%
st_transform(crs = "+proj=utm +zone=32 +ellps=GRS80 +units=m +no_defs") %>%
raster::extent()
# Define template raster
r1 <- raster::raster(nrows = 24000/grid_size,
......@@ -352,24 +353,10 @@ polygons <- read_tsv(path(resultpath, "GooseXYDump.txt"), col_types = "iiiiii")
# Tally occurences in cells and plot
raster::rasterize(y = r1, fun = "count") %>%
raster::rasterToPolygons() %>%
# st_as_sf() %>%
# st_transform(crs = "+init=epsg:4326") %>%
# ggplot() +
# geom_sf(aes(fill = layer)) +
# ylim(56.935, 57.15) +
# xlim(8.73, 9.35) +
# ggtitle("Grid cell occupancy - Pinkfoot",
# sub = paste0("Grid size = 1000m. Followinglikelyhood ", defaultval))
sp::spTransform(CRSobj = sp::CRS("+init=epsg:4326"))
# Get a map background
# fields <- readShapePoly("C:/Users/lada/Git/shiny/vejlerne/Data/fields.shp",
# proj4string = CRS("+proj=longlat +datum=WGS84 +no_defs")) %>%
# spTransform(CRSobj = CRS("+init=epsg:4326")) %>%
map <- st_read(dsn = path(o_drive, "/ALMaSS/GIS/"),
layer = "vejlerne-outline_2018-07-04",
quiet = TRUE) %>%
map <- st_read(here::here("data", "vejlerne-outline.gpkg"),
quiet = TRUE) %>%
st_transform(crs = "+init=epsg:4326") %>%
st_bbox() %>%
as.numeric() %>%
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment