exec/tbl_georch.R

In cmichaud92/UCRBtools: Useful Functions for Data Manipulation, Analysis and Visualization

library(tidyverse)
library(UCRBtools)
rch <- read_csv("~/reach-table.csv") %>%
    rename(cd_georch = cd_rch)

rch2 <- tribble(
    ~cd_georch, ~cd_rvr, ~fct_rvr, ~fct_georch, ~nhd_up, ~nhd_dn, ~nm_up, ~nm_dn,
    "UCO", "CO", "Colorado River", "Upper Colorado River", 195.4, 439.5, "Grand Lake", "Prive-Stubb Diversion",
    "GRVA", "CO", "Colorado River", "Grand Valley",  195.4, 154.6, "Price-Stubb Diversion", "Loma",
    "WWRH", "CO", "Colorado River", "Westwater-Ruby-Horsethief", 154.6, 118.0, "Loma", "above Big Hole",
    "CIS", "CO", "Colorado River", "Cisco", 118.0, 95.9, "above Big Hole", "Dewey",
    "DAL", "CO", "Colorado River", "Daily", 95.9, 75.0, "Dewey", "Takeout Beach",
    "MOB", "CO", "Colorado River", "Moab", 75.0, 47.8, "Takeout Beach", "Potash",
    "MEA", "CO", "Colorado River", "Meander Canyon", 47.8, 0, "Potash", "Green River Confluence",
    "TDV", "GR", "Green River", "Tusher Diversion", 129.5, 96.6, "Tusher Diversion", "White Wash",
    "LAB", "GR", "Green River", "Labyrinth Canyon", 96.6, 34.2, "White Wash", "Millard Canyon",
    "TPK", "GR", "Green River", "Tower Park", 34.2, 19.8, "Millard Canyon", "Deadhorse Canyon",
    "STW", "GR", "Green River", "Stillwater Canyon", 19.6, 0, "Deadhorse Canyon", "Colorado River Confluence",
    "DES", "GR", "Green River", "Desolation Canyon", 190.6, 157.9, "Jack Creek", "Wire Fence Canyon",
    "GRA", "GR", "Green River", "Gray Canyon", 157.9, 129.5, "Wire Fence Canyon", "Tusher Diversion",
    "LOD", "GR", "Green River", "Lodor Canyon", 364.8, 346.3, "Gates of Lodor", "Yampa River Confluence",
    "WPC", "GR", "Green River", "Whirlpool Canyon", 346.3, 335.9, "Yampa River Confluence", "Cove Camp",
    "IRP", "GR", "Green River", "Island and Rainbow Parks", 335.9, 328.6, "Cove Camp", "Above Moonshine Rapid",
    "SPM", "GR", "Green River", "Split Mountain Canyon", 328.6, 321.3, "Above Moonshine Rapid", "Split Mountain Boat Ramp",
    "LMGR", "GR", "Green River", "Lower Middle Green River", 240.0, 190.7, "Below White Confluence", "Jack Creek",
    "OUR", "GR", "Green River", "Ouray", 240.0, 260.0, "Ouray National Fish Hatchery", "Below White Confluence",
    "UMGR", "GR", "Green River", "Upper Middle Green River", 300.1, 260.0, "Below Jensen", "Ouray National Fish Hatchery",
    "JEN", "GR", "Green River", "Jensen", 300.1, 321.3, "Split Mountain Boat Ramp", "Below Jensen",
    "LYA", "YA", "Yampa River", "Lower Yampa Canyon", 46.9, 0, "Deer Lodge", "Green River Confluence",
    "MYA", "YA", "Yampa River", "Middle Yampa River", 189.2, 46.9, "Milner", "Deer Lodge",
    "UYA", "YA", "Yampa River", "Upper Yampa River", 222.4, 189.2, "Stagecoach Resevoir", "Milner",
    )

tbl_georch <- tbl_georch %>%
    select(-id_georch) %>%
    mutate(fct_rvr = as_factor(fct_rvr),
           fct_rvr = fct_relevel(fct_rvr, "Colorado River", "Dolores River", "Gunnison River",
                                 "Green River", "San Rafael River", "Price River",
                                 "Duchesne River", "White River", "Yampa River", "Vermillion Creek", "San Juan River",
                                 "Lake Powell")) %>%
    arrange(fct_rvr, nhd_up) %>%
    mutate(idx = row_number(),
           fct_georch = as_factor(fct_georch),
           fct_georch = fct_reorder(fct_georch, idx),
           .keep = "unused") %>%
    arrange(fct_rvr, desc(nhd_up)) %>%
    mutate(id_georch = row_number()) %>%
    select(starts_with(c("id_", "cd_", "fct_")), everything()) %>%
    select(-idx)


tbl_rch <- tbl_rch

table(tbl_georch$fct_georch)

vc <- tribble(
    ~id_georch, ~cd_georch, ~cd_rvr, ~fct_rvr, ~fct_georch, ~nhd_up, ~nhd_dn, ~nm_up, ~nm_dn,
    45, "LVC", "VC", "Vermillion Creek", "Lower Vermillion Creek", .5, 0, NA, "Green River Confluence"
)

tbl_georch <- tbl_georch %>%
    select(-c(ends_with("_rch"))) %>%
    bind_rows(vc)

vc <- tribble(
    ~id_rch, ~cd_rch, ~cd_rvr, ~fct_rvr, ~fct_rch, ~nhd_up, ~nhd_dn, ~nm_up, ~nm_dn,
    45, "LVC", "VC", "Vermillion Creek", "Lower Vermillion Creek", .5, 0, NA, "Green River Confluence"
)

tbl_rch <- tbl_rch  %>%
    bind_rows(vc) %>%
    select(-id_rch) %>%
    mutate(fct_rvr = as_factor(fct_rvr),
           fct_rvr = fct_relevel(fct_rvr, "Colorado River", "Dolores River", "Gunnison River",
                                 "Green River", "San Rafael River", "Price River",
                                 "Duchesne River", "White River", "Yampa River", "Vermillion Creek", "San Juan River",
                                 "Lake Powell")) %>%
    arrange(fct_rvr, nhd_up) %>%
    mutate(idx = row_number(),
           fct_rch = as_factor(fct_rch),
           fct_rch = fct_reorder(fct_rch, idx),
           .keep = "unused") %>%
    arrange(fct_rvr, desc(nhd_up)) %>%
    mutate(id_rch = row_number()) %>%
    select(starts_with(c("id_", "cd_", "fct_")), everything()) %>%
    select(-idx)

tbl_georch <- tbl_georch
tbl_georch$nhd_up[tbl_rch$id_rch == 16] <- 321.3
tbl_georch$nhd_dn[tbl_rch$id_rch == 16] <- 300.1

tbl_georch$nhd_up[tbl_rch$id_rch == 19] <- 260.0
tbl_georch$nhd_dn[tbl_rch$id_rch == 19] <- 240.0

tbl_georch$nhd_up[tbl_rch$id_rch == 25] <- 19.8
tbl_georch$nhd_up[tbl_rch$id_rch == 20] <- 190.7

usethis::use_data(tbl_georch)

devtools::install_github("cmichaud92/UCRBtools")

levels(rch3$fct_rvr)
levels(rch3$fct_rch)
levels(tbl_rch$fct_rch)

tbl_spp <- s
t
s$nm_com[s$cd_spp == "CH"] <- "unknown chub"
usethis::use_data(tbl_spp)

#    arrange(id_rch) %>%
select(starts_with(c("id_", "cd_", "fct_")), everything()) %>%
    select(-idx)

arrange(cd_rvr, desc(nhd_up))

# correct col names
tbl_georch <- tbl_georch %>%
    rename(fct_georch = fct_rch)