Skip to content

Repository files navigation

Martinez-Rueda polygon clipping algorithm npm version

screenshot 2016-07-26 10 54 01 screenshot 2016-07-25 18 53 44

Details

The algorithm is specifically fast and capable of working with polygons of all types: multipolygons (without cascading), polygons with holes, self-intersecting polygons and degenerate polygons with overlapping edges.

Example

Play with it by forking this Codepen

import * as martinez from 'martinez-polygon-clipping';
const gj1 = { "type": "Feature", ..., "geometry": { "type": "Polygon", "coordinates": [ [ [x, y], ... ] ]};
const gj2 = { "type": "Feature", ..., "geometry": { "type": "MultiPolygon", "coordinates": [ [ [ [x, y], ...] ] ]};

const intersection = {
  "type": "Feature",
  "properties": { ... },
  "geometry": {
    "type": "Polygon",
    "coordinates": martinez.intersection(gj1.geometry.coordinates, gj2.geometry.coordinates)
  }
};

API

  • .intersection(<Geometry>, <Geometry>) => <MultiPolygon>
  • .union(<Geometry>, <Geometry>) => <MultiPolygon>
  • .diff(<Geometry>, <Geometry>) => <MultiPolygon>
  • .xor(<Geometry>, <Geometry>) => <MultiPolygon>

<Geometry> is GeoJSON 'Polygon' or 'MultiPolygon' coordinates structure.

The result is always 'MultiPolygon' coordinates, [] if it is empty. Its rings are closed and oriented as RFC 7946 requires: exterior rings counter-clockwise, holes clockwise.

<Operation> is an enum of { INTERSECTION: 0, UNION: 1, DIFFERENCE: 2, XOR: 3 } in case you have to decide programmatically which operation do you need

Benchmarks

Operations per second, higher is better. Run npm run bench to reproduce: it benchmarks the built bundle against JSTS, polygon-clipping and polyclip-ts (the engine behind @turf/union 7). All libraries produce the same result areas on these inputs.

Benchmark Martinez JSTS 2.12 polygon-clipping 0.15 polyclip-ts 0.16
Hole_Hole union (20 vertices) 68,058 10,549 28,687 1,712
Asia union (28k vertices) 46.8 32.2 19.2 3.18
States union (2.3k vertices) 754 426 399 40.2
Asia vs Asia shifted 0.05°: union 16.2 5.29 6.06 0.97
Asia vs Asia shifted 0.05°: difference 16.0 5.25 5.11 0.72

Apple M3, Node 22.17. The "shifted" cases clip the Asia polygon against a copy of itself moved slightly east, so that nearly every edge intersects; demo/cases.html shows it with a slider for the shift.

Asia minus a copy of itself shifted 0.3° east: land in grey, the difference in red

The difference with a larger 0.3° shift, to make the slivers visible (the benchmark uses 0.05°).

Features

The algorithm of Martinez et al. was extended to work with multipolygons without cascading.

Demos

npm run dev starts the demos:

  • demo/index.html: interactive map where you can edit the polygons.
  • demo/cases.html: test case viewer. Pick any case from test/genericTestCases or test/fixtures and an operation to see the inputs, the computed result and, for generic test cases, whether it matches the expected result. Use ←/→ to step through cases, keys 1–5 to switch the operation, and F to fit the view. The selection is kept in the URL, so a case can be linked, e.g. cases.html#generic%2Fissue155/xor.

Authors

Based on

Related projects

Other JavaScript implementations of the Martinez–Rueda–Feito algorithm:

  • polygon-clipping by Mike Fogel was forked from this repository in February 2018 (see its license) and developed separately since. It snaps coordinates and intersection points to previously seen values within floating-point precision, and caps the sizes of its internal structures as a guard against infinite loops. (December 2023).
  • polyclip-ts by Luiz Barboza is a TypeScript fork of polygon-clipping, and so, indirectly, of this repository. It computes with arbitrary-precision decimals (bignumber.js) and is the engine behind @turf/union and the other Turf 7 boolean operations.

See Benchmarks for how they compare in speed on the same inputs.

License

The MIT License (MIT)

Copyright (c) 2026 Alexander Milevski

Permission is hereby granted, free of charge, to any person obtaining a copy of this software and associated documentation files (the "Software"), to deal in the Software without restriction, including without limitation the rights to use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of the Software, and to permit persons to whom the Software is furnished to do so, subject to the following conditions:

The above copyright notice and this permission notice shall be included in all copies or substantial portions of the Software.

THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.

About

Martinez-Rueda polygon clipping algorithm, does boolean operation on polygons (multipolygons, polygons with holes etc): intersection, union, difference, xor

Topics

Resources

Stars

775 stars

Watchers

17 watching

Forks

Releases

Sponsor this project

Packages

Used by

Contributors

Languages