Intersect of polygons — intersectPoly • mappER

Creates intersected polygons of catchment against Voronoi polygons

Usage

intersectPoly(voronoi, catchment, coords)

Arguments

voronoi: Voronoi/Thiessen polygon
catchment: Catchment shapefile
coords: Coordinates of rain gauges set as a polygon shapefile

Value

The calculated intersections of a polygon

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...