Creates intersected polygons of catchment against Voronoi polygons
Arguments
- voronoi
Voronoi/Thiessen polygon
- catchment
Catchment shapefile
- coords
Coordinates of rain gauges set as a polygon 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)
int
#> Simple feature collection with 4 features and 4 fields
#> Geometry type: POLYGON
#> Dimension: XY
#> Bounding box: xmin: -3.745747 ymin: 52.34419 xmax: -2.227248 ymax: 52.98128
#> Geodetic CRS: WGS 84
#> stationName WISKI Easting Northing V1
#> 1 Bettws-Y-Crwyn 441011 320360 281360 POLYGON ((-2.83667 52.53842...
#> 2 Hollies 091862 381586 322452 POLYGON ((-2.377887 52.9805...
#> 3 Bickley 445031 363130 271330 POLYGON ((-2.306061 52.4176...
#> 4 Barnhurst 091266 389910 301620 POLYGON ((-2.253267 52.6815...