HydroBASINS

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

library(dplyr)

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
library(sf)
Linking to GEOS 3.11.0, GDAL 3.5.3, PROJ 9.1.0; sf_use_s2() is TRUE
library(terra)
terra 1.7.39
library(tidyterra)

Attaching package: 'tidyterra'
The following object is masked from 'package:stats':

    filter
library(glue)

Attaching package: 'glue'
The following object is masked from 'package:terra':

    trim
library(ggplot2)
library(ggthemes)
library(stars)
Loading required package: abind
flood_polygon_data <- glue("/vsizip/vsicurl/https://data.hydrosheds.org/file/HydroBASINS/standard/hybas_na_lev06_v1c.zip") %>%
  vect() %>%
  st_as_sf(coords = c("lon","lat"))

flood_polygon_data
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

flood_polygon_data <- 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"))

flood_polygon_data
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

flood_polygon_data <- 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"))

flood_polygon_data
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 
flood_polygon_data2 <- 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"))

flood_polygon_data2
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 
flood_polygon_data3 <- 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"))

flood_polygon_data3
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

flood_polygon_data <- 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"))

flood_polygon_data
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 
flood_polygon_data2 <- 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"))

flood_polygon_data2
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 
flood_polygon_data3 <- 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"))

flood_polygon_data3
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

flood_polygon_data3 <- 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"))

flood_polygon_data3
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

flood_polygon_data3 <- 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"))

flood_polygon_data3
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

flood_polygon_data3 <- 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"))

flood_polygon_data3
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

flood_polygon_data3 <- 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"))

flood_polygon_data3
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

flood_polygon_data3 <- glue("/vsizip/vsicurl/https://data.hydrosheds.org/file/hydrosheds-associated/gloric/GloRiC_Canada_v10_shapefile.zip/GloRiC_Canada_v10_shapefile.csv") %>%
  vect()
flood_polygon_data3
ggplot() +

  geom_spatraster(data=flood_polygon_data3) +
  theme_tufte()