-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy path_targets.R
102 lines (77 loc) · 3.38 KB
/
_targets.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
# see the stan example on the reference website to investigate structure of functions
#https://github.com/wlandau/targets-stan/blob/main/R/utils.R
library(targets)
library(tarchetypes)
source("src/functions.R")
source("src/utils.R")
tar_option_set(packages = c("data.table", "leaflet", "rmarkdown", "sf", "glatos", "viridis", "leafem", "htmlwidgets", "leaflet.extras", "terra", "knitr", "rmarkdown", "geobuffer", "geosphere", "flextable"))
list(
tar_target(
raw_bathy,
"data/LH_bathy/huron_lld.tif", # raw file for lake huron bathymetry data
format = "file"
),
tar_target(
bathy,
prep_bathy(raw_bathy, xmin_out = -84.288, xmax_out = -83.540, ymin_out = 45.438, ymax_out = 45.841), # writes output to ~/Documents/cisco_stocking_survival_pilot/plan_map/output/bathy.tif, preps bathy data for use later
format = "file"
),
tar_target(
leaflet_HB,
leaflet_map(bathy, pth = "docs/index_HB.html", lines = recs_HB), #create leaflet map appropriate for webpage. Also specify bridle length classes.
format = "file"
),
tar_target(
leaflet_SB,
leaflet_map(bathy, pth = "docs/index_SB.html", lines = recs_SB, glider_lat = c(44.08313, 44.16488), glider_lon = c(-83.52757, -83.48537)), #create leaflet map appropriate for webpage. Also specify bridle length classes.
format = "file"
),
## tar_target(
## parallel_lines,
## find_parallel_lines(start_lat = 45.537, start_lon = -83.999, offset = 500, line_direction = "E", parallel_line_direction = "N", inter_pt_dist = 500, line_dist = 2000),
## format = "fst_dt"
## ),
tar_target(
parallel_lines_HB,
find_parallel_lines_split_start(start_lat = 45.5361667, start_lon = -84.00, split_dist = 250, bearing = 360-59.33614, end_lat = 45.545333, end_lon = -84.022, inter_pt_dist = 500),
format = "fst_dt"
),
tar_target(
parallel_lines_SB,
find_parallel_lines_split_start(start_lat = 44.1482, start_lon = -83.4937, split_dist = 250, end_lat = 44.16488, end_lon = -83.48537, inter_pt_dist = 500, bearing = 19.77225),
format = "fst_dt"
),
tar_target(
extract_depth_HB,
depth_extract(parallel_lines_HB, bathy, lidar = NULL),
format = "fst_dt"
),
tar_target(
extract_depth_SB,
depth_extract(parallel_lines_SB, bathy, lidar = NULL),
),
tar_target(
recs_HB,
{gear <- data.table(glatos_array = c("MBU", "VPS"), gear_config = c("VR2W-180kHz, VR2Tx-69kHz, 180kHz tag, 69 kHz tag, temp logger", "VT2W-Tx 69 khz, mid water, sync tag active"));
extract_depth_HB[gear, config := gear_config, on = "glatos_array"]},
format = "fst_dt"
),
tar_target(
recs_SB,
{gear <- data.table(glatos_array = c("MBU", "VPS"), gear_config = c("VR2W-180kHz, VR2Tx-69kHz, 180kHz tag, 69 kHz tag, temp logger", "VT2W-Tx 69 khz, mid water, sync tag active"));
extract_depth_SB[gear, config := gear_config, on = "glatos_array"]},
format = "fst_dt"
),
tar_render(report_HB, "src/rec_points_2021_HB.rmd", output_dir = "output", output_file = "output/gear_deployment_coords_HB.html"), # create html output of points),
tar_render(report_SB, "src/rec_points_2021_SB.rmd", output_dir = "output", output_file = "output/gear_deployment_coords_SB.html"),
tar_target(
gpx_SB,
make_gpx(recs_SB, pth = "output/glider_rec_SB.gpx", idcol = "site"),
format = "file"
),
tar_target(
gpx_HB,
make_gpx(recs_HB, pth = "output/glider_rec_HB.gpx", idcol = "site"),
format = "file"
)
)