smooth_chaikin: Chaikin's corner cutting algorithm

In smoothr: Smooth and Tidy Spatial Features

Chaikin's corner cutting algorithm

Description

Chaikin's corner cutting algorithm smooths a curve by iteratively replacing every point by two new points: one 1/4 of the way to the next point and one 1/4 of the way to the previous point.

Usage

smooth_chaikin(x, wrap = FALSE, refinements = 3L)

Arguments

x	numeric matrix; 2-column matrix of coordinates.
wrap	logical; whether the coordinates should be wrapped at the ends, as for polygons and closed lines, to ensure a smooth edge.
refinements	integer; number of corner cutting iterations to apply.

Details

This function works on matrices of points and is generally not called directly. Instead, use smooth() with method = "chaikin" to apply this smoothing algorithm to spatial features.

Value

A matrix with the coordinates of the smoothed curve.

References

The original reference for Chaikin's corner cutting algorithm is:

• ⁠Chaikin, G. An algorithm for high speed curve generation. Computer Graphics and Image Processing 3 (1974), 346–349⁠

This implementation was inspired by the following StackOverflow answer:

• Where to find Python implementation of Chaikin's corner cutting algorithm?

See Also

smooth()

Examples

# smooth_chaikin works on matrices of coordinates
# use the matrix of coordinates defining a polygon as an example
m <- jagged_polygons$geometry[[2]][[1]]
m_smooth <- smooth_chaikin(m, wrap = TRUE)
class(m)
class(m_smooth)
plot(m, type = "l", axes = FALSE, xlab = NA, ylab = NA)
lines(m_smooth, col = "red")

# smooth is a wrapper for smooth_chaikin that works on spatial features
library(sf)
p <- jagged_polygons$geometry[[2]]
p_smooth <- smooth(p, method = "chaikin")
class(p)
class(p_smooth)
plot(p)
plot(p_smooth, border = "red", add = TRUE)

smoothr documentation built on March 31, 2023, 11:45 p.m.