clip_multipoligon {sigugr}R Documentation

Safe Clip a Multipolygon Vector Layer

Description

Clips a 'MULTIPOLYGON' vector layer using a polygon layer, handling specific issues that might arise with geometries encoded incorrectly or containing unknown WKB types. It serves as a fallback when the 'clip_layer' function fails due to errors like 'ParseException: Unknown WKB type 12', which is associated with *MULTIPOLYGON* types.

Usage

clip_multipoligon(vector, polygon)

Arguments

vector

A 'sf' multipolygon vector layer to be clipped.

polygon

A 'sf' polygon layer used as the clipping geometry.

Details

The function ensures that the input layer is correctly encoded as 'MULTIPOLYGON' and uses GDAL utilities for re-encoding if necessary. The output is projected to the CRS of the clipping polygon.

This solution is inspired by a discussion on handling WKB type errors in R: <https://gis.stackexchange.com/questions/389814/r-st-centroid-geos-error-unknown-wkb-type-12>.

Value

A 'sf' vector layer with the clipped geometries.

See Also

Other clip functions: clip_layer(), clip_raster(), generate_bbox()

Examples

gpkg_path <- system.file("extdata", "clc.gpkg", package = "clc")

clc <- sf::st_read(gpkg_path, layer = "clc", quiet = TRUE)
lanjaron <- sf::st_read(gpkg_path, layer = "lanjaron", quiet = TRUE)

clc_clipped <- clip_multipoligon(clc, lanjaron)


[Package sigugr version 1.0.0 Index]