@loaders.gl/mvt 4.2.0-alpha.3 → 4.2.0-alpha.5

package/dist/lib/geojson-tiler/feature.js

@@ -1,35 +1,44 @@
+// loaders.gl
+// SPDX-License-Identifier: MIT
+// Copyright (c) vis.gl contributors
+// Forked from https://github.com/mapbox/geojson-vt under compatible ISC license
 export function createFeature(id, type, geom, tags) {
-  const feature = {
-    id: id == null ? null : id,
-    type,
-    geometry: geom,
-    tags,
-    minX: Infinity,
-    minY: Infinity,
-    maxX: -Infinity,
-    maxY: -Infinity
-  };
-  if (type === 'Point' || type === 'MultiPoint' || type === 'LineString') {
-    calcLineBBox(feature, geom);
-  } else if (type === 'Polygon') {
-    calcLineBBox(feature, geom[0]);
-  } else if (type === 'MultiLineString') {
-    for (const line of geom) {
-      calcLineBBox(feature, line);
+    const feature = {
+        // eslint-disable-next-line
+        id: id == null ? null : id,
+        type,
+        geometry: geom,
+        tags,
+        minX: Infinity,
+        minY: Infinity,
+        maxX: -Infinity,
+        maxY: -Infinity
+    };
+    if (type === 'Point' || type === 'MultiPoint' || type === 'LineString') {
+        calcLineBBox(feature, geom);
     }
-  } else if (type === 'MultiPolygon') {
-    for (const polygon of geom) {
-      calcLineBBox(feature, polygon[0]);
+    else if (type === 'Polygon') {
+        // the outer ring (ie [0]) contains all inner rings
+        calcLineBBox(feature, geom[0]);
     }
-  }
-  return feature;
+    else if (type === 'MultiLineString') {
+        for (const line of geom) {
+            calcLineBBox(feature, line);
+        }
+    }
+    else if (type === 'MultiPolygon') {
+        for (const polygon of geom) {
+            // the outer ring (ie [0]) contains all inner rings
+            calcLineBBox(feature, polygon[0]);
+        }
+    }
+    return feature;
 }
 function calcLineBBox(feature, geom) {
-  for (let i = 0; i < geom.length; i += 3) {
-    feature.minX = Math.min(feature.minX, geom[i]);
-    feature.minY = Math.min(feature.minY, geom[i + 1]);
-    feature.maxX = Math.max(feature.maxX, geom[i]);
-    feature.maxY = Math.max(feature.maxY, geom[i + 1]);
-  }
+    for (let i = 0; i < geom.length; i += 3) {
+        feature.minX = Math.min(feature.minX, geom[i]);
+        feature.minY = Math.min(feature.minY, geom[i + 1]);
+        feature.maxX = Math.max(feature.maxX, geom[i]);
+        feature.maxY = Math.max(feature.maxY, geom[i + 1]);
+    }
 }
-//# sourceMappingURL=feature.js.map

package/dist/lib/geojson-tiler/geojson-tiler.d.ts

@@ -1,4 +1,4 @@
-import type { GeoJSONTile, GeoJSONTileFeature } from './tile';
+import type { GeoJSONTile, GeoJSONTileFeature } from "./tile.js";
 /** Options to configure tiling */
 export type GeoJSONTilerOptions = {
     maxZoom?: number /** max zoom to preserve detail on */;

package/dist/lib/geojson-tiler/geojson-tiler.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"geojson-tiler.d.ts","sourceRoot":"","sources":["../../../src/lib/geojson-tiler/geojson-tiler.ts"],"names":[],"mappings":"AAOA,OAAO,KAAK,EAAC,WAAW,EAAE,kBAAkB,EAAC,MAAM,QAAQ,CAAC;AAQ5D,kCAAkC;AAClC,MAAM,MAAM,mBAAmB,GAAG;IAChC,OAAO,CAAC,EAAE,MAAM,CAAC,qCAAqC,CAAC;IACvD,YAAY,CAAC,EAAE,MAAM,CAAC,iCAAiC,CAAC;IACxD,cAAc,CAAC,EAAE,MAAM,CAAC,sDAAsD,CAAC;IAC/E,SAAS,CAAC,EAAE,MAAM,CAAC,sDAAsD,CAAC;IAC1E,MAAM,CAAC,EAAE,MAAM,CAAC,kBAAkB,CAAC;IACnC,MAAM,CAAC,EAAE,MAAM,CAAC,+BAA+B,CAAC;IAChD,WAAW,CAAC,EAAE,OAAO,CAAC,wCAAwC,CAAC;IAC/D,SAAS,CAAC,EAAE,MAAM,CAAC,8DAA8D,CAAC;IAClF,UAAU,CAAC,EAAE,OAAO,CAAC,qEAAqE,CAAC;IAC3F,KAAK,CAAC,EAAE,MAAM,CAAC,gCAAgC,CAAC;CACjD,CAAC;AAgBF,qBAAa,YAAY;IACvB,OAAO,EAAE,QAAQ,CAAC,mBAAmB,CAAC,CAAC;IAGvC,KAAK,EAAE,MAAM,CAAC,MAAM,EAAE,WAAW,CAAC,CAAM;IACxC,UAAU,EAAE;QAAC,CAAC,EAAE,MAAM,CAAC;QAAC,CAAC,EAAE,MAAM,CAAC;QAAC,CAAC,EAAE,MAAM,CAAA;KAAC,EAAE,CAAM;IAErD,KAAK,EAAE,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAM;IACnC,KAAK,EAAE,MAAM,CAAK;gBAEN,IAAI,KAAA,EAAE,OAAO,CAAC,EAAE,mBAAmB;IA6C/C;;;;;;OAMG;IAEH,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,WAAW,GAAG,IAAI;IAkD5D;;;;;;;OAOG;IAEH,SAAS,CACP,QAAQ,EAAE,kBAAkB,EAAE,EAC9B,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,EAAE,CAAC,EAAE,MAAM,EACX,EAAE,CAAC,EAAE,MAAM,EACX,EAAE,CAAC,EAAE,MAAM,GACV,IAAI;CAyGR"}
+{"version":3,"file":"geojson-tiler.d.ts","sourceRoot":"","sources":["../../../src/lib/geojson-tiler/geojson-tiler.ts"],"names":[],"mappings":"AAOA,OAAO,KAAK,EAAC,WAAW,EAAE,kBAAkB,EAAC,kBAAe;AAQ5D,kCAAkC;AAClC,MAAM,MAAM,mBAAmB,GAAG;IAChC,OAAO,CAAC,EAAE,MAAM,CAAC,qCAAqC,CAAC;IACvD,YAAY,CAAC,EAAE,MAAM,CAAC,iCAAiC,CAAC;IACxD,cAAc,CAAC,EAAE,MAAM,CAAC,sDAAsD,CAAC;IAC/E,SAAS,CAAC,EAAE,MAAM,CAAC,sDAAsD,CAAC;IAC1E,MAAM,CAAC,EAAE,MAAM,CAAC,kBAAkB,CAAC;IACnC,MAAM,CAAC,EAAE,MAAM,CAAC,+BAA+B,CAAC;IAChD,WAAW,CAAC,EAAE,OAAO,CAAC,wCAAwC,CAAC;IAC/D,SAAS,CAAC,EAAE,MAAM,CAAC,8DAA8D,CAAC;IAClF,UAAU,CAAC,EAAE,OAAO,CAAC,qEAAqE,CAAC;IAC3F,KAAK,CAAC,EAAE,MAAM,CAAC,gCAAgC,CAAC;CACjD,CAAC;AAgBF,qBAAa,YAAY;IACvB,OAAO,EAAE,QAAQ,CAAC,mBAAmB,CAAC,CAAC;IAGvC,KAAK,EAAE,MAAM,CAAC,MAAM,EAAE,WAAW,CAAC,CAAM;IACxC,UAAU,EAAE;QAAC,CAAC,EAAE,MAAM,CAAC;QAAC,CAAC,EAAE,MAAM,CAAC;QAAC,CAAC,EAAE,MAAM,CAAA;KAAC,EAAE,CAAM;IAErD,KAAK,EAAE,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAM;IACnC,KAAK,EAAE,MAAM,CAAK;gBAEN,IAAI,KAAA,EAAE,OAAO,CAAC,EAAE,mBAAmB;IA6C/C;;;;;;OAMG;IAEH,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,WAAW,GAAG,IAAI;IAkD5D;;;;;;;OAOG;IAEH,SAAS,CACP,QAAQ,EAAE,kBAAkB,EAAE,EAC9B,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,EAAE,CAAC,EAAE,MAAM,EACX,EAAE,CAAC,EAAE,MAAM,EACX,EAAE,CAAC,EAAE,MAAM,GACV,IAAI;CAyGR"}

package/dist/lib/geojson-tiler/geojson-tiler.js

@@ -1,170 +1,208 @@
-import { convert } from "./convert.js";
-import { clip } from "./clip.js";
-import { wrap } from "./wrap.js";
-import { transformTile } from "./transform.js";
-import { createTile } from "./tile.js";
+// loaders.gl
+// SPDX-License-Identifier: MIT
+// Copyright (c) vis.gl contributors
+// Forked from https://github.com/mapbox/geojson-vt under compatible ISC license
+import { convert } from "./convert.js"; // GeoJSON conversion and preprocessing
+import { clip } from "./clip.js"; // stripe clipping algorithm
+import { wrap } from "./wrap.js"; // date line processing
+import { transformTile } from "./transform.js"; // coordinate transformation
+import { createTile } from "./tile.js"; // final simplified tile generation
 const DEFAULT_OPTIONS = {
-  maxZoom: 14,
-  indexMaxZoom: 5,
-  indexMaxPoints: 100000,
-  tolerance: 3,
-  extent: 4096,
-  buffer: 64,
-  lineMetrics: false,
-  promoteId: undefined,
-  generateId: false,
-  debug: 0
+    maxZoom: 14, // max zoom to preserve detail on
+    indexMaxZoom: 5, // max zoom in the tile index
+    indexMaxPoints: 100000, // max number of points per tile in the tile index
+    tolerance: 3, // simplification tolerance (higher means simpler)
+    extent: 4096, // tile extent
+    buffer: 64, // tile buffer on each side
+    lineMetrics: false, // whether to calculate line metrics
+    // @ts-expect-error Ensures all these required params have defaults
+    promoteId: undefined, // name of a feature property to be promoted to feature.id
+    generateId: false, // whether to generate feature ids. Cannot be used with promoteId
+    debug: 0 // logging level (0, 1 or 2)
 };
 export class GeoJSONTiler {
-  constructor(data, options) {
-    this.options = void 0;
-    this.tiles = {};
-    this.tileCoords = [];
-    this.stats = {};
-    this.total = 0;
-    this.options = {
-      ...DEFAULT_OPTIONS,
-      ...options
-    };
-    options = this.options;
-    const debug = options.debug;
-    if (debug) console.time('preprocess data');
-    if (this.options.maxZoom < 0 || this.options.maxZoom > 24) {
-      throw new Error('maxZoom should be in the 0-24 range');
-    }
-    if (options.promoteId && this.options.generateId) {
-      throw new Error('promoteId and generateId cannot be used together.');
-    }
-    let features = convert(data, options);
-    if (debug) {
-      console.timeEnd('preprocess data');
-      console.log('index: maxZoom: %d, maxPoints: %d', options.indexMaxZoom, options.indexMaxPoints);
-      console.time('generate tiles');
-    }
-    features = wrap(features, this.options);
-    if (features.length) {
-      this.splitTile(features, 0, 0, 0);
-    }
-    if (debug) {
-      if (features.length) {
-        console.log('features: %d, points: %d', this.tiles[0].numFeatures, this.tiles[0].numPoints);
-      }
-      console.timeEnd('generate tiles');
-      console.log('tiles generated:', this.total, JSON.stringify(this.stats));
-    }
-  }
-  getTile(z, x, y) {
-    const {
-      extent,
-      debug
-    } = this.options;
-    if (z < 0 || z > 24) {
-      return null;
-    }
-    const z2 = 1 << z;
-    x = x + z2 & z2 - 1;
-    const id = toID(z, x, y);
-    if (this.tiles[id]) {
-      return transformTile(this.tiles[id], extent);
-    }
-    if (debug > 1) console.log('drilling down to z%d-%d-%d', z, x, y);
-    let z0 = z;
-    let x0 = x;
-    let y0 = y;
-    let parent;
-    while (!parent && z0 > 0) {
-      z0--;
-      x0 = x0 >> 1;
-      y0 = y0 >> 1;
-      parent = this.tiles[toID(z0, x0, y0)];
-    }
-    if (!parent || !parent.source) {
-      return null;
-    }
-    if (debug > 1) {
-      console.log('found parent tile z%d-%d-%d', z0, x0, y0);
-      console.time('drilling down');
-    }
-    this.splitTile(parent.source, z0, x0, y0, z, x, y);
-    if (debug > 1) {
-      console.timeEnd('drilling down');
+    constructor(data, options) {
+        // tiles and tileCoords are part of the public API
+        this.tiles = {};
+        this.tileCoords = [];
+        this.stats = {};
+        this.total = 0;
+        this.options = { ...DEFAULT_OPTIONS, ...options };
+        options = this.options;
+        const debug = options.debug;
+        if (debug)
+            console.time('preprocess data');
+        if (this.options.maxZoom < 0 || this.options.maxZoom > 24) {
+            throw new Error('maxZoom should be in the 0-24 range');
+        }
+        if (options.promoteId && this.options.generateId) {
+            throw new Error('promoteId and generateId cannot be used together.');
+        }
+        // projects and adds simplification info
+        let features = convert(data, options);
+        if (debug) {
+            console.timeEnd('preprocess data');
+            console.log('index: maxZoom: %d, maxPoints: %d', options.indexMaxZoom, options.indexMaxPoints);
+            console.time('generate tiles');
+        }
+        // wraps features (ie extreme west and extreme east)
+        features = wrap(features, this.options);
+        // start slicing from the top tile down
+        if (features.length) {
+            this.splitTile(features, 0, 0, 0);
+        }
+        if (debug) {
+            if (features.length) {
+                console.log('features: %d, points: %d', this.tiles[0].numFeatures, this.tiles[0].numPoints);
+            }
+            console.timeEnd('generate tiles');
+            console.log('tiles generated:', this.total, JSON.stringify(this.stats));
+        }
     }
-    return this.tiles[id] ? transformTile(this.tiles[id], extent) : null;
-  }
-  splitTile(features, z, x, y, cz, cx, cy) {
-    const stack = [features, z, x, y];
-    const options = this.options;
-    const debug = options.debug;
-    while (stack.length) {
-      y = stack.pop();
-      x = stack.pop();
-      z = stack.pop();
-      features = stack.pop();
-      const z2 = 1 << z;
-      const id = toID(z, x, y);
-      let tile = this.tiles[id];
-      if (!tile) {
+    /**
+     * Get a tile at the specified index
+     * @param z
+     * @param x
+     * @param y
+     * @returns
+     */
+    // eslint-disable-next-line complexity, max-statements
+    getTile(z, x, y) {
+        // z = +z;
+        // x = +x;
+        // y = +y;
+        const { extent, debug } = this.options;
+        if (z < 0 || z > 24) {
+            return null;
+        }
+        const z2 = 1 << z;
+        x = (x + z2) & (z2 - 1); // wrap tile x coordinate
+        const id = toID(z, x, y);
+        if (this.tiles[id]) {
+            return transformTile(this.tiles[id], extent);
+        }
+        if (debug > 1)
+            console.log('drilling down to z%d-%d-%d', z, x, y);
+        let z0 = z;
+        let x0 = x;
+        let y0 = y;
+        let parent;
+        while (!parent && z0 > 0) {
+            z0--;
+            x0 = x0 >> 1;
+            y0 = y0 >> 1;
+            parent = this.tiles[toID(z0, x0, y0)];
+        }
+        if (!parent || !parent.source) {
+            return null;
+        }
+        // if we found a parent tile containing the original geometry, we can drill down from it
         if (debug > 1) {
-          console.time('creation');
+            console.log('found parent tile z%d-%d-%d', z0, x0, y0);
+            console.time('drilling down');
         }
-        tile = this.tiles[id] = createTile(features, z, x, y, options);
-        this.tileCoords.push({
-          z,
-          x,
-          y
-        });
-        if (debug) {
-          if (debug > 1) {
-            console.log('tile z%d-%d-%d (features: %d, points: %d, simplified: %d)', z, x, y, tile.numFeatures, tile.numPoints, tile.numSimplified);
-            console.timeEnd('creation');
-          }
-          const key = `z${z}`;
-          this.stats[key] = (this.stats[key] || 0) + 1;
-          this.total++;
+        this.splitTile(parent.source, z0, x0, y0, z, x, y);
+        if (debug > 1) {
+            console.timeEnd('drilling down');
+        }
+        return this.tiles[id] ? transformTile(this.tiles[id], extent) : null;
+    }
+    /**
+     * splits features from a parent tile to sub-tiles.
+     * @param z, x, and y are the coordinates of the parent tile
+     * @param cz, cx, and cy are the coordinates of the target tile
+     *
+     * If no target tile is specified, splitting stops when we reach the maximum
+     * zoom or the number of points is low as specified in the options.
+     */
+    // eslint-disable-next-line max-params, max-statements, complexity
+    splitTile(features, z, x, y, cz, cx, cy) {
+        const stack = [features, z, x, y];
+        const options = this.options;
+        const debug = options.debug;
+        // avoid recursion by using a processing queue
+        while (stack.length) {
+            y = stack.pop();
+            x = stack.pop();
+            z = stack.pop();
+            features = stack.pop();
+            const z2 = 1 << z;
+            const id = toID(z, x, y);
+            let tile = this.tiles[id];
+            if (!tile) {
+                if (debug > 1) {
+                    console.time('creation');
+                }
+                tile = this.tiles[id] = createTile(features, z, x, y, options);
+                this.tileCoords.push({ z, x, y });
+                if (debug) {
+                    if (debug > 1) {
+                        console.log('tile z%d-%d-%d (features: %d, points: %d, simplified: %d)', z, x, y, tile.numFeatures, tile.numPoints, tile.numSimplified);
+                        console.timeEnd('creation');
+                    }
+                    const key = `z${z}`;
+                    this.stats[key] = (this.stats[key] || 0) + 1;
+                    this.total++;
+                }
+            }
+            // save reference to original geometry in tile so that we can drill down later if we stop now
+            tile.source = features;
+            // if it's the first-pass tiling
+            if (cz === undefined) {
+                // stop tiling if we reached max zoom, or if the tile is too simple
+                if (z === options.indexMaxZoom || tile.numPoints <= options.indexMaxPoints)
+                    continue;
+                // if a drilldown to a specific tile
+            }
+            else if (z === options.maxZoom || z === cz) {
+                // stop tiling if we reached base zoom or our target tile zoom
+                continue;
+            }
+            else if (cz !== undefined) {
+                // stop tiling if it's not an ancestor of the target tile
+                const zoomSteps = cz - z;
+                // @ts-expect-error TODO fix the types of cx cy
+                if (x !== cx >> zoomSteps || y !== cy >> zoomSteps)
+                    continue;
+            }
+            // if we slice further down, no need to keep source geometry
+            tile.source = null;
+            if (features.length === 0)
+                continue;
+            if (debug > 1)
+                console.time('clipping');
+            // values we'll use for clipping
+            const k1 = (0.5 * options.buffer) / options.extent;
+            const k2 = 0.5 - k1;
+            const k3 = 0.5 + k1;
+            const k4 = 1 + k1;
+            let tl = null;
+            let bl = null;
+            let tr = null;
+            let br = null;
+            let left = clip(features, z2, x - k1, x + k3, 0, tile.minX, tile.maxX, options);
+            let right = clip(features, z2, x + k2, x + k4, 0, tile.minX, tile.maxX, options);
+            // @ts-expect-error - unclear why this is needed?
+            features = null;
+            if (left) {
+                tl = clip(left, z2, y - k1, y + k3, 1, tile.minY, tile.maxY, options);
+                bl = clip(left, z2, y + k2, y + k4, 1, tile.minY, tile.maxY, options);
+                left = null;
+            }
+            if (right) {
+                tr = clip(right, z2, y - k1, y + k3, 1, tile.minY, tile.maxY, options);
+                br = clip(right, z2, y + k2, y + k4, 1, tile.minY, tile.maxY, options);
+                right = null;
+            }
+            if (debug > 1)
+                console.timeEnd('clipping');
+            stack.push(tl || [], z + 1, x * 2, y * 2);
+            stack.push(bl || [], z + 1, x * 2, y * 2 + 1);
+            stack.push(tr || [], z + 1, x * 2 + 1, y * 2);
+            stack.push(br || [], z + 1, x * 2 + 1, y * 2 + 1);
         }
-      }
-      tile.source = features;
-      if (cz === undefined) {
-        if (z === options.indexMaxZoom || tile.numPoints <= options.indexMaxPoints) continue;
-      } else if (z === options.maxZoom || z === cz) {
-        continue;
-      } else if (cz !== undefined) {
-        const zoomSteps = cz - z;
-        if (x !== cx >> zoomSteps || y !== cy >> zoomSteps) continue;
-      }
-      tile.source = null;
-      if (features.length === 0) continue;
-      if (debug > 1) console.time('clipping');
-      const k1 = 0.5 * options.buffer / options.extent;
-      const k2 = 0.5 - k1;
-      const k3 = 0.5 + k1;
-      const k4 = 1 + k1;
-      let tl = null;
-      let bl = null;
-      let tr = null;
-      let br = null;
-      let left = clip(features, z2, x - k1, x + k3, 0, tile.minX, tile.maxX, options);
-      let right = clip(features, z2, x + k2, x + k4, 0, tile.minX, tile.maxX, options);
-      features = null;
-      if (left) {
-        tl = clip(left, z2, y - k1, y + k3, 1, tile.minY, tile.maxY, options);
-        bl = clip(left, z2, y + k2, y + k4, 1, tile.minY, tile.maxY, options);
-        left = null;
-      }
-      if (right) {
-        tr = clip(right, z2, y - k1, y + k3, 1, tile.minY, tile.maxY, options);
-        br = clip(right, z2, y + k2, y + k4, 1, tile.minY, tile.maxY, options);
-        right = null;
-      }
-      if (debug > 1) console.timeEnd('clipping');
-      stack.push(tl || [], z + 1, x * 2, y * 2);
-      stack.push(bl || [], z + 1, x * 2, y * 2 + 1);
-      stack.push(tr || [], z + 1, x * 2 + 1, y * 2);
-      stack.push(br || [], z + 1, x * 2 + 1, y * 2 + 1);
     }
-  }
 }
 function toID(z, x, y) {
-  return ((1 << z) * y + x) * 32 + z;
+    return ((1 << z) * y + x) * 32 + z;
 }
-//# sourceMappingURL=geojson-tiler.js.map

package/dist/lib/geojson-tiler/simplify.js

@@ -1,46 +1,66 @@
+// loaders.gl
+// SPDX-License-Identifier: MIT
+// Copyright (c) vis.gl contributors
+// Forked from https://github.com/mapbox/geojson-vt under compatible ISC license
+/**
+ * Calculate simplification data using optimized Douglas-Peucker algorithm
+ *
+ * @param coords contiguous list of coordinates
+ * @param first first coord to simplify
+ * @param last last coord to simplify
+ * @param sqTolerance tolerance (square distance)
+ */
 export function simplify(coords, first, last, sqTolerance) {
-  let maxSqDist = sqTolerance;
-  const mid = last - first >> 1;
-  let minPosToMid = last - first;
-  let index;
-  const ax = coords[first];
-  const ay = coords[first + 1];
-  const bx = coords[last];
-  const by = coords[last + 1];
-  for (let i = first + 3; i < last; i += 3) {
-    const d = getSqSegDist(coords[i], coords[i + 1], ax, ay, bx, by);
-    if (d > maxSqDist) {
-      index = i;
-      maxSqDist = d;
-    } else if (d === maxSqDist) {
-      const posToMid = Math.abs(i - mid);
-      if (posToMid < minPosToMid) {
-        index = i;
-        minPosToMid = posToMid;
-      }
+    let maxSqDist = sqTolerance;
+    const mid = (last - first) >> 1;
+    let minPosToMid = last - first;
+    let index;
+    const ax = coords[first];
+    const ay = coords[first + 1];
+    const bx = coords[last];
+    const by = coords[last + 1];
+    for (let i = first + 3; i < last; i += 3) {
+        const d = getSqSegDist(coords[i], coords[i + 1], ax, ay, bx, by);
+        if (d > maxSqDist) {
+            index = i;
+            maxSqDist = d;
+        }
+        else if (d === maxSqDist) {
+            // a workaround to ensure we choose a pivot close to the middle of the list,
+            // reducing recursion depth, for certain degenerate inputs
+            // https://github.com/mapbox/geojson-vt/issues/104
+            const posToMid = Math.abs(i - mid);
+            if (posToMid < minPosToMid) {
+                index = i;
+                minPosToMid = posToMid;
+            }
+        }
+    }
+    if (maxSqDist > sqTolerance) {
+        if (index - first > 3)
+            simplify(coords, first, index, sqTolerance);
+        coords[index + 2] = maxSqDist;
+        if (last - index > 3)
+            simplify(coords, index, last, sqTolerance);
     }
-  }
-  if (maxSqDist > sqTolerance) {
-    if (index - first > 3) simplify(coords, first, index, sqTolerance);
-    coords[index + 2] = maxSqDist;
-    if (last - index > 3) simplify(coords, index, last, sqTolerance);
-  }
 }
+/** square distance from a point to a segment */
+// eslint-disable-next-line max-params
 function getSqSegDist(px, py, x, y, bx, by) {
-  let dx = bx - x;
-  let dy = by - y;
-  if (dx !== 0 || dy !== 0) {
-    const t = ((px - x) * dx + (py - y) * dy) / (dx * dx + dy * dy);
-    if (t > 1) {
-      x = bx;
-      y = by;
-    } else if (t > 0) {
-      x += dx * t;
-      y += dy * t;
+    let dx = bx - x;
+    let dy = by - y;
+    if (dx !== 0 || dy !== 0) {
+        const t = ((px - x) * dx + (py - y) * dy) / (dx * dx + dy * dy);
+        if (t > 1) {
+            x = bx;
+            y = by;
+        }
+        else if (t > 0) {
+            x += dx * t;
+            y += dy * t;
+        }
     }
-  }
-  dx = px - x;
-  dy = py - y;
-  return dx * dx + dy * dy;
+    dx = px - x;
+    dy = py - y;
+    return dx * dx + dy * dy;
 }
-//# sourceMappingURL=simplify.js.map