# alphahull: an R Package for Alpha-Convex Hull

[This article was first published on

Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

A new paper on the α-convex hull appeared in the Journal of Statistical Software today (http://www.jstatsoft.org/v34/i05/paper). The α-convex hull is an interesting problem which caught my attention long time ago but I didn’t know a solution then. R has a function **Statistics, R, Graphics and Fun » R Language**, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

`chull()`

which can generate (indices of) the convex hull for a series of points. Now we can use the R package `alphahull`

to compute the α-convex hull. For those who are not familiar with the α-convex hull, the animation below might be a good illustration for the difference between a convex hull and an α-convex hull. Note how the parameter α affects the shape of the hull:
The above animation can be reproduced with the code below (uncomment the lines to create a GIF animation with the `animation`

package):

set.seed(123) theta = runif(n <- 300, 0, 2 * pi) r = sqrt(runif(n, 0.25^2, 0.5^2)) x = cbind(0.5 + r * cos(theta), 0.5 + r * sin(theta)) ## library(animation) ## saveMovie({ library(alphahull) par(mar = rep(0, 4), xaxt = "n", yaxt = "n", bg = "black", col = "white") for (alpha in seq(0.25, 0, -0.01)) { plot(ahull(x, alpha = alpha), pch = 20, col = "white", panel.last = text(0.5, 0.5, sprintf("alpha = %.2f", alpha), cex = 2)) } ## }, moviename = "alpha-convex-hull")

Again, you may interactively play with the convex hull using the `gWidgets`

package:

## install.packages('gWidgetsRGtk2') first if not installed library(gWidgetsRGtk2) options(guiToolkit = "RGtk2") g = glayout(container = gwindow("alpha-convex hull demo")) g[1, 1:2, expand = TRUE] = ggraphics(container = g) g[2, 1] = "alpha" g[2, 2, expand = TRUE] = gslider(from = 0, to = 0.3, value = 0.2, by = 0.01, container = g, handler = function(h, ...) { par(mar = rep(0, 4), xaxt = "n", yaxt = "n") plot(ahull(x, alpha = svalue(h$obj)), pch = 20, ann = FALSE) })

## Related Posts

To

**leave a comment**for the author, please follow the link and comment on their blog:**Statistics, R, Graphics and Fun » R Language**.R-bloggers.com offers

**daily e-mail updates**about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.

Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.