forked from USGS-R/delaware-model-prep
-
Notifications
You must be signed in to change notification settings - Fork 0
/
8_visualize.yml
139 lines (119 loc) · 4.92 KB
/
8_visualize.yml
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
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
target_default: 8_visualize
packages:
- sf
- tidyverse
- lwgeom
- smoothr
- igraph
- purrr
- lubridate
sources:
- 8_visualize/src/visualization_functions.R
- 8_visualize/src/freq_dist_functions.R
- 8_visualize/src/visualize_sntemp_predictions.R
- 8_visualize/src/plot_rmse_map.R
targets:
8_visualize:
depends:
- 8_visualize/out/map_network_subset.png
- 8_visualize/out/map_distance_down_from_579.png
- 8_visualize/out/map_distance_up_from_579.png
- 8_visualize/out/distance_heatmap_full_up.png
- 8_visualize/out/distance_heatmap_full_down.png
- 8_visualize/out/distance_heatmap_full_updown.png
- 8_visualize/out/site_obscount_freq_dist.rds
- 8_visualize/out/site_obscount_freq_dist.png
- 8_visualize/out/nreaches_nobs_nyears.png
- 8_visualize/out/pred_obs_1578.png
- 8_visualize/out/rgcn_mods_comparison.png
- 8_visualize/out/rgcn2_full_rmse_on_network.png
- 8_visualize/out/sntemp_rmse_on_network.png
# map of network and subset
8_visualize/out/map_network_subset.png:
command: plot_subnet(
out_file = target_name,
subnet_ind = '9_collaborator_data/umn/network_subset.rds.ind',
network_ind = '1_network/out/network.rds.ind',
summary_ind = '9_collaborator_data/mkrsp/delaware_sites_summary.rds.ind',
sites_ind = '2_observations/out/drb_filtered_sites.rds.ind')
# plot distances from a particular reach
8_visualize/out/map_distance_down_from_579.png:
command: plot_dists(
out_file = target_name,
from_reach = I('579_1'),
dist_mat_ind = '9_collaborator_data/umn/distance_matrix_subset.rds.ind',
direction = I('downstream'),
labels = I('subseg_id'),
network_ind = '9_collaborator_data/umn/network_subset.rds.ind',
title = I('Subnetwork - Downstream'))
8_visualize/out/map_distance_up_from_579.png:
command: plot_dists(
out_file = target_name,
from_reach = I('579_1'),
dist_mat_ind = '9_collaborator_data/umn/distance_matrix_subset.rds.ind',
direction = I('upstream'),
label = I('subseg_id'),
network_ind = '9_collaborator_data/umn/network_subset.rds.ind',
title = I('Subnetwork - Upstream'))
# last plot is distance heatmaps that uses dist_heatmap
8_visualize/out/distance_heatmap_full_up.png:
command: dist_heatmap(
out_file = target_name,
dist_ind = '1_network/out/subseg_distance_matrix.rds.ind',
labels = I('seg_id_nat'),
title = I('Full Network - upstream'),
direction = I('upstream'))
8_visualize/out/distance_heatmap_full_down.png:
command: dist_heatmap(
out_file = target_name,
dist_ind = '1_network/out/subseg_distance_matrix.rds.ind',
labels = I('seg_id_nat'),
title = I('Full Network - downstream'),
direction = I('downstream'))
8_visualize/out/distance_heatmap_full_updown.png:
command: dist_heatmap(
out_file = target_name,
dist_ind = '1_network/out/subseg_distance_matrix.rds.ind',
labels = I('seg_id_nat'),
title = I('Full Network - Upstream & Downstream'),
direction = I('updown'))
# plot how many sites have x number of observation days
8_visualize/out/site_obscount_freq_dist.rds:
command: get_freq_dist(
out_file = target_name,
dat_ind = '2_observations/out/obs_temp_drb.rds.ind')
8_visualize/out/site_obscount_freq_dist.png:
command: plot_freq_dist(
out_file = target_name,
freq_dat = '8_visualize/out/site_obscount_freq_dist.rds',
cutoffs = I(c(1, 10, 100, 1000, 10000)))
8_visualize/out/nreaches_nobs_nyears.png:
command: plot_year_obs_tradeoff(
dat_ind = '2_observations/out/obs_temp_drb.rds.ind',
out_file = target_name)
# plot SNTemp predictions vs observations for a given reach
8_visualize/out/pred_obs_1578.png:
command: plot_pred_vs_obs(
out_file = target_name,
site = I('1578'),
start_date = I('2010-01-01'),
end_date = I('2015-12-31'),
dat_ind = '2_observations/out/obs_temp_drb.rds.ind',
pred_ind = '3_predictions/out/uncal_sntemp_preds.csv.ind')
# compare Jeff and Xiaowei's model output
8_visualize/out/rgcn_mods_comparison.png:
command: plot_rgcns(compare_ind = '3_predictions/out/compare_predictions_obs.feather.ind', out_file = target_name)
# plot Jeff's model RMSE on the network
# can plot any model that's in 3_predictions/out/compare_predictions_obs.feather
8_visualize/out/rgcn2_full_rmse_on_network.png:
command: plot_rmse_map(
out_file = target_name,
network_ind = '1_network/out/network.rds.ind',
compare_ind = '3_predictions/out/compare_predictions_obs.feather.ind',
plot_model = I('rgcn2_full_temp_c'))
8_visualize/out/sntemp_rmse_on_network.png:
command: plot_rmse_map(
out_file = target_name,
network_ind = '1_network/out/network.rds.ind',
compare_ind = '3_predictions/out/compare_predictions_obs.feather.ind',
plot_model = I('sntemp_temp_c'))