smoothify 0.3.2__tar.gz → 0.3.3__tar.gz

{smoothify-0.3.2 → smoothify-0.3.3}/.github/workflows/ci.yml

@@ -33,8 +33,11 @@ jobs:
         run: uv run mypy smoothify/
 
       - name: Test
-        run: uv run pytest tests/ -x -q
+        # -n 2 overrides the local `-n auto` default: GitHub runners have 4
+        # vCPUs and some tests spawn smoothify's own worker pool, so a smaller
+        # fan-out avoids oversubscribing the runner.
+        run: uv run pytest tests/ -n 2 -x -q
 
       - name: Notebook smoke test
         if: matrix.python-version == '3.11'
-        run: uv run pytest --nbmake examples/*.ipynb -q
+        run: uv run pytest --nbmake examples/*.ipynb -n 2 -q

{smoothify-0.3.2 → smoothify-0.3.3}/CHANGELOG.md

@@ -4,6 +4,15 @@ All notable changes to this project will be documented in this file.
 
 ## [Unreleased]
 
+## [0.3.3] - 2026-07-02
+
+### Changed
+- Start-point variants are now merged by a per-point median consensus instead of a geometric union. Douglas-Peucker is start-vertex dependent, so the four rotated variants disagree slightly about where a feature's vertices land; the union superimposed those disagreements, whereas the median resamples the variants to a common point count, brings them into phase via an FFT cross-correlation, and takes the coordinate-wise median — a start-invariant consensus that outvotes a single bad variant instead of averaging toward it. This is also faster than the union it replaces (~1.3x end-to-end at the default 3 iterations, ~1.4x at 5, on `examples/Water.gpkg`). The median sits at the consensus interior rather than the outer envelope, so the merged ring is marginally smaller before area preservation restores it (with `preserve_area=True`); with `preserve_area=False` the small reduction is left as-is.
+- Faster smoothing with no meaningful change to output. Start-point variants now begin at evenly spaced arc-length positions, computed by interpolating a single new start vertex, rather than rotating to a vertex index on a fully densified ring. The old path segmentized the whole ring only so that index-based rotation would land on evenly spaced points — but every vertex it added was collinear and immediately stripped by the following Douglas-Peucker. Feeding DP just the original vertices plus one cuts its input roughly 3.4x and the simplify step ~3.5x on `examples/Water.gpkg` (total single-core CPU work ~4.4s → ~2.7s, about 1.25x faster end-to-end at the default 3 iterations). The per-point median phase-alignment and Chaikin corner cutting also replace `numpy.roll` with a slice-based cyclic shift (bit-identical, fewer allocations). Output is unchanged within tolerance (per-polygon area drift ≤ 0.01%, worst concave turn unchanged).
+
+### Fixed
+- Eliminated the residual "dimple" fold that appeared where the start-point variants disagreed about a shallow feature near the simplify tolerance. Their union superimposed the disagreement as a double peak with a valley that the area-preservation buffer then sharpened into a fold slipping under the concave-turn seal threshold. The per-point median consensus (see Changed) resolves the variants to one boundary instead of overlaying them, removing the artifact (raster-noise fuzz harness: 4 → 0 residual folds; `examples/Water.gpkg`: 8 → 0).
+
 ## [0.3.2] - 2026-06-30
 
 ### Changed

{smoothify-0.3.2/smoothify.egg-info → smoothify-0.3.3}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: smoothify
-Version: 0.3.2
+Version: 0.3.3
 Summary: Transform pixelated geometries from raster data into smooth natural looking features
 Author-email: Nick Wright <nicholas.wright@dpird.wa.gov.au>
 License-Expression: MIT
@@ -229,22 +229,21 @@ smoothed = smoothify(
 
 ## How It Works
 
-Smoothify uses an advanced multi-step smoothing pipeline:
+Smoothify uses an advanced multi-step smoothing pipeline. The numbered steps below correspond to the panels in the figure:
 
 <p align="left">
   <img src="https://raw.githubusercontent.com/DPIRD-DMA/Smoothify/main/images/pipeline_steps.png" alt="Smoothify pipeline steps" width="800">
 </p>
 
 
-1. Joins touching holes (for Polygons, when `merge_holes=True`) so they smooth as one opening
-2. Adds intermediate vertices along line segments (segmentize)
-3. Generates multiple rotated variants (for Polygons) to avoid artifacts
-4. Simplifies each variant to remove noise
-5. Applies Chaikin corner cutting to smooth
-6. Merges all variants via union to eliminate start-point artifacts
-7. Applies final smoothing pass
-8. Optionally restores original area via buffering (for Polygons)
-9. Detects and repairs any sharp folds left by features near the smoothing scale (e.g. one-pixel-wide arms), using a small morphological opening/closing bounded at `segment_length / 4`
+1. **Pixelated input** — a polygon straight from raster-to-vector conversion, with a stair-stepped boundary
+2. **Multiple variants** (for Polygons) that start at evenly spaced arc-length positions, each simplified to strip staircase noise, so no artifact is tied to a fixed start vertex
+3. **Chaikin corner cutting** applied to each variant
+4. **Per-point median merge** — a start-invariant consensus that resolves the variants' disagreements
+5. **Final smoothing pass** on the merged result
+6. **Restore original area** via buffering (for Polygons, when `preserve_area=True`)
+
+Two steps are not shown in the figure: before step 2, touching holes are joined (for Polygons, when `merge_holes=True`) so they smooth as one opening; and after step 6, the result is checked for sharp concave folds left by features near the smoothing scale (e.g. one-pixel-wide arms) and repaired with a small morphological opening/closing bounded at `segment_length / 4`.
 
 ## Invalid Geometries
 
@@ -278,17 +277,17 @@ Smoothify uses [pytest](https://pytest.org/). After cloning the repository, inst
 # Install dependencies (including the dev group)
 uv sync
 
-# Run all tests
+# Run all tests (parallelised across CPU cores by default via pytest-xdist)
 uv run pytest tests/
 
 # Run with coverage
 uv run pytest tests/ --cov=smoothify --cov-report=html
 
-# Run a single test
-uv run pytest tests/test_chaikin.py::TestChaikinCornerCutting::test_simple_square_polygon
+# Run a single test (add `-n 0` to disable parallelism for clearer output)
+uv run pytest tests/test_chaikin.py::TestChaikinCornerCutting::test_simple_square_polygon -n 0
 ```
 
-If you prefer not to use uv, install the dev dependencies into your environment and run `pytest tests/` directly.
+The suite runs in parallel by default (`-n auto` in `pytest.ini`); pass `-n 0` to run serially when debugging. If you prefer not to use uv, install the dev dependencies into your environment and run `pytest tests/` directly.
 
 ## Contributing
 

{smoothify-0.3.2 → smoothify-0.3.3}/README.md

@@ -196,22 +196,21 @@ smoothed = smoothify(
 
 ## How It Works
 
-Smoothify uses an advanced multi-step smoothing pipeline:
+Smoothify uses an advanced multi-step smoothing pipeline. The numbered steps below correspond to the panels in the figure:
 
 <p align="left">
   <img src="https://raw.githubusercontent.com/DPIRD-DMA/Smoothify/main/images/pipeline_steps.png" alt="Smoothify pipeline steps" width="800">
 </p>
 
 
-1. Joins touching holes (for Polygons, when `merge_holes=True`) so they smooth as one opening
-2. Adds intermediate vertices along line segments (segmentize)
-3. Generates multiple rotated variants (for Polygons) to avoid artifacts
-4. Simplifies each variant to remove noise
-5. Applies Chaikin corner cutting to smooth
-6. Merges all variants via union to eliminate start-point artifacts
-7. Applies final smoothing pass
-8. Optionally restores original area via buffering (for Polygons)
-9. Detects and repairs any sharp folds left by features near the smoothing scale (e.g. one-pixel-wide arms), using a small morphological opening/closing bounded at `segment_length / 4`
+1. **Pixelated input** — a polygon straight from raster-to-vector conversion, with a stair-stepped boundary
+2. **Multiple variants** (for Polygons) that start at evenly spaced arc-length positions, each simplified to strip staircase noise, so no artifact is tied to a fixed start vertex
+3. **Chaikin corner cutting** applied to each variant
+4. **Per-point median merge** — a start-invariant consensus that resolves the variants' disagreements
+5. **Final smoothing pass** on the merged result
+6. **Restore original area** via buffering (for Polygons, when `preserve_area=True`)
+
+Two steps are not shown in the figure: before step 2, touching holes are joined (for Polygons, when `merge_holes=True`) so they smooth as one opening; and after step 6, the result is checked for sharp concave folds left by features near the smoothing scale (e.g. one-pixel-wide arms) and repaired with a small morphological opening/closing bounded at `segment_length / 4`.
 
 ## Invalid Geometries
 
@@ -245,17 +244,17 @@ Smoothify uses [pytest](https://pytest.org/). After cloning the repository, inst
 # Install dependencies (including the dev group)
 uv sync
 
-# Run all tests
+# Run all tests (parallelised across CPU cores by default via pytest-xdist)
 uv run pytest tests/
 
 # Run with coverage
 uv run pytest tests/ --cov=smoothify --cov-report=html
 
-# Run a single test
-uv run pytest tests/test_chaikin.py::TestChaikinCornerCutting::test_simple_square_polygon
+# Run a single test (add `-n 0` to disable parallelism for clearer output)
+uv run pytest tests/test_chaikin.py::TestChaikinCornerCutting::test_simple_square_polygon -n 0
 ```
 
-If you prefer not to use uv, install the dev dependencies into your environment and run `pytest tests/` directly.
+The suite runs in parallel by default (`-n auto` in `pytest.ini`); pass `-n 0` to run serially when debugging. If you prefer not to use uv, install the dev dependencies into your environment and run `pytest tests/` directly.
 
 ## Contributing