.Rhistory

longitude=as.numeric(longitude),
start_date=as.numeric(start_date),
end_date=as.numeric(end_date)) %>%
mutate(rec_len = end_date - start_date) %>%
filter(!is.na(latitude) | !is.na(longitude)) %>%
filter(network=='USGS')  %>%
distinct()
glimpse(gwwld)
glimpse(usgs_floodperc)
mutate(site_number=as.numeric(site_number),
latitude=as.numeric(latitude),
read.csv('../data/gwwld_scope_2023.csv') %>%
mutate(site_number=as.numeric(site_number),
gwwld <-
read.csv('../data/gwwld_scope_2023.csv') %>%
select(network, network_id, latitude, longitude, start_date, end_date) %>%
mutate(site_number=as.numeric(site_number),
latitude=as.numeric(latitude),
longitude=as.numeric(longitude),
start_date=as.numeric(start_date),
end_date=as.numeric(end_date)) %>%
mutate(rec_len = end_date - start_date) %>%
filter(!is.na(latitude) | !is.na(longitude)) %>%
filter(network=='USGS')  %>%
distinct()
qqq
q
q
q
q
q
gwwld <-
read.csv('../data/gwwld_scope_2023.csv') %>%
select(network, network_id, latitude, longitude, start_date, end_date) %>%
mutate(site_number=as.numeric(site_number),
latitude=as.numeric(latitude),
longitude=as.numeric(longitude),
start_date=as.numeric(start_date),
end_date=as.numeric(end_date)) %>%
mutate(rec_len = end_date - start_date) %>%
filter(!is.na(latitude) | !is.na(longitude)) %>%
filter(network=='USGS')  %>%
distinct()
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
filter(percent_flooded > 0.01)
glimpse(usgs_floodperc)
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv')
glimpse(usgs_floodperc)
sf::st_as_sf( usgs_floodperc, wkt = "geometry" )
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
filter(percent_flooded > 0.01) %>%
mutate(keyword = strsplit(keyword, ",")) %>%
unnest(keyword) %>%
mutate(value = 1)
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
filter(percent_flooded > 0.01) %>%
mutate(geometry = strsplit(geometry, ",")) %>%
unnest(geometry) #%>%
glimpse(usgs_floodperc)
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
filter(percent_flooded > 0.01) %>%
mutate(geometry = strsplit(geometry, ","))
glimpse(usgs_floodperc)
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
filter(percent_flooded > 0.01) %>%
mutate(geometry = strsplit(geometry, ",")) %>%
unnest(geometry) #%>%
glimpse(usgs_floodperc)
gwwld <-
read.csv('../data/gwwld_scope_2023.csv') %>%
select(network, network_id, latitude, longitude, start_date, end_date) %>%
mutate(site_number=as.numeric(site_number),
latitude=as.numeric(latitude),
longitude=as.numeric(longitude)) %>%
# start_date=as.numeric(start_date),
# end_date=as.numeric(end_date)) %>%
mutate(rec_len = end_date - start_date) %>%
filter(!is.na(latitude) | !is.na(longitude)) %>%
filter(network=='USGS')  %>%
distinct()
gwwld <-
read.csv('../data/gwwld_scope_2023.csv') %>%
select(network, network_id, latitude, longitude, start_date, end_date) %>%
mutate(#site_number=as.numeric(site_number),
latitude=as.numeric(latitude),
longitude=as.numeric(longitude)) %>%
# start_date=as.numeric(start_date),
# end_date=as.numeric(end_date)) %>%
mutate(rec_len = end_date - start_date) %>%
filter(!is.na(latitude) | !is.na(longitude)) %>%
filter(network=='USGS')  %>%
distinct()
glimpse(gwwld)
gwwld <-
read.csv('../data/gwwld_scope_2023.csv') %>%
select(network, network_id, latitude, longitude, start_date, end_date) %>%
mutate(#site_number=as.numeric(site_number),
latitude=as.numeric(latitude),
longitude=as.numeric(longitude)) %>%
# start_date=as.numeric(start_date),
# end_date=as.numeric(end_date)) %>%
# mutate(rec_len = end_date - start_date) %>%
# filter(!is.na(latitude) | !is.na(longitude)) %>%
filter(network=='USGS')  %>%
distinct()
glimpse(gwwld)
gwwld <-
read.csv('../data/gwwld_scope_2023.csv') %>%
select(network, network_id, latitude, longitude, start_date, end_date) %>%
mutate(#site_number=as.numeric(site_number),
latitude=as.numeric(latitude),
longitude=as.numeric(longitude)) %>%
# start_date=as.numeric(start_date),
# end_date=as.numeric(end_date)) %>%
# mutate(rec_len = end_date - start_date) %>%
filter(!is.na(latitude) | !is.na(longitude)) %>%
filter(network=='USGS')  %>%
distinct()
glimpse(gwwld)
gwwld <-
read.csv('../data/gwwld_scope_2023.csv') %>%
select(network, network_id, latitude, longitude, start_date, end_date) %>%
mutate(#site_number=as.numeric(site_number),
latitude=as.numeric(latitude),
longitude=as.numeric(longitude)) %>%
# start_date=as.numeric(start_date),
# end_date=as.numeric(end_date)) %>%
# mutate(rec_len = end_date - start_date) %>%
filter(!is.na(latitude) | !is.na(longitude) | !is.na(start_date)) %>%
filter(network=='USGS')  %>%
distinct()
glimpse(gwwld)
gwwld <-
read.csv('../data/gwwld_scope_2023.csv') %>%
select(network, network_id, latitude, longitude, start_date, end_date) %>%
mutate(network_id=as.numeric(network_id),
latitude=as.numeric(latitude),
longitude=as.numeric(longitude)) %>%
# start_date=as.numeric(start_date),
# end_date=as.numeric(end_date)) %>%
# mutate(rec_len = end_date - start_date) %>%
filter(!is.na(latitude) | !is.na(longitude) | !is.na(start_date)) %>%
filter(network=='USGS')  %>%
distinct()
glimpse(gwwld)
glimpse(usgs_floodperc)
gwwld <-
read.csv('../data/gwwld_scope_2023.csv') %>%
select(network, network_id, latitude, longitude, start_date, end_date) %>%
mutate(network_id=as.numeric(network_id),
latitude=as.numeric(latitude),
longitude=as.numeric(longitude)) %>%
# start_date=as.numeric(start_date),
# end_date=as.numeric(end_date)) %>%
# mutate(rec_len = end_date - start_date) %>%
filter(!is.na(latitude) | !is.na(longitude) | !is.na(start_date)) %>%
filter(network=='USGS')  %>%
distinct()  %>%
filter(network_id %in% usgs_floodperc$site_name)
glimpse(gwwld)
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
filter(percent_flooded > 0.01) %>%
mutate(geometry = strsplit(geometry, ",")) %>%
unnest(geometry) #%>%
glimpse(usgs_floodperc)
strsplit
?strsplit
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
filter(percent_flooded > 0.01) %>%
mutate(geometry = strsplit(geometry, ",")) %>%
unlist(geometry) #%>%
glimpse(usgs_floodperc)
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
filter(percent_flooded > 0.01) %>%
mutate(geometry = strsplit(geometry, ",")) %>%
unnest(geometry) #%>%
glimpse(usgs_floodperc)
%>% unnest(geometry) #%>%
glimpse(usgs_floodperc)
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
filter(percent_flooded > 0.01) %>%
mutate(geometry = strsplit(geometry, ","))
glimpse(usgs_floodperc)
View(usgs_floodperc)
glimpse(usgs_floodperc$geometry[1])
usgs_floodperc$geometry[[1]]
usgs_floodperc$geometry[1]
usgs_floodperc$geometry
data.frame(usgs_floodperc$geometry)
a <- data.frame(usgs_floodperc$geometry)
a
View(a)
a <- data.frame(usgs_floodperc$geometry) %>% pivot_longer()
a <- data.frame(usgs_floodperc$geometry) %>% pivot_longer(cols = everything())
glimpse(usgs_floodperc$geometry[1])
View(a)
data.frame(usgs_floodperc$geometry)
View(a)
?strsplit
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
filter(percent_flooded > 0.01) %>%
mutate(geometry = strsplit(geometry, "[,]"))
glimpse(usgs_floodperc$geometry[1])
strsplit(usgs_floodperc$geometry, "[,]")
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
filter(percent_flooded > 0.01)
glimpse(usgs_floodperc$geometry)
strsplit(usgs_floodperc$geometry, "[,]")
strsplit(usgs_floodperc$geometry, "[,]") %>% unnest()
strsplit(usgs_floodperc$geometry, "[,]") %>% unlist()
separate_rows(usgs_floodperc, geometry, sep = ",")
a <- separate_rows(usgs_floodperc, geometry, sep = ",")
View(a)
usgs_floodperc %>% separate(geometry, c('lat', 'lon'))
?separate
usgs_floodperc %>% separate(geometry, c('lat', 'lon'), sep=',')
usgs_floodperc %>% separate(geometry, c('lat', 'lon'), sep=',') %>%
mutate(lat=as.numeric(lat),
lon=as.numeric(lon))
usgs_floodperc %>% separate(geometry, c('lat', 'lon'), sep=',') %>%
mutate(lat=as.numeric(lat),
lon=as.numeric(lon))
usgs_floodperc %>% separate(geometry, c('lat', 'lon'), sep=',')
usgs_floodperc %>% separate(geometry, c('lat', 'lon'), sep=',') %>%
mutate(across(c("lat", "lon"), readr::parse_number))
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
filter(percent_flooded > 0.01)
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
filter(percent_flooded > 0.01)
usgs_floodperc <-
usgs_floodperc %>% separate(geometry, c('lat', 'lon'), sep=',') %>%
mutate(across(c("lat", "lon"), readr::parse_number))
usgs_floodperc
glimpse(usgs_floodperc)
doe_sites <- data.frame(c(network='COMPASS-FME',
site_id="SERC",
Longitude= -112.6246417,
Latitude= 33.3396861),
c(network='COMPASS-FME',
site_id="SERC",
Longitude= -112.6246417,
Latitude= 33.3396861))
doe_sites
View(doe_sites)
glimpse(swot_proposal_wells)
swot_proposal_wells <-
usgs_floodperc %>%
bind_rows(., serc) %>%
bind_rows(., winous)
serc <- data.frame(c(network='COMPASS-FME',
site_id="SERC",
Longitude=-76.556279,
Latitude= 38.901403))
winous <- data.frame(c(network='COMPASS-FME',
site_id="Winous Point",
Longitude= -82.999265,
Latitude= 41.463062))
swot_proposal_wells <-
usgs_floodperc %>%
bind_rows(., serc) %>%
bind_rows(., winous)
glimpse(swot_proposal_wells)
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
filter(percent_flooded > 0.01)
usgs_floodperc <-
usgs_floodperc %>% separate(geometry, c('lat', 'lon'), sep=',') %>%
mutate(across(c("Longitude", "Latitude"), readr::parse_number))
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
filter(percent_flooded > 0.01)
usgs_floodperc <-
usgs_floodperc %>% separate(geometry, c('lat', 'lon'), sep=',') %>%
mutate(across(c("Longitude", "Latitude"), readr::parse_number))
serc <- data.frame(network='COMPASS-FME',
site_id="SERC",
Longitude=-76.556279,
Latitude= 38.901403)
winous <- data.frame(network='COMPASS-FME',
site_id="Winous Point",
Longitude= -82.999265,
Latitude= 41.463062)
swot_proposal_wells <-
usgs_floodperc %>%
bind_rows(., serc) %>%
bind_rows(., winous)
glimpse(swot_proposal_wells)
View(swot_proposal_wells)
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
filter(percent_flooded > 0.01)
usgs_floodperc <-
usgs_floodperc %>% separate(geometry, c('lat', 'lon'), sep=',') %>%
mutate(across(c("Longitude", "Latitude"), readr::parse_number))
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
filter(percent_flooded > 0.01)
usgs_floodperc <-
usgs_floodperc %>% separate(geometry, c('lat', 'lon'), sep=',') %>%
mutate(across(c("Longitude", "Latitude"), readr::parse_number))
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
filter(percent_flooded > 0.01)
usgs_floodperc <-
usgs_floodperc %>% separate(geometry, c('lat', 'lon'), sep=',') %>%
mutate(across(c("Longitude", "Latitude"), readr::parse_number))
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
filter(percent_flooded > 0.01)
usgs_floodperc <-
usgs_floodperc %>% separate(geometry, c('Longitude', 'Latitude'), sep=',') %>%
mutate(across(c("Longitude", "Latitude"), readr::parse_number))
serc <- data.frame(network='COMPASS-FME',
site_name="SERC",
Longitude=-76.556279,
Latitude= 38.901403)
winous <- data.frame(network='COMPASS-FME',
site_name="Winous Point",
Longitude= -82.999265,
Latitude= 41.463062)
swot_proposal_wells <-
usgs_floodperc %>%
bind_rows(., serc) %>%
bind_rows(., winous)
winous
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
filter(percent_flooded > 0.01)
usgs_floodperc <-
usgs_floodperc %>% separate(geometry, c('longitude', 'latitude'), sep=',') %>%
mutate(across(c("longitude", "latitude"), readr::parse_number))
serc <- data.frame(network='COMPASS-FME',
site_name="SERC",
longitude=-76.556279,
latitude= 38.901403)
winous <- data.frame(network='COMPASS-FME',
site_name="Winous Point",
longitude= -82.999265,
latitude= 41.463062)
swot_proposal_wells <-
usgs_floodperc %>%
bind_rows(., serc) %>%
bind_rows(., winous)
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
mutate(site_name =as.character(site_name)) %>%
filter(percent_flooded > 0.01)
usgs_floodperc <-
usgs_floodperc %>% separate(geometry, c('longitude', 'latitude'), sep=',') %>%
mutate(across(c("longitude", "latitude"), readr::parse_number))
serc <- data.frame(network='COMPASS-FME',
site_name="SERC",
longitude=-76.556279,
latitude= 38.901403)
winous <- data.frame(network='COMPASS-FME',
site_name="Winous Point",
longitude= -82.999265,
latitude= 41.463062)
swot_proposal_wells <-
usgs_floodperc %>%
bind_rows(., serc) %>%
bind_rows(., winous)
glimpse(swot_proposal_wells)
View(swot_proposal_wells)
swot_proposal_wells <-
usgs_floodperc %>%
bind_rows(., serc) %>%
bind_rows(., winous) %>%
select(network, site_name, begin_date, end_date, percent_flooded, longitude, latitude)
glimpse(swot_proposal_wells)
swot_proposal_wells_sf = st_as_sf(swot_proposal_wells, coords = c("longitude", "latitude"), crs = 4326)#, relation_to_geometry = "field")
swot_proposal_wells_sf
getwd()
st_write(swot_proposal_wells_sf, "../../swot_proposal_wells.shp")
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
mutate(site_name =as.character(site_name)) %>%
filter(percent_flooded > 0.000001)
usgs_floodperc <-
usgs_floodperc %>% separate(geometry, c('longitude', 'latitude'), sep=',') %>%
mutate(across(c("longitude", "latitude"), readr::parse_number))
serc <- data.frame(network='COMPASS-FME',
site_name="SERC",
longitude=-76.556279,
latitude= 38.901403)
winous <- data.frame(network='COMPASS-FME',
site_name="Winous Point",
longitude= -82.999265,
latitude= 41.463062)
swot_proposal_wells <-
usgs_floodperc %>%
bind_rows(., serc) %>%
bind_rows(., winous) %>%
select(network, site_name, begin_date, end_date, percent_flooded, longitude, latitude)
glimpse(swot_proposal_wells)
usgs_floodperc <- read.csv('../data/sources/usgs/percent_flooded_by_ecoregion.csv') %>%
select(site_number, begin_date, end_date, percent_flooded, geometry) %>%
mutate(network='USGS') %>%
rename(site_name= site_number) %>%
mutate(site_name =as.character(site_name)) %>%
filter(percent_flooded > 0)
usgs_floodperc <-
usgs_floodperc %>% separate(geometry, c('longitude', 'latitude'), sep=',') %>%
mutate(across(c("longitude", "latitude"), readr::parse_number))
serc <- data.frame(network='COMPASS-FME',
site_name="SERC",
longitude=-76.556279,
latitude= 38.901403)
winous <- data.frame(network='COMPASS-FME',
site_name="Winous Point",
longitude= -82.999265,
latitude= 41.463062)
swot_proposal_wells <-
usgs_floodperc %>%
bind_rows(., serc) %>%
bind_rows(., winous) %>%
select(network, site_name, begin_date, end_date, percent_flooded, longitude, latitude)
glimpse(swot_proposal_wells)
swot_proposal_wells_sf = st_as_sf(swot_proposal_wells, coords = c("longitude", "latitude"), crs = 4326)#, relation_to_geometry = "field")
swot_proposal_wells_sf
st_write(swot_proposal_wells_sf, "../../swot_proposal_wells.shp")
st_write(swot_proposal_wells_sf, "../swot_proposal_wells.shp")
gw <- st_read("../combined/gw_wells_combined_pts_v2.shp")
gw
swot_proposal_wells <-
bind_rows(serc, winous) #%>%
swot_proposal_wells <-
bind_rows(serc, winous) %>%
st_as_sf(., coords = c("longitude", "latitude"), crs = 4326)
swot_proposal_wells
swot_proposal_wells
swot_proposal_wells <-
bind_rows(serc, winous) %>%
st_as_sf(., coords = c("longitude", "latitude"), crs = 4326)
gw <- st_read("../combined/gw_wells_combined_pts_v2.shp") %>%
mutate(site_name = site_no) %>%
mutate(network= source)
swot_proposal_wells <- bind_rows(gw, swot_proposal_wells)
gw <- st_read("../combined/gw_wells_combined_pts_v2.shp") %>%
mutate(site_name = as.character(site_no)) %>%
mutate(network= source)
swot_proposal_wells <- bind_rows(gw, swot_proposal_wells)
swot_proposal_wells
View(swot_proposal_wells)
swot_proposal_wells <-
bind_rows(serc, winous) %>%
st_as_sf(., coords = c("longitude", "latitude"), crs = 4326)
#%>%
# select(network, site_name, begin_date, end_date, percent_flooded, longitude, latitude)
gw <- st_read("../combined/gw_wells_combined_pts_v2.shp") %>%
mutate(site_name = as.character(site_no)) %>%
mutate(network= source)
swot_proposal_wells <- bind_rows(gw, swot_proposal_wells) %>%
select(network, site_name, geometry)
st_write(swot_proposal_wells, "../swot_proposal_gw_wells.shp")