Skip to contents

Safe Clip a Multipolygon Vector Layer

Source: R/clip_layer.R
clip_multipoligon.Rd

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.

Value

A `sf` vector layer with the clipped geometries.

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

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)