mapshaper 0.7.58 → 0.7.60

package/mapshaper.js

@@ -18168,10 +18168,29 @@
     parseSizeParam: parseSizeParam
   });
 
-  // Default gap-width as a fraction of the layer's median polygon-ring segment
-  // length. Shared by interior fill and exterior close-outer-gaps.
+  // Conservative width as a fraction of the layer's median polygon-ring segment
+  // length. Used for partition min-width and duplicate-boundary search — not for
+  // the default gap-width=auto fill threshold (see getDefaultGapWidth).
   var GAP_WIDTH_SEGMENT_FRACTION = 0.01;
 
+  // Inject gap_width=auto unless the caller already chose a width or a legacy
+  // area/sliver option. Used by -clean, -filter-slivers, and gap partition.
+  function applyDefaultGapWidthOpts(opts) {
+    opts = Object.assign({}, opts);
+    if (opts.gap_width != null) return opts;
+    // Legacy area/sliver options: keep the historical gap_fill_area=auto default.
+    if (opts.gap_fill_area != null || opts.min_gap_area != null ||
+        opts.min_area != null || opts.sliver_control != null) {
+      if (opts.gap_fill_area == null && opts.min_gap_area == null &&
+          opts.min_area == null) {
+        opts.gap_fill_area = 'auto';
+      }
+      return opts;
+    }
+    opts.gap_width = 'auto';
+    return opts;
+  }
+
   // Used by -clean -dissolve -filter-slivers -filter-islands to generate filters
   // for removing small polygon rings / filling mosaic gaps.
   // Prefers gap_width= (including 'auto'). Legacy gap_fill_area / sliver_control
@@ -18317,9 +18336,10 @@
     };
   }
 
-  // Median polygon-ring segment length * GAP_WIDTH_SEGMENT_FRACTION.
+  // Width equivalent of the legacy gap-fill-area=auto + sliver-control=1 test:
+  // W = sqrt(A0 / π). The two filters accept the same rings (see getGapWidthTest).
   function getDefaultGapWidth(lyr, arcs) {
-    return getMedianPolygonSegmentLength(lyr, arcs) * GAP_WIDTH_SEGMENT_FRACTION;
+    return Math.sqrt(getDefaultSliverThreshold(lyr, arcs) / Math.PI);
   }
 
   function getMedianPolygonSegmentLength(lyr, arcs) {
@@ -18400,6 +18420,7 @@
   var Slivers = /*#__PURE__*/Object.freeze({
     __proto__: null,
     GAP_WIDTH_SEGMENT_FRACTION: GAP_WIDTH_SEGMENT_FRACTION,
+    applyDefaultGapWidthOpts: applyDefaultGapWidthOpts,
     calcMaxSliverArea: calcMaxSliverArea,
     getDefaultGapWidth: getDefaultGapWidth,
     getDefaultSliverThreshold: getDefaultSliverThreshold,
@@ -23467,8 +23488,7 @@
     var spherical = isLatLngCRS(crs);
     var medianSeg = measureMedianSegment(lyr, arcs);
     // Same width as interior gap filling, from the same automatic default.
-    var distance = resolveCloseDistance(opts, crs,
-      medianSeg * GAP_WIDTH_SEGMENT_FRACTION);
+    var distance = resolveCloseDistance(opts, crs, getDefaultGapWidth(lyr, arcs));
     if (!(distance > 0)) return 0;
 
     profileStart('cg.outsideFacingArcs');
@@ -23518,7 +23538,7 @@
     return {seams: seams, seeds: seeds, arcsById: arcsById, footpoints: footpoints};
   }
 
-  // Same segment median as interior gap-width=auto (all polygon-ring segments).
+  // Median of every polygon-ring segment; used for min-seam length.
   function measureMedianSegment(lyr, arcs) {
     return profileWrap('cg.medianSegment', function() {
       return getMedianPolygonSegmentLength(lyr, arcs);
@@ -24154,18 +24174,7 @@
   function partitionPolygonMosaicGaps(lyr, dataset, nodes, opts) {
     if (!lyr.shapes || lyr.shapes.length < MIN_GAP_OWNERS) return false;
     var mosaicIndex = new MosaicIndex(lyr, nodes, {flat: true});
-    var sliverOpts = utils.extend({}, opts);
-    if (sliverOpts.gap_width == null) {
-      if (sliverOpts.gap_fill_area != null || sliverOpts.min_gap_area != null ||
-          sliverOpts.min_area != null || sliverOpts.sliver_control != null) {
-        if (sliverOpts.gap_fill_area == null && sliverOpts.min_gap_area == null &&
-            sliverOpts.min_area == null) {
-          sliverOpts.gap_fill_area = 'auto';
-        }
-      } else {
-        sliverOpts.gap_width = 'auto';
-      }
-    }
+    var sliverOpts = applyDefaultGapWidthOpts(opts);
     var filter = getSliverFilter(lyr, dataset, sliverOpts).filter;
     var gaps = mosaicIndex.getUnusedTileData(filter).filter(function(gap) {
       return countOwners(gap.boundary) >= MIN_GAP_OWNERS;
@@ -24601,10 +24610,11 @@
   // several features with hairline slivers, which is worse for having more of them.
   //
   // A layer's median segment length is the scale at which its boundaries were
-  // drawn, so the same fraction of it that serves as the default gap-width -- what
-  // mapshaper already treats as too narrow to be intended -- marks where dividing
-  // stops being worth it. This is deliberately independent of the gap-width in
-  // effect: raising that to fill wider gaps should not coarsen the data's own scale.
+  // drawn. A small fraction of that scale marks where dividing a gap stops being
+  // worth it: the displacement would be smaller than the data's own detail, and
+  // the pieces would be hairline slivers. This is deliberately independent of
+  // the gap-width in effect: raising that to fill wider gaps should not coarsen
+  // the data's own scale.
   function getMinDivisibleGapWidth(lyr, arcs) {
     var width = getMedianPolygonSegmentLength(lyr, arcs) *
       GAP_WIDTH_SEGMENT_FRACTION;
@@ -24741,7 +24751,7 @@
 
   function cleanPolygonLayerGeometry(lyr, dataset, opts) {
     // clean polygons by apply the 'dissolve2' function to each feature
-    opts = withDefaultGapWidth(opts);
+    opts = applyDefaultGapWidthOpts(opts);
     var groups = lyr.shapes.map(function(shp, i) {
       return [i];
     });
@@ -24750,22 +24760,6 @@
     markLayerChanged(lyr, {operation: 'cleanPolygonLayerGeometry', unit: 'shapes'});
   }
 
-  function withDefaultGapWidth(opts) {
-    opts = Object.assign({}, opts);
-    if (opts.gap_width != null) return opts;
-    // Legacy area/sliver options: keep the historical gap_fill_area=auto default.
-    if (opts.gap_fill_area != null || opts.min_gap_area != null ||
-        opts.min_area != null || opts.sliver_control != null) {
-      if (opts.gap_fill_area == null && opts.min_gap_area == null &&
-          opts.min_area == null) {
-        opts.gap_fill_area = 'auto';
-      }
-      return opts;
-    }
-    opts.gap_width = 'auto';
-    return opts;
-  }
-
   function collectArcIds(shapes) {
     var ids = {};
     traversePaths(shapes, function(o) {
@@ -40329,12 +40323,13 @@ ${svg}
 
     parser.command('filter-slivers')
       .describe('remove small polygon rings')
+      .option('gap-width', gapWidthOpt)
       .option('min-area', {
         type: 'area',
-        describe: 'area threshold (e.g. 2sqkm)'
+        describe: '(deprecated) use gap-width= instead'
       })
       .option('sliver-control', {
-        describe: 'boost area threshold of slivers (0-1, default is 1)',
+        describe: '(deprecated) use gap-width= instead',
         type: 'number'
       })
       .option('weighted', {
@@ -40383,9 +40378,9 @@ ${svg}
       .option('no-replace', noReplaceOpt);
 
     parser.command('grid')
-      .describe('create a grid of square, hexagonal, rhombus or triangle polygons')
+      .describe('create a grid of square, hexagonal, rhombus, triangle or cairo polygons')
       .option('type', {
-        describe: 'square, square2, hex, hex2, rhombus, rhombus2, triangle or triangle2 (default is square)'
+        describe: 'square, square2, hex, hex2, rhombus, rhombus2, triangle, triangle2 or cairo (default is square)'
       })
       .option('interval', {
         describe: 'side length (e.g. 500m, 12km)',
@@ -40407,6 +40402,10 @@ ${svg}
         describe: 'scale factor for cells, between 0 and 2',
         type: 'number'
       })
+      .option('rotate', {
+        describe: 'rotation angle in degrees (default is 0)',
+        type: 'number'
+      })
       // .option('bbox', {
       //   type: 'bbox',
       //   describe: 'xmin,ymin,xmax,ymax (default is bbox of data)'
@@ -56569,7 +56568,13 @@ ${svg}
   };
 
   function filterSlivers(lyr, dataset, optsArg) {
-    var opts = utils.extend({sliver_control: 1}, optsArg);
+    var opts = applyDefaultGapWidthOpts(optsArg);
+    // Legacy area path only: keep the historical sliver-control=1 default.
+    // Do not inject it before applyDefaultGapWidthOpts, or the implicit default
+    // would lock the command onto the area path.
+    if (opts.gap_width == null && opts.sliver_control == null) {
+      opts.sliver_control = 1;
+    }
     var filterData = getSliverFilter(lyr, dataset, opts);
     var ringTest = filterData.filter;
     var removed = 0;
@@ -74555,7 +74560,7 @@ ${svg}
   cmd.polygonGrid = function(targetLayers, targetDataset, opts) {
     requireProjectedDataset(targetDataset);
     var params = getGridParams(targetLayers, targetDataset, opts);
-    var gridDataset = makeGridDataset(params); // grid is a new dataset
+    var gridDataset = makeGridDataset(params, opts); // grid is a new dataset
     gridDataset.info = copyDatasetInfo(targetDataset.info);
     setOutputLayerName(gridDataset.layers[0], null, 'grid', opts);
     if (opts.debug) gridDataset.layers.push(cmd.pointGrid2(targetLayers, targetDataset, opts));
@@ -74619,35 +74624,67 @@ ${svg}
 
   function makeGridDataset(params, opts) {
     var geojson, dataset;
+    var rotation = opts.rotate || 0;
+    var genParams = rotation ? getRotatedGridParams(params, rotation) : params;
     if (params.type == 'square') {
-      geojson = getSquareGridGeoJSON(getSquareGridCoordinates(params));
+      geojson = getSquareGridGeoJSON(getSquareGridCoordinates(genParams));
     } else if (params.type == 'square2') {
-      geojson = getRotatedSquareGridGeoJSON(params);
+      geojson = getRotatedSquareGridGeoJSON(genParams);
     } else if (params.type == 'hex') {
-      geojson = getHexGridGeoJSON(getHexGridCoordinates(params));
+      geojson = getHexGridGeoJSON(getHexGridCoordinates(genParams));
     } else if (params.type == 'hex2') {
       // use rotated grid
-      geojson = getHexGridGeoJSON(getHexGridCoordinates(swapGridParams(params)));
+      geojson = getHexGridGeoJSON(getHexGridCoordinates(swapGridParams(genParams)));
       swapPolygonCoords(geojson);
     } else if (params.type == 'rhombus') {
-      geojson = getRhombusGridGeoJSON(params, false);
+      geojson = getRhombusGridGeoJSON(genParams, false);
     } else if (params.type == 'rhombus2') {
-      geojson = getRhombusGridGeoJSON(params, true);
+      geojson = getRhombusGridGeoJSON(genParams, true);
     } else if (params.type == 'triangle') {
-      geojson = getTriangleGridGeoJSON(params, false);
+      geojson = getTriangleGridGeoJSON(genParams, false);
     } else if (params.type == 'triangle2') {
-      geojson = getTriangleGridGeoJSON(params, true);
+      geojson = getTriangleGridGeoJSON(genParams, true);
+    } else if (params.type == 'cairo') {
+      geojson = getCairoGridGeoJSON(genParams);
     } else {
       stop$1('Unsupported grid type');
     }
     scaleGridCells(geojson, params.cellScale);
     alignGridToBounds(geojson, params.bbox);
+    if (rotation) {
+      rotateGridCells(geojson, rotation, params.bbox);
+    }
     cullGridCells(geojson, params.bbox);
     dataset = importGeoJSON(geojson, {});
     buildTopology(dataset);
     return dataset;
   }
 
+  // Expand the generation rectangle so that, after rotation around the shared
+  // center, the grid still covers the original target bbox.
+  function getRotatedGridParams(params, rotation) {
+    var angle = rotation * Math.PI / 180;
+    var cos = Math.abs(Math.cos(angle));
+    var sin = Math.abs(Math.sin(angle));
+    return utils.defaults({
+      width: params.width * cos + params.height * sin,
+      height: params.width * sin + params.height * cos
+    }, params);
+  }
+
+  function rotateGridCells(geojson, rotation, bbox) {
+    var cx = (bbox[0] + bbox[2]) / 2;
+    var cy = (bbox[1] + bbox[3]) / 2;
+    var transform = getAffineTransform(rotation, 1, [0, 0], [cx, cy]);
+    geojson.geometries.forEach(function(geom) {
+      if (geom.type == 'Polygon') {
+        geom.coordinates[0] = geom.coordinates[0].map(function(xy) {
+          return transform(xy[0], xy[1]);
+        });
+      }
+    });
+  }
+
   function swapGridParams(params) {
     var bbox = params.bbox;
     return utils.defaults({
@@ -74729,6 +74766,7 @@ ${svg}
     if (type == 'hex' || type == 'hex2') return 3 * Math.sqrt(3) / 2;
     if (type == 'rhombus' || type == 'rhombus2') return Math.sqrt(3) / 2;
     if (type == 'triangle' || type == 'triangle2') return Math.sqrt(3) / 4;
+    if (type == 'cairo') return (4 + Math.sqrt(7)) / 4;
     stop$1('Unsupported grid type');
   }
 
@@ -74882,6 +74920,95 @@ ${svg}
     return geojson;
   }
 
+  // Equilateral type-4 Cairo tessellation. Four-valent (90°) vertices sit on a
+  // square lattice of spacing S = interval * (√7+1)/2; each lattice square
+  // contains a central bar of length `interval`, alternating horizontal and
+  // vertical. Two pentagons share each bar.
+  //
+  // Vertices are computed once and reused. Algebraically equal expressions
+  // such as i*S and (i-1)*S + S are not bit-identical, and independently
+  // computed copies of the same 4-way junction become T-intersections.
+  function getCairoGridGeoJSON(params) {
+    var side = params.interval;
+    var spacing = side * (Math.sqrt(7) + 1) / 2;
+    var halfBar = side / 2;
+    var iMin = -1;
+    var jMin = -1;
+    var iMax = Math.ceil(params.width / spacing) + 1;
+    var jMax = Math.ceil(params.height / spacing) + 1;
+    var lattice = [];
+    var bars = [];
+    var geometries = [];
+    var i, j, y, cx, cy;
+
+    for (j = jMin - 1; j <= jMax + 1; j++) {
+      y = j * spacing;
+      lattice[j] = [];
+      for (i = iMin - 1; i <= iMax + 1; i++) {
+        lattice[j][i] = [i * spacing, y];
+      }
+    }
+
+    for (j = jMin - 1; j < jMax + 1; j++) {
+      bars[j] = [];
+      for (i = iMin - 1; i < iMax + 1; i++) {
+        cx = (i + 0.5) * spacing;
+        cy = (j + 0.5) * spacing;
+        bars[j][i] = (i + j) % 2 === 0 ? {
+          a: [cx - halfBar, cy],
+          b: [cx + halfBar, cy]
+        } : {
+          a: [cx, cy - halfBar],
+          b: [cx, cy + halfBar]
+        };
+      }
+    }
+
+    for (j = jMin; j < jMax; j++) {
+      for (i = iMin; i < iMax; i++) {
+        if ((i + j) % 2 === 0) {
+          appendCairoPentagon(geometries, [
+            lattice[j][i],
+            bars[j - 1][i].b,
+            lattice[j][i + 1],
+            bars[j][i].b,
+            bars[j][i].a
+          ]);
+          appendCairoPentagon(geometries, [
+            lattice[j + 1][i],
+            bars[j][i].a,
+            bars[j][i].b,
+            lattice[j + 1][i + 1],
+            bars[j + 1][i].a
+          ]);
+        } else {
+          appendCairoPentagon(geometries, [
+            lattice[j][i],
+            bars[j][i].a,
+            bars[j][i].b,
+            lattice[j + 1][i],
+            bars[j][i - 1].b
+          ]);
+          appendCairoPentagon(geometries, [
+            lattice[j][i + 1],
+            bars[j][i + 1].a,
+            lattice[j + 1][i + 1],
+            bars[j][i].b,
+            bars[j][i].a
+          ]);
+        }
+      }
+    }
+    return {type: 'GeometryCollection', geometries: geometries};
+  }
+
+  function appendCairoPentagon(geometries, verts) {
+    geometries.push({
+      type: 'Polygon',
+      coordinates: [verts.concat([verts[0]])]
+    });
+  }
+
   function subdivideHexagon(coords) {
     var center = getPolygonRingCenter(coords);
     return [
@@ -80084,7 +80211,7 @@ ${svg}
     return name == 'rectangle' || name == 'rectangles' || name == 'filter' && opts.cleanup;
   }
 
-  var version = "0.7.58";
+  var version = "0.7.60";
 
   // Parse command line args into commands and run them
   // Function takes an optional Node-style callback. A Promise is returned if no callback is given.

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "mapshaper",
-  "version": "0.7.58",
+  "version": "0.7.60",
   "description": "A tool for editing geospatial data for mapping and GIS.",
   "keywords": [
     "shapefile",

package/www/basemap.js

@@ -3,7 +3,7 @@ window.mapboxParams = {
   // css: 'https://api.mapbox.com/mapbox-gl-js/v2.11.1/mapbox-gl.css',
   js: 'https://api.mapbox.com/mapbox-gl-js/v3.15.0/mapbox-gl.js',
   css: 'https://api.mapbox.com/mapbox-gl-js/v3.15.0/mapbox-gl.css',
-  production_key: 'pk.eyJ1IjoiZ3JhbW1hdGEiLCJhIjoiY21jZ2RlbGk1MGhvNzJrcHV2dTYzZWE2byJ9.btkRnRrPW75c0AlQ5lH2cg',
+  production_key: 'pk.eyJ1IjoiZ3JhbW1hdGEiLCJhIjoiY210cjg2NTk5MDBiMTJ4b3dvbDNwMjdudSJ9.EohCQlad59ntK8CqWlopKw',
   localhost_key: 'pk.eyJ1IjoiZ3JhbW1hdGEiLCJhIjoiY21jZ2Rna2hrMGo2ZjJqcHFxMHAycDlyMSJ9.q8FchhB-NAN_FGBHN0AJeg',
   styles: [{
     name: 'Reference map',

package/www/mapshaper-gui.js

@@ -11109,6 +11109,8 @@
 
   var openMenu;
   var openMenuId;
+  // One menu element per parent element, created on demand.
+  var menus = new WeakMap();
   // Prefixes shown to the left of a menu item. The checkmark is what a copyable
   // item shows after it has been copied.
   var BULLET = '•  ';
@@ -11132,14 +11134,26 @@
     }
   }
 
+  // Menus are re-used, rather than created per click: a discarded menu stays in
+  // the DOM, and an empty one is a 10px sliver at the bottom of the page, which
+  // is enough to make the page scrollable (see the note in ContextMenu()).
+  function getContextMenu(parentArg) {
+    var parent = parentArg || document.querySelector('body');
+    var menu = menus.get(parent);
+    if (!menu) {
+      menu = new ContextMenu(parent);
+      menus.set(parent, menu);
+    }
+    return menu;
+  }
+
   function openContextMenu(e, lyr, parent) {
-    var menu = new ContextMenu(parent);
     if (e.contextMenuId && e.contextMenuId == openMenuId) {
       closeOpenMenu(true);
       return;
     }
     closeOpenMenu(true);
-    menu.open(e, lyr);
+    getContextMenu(parent).open(e, lyr);
     openMenuId = e.contextMenuId || null;
   }
 
@@ -11151,6 +11165,14 @@
     var _open = false;
     var _openCount = 0;
 
+    // The menu is hidden until it is opened. A visible empty menu is 10px of
+    // padding and border at the end of the page, and the page is exactly as tall
+    // as the window, so those 10px make the document scrollable. The body is
+    // overflow:hidden, but that only stops the user from scrolling: the browser
+    // still scrolls the page to reveal a focused element or a caret, and once it
+    // does, the whole UI stays shifted up (mbloch/mapshaper#702).
+    menu.hide();
+
     this.isOpen = function() {
       return _open;
     };
@@ -11292,10 +11314,23 @@
       }
       menu.css('top', (e.pageY - offs.top - 15) + 'px');
       menu.show();
+      keepInsideWindow();
 
       _open = true;
       _openCount++;
 
+      // A menu opened near the bottom of the window would otherwise hang off the
+      // screen, out of reach and adding to the page's scrollable overflow.
+      function keepInsideWindow() {
+        var pad = 6;
+        var rect = menu.node().getBoundingClientRect();
+        var overflow = rect.bottom + pad - document.documentElement.clientHeight;
+        var shift = Math.min(overflow, rect.top - pad);
+        if (shift > 0) {
+          menu.css('top', (parseFloat(menu.node().style.top) - shift) + 'px');
+        }
+      }
+
       function getParentOffset() { // crossbrowser version
         if (parent == body) {
           return {top: 0, left: 0};
@@ -25418,7 +25453,7 @@ GUI and setting the size and crop of SVG output.</p><div><input type="text" clas
     gui.messages = new MessageControl(gui);
     gui.basemap = new Basemap(gui);
     gui.session = new SessionHistory(gui);
-    gui.contextMenu = new ContextMenu();
+    gui.contextMenu = getContextMenu();
     gui.undo = new Undo(gui);
     gui.map = new MshpMap(gui);
     // gui.project = new ProjectOptions(gui);

package/www/mapshaper.js

@@ -18168,10 +18168,29 @@
     parseSizeParam: parseSizeParam
   });
 
-  // Default gap-width as a fraction of the layer's median polygon-ring segment
-  // length. Shared by interior fill and exterior close-outer-gaps.
+  // Conservative width as a fraction of the layer's median polygon-ring segment
+  // length. Used for partition min-width and duplicate-boundary search — not for
+  // the default gap-width=auto fill threshold (see getDefaultGapWidth).
   var GAP_WIDTH_SEGMENT_FRACTION = 0.01;
 
+  // Inject gap_width=auto unless the caller already chose a width or a legacy
+  // area/sliver option. Used by -clean, -filter-slivers, and gap partition.
+  function applyDefaultGapWidthOpts(opts) {
+    opts = Object.assign({}, opts);
+    if (opts.gap_width != null) return opts;
+    // Legacy area/sliver options: keep the historical gap_fill_area=auto default.
+    if (opts.gap_fill_area != null || opts.min_gap_area != null ||
+        opts.min_area != null || opts.sliver_control != null) {
+      if (opts.gap_fill_area == null && opts.min_gap_area == null &&
+          opts.min_area == null) {
+        opts.gap_fill_area = 'auto';
+      }
+      return opts;
+    }
+    opts.gap_width = 'auto';
+    return opts;
+  }
+
   // Used by -clean -dissolve -filter-slivers -filter-islands to generate filters
   // for removing small polygon rings / filling mosaic gaps.
   // Prefers gap_width= (including 'auto'). Legacy gap_fill_area / sliver_control
@@ -18317,9 +18336,10 @@
     };
   }
 
-  // Median polygon-ring segment length * GAP_WIDTH_SEGMENT_FRACTION.
+  // Width equivalent of the legacy gap-fill-area=auto + sliver-control=1 test:
+  // W = sqrt(A0 / π). The two filters accept the same rings (see getGapWidthTest).
   function getDefaultGapWidth(lyr, arcs) {
-    return getMedianPolygonSegmentLength(lyr, arcs) * GAP_WIDTH_SEGMENT_FRACTION;
+    return Math.sqrt(getDefaultSliverThreshold(lyr, arcs) / Math.PI);
   }
 
   function getMedianPolygonSegmentLength(lyr, arcs) {
@@ -18400,6 +18420,7 @@
   var Slivers = /*#__PURE__*/Object.freeze({
     __proto__: null,
     GAP_WIDTH_SEGMENT_FRACTION: GAP_WIDTH_SEGMENT_FRACTION,
+    applyDefaultGapWidthOpts: applyDefaultGapWidthOpts,
     calcMaxSliverArea: calcMaxSliverArea,
     getDefaultGapWidth: getDefaultGapWidth,
     getDefaultSliverThreshold: getDefaultSliverThreshold,
@@ -23467,8 +23488,7 @@
     var spherical = isLatLngCRS(crs);
     var medianSeg = measureMedianSegment(lyr, arcs);
     // Same width as interior gap filling, from the same automatic default.
-    var distance = resolveCloseDistance(opts, crs,
-      medianSeg * GAP_WIDTH_SEGMENT_FRACTION);
+    var distance = resolveCloseDistance(opts, crs, getDefaultGapWidth(lyr, arcs));
     if (!(distance > 0)) return 0;
 
     profileStart('cg.outsideFacingArcs');
@@ -23518,7 +23538,7 @@
     return {seams: seams, seeds: seeds, arcsById: arcsById, footpoints: footpoints};
   }
 
-  // Same segment median as interior gap-width=auto (all polygon-ring segments).
+  // Median of every polygon-ring segment; used for min-seam length.
   function measureMedianSegment(lyr, arcs) {
     return profileWrap('cg.medianSegment', function() {
       return getMedianPolygonSegmentLength(lyr, arcs);
@@ -24154,18 +24174,7 @@
   function partitionPolygonMosaicGaps(lyr, dataset, nodes, opts) {
     if (!lyr.shapes || lyr.shapes.length < MIN_GAP_OWNERS) return false;
     var mosaicIndex = new MosaicIndex(lyr, nodes, {flat: true});
-    var sliverOpts = utils.extend({}, opts);
-    if (sliverOpts.gap_width == null) {
-      if (sliverOpts.gap_fill_area != null || sliverOpts.min_gap_area != null ||
-          sliverOpts.min_area != null || sliverOpts.sliver_control != null) {
-        if (sliverOpts.gap_fill_area == null && sliverOpts.min_gap_area == null &&
-            sliverOpts.min_area == null) {
-          sliverOpts.gap_fill_area = 'auto';
-        }
-      } else {
-        sliverOpts.gap_width = 'auto';
-      }
-    }
+    var sliverOpts = applyDefaultGapWidthOpts(opts);
     var filter = getSliverFilter(lyr, dataset, sliverOpts).filter;
     var gaps = mosaicIndex.getUnusedTileData(filter).filter(function(gap) {
       return countOwners(gap.boundary) >= MIN_GAP_OWNERS;
@@ -24601,10 +24610,11 @@
   // several features with hairline slivers, which is worse for having more of them.
   //
   // A layer's median segment length is the scale at which its boundaries were
-  // drawn, so the same fraction of it that serves as the default gap-width -- what
-  // mapshaper already treats as too narrow to be intended -- marks where dividing
-  // stops being worth it. This is deliberately independent of the gap-width in
-  // effect: raising that to fill wider gaps should not coarsen the data's own scale.
+  // drawn. A small fraction of that scale marks where dividing a gap stops being
+  // worth it: the displacement would be smaller than the data's own detail, and
+  // the pieces would be hairline slivers. This is deliberately independent of
+  // the gap-width in effect: raising that to fill wider gaps should not coarsen
+  // the data's own scale.
   function getMinDivisibleGapWidth(lyr, arcs) {
     var width = getMedianPolygonSegmentLength(lyr, arcs) *
       GAP_WIDTH_SEGMENT_FRACTION;
@@ -24741,7 +24751,7 @@
 
   function cleanPolygonLayerGeometry(lyr, dataset, opts) {
     // clean polygons by apply the 'dissolve2' function to each feature
-    opts = withDefaultGapWidth(opts);
+    opts = applyDefaultGapWidthOpts(opts);
     var groups = lyr.shapes.map(function(shp, i) {
       return [i];
     });
@@ -24750,22 +24760,6 @@
     markLayerChanged(lyr, {operation: 'cleanPolygonLayerGeometry', unit: 'shapes'});
   }
 
-  function withDefaultGapWidth(opts) {
-    opts = Object.assign({}, opts);
-    if (opts.gap_width != null) return opts;
-    // Legacy area/sliver options: keep the historical gap_fill_area=auto default.
-    if (opts.gap_fill_area != null || opts.min_gap_area != null ||
-        opts.min_area != null || opts.sliver_control != null) {
-      if (opts.gap_fill_area == null && opts.min_gap_area == null &&
-          opts.min_area == null) {
-        opts.gap_fill_area = 'auto';
-      }
-      return opts;
-    }
-    opts.gap_width = 'auto';
-    return opts;
-  }
-
   function collectArcIds(shapes) {
     var ids = {};
     traversePaths(shapes, function(o) {
@@ -40329,12 +40323,13 @@ ${svg}
 
     parser.command('filter-slivers')
       .describe('remove small polygon rings')
+      .option('gap-width', gapWidthOpt)
       .option('min-area', {
         type: 'area',
-        describe: 'area threshold (e.g. 2sqkm)'
+        describe: '(deprecated) use gap-width= instead'
       })
       .option('sliver-control', {
-        describe: 'boost area threshold of slivers (0-1, default is 1)',
+        describe: '(deprecated) use gap-width= instead',
         type: 'number'
       })
       .option('weighted', {
@@ -40383,9 +40378,9 @@ ${svg}
       .option('no-replace', noReplaceOpt);
 
     parser.command('grid')
-      .describe('create a grid of square, hexagonal, rhombus or triangle polygons')
+      .describe('create a grid of square, hexagonal, rhombus, triangle or cairo polygons')
       .option('type', {
-        describe: 'square, square2, hex, hex2, rhombus, rhombus2, triangle or triangle2 (default is square)'
+        describe: 'square, square2, hex, hex2, rhombus, rhombus2, triangle, triangle2 or cairo (default is square)'
       })
       .option('interval', {
         describe: 'side length (e.g. 500m, 12km)',
@@ -40407,6 +40402,10 @@ ${svg}
         describe: 'scale factor for cells, between 0 and 2',
         type: 'number'
       })
+      .option('rotate', {
+        describe: 'rotation angle in degrees (default is 0)',
+        type: 'number'
+      })
       // .option('bbox', {
       //   type: 'bbox',
       //   describe: 'xmin,ymin,xmax,ymax (default is bbox of data)'
@@ -56569,7 +56568,13 @@ ${svg}
   };
 
   function filterSlivers(lyr, dataset, optsArg) {
-    var opts = utils.extend({sliver_control: 1}, optsArg);
+    var opts = applyDefaultGapWidthOpts(optsArg);
+    // Legacy area path only: keep the historical sliver-control=1 default.
+    // Do not inject it before applyDefaultGapWidthOpts, or the implicit default
+    // would lock the command onto the area path.
+    if (opts.gap_width == null && opts.sliver_control == null) {
+      opts.sliver_control = 1;
+    }
     var filterData = getSliverFilter(lyr, dataset, opts);
     var ringTest = filterData.filter;
     var removed = 0;
@@ -74555,7 +74560,7 @@ ${svg}
   cmd.polygonGrid = function(targetLayers, targetDataset, opts) {
     requireProjectedDataset(targetDataset);
     var params = getGridParams(targetLayers, targetDataset, opts);
-    var gridDataset = makeGridDataset(params); // grid is a new dataset
+    var gridDataset = makeGridDataset(params, opts); // grid is a new dataset
     gridDataset.info = copyDatasetInfo(targetDataset.info);
     setOutputLayerName(gridDataset.layers[0], null, 'grid', opts);
     if (opts.debug) gridDataset.layers.push(cmd.pointGrid2(targetLayers, targetDataset, opts));
@@ -74619,35 +74624,67 @@ ${svg}
 
   function makeGridDataset(params, opts) {
     var geojson, dataset;
+    var rotation = opts.rotate || 0;
+    var genParams = rotation ? getRotatedGridParams(params, rotation) : params;
     if (params.type == 'square') {
-      geojson = getSquareGridGeoJSON(getSquareGridCoordinates(params));
+      geojson = getSquareGridGeoJSON(getSquareGridCoordinates(genParams));
     } else if (params.type == 'square2') {
-      geojson = getRotatedSquareGridGeoJSON(params);
+      geojson = getRotatedSquareGridGeoJSON(genParams);
     } else if (params.type == 'hex') {
-      geojson = getHexGridGeoJSON(getHexGridCoordinates(params));
+      geojson = getHexGridGeoJSON(getHexGridCoordinates(genParams));
     } else if (params.type == 'hex2') {
       // use rotated grid
-      geojson = getHexGridGeoJSON(getHexGridCoordinates(swapGridParams(params)));
+      geojson = getHexGridGeoJSON(getHexGridCoordinates(swapGridParams(genParams)));
       swapPolygonCoords(geojson);
     } else if (params.type == 'rhombus') {
-      geojson = getRhombusGridGeoJSON(params, false);
+      geojson = getRhombusGridGeoJSON(genParams, false);
     } else if (params.type == 'rhombus2') {
-      geojson = getRhombusGridGeoJSON(params, true);
+      geojson = getRhombusGridGeoJSON(genParams, true);
     } else if (params.type == 'triangle') {
-      geojson = getTriangleGridGeoJSON(params, false);
+      geojson = getTriangleGridGeoJSON(genParams, false);
     } else if (params.type == 'triangle2') {
-      geojson = getTriangleGridGeoJSON(params, true);
+      geojson = getTriangleGridGeoJSON(genParams, true);
+    } else if (params.type == 'cairo') {
+      geojson = getCairoGridGeoJSON(genParams);
     } else {
       stop$1('Unsupported grid type');
     }
     scaleGridCells(geojson, params.cellScale);
     alignGridToBounds(geojson, params.bbox);
+    if (rotation) {
+      rotateGridCells(geojson, rotation, params.bbox);
+    }
     cullGridCells(geojson, params.bbox);
     dataset = importGeoJSON(geojson, {});
     buildTopology(dataset);
     return dataset;
   }
 
+  // Expand the generation rectangle so that, after rotation around the shared
+  // center, the grid still covers the original target bbox.
+  function getRotatedGridParams(params, rotation) {
+    var angle = rotation * Math.PI / 180;
+    var cos = Math.abs(Math.cos(angle));
+    var sin = Math.abs(Math.sin(angle));
+    return utils.defaults({
+      width: params.width * cos + params.height * sin,
+      height: params.width * sin + params.height * cos
+    }, params);
+  }
+
+  function rotateGridCells(geojson, rotation, bbox) {
+    var cx = (bbox[0] + bbox[2]) / 2;
+    var cy = (bbox[1] + bbox[3]) / 2;
+    var transform = getAffineTransform(rotation, 1, [0, 0], [cx, cy]);
+    geojson.geometries.forEach(function(geom) {
+      if (geom.type == 'Polygon') {
+        geom.coordinates[0] = geom.coordinates[0].map(function(xy) {
+          return transform(xy[0], xy[1]);
+        });
+      }
+    });
+  }
+
   function swapGridParams(params) {
     var bbox = params.bbox;
     return utils.defaults({
@@ -74729,6 +74766,7 @@ ${svg}
     if (type == 'hex' || type == 'hex2') return 3 * Math.sqrt(3) / 2;
     if (type == 'rhombus' || type == 'rhombus2') return Math.sqrt(3) / 2;
     if (type == 'triangle' || type == 'triangle2') return Math.sqrt(3) / 4;
+    if (type == 'cairo') return (4 + Math.sqrt(7)) / 4;
     stop$1('Unsupported grid type');
   }
 
@@ -74882,6 +74920,95 @@ ${svg}
     return geojson;
   }
 
+  // Equilateral type-4 Cairo tessellation. Four-valent (90°) vertices sit on a
+  // square lattice of spacing S = interval * (√7+1)/2; each lattice square
+  // contains a central bar of length `interval`, alternating horizontal and
+  // vertical. Two pentagons share each bar.
+  //
+  // Vertices are computed once and reused. Algebraically equal expressions
+  // such as i*S and (i-1)*S + S are not bit-identical, and independently
+  // computed copies of the same 4-way junction become T-intersections.
+  function getCairoGridGeoJSON(params) {
+    var side = params.interval;
+    var spacing = side * (Math.sqrt(7) + 1) / 2;
+    var halfBar = side / 2;
+    var iMin = -1;
+    var jMin = -1;
+    var iMax = Math.ceil(params.width / spacing) + 1;
+    var jMax = Math.ceil(params.height / spacing) + 1;
+    var lattice = [];
+    var bars = [];
+    var geometries = [];
+    var i, j, y, cx, cy;
+
+    for (j = jMin - 1; j <= jMax + 1; j++) {
+      y = j * spacing;
+      lattice[j] = [];
+      for (i = iMin - 1; i <= iMax + 1; i++) {
+        lattice[j][i] = [i * spacing, y];
+      }
+    }
+
+    for (j = jMin - 1; j < jMax + 1; j++) {
+      bars[j] = [];
+      for (i = iMin - 1; i < iMax + 1; i++) {
+        cx = (i + 0.5) * spacing;
+        cy = (j + 0.5) * spacing;
+        bars[j][i] = (i + j) % 2 === 0 ? {
+          a: [cx - halfBar, cy],
+          b: [cx + halfBar, cy]
+        } : {
+          a: [cx, cy - halfBar],
+          b: [cx, cy + halfBar]
+        };
+      }
+    }
+
+    for (j = jMin; j < jMax; j++) {
+      for (i = iMin; i < iMax; i++) {
+        if ((i + j) % 2 === 0) {
+          appendCairoPentagon(geometries, [
+            lattice[j][i],
+            bars[j - 1][i].b,
+            lattice[j][i + 1],
+            bars[j][i].b,
+            bars[j][i].a
+          ]);
+          appendCairoPentagon(geometries, [
+            lattice[j + 1][i],
+            bars[j][i].a,
+            bars[j][i].b,
+            lattice[j + 1][i + 1],
+            bars[j + 1][i].a
+          ]);
+        } else {
+          appendCairoPentagon(geometries, [
+            lattice[j][i],
+            bars[j][i].a,
+            bars[j][i].b,
+            lattice[j + 1][i],
+            bars[j][i - 1].b
+          ]);
+          appendCairoPentagon(geometries, [
+            lattice[j][i + 1],
+            bars[j][i + 1].a,
+            lattice[j + 1][i + 1],
+            bars[j][i].b,
+            bars[j][i].a
+          ]);
+        }
+      }
+    }
+    return {type: 'GeometryCollection', geometries: geometries};
+  }
+
+  function appendCairoPentagon(geometries, verts) {
+    geometries.push({
+      type: 'Polygon',
+      coordinates: [verts.concat([verts[0]])]
+    });
+  }
+
   function subdivideHexagon(coords) {
     var center = getPolygonRingCenter(coords);
     return [
@@ -80084,7 +80211,7 @@ ${svg}
     return name == 'rectangle' || name == 'rectangles' || name == 'filter' && opts.cleanup;
   }
 
-  var version = "0.7.58";
+  var version = "0.7.60";
 
   // Parse command line args into commands and run them
   // Function takes an optional Node-style callback. A Promise is returned if no callback is given.