R: Generate 2D Quality meshes and constrained Delaunay...

RTriangle-package {RTriangle} | R Documentation

## Generate 2D Quality meshes and constrained Delaunay triangulations

### Description

Generate 2D Quality meshes and constrained Delaunay triangulations

### Details

This package is a wrapper of Jonathan Richard Shewchuk's Triangle package. `triangulate` triangulates a Planar Straight Line Graph (PSLG), a collection of vertices and segments created with `pslg`. A mesh in the can be created within an arbitrary closed outline and the maximum area and minimum angle of the triangles in the mesh can be specified.

### Author(s)

David C. Sterratt david.c.sterratt@ed.ac.uk

### References

- http://www.cs.cmu.edu/~quake/triangle.html

- Jonathan Richard Shewchuk, Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator, in “Applied Computational Geometry: Towards Geometric Engineering” (Ming C. Lin and Dinesh Manocha, editors), volume 1148 of Lecture Notes in Computer Science, pages 203-222, Springer-Verlag, Berlin, May 1996. (From the First ACM Workshop on Applied Computational Geometry.)

- Jonathan Richard Shewchuk, Delaunay Refinement Algorithms for Triangular Mesh Generation, Computational Geometry: Theory and Applications 22(1-3):21-74, May 2002.

### See Also

`triangulate`

### Examples

    ## Create an object with a concavity
    p <- pslg(P=rbind(c(0, 0), c(0, 1), c(0.5, 0.5), c(1, 1), c(1, 0)),
              S=rbind(c(1, 2), c(2, 3), c(3, 4), c(4, 5), c(5, 1)))
    ## Plot it
    plot(p)
    ## Triangulate it
    tp <- triangulate(p)
    plot(tp)
    ## Triangulate it subject to minimum area constraint
    tp <- triangulate(p, a=0.01)
    plot(tp)
    ## Load a data set containing a hole
    A <- read.pslg(file.path(system.file(package = "RTriangle"), "extdata", "A.poly"))
    plot(A)
    ## Triangulate the PSLG
    tA <- triangulate(A)
    plot(tA)
    ## Triangulate the PSLG with triangles in which no angle
    ## is smaller than 20 degrees
    tA <- triangulate(A, q=20)
    plot(tA)
    ## Triangulate the PSLG with triangles in which no triangle has
    ## area greater than 0.001
    tA <- triangulate(A, a=0.001)
    plot(tA)

[Package RTriangle version 1.6-0.10 Index]