Skip to content

R 📦 for fast and scalable wind fetch calculations

License

Unknown, MIT licenses found

Licenses found

Unknown
LICENSE
MIT
LICENSE.md
Notifications You must be signed in to change notification settings

anguswg-ucsb/fetchr

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

23 Commits
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

fetchr

Dependencies License: MIT


The goal of fetchr is to provide a fast and efficient raster based method for calculating fetch lengths across thousands of water grid cells. Calculating fetch lengths, the distance that wind can blow in a constant direction over a body of water without interruption, can be a slow and memory intensive process when done on thousands of points of interest in multiple directions. fetchr attempts to fix this problem and allows for thousands of fetch calculations to be performed in a fraction of the time that other methods take.


Installation

You can install the development version of fetchr from GitHub with:

# install.packages("devtools")
devtools::install_github("anguswg-ucsb/fetchr")

Coastline raster

If we start with a raster representing the Southern California coast near Santa Barbara, CA.

library(fetchr)

# land raster
land_rast <- terra::rast(fetchr::land)

terra::plot(land_rast, col = "#2e8b57")


Make a binary land water rasters

We can take this land raster, indicate which cells are water cells, and create a binary land water raster (land = 0, water = 1).

# binary land water raster 
landwater <- fetchr::get_landwater(
  r           = land_rast,                 # land raster
  water_value = NA,                        # cells with a value of NA are classified as water, all other cells are land
  res         = terra::res(land_rast)[1]   # return raster with the same cell resolution as input raster
  )
#> Creating landwater raster...

terra::plot(landwater, col = c("#2e8b57", "#add8e6"))

This raster now meets all the specification for using get_fetch():

  • Binary cell values (land cells = 0 and water cells = 1)
  • Projected Coordinate Reference System
  • Regular grid cell size (same x and y cell resolution)

Calculate fetch length

Internally, get_fetch() will coerce polygon/multipolygon geometries and rasters into the required binary landwater raster. It is recommended to provide either an sf/terra polygon or a raster/terra raster with a single value for land cells and NA values for water cells. To calculate fetch distances, we can simply provide an sf/terra polygon or raster/terra raster to get_fetch()

system.time(
  
  fetch <- fetchr::get_fetch(
    r           = land_rast,     # binary land water raster
    max_dist    = 200000,        # maximum distance to calculate fetch in meters (200km)
    in_parallel = TRUE,          # run calculations in parallel
    verbose     = TRUE
    )
  
)
#> Calculating Fetch...
#>    user  system elapsed 
#>   0.037   0.012  12.351

plot(fetch)

In this example here, my machine calculated fetch distances for > 24,000 water cells in about ~3 seconds, or ~0.00014 seconds per point. That is a ~99.99% reduction in computation time compared to various other polygon based methods out there!


Multidirectional fetch calculations

The get_fetch_directions() function allows you to get get individual fetch distances in 8 directions:

multi_fetch <- fetchr::get_fetch_directions(
  r           = land_rast,     # binary land water raster
  max_dist    = 200000,        # maximum distance to calculate fetch in meters (200km)
  in_parallel = TRUE,          # run calculations in parallel
  verbose     = TRUE
  )
#> Calculating Fetch...

plot(multi_fetch)

About

R 📦 for fast and scalable wind fetch calculations

Topics

Resources

License

Unknown, MIT licenses found

Licenses found

Unknown
LICENSE
MIT
LICENSE.md

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages