Lehner, B., Grill G. (2013). Global river hydrography and network routing: baseline data and new approaches to study the world’s large river systems. Hydrological Processes, 27(15): 2171–2186. https://doi.org/10.1002/hyp.9740
https://www.hydrosheds.org/products/hydrobasins
Attaching package: 'dplyr'
The following objects are masked from 'package:stats':
filter, lag
The following objects are masked from 'package:base':
intersect, setdiff, setequal, union
Linking to GEOS 3.11.0, GDAL 3.5.3, PROJ 9.1.0; sf_use_s2() is TRUE
Attaching package: 'tidyterra'
The following object is masked from 'package:stats':
filter
Attaching package: 'glue'
The following object is masked from 'package:terra':
trim
library (ggplot2)library (ggthemes)library (stars)
Loading required package: abind
<- glue ("/vsizip/vsicurl/https://data.hydrosheds.org/file/HydroBASINS/standard/hybas_na_lev06_v1c.zip" ) %>% vect () %>% st_as_sf (coords = c ("lon" ,"lat" ))
Simple feature collection with 2043 features and 13 fields
Geometry type: GEOMETRY
Dimension: XY
Bounding box: xmin: -137.9625 ymin: 5.495833 xmax: -52.61605 ymax: 62.74232
Geodetic CRS: WGS 84
First 10 features:
HYBAS_ID NEXT_DOWN NEXT_SINK MAIN_BAS DIST_SINK DIST_MAIN SUB_AREA
1 -1529934582 0 -1529934582 -1529934582 0 0 5285.0
2 -1529934352 0 -1529934352 -1529934352 0 0 3963.9
3 -1529934342 0 -1529934342 -1529934342 0 0 25440.1
4 -1529933372 0 -1529933372 -1529933372 0 0 5014.5
5 -1529933362 0 -1529933362 -1529933362 0 0 10703.8
6 -1529932942 0 -1529932942 -1529932942 0 0 3269.3
7 -1529932932 0 -1529932932 -1529932932 0 0 11559.6
8 -1529932252 0 -1529932252 -1529932252 0 0 3698.1
9 -1529932242 0 -1529932242 -1529932242 0 0 2494.1
10 -1529932192 0 -1529932192 -1529932192 0 0 8088.1
UP_AREA PFAF_ID ENDO COAST ORDER SORT geometry
1 5285.0 771111 0 1 0 1 MULTIPOLYGON (((-79.51389 8...
2 3963.9 771112 0 0 1 2 POLYGON ((-81.02917 8.19166...
3 25440.1 771113 0 1 0 3 MULTIPOLYGON (((-81.75972 8...
4 5014.6 771114 0 0 1 4 POLYGON ((-83.61667 9.55833...
5 10703.8 771115 0 1 0 5 POLYGON ((-84.64167 9.7875,...
6 3269.7 771116 0 0 1 6 POLYGON ((-85.60833 10.8458...
7 11559.6 771117 0 1 0 7 MULTIPOLYGON (((-84.95972 9...
8 3698.3 771118 0 0 1 8 POLYGON ((-87.27917 12.8125...
9 2494.1 771119 0 1 0 9 POLYGON ((-87.05 13.24167, ...
10 8088.1 771120 0 0 1 10 POLYGON ((-86.57083 13.7875...
ggplot (data= flood_polygon_data, aes (fill= ORDER)) + geom_sf () + theme_tufte ()
https://figshare.com/ndownloader/files/20087237
<- glue ("/vsizip/vsicurl/https://data.hydrosheds.org/file/hydrosheds-v1-dem/hyd_na_dem_15s.zip/hyd_na_dem_15s.tif" ) %>% rast () #%>% #st_as_sf(coords = c("lon","lat"))
class : SpatRaster
dimensions : 13920, 20640, 1 (nrow, ncol, nlyr)
resolution : 0.004166667, 0.004166667 (x, y)
extent : -138, -52, 5, 63 (xmin, xmax, ymin, ymax)
coord. ref. : lon/lat WGS 84 (EPSG:4326)
source : hyd_na_dem_15s.tif
name : Band_1
ggplot () + geom_spatraster (data= flood_polygon_data) + theme_tufte ()
SpatRaster resampled to ncells = 500822
<- glue ("/vsizip/vsicurl/https://data.hydrosheds.org/file/hydrosheds-v1-con/na_con_3s/n20w110_con.zip/n20w110_con.tif" ) %>% rast () #%>% #st_as_sf(coords = c("lon","lat"))
class : SpatRaster
dimensions : 12000, 12000, 1 (nrow, ncol, nlyr)
resolution : 0.0008333333, 0.0008333333 (x, y)
extent : -110, -100, 20, 30 (xmin, xmax, ymin, ymax)
coord. ref. : lon/lat WGS 84 (EPSG:4326)
source : n20w110_con.tif
name : n25w105_con
<- glue ("/vsizip/vsicurl/https://data.hydrosheds.org/file/hydrosheds-v1-con/na_con_3s/n20w120_con.zip/n20w120_con.tif" ) %>% rast () #%>% #st_as_sf(coords = c("lon","lat"))
class : SpatRaster
dimensions : 12000, 12000, 1 (nrow, ncol, nlyr)
resolution : 0.0008333333, 0.0008333333 (x, y)
extent : -120, -110, 20, 30 (xmin, xmax, ymin, ymax)
coord. ref. : lon/lat WGS 84 (EPSG:4326)
source : n20w120_con.tif
name : n25w115_con
<- glue ("/vsizip/vsicurl/https://data.hydrosheds.org/file/hydrosheds-v1-con/na_con_3s/n40w080_con.zip/n40w080_con.tif" ) %>% rast () #%>% #st_as_sf(coords = c("lon","lat"))
class : SpatRaster
dimensions : 12000, 12000, 1 (nrow, ncol, nlyr)
resolution : 0.0008333333, 0.0008333333 (x, y)
extent : -80, -70, 40, 50 (xmin, xmax, ymin, ymax)
coord. ref. : lon/lat WGS 84 (EPSG:4326)
source : n40w080_con.tif
name : n45w075_con
ggplot () + geom_spatraster (data= flood_polygon_data) + geom_spatraster (data= flood_polygon_data2) + geom_spatraster (data= flood_polygon_data3) + theme_tufte ()
SpatRaster resampled to ncells = 501264
SpatRaster resampled to ncells = 501264
SpatRaster resampled to ncells = 501264
<- glue ("/vsizip/vsicurl/https://data.hydrosheds.org/file/hydrosheds-v1-dir/na_dir_3s/n50w100_dir.zip/n50w100_dir.tif" ) %>% rast () #%>% #st_as_sf(coords = c("lon","lat"))
class : SpatRaster
dimensions : 12000, 12000, 1 (nrow, ncol, nlyr)
resolution : 0.0008333333, 0.0008333333 (x, y)
extent : -100, -90, 50, 60 (xmin, xmax, ymin, ymax)
coord. ref. : lon/lat WGS 84 (EPSG:4326)
source : n50w100_dir.tif
name : n55w095_dir
<- glue ("/vsizip/vsicurl/https://data.hydrosheds.org/file/hydrosheds-v1-dir/na_dir_3s/n50w110_dir.zip/n50w110_dir.tif" ) %>% rast () #%>% #st_as_sf(coords = c("lon","lat"))
class : SpatRaster
dimensions : 12000, 12000, 1 (nrow, ncol, nlyr)
resolution : 0.0008333333, 0.0008333333 (x, y)
extent : -110, -100, 50, 60 (xmin, xmax, ymin, ymax)
coord. ref. : lon/lat WGS 84 (EPSG:4326)
source : n50w110_dir.tif
name : n55w105_dir
<- glue ("/vsizip/vsicurl/https://data.hydrosheds.org/file/hydrosheds-v1-dir/na_dir_3s/n50w120_dir.zip/n50w120_dir.tif" ) %>% rast () #%>% #st_as_sf(coords = c("lon","lat"))
class : SpatRaster
dimensions : 12000, 12000, 1 (nrow, ncol, nlyr)
resolution : 0.0008333333, 0.0008333333 (x, y)
extent : -120, -110, 50, 60 (xmin, xmax, ymin, ymax)
coord. ref. : lon/lat WGS 84 (EPSG:4326)
source : n50w120_dir.tif
name : n55w115_dir
ggplot () + geom_spatraster (data= flood_polygon_data) + geom_spatraster (data= flood_polygon_data2) + geom_spatraster (data= flood_polygon_data3) + theme_tufte ()
SpatRaster resampled to ncells = 501264
SpatRaster resampled to ncells = 501264
SpatRaster resampled to ncells = 501264
<- glue ("/vsizip/vsicurl/https://data.hydrosheds.org/file/hydrosheds-v1-aca/hyd_na_aca_15s.zip/hyd_na_aca_15s.tif" ) %>% rast () #%>% #st_as_sf(coords = c("lon","lat"))
class : SpatRaster
dimensions : 13920, 20640, 1 (nrow, ncol, nlyr)
resolution : 0.004166667, 0.004166667 (x, y)
extent : -138, -52, 5, 63 (xmin, xmax, ymin, ymax)
coord. ref. : lon/lat WGS 84 (EPSG:4326)
source : hyd_na_aca_15s.tif
name : Band_1
ggplot () + geom_spatraster (data= flood_polygon_data3) + theme_tufte ()
SpatRaster resampled to ncells = 500822
<- glue ("/vsizip/vsicurl/https://data.hydrosheds.org/file/hydrosheds-v1-ldn/hyd_na_ldn_15s.zip/hyd_na_ldn_15s.tif" ) %>% rast () #%>% #st_as_sf(coords = c("lon","lat"))
class : SpatRaster
dimensions : 13920, 20640, 1 (nrow, ncol, nlyr)
resolution : 0.004166667, 0.004166667 (x, y)
extent : -138, -52, 5, 63 (xmin, xmax, ymin, ymax)
coord. ref. : lon/lat WGS 84 (EPSG:4326)
source : hyd_na_ldn_15s.tif
name : af_len_dn_15s
ggplot () + geom_spatraster (data= flood_polygon_data3) + theme_tufte ()
SpatRaster resampled to ncells = 500822
<- glue ("/vsizip/vsicurl/https://data.hydrosheds.org/file/hydrosheds-v1-msk/hyd_na_msk_15s.zip/hyd_na_msk_15s.tif" ) %>% rast () #%>% #st_as_sf(coords = c("lon","lat"))
class : SpatRaster
dimensions : 13920, 20640, 1 (nrow, ncol, nlyr)
resolution : 0.004166667, 0.004166667 (x, y)
extent : -138, -52, 5, 63 (xmin, xmax, ymin, ymax)
coord. ref. : lon/lat WGS 84 (EPSG:4326)
source : hyd_na_msk_15s.tif
name : Band_1
ggplot () + geom_spatraster (data= flood_polygon_data3) + theme_tufte ()
SpatRaster resampled to ncells = 500822
<- glue ("/vsizip/vsicurl/https://data.hydrosheds.org/file/hydrosheds-associated/gloric/GloRiC_v10_shapefile.zip/GloRiC_v10_shapefile/GloRiC_v10.shp" ) %>% st_read () %>% st_as_sf (coords = c ("lon" ,"lat" ))
ggplot () + geom_spatraster (data= flood_polygon_data3) + theme_tufte ()
https://figshare.com/ndownloader/files/31910714
https://data.hydrosheds.org/file/hydrosheds-associated/gloric/GloRiC_Canada_v10_shapefile.zip
<- glue ("/vsizip/vsicurl/https://data.hydrosheds.org/file/hydrosheds-associated/gloric/GloRiC_Canada_v10_shapefile.zip/GloRiC_Canada_v10_shapefile.csv" ) %>% vect ()
ggplot () + geom_spatraster (data= flood_polygon_data3) + theme_tufte ()
