-
Notifications
You must be signed in to change notification settings - Fork 7
/
Copy path05.3_alltrips_scenario_building_bodge.R
52 lines (44 loc) · 3.2 KB
/
05.3_alltrips_scenario_building_bodge.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
#SET UP
rm(list = ls())
source("00_setup_and_funs.R")
memory.limit(size=1000000)
purpose <- "alltrips"
geography <- "msoa"
#########################
### TEMP BODGE: DOWNLOAD EXISTING ALL TRIPS DATA AND RENAME VARIABLES TO WORK IN LATEST VERSION
#########################
region <- "greater-manchester"
z_codebook <- read_csv(file.path(path_outputs_regional_R, purpose, geography, region, "alltrips_codebooks_temp/z_codebook.csv"))
c_codebook <- read_csv(file.path(path_outputs_regional_R, purpose, geography, region, "alltrips_codebooks_temp/c_codebook.csv"))
od_l_rf_codebook <- read_csv(file.path(path_outputs_regional_R, purpose, geography, region, "alltrips_codebooks_temp/od_l_rf_codebook.csv"))
rq_codebook <- read_csv(file.path(path_outputs_regional_R, purpose, geography, region, "alltrips_codebooks_temp/rq_codebook.csv"))
rnet_codebook <- read_csv(file.path(path_outputs_regional_R, purpose, geography, region, "alltrips_codebooks_temp/rnet_codebook.csv"))
z <- readRDS(file.path(path_outputs_regional_R, purpose, geography, region, "github_originals/z.Rds"))
z@data <- dplyr::rename(z@data, geo_name = `geo_label`)
z@data <- z@data[z_codebook$`Variable name`]
saveRDS(z, file.path(path_outputs_regional_R, purpose, geography, region, "z.Rds"), version = 2)
c <- readRDS(file.path(path_outputs_regional_R, purpose, geography, region, "github_originals/c.Rds"))
c@data <- dplyr::rename(c@data, geo_name = `geo_label`)
c@data <- c@data[c_codebook$`Variable name`]
saveRDS(c, file.path(path_outputs_regional_R, purpose, geography, region, "c.Rds"), version = 2)
l <- readRDS(file.path(path_outputs_regional_R, purpose, geography, region, "github_originals/l.Rds"))
l@data <- dplyr::rename(l@data, geo_code1 = `msoa1`, geo_code2 = `msoa2`, geo_name1 = `geo_label1`,
geo_name2 = `geo_label2`, e_dist_km = `dist`, rf_dist_km = `dist_fast`,
rq_dist_km = `dist_quiet`, dist_rf_e = `cirquity`,
dist_rq_rf = `distq_f`, rf_avslope_perc = `avslope`,
rq_avslope_perc = `avslope_q`, rf_time_min = `time_fast`, rq_time_min = `time_quiet`)
l@data <- l@data[od_l_rf_codebook$`Variable name`]
saveRDS(l, file.path(path_outputs_regional_R, purpose, geography, region, "l.Rds"), version = 2)
rf <- readRDS(file.path(path_outputs_regional_R, purpose, geography, region, "github_originals/rf.Rds"))
rf@data <- left_join(rf@data, l@data, by="id")
rf@data <- rf@data[od_l_rf_codebook$`Variable name`]
saveRDS(rf, file.path(path_outputs_regional_R, purpose, geography, region, "rf.Rds"), version = 2)
rq <- readRDS(file.path(path_outputs_regional_R, purpose, geography, region, "github_originals/rq.Rds"))
rq@data <- left_join(rq@data, l@data, by="id")
rq@data <- rq@data[rq_codebook$`Variable name`]
saveRDS(rq, file.path(path_outputs_regional_R, purpose, geography, region, "rq.Rds"), version = 2)
rnet <- readRDS(file.path(path_outputs_regional_R, purpose, geography, region, "github_originals/rnet.Rds"))
rnet@data <- dplyr::rename(rnet@data, singlezone = `Singlezone`)
rnet@data$local_id <- 1:nrow(rnet)
rnet@data <- rnet@data[rnet_codebook$`Variable name`]
saveRDS(rnet, file.path(path_outputs_regional_R, purpose, geography, region, "rnet.Rds"), version = 2)