cdt2d

A robust 2D constrained Delaunay triangulation library written in JavaScript.

WORK IN PROGRESS

Demo

To test out this module, you can open up a demo in your browser with the following link:

cdt2d demo

• Click to add points
• Click on a point to remove it
• Drag one point onto another to add an edge constraint
• Click on a green edge to remove a constraint
• Toggle options by clicking on the checkboxes on the left
• Click reset to clear all points

Examples

Simple example

Here is a simple example showing how to invoke cdt2d:

//First we need to reqire the module 
var cdt2d = require('cdt2d')
 
//Then we define a list of points, represented as pairs of x,y coordinates 
var points = [
  [-2,-2],
  [-2, 2],
  [ 2, 2],
  [ 2,-2],
  [ 1, 0],
  [ 0, 1],
  [-1, 0],
  [ 0,-1]
]
 
//Next we can optionally define some edge constraints 
// This set of edges determines a pair of loops 
var edges = [
 //Outer loop 
 [0, 1],
 [1, 2],
 [2, 3],
 [3, 0],
 
 //Inner loop 
 [4, 5],
 [5, 6],
 [6, 7],
 [7, 4]
]
 
//Finally we call cdt2d with the points and edges 
// The flag {exterior: false} tells  it to remove exterior faces 
console.log(cdt2d(points, edges, {exterior: false}))

Output

The above program will output the following triangles:

[ [ 0, 3, 7 ],
  [ 0, 6, 1 ],
  [ 0, 7, 6 ],
  [ 1, 5, 2 ],
  [ 1, 6, 5 ],
  [ 2, 4, 3 ],
  [ 2, 5, 4 ],
  [ 3, 4, 7 ] ]

Each triangle is represented as an array of 3 indices of points. We can visualize this data in the following figure:

Messy graphs

If your input doesn't satisfy the validity invariants (ie no self intersections, duplicate vertices or t-junctions), then you will need to preprocess it to clean it up. One way to do this is with the clean-pslg module. Here is an example showing how to do this:

var cleanPSLG = require('clean-pslg')
var cdt2d = require('cdt2d')
 
var points = [
  [-1, 0],
  [ 1, 0],
  [ 0,-1],
  [ 0, 1]
]
 
var edges = [
  [0, 1],
  [1, 2]
]
 
//This updates points/edges so that they now form a valid PSLG 
cleanPSLG(points, edges)
 
//Generate the triangulation 
console.log({
  points: points,
  edges: edges,
  triangles: cdt2d(points, edges)
})

Output

TODO

Polygon example

It is also pretty easy to use this module with polygons, as one would get from a GeoJSON file. To do this, it is first necessary to convert them into a planar straight line graph. This can be done using the poly-to-pslg module:

var toPSLG = require('poly-to-pslg')
var cdt2d = require('cdt2d')
 

TODO

Polygon with holes example

The above procedure even works if the polygons have holes:

var toPSLG = require('poly-to-pslg')
var cdt2d = require('cdt2d')
 

TODO

Delaunay triangulation

You can also use cdt2d to generate Delaunay triangulations of arbitrary point sets in the plane:

TODO

Install

This module works in any modern CommonJS environment. You can install it using npm with the following command:

npm i cdt2d

You should be able to then use it in node or on the web with browserify.

API

var cells = require('cdt2d')(points[, edges, options])

Constructs a constrained Delaunay triangulation of a planar straight-line graph.

• points are the vertices of the triangulation, represented by pairs of numbers.
• edges is an optional list of edge constraints which must occur within the triangulation. These constraints are given by pairs of indices of points. If not specified, then no constraints are used.
• options is an object which takes some optional parameters.
  • delaunay if this flag is set to true, then the resulting triangulation is converted to a Delaunay triangulation by edge flipping. Otherwise if it is false, then an arbitrary triangulation is returned. (Default true)
  • interior if set, only return interior faces. See note. (Default true)
  • exterior if set, only return exterior faces. See note. (Default true)
  • infinity if set, then the triangulation is augmented with a point at infinity represented by the index -1. (Default false)

Returns A list of all triangles represented as triples of indices of vertices

Note on interior/exterior classification Interior/exterior faces are classified by treating the constraint edges as the boundary and traversing the triangulation. The point at infinity is in the exterior of the set, and other faces are classified by the parity of the path with fewest crossings from the face to the point at infinity.

Assumptions This module makes the following assumptions about the points and edge constraints:

• No point in the input is duplicated
• No pair of edge constraints cross in their relative interior
• No point is contained in the relative interior of an edge (ie no T-junctions)

If your input does not satisfy these conditions, you will need to preprocess it first (using clean-pslg for example) otherwise cdt2d may return incorrect results.

Limitations Currently there is no way to specify that only some edge constraints are to be included in the boundary. It is also not possible to add a constraint from a vertex to the point at infinity. If there is enough demand I may add these features or perhaps create a separate module.

Benchmarks and comparisons

Assertion: cdt2d is the only non-broken triangulation library in JavaScript.

• TODO Catalogue failing cases for other libraries
• TODO Need to measure performance and finetune

Libraries to compare against:

• earcut
• poly2tri
• pnltri
• libtess.js

License

(c) 2015 Mikola Lysenko. MIT License