@math.gl/geometry 5.0.0-alpha.2 → 5.0.0-alpha.20

package/dist/polygon/cut-by-mercator-bounds.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"cut-by-mercator-bounds.d.ts","sourceRoot":"","sources":["../../src/polygon/cut-by-mercator-bounds.ts"],"names":[],"mappings":"AAMA,OAAO,KAAK,EAAC,OAAO,EAAC,yBAAsB;AAC3C,OAAO,KAAK,EAAC,YAAY,EAAC,MAAM,eAAe,CAAC;AAKhD,0GAA0G;AAC1G,wBAAgB,2BAA2B,CACzC,SAAS,EAAE,QAAQ,CAAC,YAAY,CAAC,EACjC,OAAO,CAAC,EAAE;IACR,IAAI,CAAC,EAAE,MAAM,CAAC;IACd,UAAU,CAAC,EAAE,MAAM,CAAC;IACpB,QAAQ,CAAC,EAAE,MAAM,CAAC;IAClB,SAAS,CAAC,EAAE,OAAO,CAAC;CACrB,GACA,MAAM,EAAE,EAAE,CAsBZ;AAED,0GAA0G;AAC1G,wBAAgB,0BAA0B,CACxC,SAAS,EAAE,QAAQ,CAAC,YAAY,CAAC,EACjC,WAAW,GAAE,QAAQ,CAAC,YAAY,CAAC,GAAG,IAAW,EACjD,OAAO,CAAC,EAAE;IACR,IAAI,CAAC,EAAE,MAAM,CAAC;IACd,SAAS,CAAC,EAAE,OAAO,CAAC;IACpB,WAAW,CAAC,EAAE,MAAM,CAAC;IACrB,SAAS,CAAC,EAAE,OAAO,CAAC;CACrB,GACA,OAAO,EAAE,CAoDX"}

package/dist/polygon/cut-by-mercator-bounds.js

@@ -0,0 +1,149 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+import { cutPolylineByGrid, cutPolygonByGrid } from "./cut-by-grid.js";
+import { getPointAtIndex, push } from "./utils.js";
+// https://en.wikipedia.org/wiki/Web_Mercator_projection
+const DEFAULT_MAX_LATITUDE = 85.051129;
+/** https://user-images.githubusercontent.com/2059298/78465769-938b7a00-76ae-11ea-9b95-1f4c26425ab9.png */
+export function cutPolylineByMercatorBounds(positions, options) {
+    const { size = 2, startIndex = 0, endIndex = positions.length, normalize = true } = options || {};
+    // Remap longitudes so that each segment takes the shorter path
+    const newPositions = positions.slice(startIndex, endIndex);
+    wrapLongitudesForShortestPath(newPositions, size, 0, endIndex - startIndex);
+    const parts = cutPolylineByGrid(newPositions, {
+        size,
+        broken: true,
+        gridResolution: 360,
+        gridOffset: [-180, -180]
+    });
+    if (normalize) {
+        // Each part is guaranteed to be in a single copy of the world
+        // Map longitudes back to [-180, 180]
+        for (const part of parts) {
+            shiftLongitudesIntoRange(part, size);
+        }
+    }
+    return parts;
+}
+/** https://user-images.githubusercontent.com/2059298/78465770-94241080-76ae-11ea-809a-6a8534dac1d9.png */
+export function cutPolygonByMercatorBounds(positions, holeIndices = null, options) {
+    const { size = 2, normalize = true, edgeTypes = false } = options || {};
+    holeIndices = holeIndices || [];
+    const newPositions = [];
+    const newHoleIndices = [];
+    let srcStartIndex = 0;
+    let targetIndex = 0;
+    for (let ringIndex = 0; ringIndex <= holeIndices.length; ringIndex++) {
+        // srcStartIndex/srcEndIndex define the ring in the original positions
+        const srcEndIndex = holeIndices[ringIndex] || positions.length;
+        // targetStartIndex/targetIndex define the ring in newPositions
+        const targetStartIndex = targetIndex;
+        // In case the ring contains a pole (e.g. Antarctica), we'll have to insert vertices
+        // The insertion point is defined by the vertex closest to the pole
+        // Split the the ring by the insertion point when copying to newPositions
+        const splitIndex = findSplitIndex(positions, size, srcStartIndex, srcEndIndex);
+        for (let i = splitIndex; i < srcEndIndex; i++) {
+            newPositions[targetIndex++] = positions[i];
+        }
+        for (let i = srcStartIndex; i < splitIndex; i++) {
+            newPositions[targetIndex++] = positions[i];
+        }
+        // Remap longitudes so that each segment takes the shorter path
+        wrapLongitudesForShortestPath(newPositions, size, targetStartIndex, targetIndex);
+        // Handle the case when the ring contains a pole
+        insertPoleVertices(newPositions, size, targetStartIndex, targetIndex, options?.maxLatitude);
+        srcStartIndex = srcEndIndex;
+        newHoleIndices[ringIndex] = targetIndex;
+    }
+    newHoleIndices.pop();
+    const parts = cutPolygonByGrid(newPositions, newHoleIndices, {
+        size,
+        gridResolution: 360,
+        gridOffset: [-180, -180],
+        edgeTypes
+    });
+    if (normalize) {
+        // Each part is guaranteed to be in a single copy of the world
+        // Map longitudes back to [-180, 180]
+        for (const part of parts) {
+            // May mutate newPositions, which is created by us.
+            shiftLongitudesIntoRange(part.positions, size);
+        }
+    }
+    return parts;
+}
+/* Helpers */
+// See comments for insertPoleVertices
+function findSplitIndex(positions, size, startIndex, endIndex) {
+    let maxLat = -1;
+    let pointIndex = -1;
+    for (let i = startIndex + 1; i < endIndex; i += size) {
+        const lat = Math.abs(positions[i]);
+        if (lat > maxLat) {
+            maxLat = lat;
+            pointIndex = i - 1;
+        }
+    }
+    return pointIndex;
+}
+// https://user-images.githubusercontent.com/2059298/78857483-5987e400-79de-11ea-98fc-0631287a8431.png
+//
+// If the polygon contains a pole, to tesselate it correctly, we need to insert the edge
+// of map into the polygon. This requires adding two vertices that represent the pole, by
+// drawing a perpendicular line to the Mercator map edge from a selected vertex on the ring.
+//
+// We select the insertion position carefully so that the inserted line segments do not
+// intersect with the ring itself. This is ensured by findSplitIndex, which returns the
+// vertex closest to the pole.
+function insertPoleVertices(positions, size, startIndex, endIndex, maxLatitude = DEFAULT_MAX_LATITUDE) {
+    // Check if the ring contains a pole
+    const firstLng = positions[startIndex];
+    const lastLng = positions[endIndex - size];
+    if (Math.abs(firstLng - lastLng) > 180) {
+        // The ring does not make a round trip
+        // Add the nearest pole to the vertices so that the polygon tesselates correctly
+        const p = getPointAtIndex(positions, 0, size, startIndex);
+        // Copy the first vertex to the world of the last vertex
+        p[0] += Math.round((lastLng - firstLng) / 360) * 360;
+        push(positions, p);
+        // Project the copied vertex to the edge of the map
+        p[1] = Math.sign(p[1]) * maxLatitude;
+        push(positions, p);
+        // Project the first vertex to the edge of the map
+        p[0] = firstLng;
+        push(positions, p);
+    }
+}
+function wrapLongitudesForShortestPath(positions, size, startIndex, endIndex) {
+    let prevLng = positions[0];
+    let lng;
+    for (let i = startIndex; i < endIndex; i += size) {
+        lng = positions[i];
+        const delta = lng - prevLng;
+        if (delta > 180 || delta < -180) {
+            lng -= Math.round(delta / 360) * 360;
+        }
+        positions[i] = prevLng = lng;
+    }
+}
+function shiftLongitudesIntoRange(positions, size) {
+    let refLng;
+    const pointCount = positions.length / size;
+    // Find a longitude that is not on the edge of a world
+    // Which we will use to determine which world copy it is
+    for (let i = 0; i < pointCount; i++) {
+        refLng = positions[i * size];
+        if ((refLng + 180) % 360 !== 0) {
+            break;
+        }
+    }
+    const delta = -Math.round(refLng / 360) * 360;
+    if (delta === 0) {
+        return;
+    }
+    for (let i = 0; i < pointCount; i++) {
+        positions[i * size] += delta;
+    }
+}
+//# sourceMappingURL=cut-by-mercator-bounds.js.map

package/dist/polygon/cut-by-mercator-bounds.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"cut-by-mercator-bounds.js","sourceRoot":"","sources":["../../src/polygon/cut-by-mercator-bounds.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAE5D,OAAO,EAAC,iBAAiB,EAAE,gBAAgB,EAAC,yBAAsB;AAClE,OAAO,EAAC,eAAe,EAAE,IAAI,EAAC,mBAAgB;AAI9C,wDAAwD;AACxD,MAAM,oBAAoB,GAAG,SAAS,CAAC;AAEvC,0GAA0G;AAC1G,MAAM,UAAU,2BAA2B,CACzC,SAAiC,EACjC,OAKC;IAED,MAAM,EAAC,IAAI,GAAG,CAAC,EAAE,UAAU,GAAG,CAAC,EAAE,QAAQ,GAAG,SAAS,CAAC,MAAM,EAAE,SAAS,GAAG,IAAI,EAAC,GAAG,OAAO,IAAI,EAAE,CAAC;IAEhG,+DAA+D;IAC/D,MAAM,YAAY,GAAG,SAAS,CAAC,KAAK,CAAC,UAAU,EAAE,QAAQ,CAAC,CAAC;IAC3D,6BAA6B,CAAC,YAAY,EAAE,IAAI,EAAE,CAAC,EAAE,QAAQ,GAAG,UAAU,CAAC,CAAC;IAE5E,MAAM,KAAK,GAAG,iBAAiB,CAAC,YAAY,EAAE;QAC5C,IAAI;QACJ,MAAM,EAAE,IAAI;QACZ,cAAc,EAAE,GAAG;QACnB,UAAU,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC;KACzB,CAAe,CAAC;IAEjB,IAAI,SAAS,EAAE,CAAC;QACd,8DAA8D;QAC9D,qCAAqC;QACrC,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;YACzB,wBAAwB,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;QACvC,CAAC;IACH,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,0GAA0G;AAC1G,MAAM,UAAU,0BAA0B,CACxC,SAAiC,EACjC,cAA6C,IAAI,EACjD,OAKC;IAED,MAAM,EAAC,IAAI,GAAG,CAAC,EAAE,SAAS,GAAG,IAAI,EAAE,SAAS,GAAG,KAAK,EAAC,GAAG,OAAO,IAAI,EAAE,CAAC;IACtE,WAAW,GAAG,WAAW,IAAI,EAAE,CAAC;IAChC,MAAM,YAAY,GAAa,EAAE,CAAC;IAClC,MAAM,cAAc,GAAa,EAAE,CAAC;IACpC,IAAI,aAAa,GAAG,CAAC,CAAC;IACtB,IAAI,WAAW,GAAG,CAAC,CAAC;IAEpB,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,IAAI,WAAW,CAAC,MAAM,EAAE,SAAS,EAAE,EAAE,CAAC;QACrE,sEAAsE;QACtE,MAAM,WAAW,GAAG,WAAW,CAAC,SAAS,CAAC,IAAI,SAAS,CAAC,MAAM,CAAC;QAC/D,+DAA+D;QAC/D,MAAM,gBAAgB,GAAG,WAAW,CAAC;QAErC,oFAAoF;QACpF,mEAAmE;QACnE,yEAAyE;QACzE,MAAM,UAAU,GAAG,cAAc,CAAC,SAAS,EAAE,IAAI,EAAE,aAAa,EAAE,WAAW,CAAC,CAAC;QAC/E,KAAK,IAAI,CAAC,GAAG,UAAU,EAAE,CAAC,GAAG,WAAW,EAAE,CAAC,EAAE,EAAE,CAAC;YAC9C,YAAY,CAAC,WAAW,EAAE,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;QAC7C,CAAC;QACD,KAAK,IAAI,CAAC,GAAG,aAAa,EAAE,CAAC,GAAG,UAAU,EAAE,CAAC,EAAE,EAAE,CAAC;YAChD,YAAY,CAAC,WAAW,EAAE,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;QAC7C,CAAC;QAED,+DAA+D;QAC/D,6BAA6B,CAAC,YAAY,EAAE,IAAI,EAAE,gBAAgB,EAAE,WAAW,CAAC,CAAC;QAEjF,gDAAgD;QAChD,kBAAkB,CAAC,YAAY,EAAE,IAAI,EAAE,gBAAgB,EAAE,WAAW,EAAE,OAAO,EAAE,WAAW,CAAC,CAAC;QAE5F,aAAa,GAAG,WAAW,CAAC;QAC5B,cAAc,CAAC,SAAS,CAAC,GAAG,WAAW,CAAC;IAC1C,CAAC;IACD,cAAc,CAAC,GAAG,EAAE,CAAC;IAErB,MAAM,KAAK,GAAG,gBAAgB,CAAC,YAAY,EAAE,cAAc,EAAE;QAC3D,IAAI;QACJ,cAAc,EAAE,GAAG;QACnB,UAAU,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC;QACxB,SAAS;KACV,CAAC,CAAC;IAEH,IAAI,SAAS,EAAE,CAAC;QACd,8DAA8D;QAC9D,qCAAqC;QACrC,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;YACzB,mDAAmD;YACnD,wBAAwB,CAAC,IAAI,CAAC,SAAyB,EAAE,IAAI,CAAC,CAAC;QACjE,CAAC;IACH,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,aAAa;AAEb,sCAAsC;AACtC,SAAS,cAAc,CACrB,SAAiC,EACjC,IAAY,EACZ,UAAkB,EAClB,QAAgB;IAEhB,IAAI,MAAM,GAAG,CAAC,CAAC,CAAC;IAChB,IAAI,UAAU,GAAG,CAAC,CAAC,CAAC;IACpB,KAAK,IAAI,CAAC,GAAG,UAAU,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,IAAI,IAAI,EAAE,CAAC;QACrD,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC;QACnC,IAAI,GAAG,GAAG,MAAM,EAAE,CAAC;YACjB,MAAM,GAAG,GAAG,CAAC;YACb,UAAU,GAAG,CAAC,GAAG,CAAC,CAAC;QACrB,CAAC;IACH,CAAC;IACD,OAAO,UAAU,CAAC;AACpB,CAAC;AAED,sGAAsG;AACtG,EAAE;AACF,wFAAwF;AACxF,yFAAyF;AACzF,4FAA4F;AAC5F,EAAE;AACF,uFAAuF;AACvF,uFAAuF;AACvF,8BAA8B;AAC9B,SAAS,kBAAkB,CACzB,SAAmB,EACnB,IAAY,EACZ,UAAkB,EAClB,QAAgB,EAChB,cAAsB,oBAAoB;IAE1C,oCAAoC;IACpC,MAAM,QAAQ,GAAG,SAAS,CAAC,UAAU,CAAC,CAAC;IACvC,MAAM,OAAO,GAAG,SAAS,CAAC,QAAQ,GAAG,IAAI,CAAC,CAAC;IAC3C,IAAI,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,OAAO,CAAC,GAAG,GAAG,EAAE,CAAC;QACvC,sCAAsC;QACtC,gFAAgF;QAChF,MAAM,CAAC,GAAG,eAAe,CAAC,SAAS,EAAE,CAAC,EAAE,IAAI,EAAE,UAAU,CAAC,CAAC;QAC1D,wDAAwD;QACxD,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,CAAC,OAAO,GAAG,QAAQ,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;QACrD,IAAI,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;QACnB,mDAAmD;QACnD,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC;QACrC,IAAI,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;QACnB,kDAAkD;QAClD,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC;QAChB,IAAI,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;IACrB,CAAC;AACH,CAAC;AAED,SAAS,6BAA6B,CACpC,SAAuB,EACvB,IAAY,EACZ,UAAkB,EAClB,QAAgB;IAEhB,IAAI,OAAO,GAAW,SAAS,CAAC,CAAC,CAAC,CAAC;IACnC,IAAI,GAAW,CAAC;IAChB,KAAK,IAAI,CAAC,GAAG,UAAU,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,IAAI,IAAI,EAAE,CAAC;QACjD,GAAG,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;QACnB,MAAM,KAAK,GAAG,GAAG,GAAG,OAAO,CAAC;QAC5B,IAAI,KAAK,GAAG,GAAG,IAAI,KAAK,GAAG,CAAC,GAAG,EAAE,CAAC;YAChC,GAAG,IAAI,IAAI,CAAC,KAAK,CAAC,KAAK,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;QACvC,CAAC;QACD,SAAS,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,GAAG,CAAC;IAC/B,CAAC;AACH,CAAC;AAED,SAAS,wBAAwB,CAAC,SAAuB,EAAE,IAAY;IACrE,IAAI,MAAc,CAAC;IACnB,MAAM,UAAU,GAAG,SAAS,CAAC,MAAM,GAAG,IAAI,CAAC;IAE3C,sDAAsD;IACtD,wDAAwD;IACxD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,EAAE,CAAC,EAAE,EAAE,CAAC;QACpC,MAAM,GAAG,SAAS,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;QAC7B,IAAI,CAAC,MAAM,GAAG,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,EAAE,CAAC;YAC/B,MAAM;QACR,CAAC;IACH,CAAC;IAED,MAAM,KAAK,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,MAAM,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;IAC9C,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;QAChB,OAAO;IACT,CAAC;IACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,EAAE,CAAC,EAAE,EAAE,CAAC;QACpC,SAAS,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,KAAK,CAAC;IAC/B,CAAC;AACH,CAAC"}

package/dist/polygon/earcut.d.ts

@@ -0,0 +1,13 @@
+import type { NumericArray } from '@math.gl/core';
+import { Plane2D } from "./polygon-utils.js";
+/**
+ * Computes a triangulation of a polygon
+ * @param positions a flat array of the vertex positions that define the polygon.
+ * @param holeIndices an array of hole indices if any (e.g. [5, 8] for a 12-vertex input would mean one hole with vertices 5–7 and another with 8–11).
+ * @param dim the number of elements in each vertex. Size `2` will interpret `positions` as `[x0, y0, x1, y1, ...]` and size `3` will interpret `positions` as `[x0, y0, z0, x1, y1, z1, ...]`. Default `2`.
+ * @param areas areas of outer polygon and holes as computed by `getPolygonSignedArea()`. Can be optionally supplied to speed up triangulation
+ * @returns array of indices into the `positions` array that describes the triangulation of the polygon
+ * Adapted from https://github.com/mapbox/earcut
+ */
+export declare function earcut(positions: NumericArray, holeIndices?: NumericArray, dim?: number, areas?: NumericArray, plane?: Plane2D): number[];
+//# sourceMappingURL=earcut.d.ts.map

package/dist/polygon/earcut.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"earcut.d.ts","sourceRoot":"","sources":["../../src/polygon/earcut.ts"],"names":[],"mappings":"AAkCA,OAAO,KAAK,EAAC,YAAY,EAAC,MAAM,eAAe,CAAC;AAChD,OAAO,EAAiC,OAAO,EAAC,2BAAwB;AAExE;;;;;;;;GAQG;AACH,wBAAgB,MAAM,CACpB,SAAS,EAAE,YAAY,EACvB,WAAW,CAAC,EAAE,YAAY,EAC1B,GAAG,GAAE,MAAU,EACf,KAAK,CAAC,EAAE,YAAY,EACpB,KAAK,GAAE,OAAc,GACpB,MAAM,EAAE,CAwCV"}