Calculates the the proportion of catchment that the gauge occupies using Thiessen / Voronoi polygons
Arguments
- intersectPoly
Shapefile of the intersected catchment polygon
- catchment
Catchment shapefile
Examples
library(riskyData)
data(crowle); data(bickley); data(barnhurst); data(hollies); data(ledbury);
data(bettwsYCrwyn); data(bewdCatch)
# Obtain gauge coordinates
gcs <- getCoords(crowle,
bickley,
barnhurst,
hollies,
ledbury,
bettwsYCrwyn)
bewdTeeSun <- teeSun(gaugeCoords = gcs, catchment = bewdCatch)
int <- intersectPoly(coords = gcs,
voronoi = bewdTeeSun,
catchment = bewdCatch)
gaugeProp(intersectPoly = int)
#> Station WISKI Area Proportion
#> 1: Bettws-Y-Crwyn 441011 2465.18 57.127827
#> 2: Hollies 091862 1050.89 24.353217
#> 3: Bickley 445031 632.26 14.651928
#> 4: Barnhurst 091266 166.87 3.867028