Skip to contents

Calculates the the proportion of catchment that the gauge occupies using Thiessen / Voronoi polygons

Usage

gaugeProp(intersectPoly)

Arguments

intersectPoly

Shapefile of the intersected catchment polygon

catchment

Catchment shapefile

Value

Rainfall proportions for each rain gauge

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