data-raw/snow-station.R

In bcgov/cccharts: Plot Climate Change Indicator Trends for British Columbia

# Copyright 2016 Province of British Columbia
#
# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at
#
# http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing, software distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and limitations under the License.

source("data-raw/header.R")

## Load CSV data file from BC Data Catalogue. Data licensed under the Open Data License-BC
## See metadata record in BC Data Catalogue for details on the data set.
snow_station <- read_csv("https://catalogue.data.gov.bc.ca/dataset/86526746-40dd-41d2-82c0-fbee3a2e93a2/resource/f1c6257b-6596-483e-82ac-59ccdc8df898/download/bcsnowdepthswetrendsbystation1950-2014.csv")

snow_station_observed <- read_csv("data-raw/raw_snow_data/swe_depth_anomalies_1950_2014.csv")

snow_station$StartYear <- 1950L
snow_station$EndYear <- 2014L

snow_station %<>% rename(Station = station_id,
                         Latitude = latitude,
                         Longitude = longitude,
                         Intercept = intercept,
                         Significant = sigstat)

snow_station %<>% get_ecoprovince()
snow_station %<>% arrange(Ecoprovince, Longitude, Latitude)

snow_station$Indicator <- NA
snow_station$Indicator[snow_station$measure == "depth"] <- "Snow Depth"
snow_station$Indicator[snow_station$measure == "swe"] <- "Snow Water Equivalent"

snow_station$Units <- "Percent"
snow_station$Period <- 1L

snow_station$Ecoprovince %<>%  factor(levels = ecoprovince)


snow_station %<>% mutate(Estimate = slope_percentperyear,
                         Lower = slope_lci,
                         Upper = slope_uci)

snow_station %<>% select(
  Indicator, Units, Period, StartYear, EndYear, Ecoprovince, Station, Latitude, Longitude,
  Estimate, Lower, Upper, Intercept, Significant)

snow_station %<>% arrange(Indicator, Ecoprovince, Station, StartYear, EndYear)

snow_station_observed %<>% select(Station = STATIONID, Year = YEAR, `Snow Depth` = depth_anom, `Snow Water Equivalent` = swe_anom)

snow_station_observed %<>% gather(Indicator, Value, -Station, -Year)

snow_station_observed %<>% inner_join(snow_station, by = c("Indicator", "Station"))
snow_station_observed %<>% filter(!is.na(Value))
snow_station_observed %<>% filter(Year >= StartYear & Year <= EndYear)

snow_station %<>% select(
  Indicator, Units, Period, StartYear, EndYear, Ecoprovince, Station, Latitude, Longitude,
  Estimate, Lower, Upper, Intercept, Significant)

snow_station_observed$Units <- "Percent"

snow_station_observed %<>% select(Indicator, Ecoprovince, Station, Year, Value, Units)

snow_station$Station %<>% factor(unique(snow_station$Station))
snow_station_observed$Station %<>% factor(levels = levels(snow_station$Station))

snow_station_observed$Value %<>% as.numeric()

snow_station %<>% mutate(Intercept = Intercept + Estimate * StartYear)

snow_station %<>% cccharts::change_period(10L)

use_data(snow_station, overwrite = TRUE)
use_data(snow_station_observed, overwrite = TRUE)