Is this package useful?

[bevingtona]

Z values are already in the FWA, see example below.

Is this package still useful?

Some options:

• User can select which elevation data to use: fwa or cded or srtm etc. ?
• With a DEM, user can query elevation of a buffer region

To do:

• Compare different methods

library(stars)
#> Loading required package: abind
#> Loading required package: sf
#> Linking to GEOS 3.9.1, GDAL 3.2.1, PROJ 7.2.1
library(bcmaps)
library(bcdata)
#> 
#> Attaching package: 'bcdata'
#> The following object is masked from 'package:stats':
#> 
#>     filter
library(sf)
library(tidyverse)

my_river <- bcdc_query_geodata("freshwater-atlas-rivers") %>%
  filter(GNIS_NAME_1 == "Parsnip River") %>%
  collect()

my_river_code <- unique(my_river$FWA_WATERSHED_CODE)

my_stream_network <- bcdc_query_geodata("freshwater-atlas-stream-network") %>%
  filter(FWA_WATERSHED_CODE == my_river_code) %>%
  collect()

my_stream_network <-
  my_stream_network %>%
  filter(BLUE_LINE_KEY == unique(my_stream_network$WATERSHED_KEY)) %>% st_as_sf()

p <- st_cast(my_stream_network, to = "POINT") %>% 
  select(geometry)
  
d <- p %>% 
  mutate(coords = st_coordinates(p)) %>% 
  mutate(X = coords[,"X"],
         Y = coords[,"Y"],
         Z = coords[,"Z"]) %>% 
  select(-coords) %>% 
  arrange(Z) %>% 
  mutate(dist_seg_m = replace_na(as.numeric(st_distance(geometry, lag(geometry), by_element = TRUE)),0),
         dist_tot_m = cumsum(dist_seg_m))

d %>% 
  ggplot() + 
  geom_point(aes(dist_tot_m/1000,Z,color=Z))

^{Created on 2022-01-23 by the reprex package (v2.0.1)}