mesh2cad 0.3.0__tar.gz

mesh2cad-0.3.0/.github/workflows/ci.yml

@@ -0,0 +1,12 @@
+name: ci
+on: [push, pull_request]
+jobs:
+  test:
+    runs-on: ubuntu-latest
+    steps:
+      - uses: actions/checkout@v4
+      - uses: actions/setup-python@v5
+        with:
+          python-version: "3.12"
+      - run: pip install -e .[dev]
+      - run: pytest -q

mesh2cad-0.3.0/.github/workflows/publish.yml

@@ -0,0 +1,24 @@
+name: publish
+# Manual only: publishes to PyPI via Trusted Publishing (no token stored).
+on:
+  workflow_dispatch:
+jobs:
+  build:
+    runs-on: ubuntu-latest
+    steps:
+      - uses: actions/checkout@v4
+      - uses: actions/setup-python@v5
+        with: {python-version: "3.12"}
+      - run: pip install build twine && python -m build && twine check dist/*
+      - uses: actions/upload-artifact@v4
+        with: {name: dist, path: dist}
+  pypi:
+    needs: build
+    runs-on: ubuntu-latest
+    environment: pypi
+    permissions:
+      id-token: write
+    steps:
+      - uses: actions/download-artifact@v4
+        with: {name: dist, path: dist}
+      - uses: pypa/gh-action-pypi-publish@release/v1

mesh2cad-0.3.0/.gitignore

@@ -0,0 +1,10 @@
+*.egg-info/
+*.py[cod]
+*.step
+.pytest_cache/
+.venv/
+__pycache__/
+build/
+dist/
+out/
+docs/__pycache__/

mesh2cad-0.3.0/CHANGELOG.md

@@ -0,0 +1,29 @@
+# Changelog
+
+## Unreleased
+- Benchmark (`benchmarks/run.py`) now includes the cone parts (countersunk plate, tapered boss, chamfered boss) and a cone face column; README table refreshed.
+- Docs: Limits and `server.json` description now list cones; Glama score badge and MCP directory links in the README.
+
+## 0.3.0 — 2026-10-15
+- Cone primitive: `Cone`, `fit_cone` (constant normal-to-axis angle + tangent planes through the apex), `refine_cone`.
+- Segmentation fits cones for chamfers, countersinks and tapered bosses; a coaxial cylinder + cone component is split by normal-to-axis angle before fitting. Cones win over sphere/cylinder fits only when clearly better.
+- Exact B-rep conical frustum faces between two planes perpendicular to the axis, with plane ∩ cone circles on the planar faces. Hybrid fallback snaps vertices onto the cone.
+- Fidelity report counts `cone` faces. CLI `--no-cones` and MCP `cones` flag (default on).
+- Synthetic roundtrip: countersunk plate, 10° tapered boss, chamfered boss.
+
+## 0.2.0 — 2026-10-14
+- Sphere primitive: `Sphere`, `fit_sphere` (normal rays + algebraic fit), `refine_sphere`.
+- Segmentation emits `Region(surface=Sphere)` for leftover caps and for clusters of small patches that tile a sphere (noisy or faceted ball ends).
+- Exact B-rep spherical faces when the region is bounded by plane ∩ sphere circles (cap or zone). Hybrid fallback snaps vertices onto the sphere.
+- Fidelity report counts `sphere` faces. CLI `--no-spheres` and MCP `spheres` flag (default on).
+- Synthetic roundtrip: hemispherical dome, spherical pocket, and a sphere clipped to a zone between two planes.
+
+## 0.1.0 — 2026-10-14
+First release.
+- Mesh loading/cleaning for STL, OBJ, PLY (trimesh).
+- Segmentation into planes and cylinders (region growing, strip-chain cylinder detection, robust cylinder refinement, region merge/absorb).
+- Exact B-rep STEP export via OpenCascade for planar parts with through-holes and bosses; hybrid/faceted fallbacks.
+- Fidelity report (Hausdorff, RMS, volume error, face types, validity) as JSON/Markdown.
+- CLI: `mesh2cad convert`, `mesh2cad report`, `mesh2cad mcp`.
+- MCP server with `convert`, `analyze`, `report` tools.
+- Synthetic self-test benchmark (6 parts × 3 noise levels × 3 seeds, all exact) in CI.

mesh2cad-0.3.0/Dockerfile

@@ -0,0 +1,10 @@
+# Stdio MCP server for Glama / local Docker use.
+FROM python:3.12-slim
+
+WORKDIR /app
+COPY pyproject.toml README.md LICENSE ./
+COPY src ./src
+
+RUN pip install --no-cache-dir ".[mcp]"
+
+ENTRYPOINT ["mesh2cad", "mcp"]

mesh2cad-0.3.0/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2026 Emir Mrkaljevic
+
+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.

mesh2cad-0.3.0/PKG-INFO

@@ -0,0 +1,117 @@
+Metadata-Version: 2.5
+Name: mesh2cad
+Version: 0.3.0
+Summary: Convert STL/OBJ/PLY meshes into B-rep CAD (STEP) with primitive fitting
+License-Expression: MIT
+License-File: LICENSE
+Requires-Python: >=3.10
+Requires-Dist: cadquery-ocp
+Requires-Dist: networkx
+Requires-Dist: numpy
+Requires-Dist: rtree
+Requires-Dist: scipy
+Requires-Dist: trimesh
+Provides-Extra: dev
+Requires-Dist: cadquery; extra == 'dev'
+Requires-Dist: mcp; extra == 'dev'
+Requires-Dist: pytest; extra == 'dev'
+Provides-Extra: mcp
+Requires-Dist: mcp; extra == 'mcp'
+Description-Content-Type: text/markdown
+
+# mesh2cad
+[![ci](https://github.com/emirmrkaljevic/mesh2cad/actions/workflows/ci.yml/badge.svg)](https://github.com/emirmrkaljevic/mesh2cad/actions/workflows/ci.yml)
+[![emirmrkaljevic/mesh2cad MCP server](https://glama.ai/mcp/servers/emirmrkaljevic/mesh2cad/badges/score.svg)](https://glama.ai/mcp/servers/emirmrkaljevic/mesh2cad)
+
+Turn 3D meshes (STL/OBJ/PLY) into **real, editable CAD (STEP)** — fitted planes, cylinders, spheres, and cones instead of triangle soup. CLI + MCP server for AI agents. MIT, pure Python on top of OpenCascade.
+
+![demo](https://raw.githubusercontent.com/emirmrkaljevic/mesh2cad/main/docs/demo.gif)
+
+![before/after](https://raw.githubusercontent.com/emirmrkaljevic/mesh2cad/main/docs/before_after.png)
+
+*Left: 21k-triangle mesh with 0.1 mm noise. Right: the STEP mesh2cad writes — 10 planes + 2 cylinders, max deviation 0.05 mm from the true part.*
+
+## Install
+```
+pip install "mesh2cad[mcp] @ git+https://github.com/emirmrkaljevic/mesh2cad"
+```
+(PyPI release coming.)
+
+## Use
+```
+mesh2cad convert part.stl -o part.step --tol 0.1 --report report.md
+mesh2cad report part.stl --json
+```
+`--tol` ≈ the noise / accuracy of your mesh in model units (default 0.1 % of the bounding-box diagonal).
+
+Output: STEP + fidelity report — method (`exact` / `faceted-hybrid` / `faceted`), face types, solid validity, two-sided Hausdorff / RMS deviation, volume error, % area per surface type.
+
+### MCP server (Claude Desktop, Cursor, …)
+```json
+{ "mcpServers": { "mesh2cad": { "command": "mesh2cad", "args": ["mcp"] } } }
+```
+Tools: `convert(input_path, output_path?, tol?, cylinders?, spheres?, cones?)`, `analyze(input_path, tol?, cylinders?, spheres?, cones?)`, `report(input_path, tol?, cylinders?, spheres?, cones?)`. `spheres` and `cones` default on.
+
+Listed on [Glama](https://glama.ai/mcp/servers/emirmrkaljevic/mesh2cad) and the [TensorBlock MCP catalog](https://tensorblock.co/mcp/servers/github-emirmrkaljevic-mesh2cad-0f1c0e66).
+
+## How it works
+1. Load + clean mesh (trimesh); Taubin-smoothed copy for normals.
+2. Region growing into planes; strip chains and leftovers are fitted as cylinders, spheres, or cones (normal-ray / algebraic / constant-normal-angle fits, robust nonlinear refinement). A coaxial cylinder + cone component (countersink, chamfered pin) is split by normal-to-axis angle first. Small noisy patches that tile a sphere are merged into one spherical region. Split/redundant regions are merged and absorbed.
+3. Exact B-rep (OCP / OpenCascade): plane polygons from triple-plane intersections, circles from plane ∩ cylinder and plane ∩ sphere, plane ∩ cone, cylindrical faces, spherical caps and zones, conical frustums, sewn into a solid and validated.
+4. Fallbacks: hybrid faceted solid (vertices snapped to fitted surfaces, coplanar triangles merged), then raw faceted solid.
+5. Fidelity report against the input mesh.
+
+## Self-test benchmark
+Ground-truth cadquery parts → tessellated + subdivided (≤ 2 mm edges) → Gaussian vertex noise → mesh2cad → compared with the truth. 3 seeds, worst case shown (`python benchmarks/run.py`); the same cases run in CI.
+
+| part | noise σ (mm) | exact | faces (plane/cyl[/sph[/cone]]) | Hausdorff (mm) | RMS (mm) | vol err % | time s |
+|---|---|---|---|---|---|---|---|
+| cube | 0.0 | 3/3 | 6/0 | 0.000 | 0.000 | 0.00 | 0.7 |
+| cube | 0.05 | 3/3 | 6/0 | 0.019 | 0.006 | 0.07 | 0.6 |
+| cube | 0.2 | 3/3 | 6/0 | 0.075 | 0.026 | 0.34 | 0.6 |
+| l_bracket | 0.0 | 3/3 | 8/0 | 0.000 | 0.000 | 0.00 | 0.9 |
+| l_bracket | 0.05 | 3/3 | 8/0 | 0.040 | 0.006 | 0.03 | 1.0 |
+| l_bracket | 0.2 | 3/3 | 8/0 | 0.216 | 0.026 | 0.15 | 1.0 |
+| stepped_block | 0.0 | 3/3 | 11/0 | 0.000 | 0.000 | 0.00 | 2.1 |
+| stepped_block | 0.05 | 3/3 | 11/0 | 0.032 | 0.006 | 0.02 | 2.3 |
+| stepped_block | 0.2 | 3/3 | 11/0 | 0.136 | 0.022 | 0.07 | 2.1 |
+| plate_with_hole | 0.0 | 3/3 | 6/1 | 0.015 | 0.002 | 0.00 | 1.9 |
+| plate_with_hole | 0.05 | 3/3 | 6/1 | 0.028 | 0.006 | 0.05 | 1.9 |
+| plate_with_hole | 0.2 | 3/3 | 6/1 | 0.078 | 0.023 | 0.18 | 1.9 |
+| boss | 0.0 | 3/3 | 7/1 | 0.027 | 0.006 | 0.01 | 1.9 |
+| boss | 0.05 | 3/3 | 7/1 | 0.031 | 0.007 | 0.04 | 2.5 |
+| boss | 0.2 | 3/3 | 7/1 | 0.090 | 0.016 | 0.07 | 3.5 |
+| ribbed_plate | 0.0 | 3/3 | 10/2 | 0.013 | 0.002 | 0.00 | 2.3 |
+| ribbed_plate | 0.05 | 3/3 | 10/2 | 0.028 | 0.005 | 0.03 | 2.4 |
+| ribbed_plate | 0.2 | 3/3 | 10/2 | 0.108 | 0.020 | 0.12 | 2.3 |
+| spherical_dome | 0.0 | 3/3 | 6/0/1 | 0.065 | 0.007 | 0.04 | 1.9 |
+| spherical_dome | 0.05 | 3/3 | 6/0/1 | 0.067 | 0.008 | 0.06 | 2.0 |
+| spherical_dome | 0.2 | 3/3 | 6/0/1 | 0.241 | 0.028 | 0.25 | 1.9 |
+| spherical_pocket | 0.0 | 3/3 | 6/0/1 | 0.046 | 0.004 | 0.02 | 2.2 |
+| spherical_pocket | 0.05 | 3/3 | 6/0/1 | 0.057 | 0.005 | 0.03 | 2.1 |
+| spherical_pocket | 0.2 | 3/3 | 6/0/1 | 0.075 | 0.015 | 0.13 | 2.1 |
+| countersunk_plate | 0.0 | 3/3 | 6/1/0/1 | 0.039 | 0.003 | 0.01 | 2.2 |
+| countersunk_plate | 0.05 | 3/3 | 6/1/0/1 | 0.055 | 0.009 | 0.07 | 2.1 |
+| countersunk_plate | 0.2 | 0/3 | 16322/0, 16324/0 | 0.745 | 0.130 | 0.17 | 41.6 |
+| tapered_boss | 0.0 | 3/3 | 7/0/0/1 | 0.086 | 0.011 | 0.05 | 1.8 |
+| tapered_boss | 0.05 | 3/3 | 7/0/0/1 | 0.095 | 0.012 | 0.05 | 1.7 |
+| tapered_boss | 0.2 | 3/3 | 7/0/0/1 | 0.123 | 0.025 | 0.22 | 1.7 |
+| chamfered_boss | 0.0 | 3/3 | 7/1/0/1 | 0.053 | 0.005 | 0.02 | 2.9 |
+| chamfered_boss | 0.05 | 1/3 | 16461/0, 16473/0, 7/1/0/1 | 0.160 | 0.035 | 0.08 | 31.6 |
+| chamfered_boss | 0.2 | 0/3 | 16565/0, 16569/0, 16570/0 | 0.703 | 0.133 | 0.26 | 34.8 |
+
+"exact" = valid B-rep with the true face count. All rows are the worst of 3 seeds from `python benchmarks/run.py` (v0.3.0). Rows with fewer than 3/3 exact fell back to a faceted solid (face count = triangles): noisy narrow countersink/chamfer bands don't constrain a cone yet, so CI checks those parts only where they are exact.
+
+## Limits (v0.3)
+Planes, cylinders, spheres, and cones. Cylinders must end on planes perpendicular to their axis (through-holes, bosses, pins). Spheres are exact when each boundary loop is a single plane ∩ sphere circle (domes, spherical bosses, spherical pockets, and zones between two parallel planes). A sphere that runs into an edge or a cylinder falls back to the hybrid solid — vertices are still snapped onto the fitted sphere, but the face stays faceted. Cones are exact when they run between two planes perpendicular to their axis: countersinks, chamfered holes and bosses, drafted/tapered bosses. Cones may also meet a coaxial cylinder (countersunk hole, chamfered pin). A cone that ends in its apex falls back to the hybrid solid, and countersinks are reliable up to ~0.05 mm noise, narrow chamfers on clean meshes only for now (wide drafted walls up to 0.2 mm). No fillets or free-form yet.
+
+## Development
+```
+pip install -e .[dev]
+pytest
+python docs/make_demo.py   # regenerate images (needs matplotlib, pillow)
+```
+
+License: MIT
+
+<!-- mcp-name: io.github.emirmrkaljevic/mesh2cad -->

mesh2cad-0.3.0/README.md

@@ -0,0 +1,96 @@
+# mesh2cad
+[![ci](https://github.com/emirmrkaljevic/mesh2cad/actions/workflows/ci.yml/badge.svg)](https://github.com/emirmrkaljevic/mesh2cad/actions/workflows/ci.yml)
+[![emirmrkaljevic/mesh2cad MCP server](https://glama.ai/mcp/servers/emirmrkaljevic/mesh2cad/badges/score.svg)](https://glama.ai/mcp/servers/emirmrkaljevic/mesh2cad)
+
+Turn 3D meshes (STL/OBJ/PLY) into **real, editable CAD (STEP)** — fitted planes, cylinders, spheres, and cones instead of triangle soup. CLI + MCP server for AI agents. MIT, pure Python on top of OpenCascade.
+
+![demo](https://raw.githubusercontent.com/emirmrkaljevic/mesh2cad/main/docs/demo.gif)
+
+![before/after](https://raw.githubusercontent.com/emirmrkaljevic/mesh2cad/main/docs/before_after.png)
+
+*Left: 21k-triangle mesh with 0.1 mm noise. Right: the STEP mesh2cad writes — 10 planes + 2 cylinders, max deviation 0.05 mm from the true part.*
+
+## Install
+```
+pip install "mesh2cad[mcp] @ git+https://github.com/emirmrkaljevic/mesh2cad"
+```
+(PyPI release coming.)
+
+## Use
+```
+mesh2cad convert part.stl -o part.step --tol 0.1 --report report.md
+mesh2cad report part.stl --json
+```
+`--tol` ≈ the noise / accuracy of your mesh in model units (default 0.1 % of the bounding-box diagonal).
+
+Output: STEP + fidelity report — method (`exact` / `faceted-hybrid` / `faceted`), face types, solid validity, two-sided Hausdorff / RMS deviation, volume error, % area per surface type.
+
+### MCP server (Claude Desktop, Cursor, …)
+```json
+{ "mcpServers": { "mesh2cad": { "command": "mesh2cad", "args": ["mcp"] } } }
+```
+Tools: `convert(input_path, output_path?, tol?, cylinders?, spheres?, cones?)`, `analyze(input_path, tol?, cylinders?, spheres?, cones?)`, `report(input_path, tol?, cylinders?, spheres?, cones?)`. `spheres` and `cones` default on.
+
+Listed on [Glama](https://glama.ai/mcp/servers/emirmrkaljevic/mesh2cad) and the [TensorBlock MCP catalog](https://tensorblock.co/mcp/servers/github-emirmrkaljevic-mesh2cad-0f1c0e66).
+
+## How it works
+1. Load + clean mesh (trimesh); Taubin-smoothed copy for normals.
+2. Region growing into planes; strip chains and leftovers are fitted as cylinders, spheres, or cones (normal-ray / algebraic / constant-normal-angle fits, robust nonlinear refinement). A coaxial cylinder + cone component (countersink, chamfered pin) is split by normal-to-axis angle first. Small noisy patches that tile a sphere are merged into one spherical region. Split/redundant regions are merged and absorbed.
+3. Exact B-rep (OCP / OpenCascade): plane polygons from triple-plane intersections, circles from plane ∩ cylinder and plane ∩ sphere, plane ∩ cone, cylindrical faces, spherical caps and zones, conical frustums, sewn into a solid and validated.
+4. Fallbacks: hybrid faceted solid (vertices snapped to fitted surfaces, coplanar triangles merged), then raw faceted solid.
+5. Fidelity report against the input mesh.
+
+## Self-test benchmark
+Ground-truth cadquery parts → tessellated + subdivided (≤ 2 mm edges) → Gaussian vertex noise → mesh2cad → compared with the truth. 3 seeds, worst case shown (`python benchmarks/run.py`); the same cases run in CI.
+
+| part | noise σ (mm) | exact | faces (plane/cyl[/sph[/cone]]) | Hausdorff (mm) | RMS (mm) | vol err % | time s |
+|---|---|---|---|---|---|---|---|
+| cube | 0.0 | 3/3 | 6/0 | 0.000 | 0.000 | 0.00 | 0.7 |
+| cube | 0.05 | 3/3 | 6/0 | 0.019 | 0.006 | 0.07 | 0.6 |
+| cube | 0.2 | 3/3 | 6/0 | 0.075 | 0.026 | 0.34 | 0.6 |
+| l_bracket | 0.0 | 3/3 | 8/0 | 0.000 | 0.000 | 0.00 | 0.9 |
+| l_bracket | 0.05 | 3/3 | 8/0 | 0.040 | 0.006 | 0.03 | 1.0 |
+| l_bracket | 0.2 | 3/3 | 8/0 | 0.216 | 0.026 | 0.15 | 1.0 |
+| stepped_block | 0.0 | 3/3 | 11/0 | 0.000 | 0.000 | 0.00 | 2.1 |
+| stepped_block | 0.05 | 3/3 | 11/0 | 0.032 | 0.006 | 0.02 | 2.3 |
+| stepped_block | 0.2 | 3/3 | 11/0 | 0.136 | 0.022 | 0.07 | 2.1 |
+| plate_with_hole | 0.0 | 3/3 | 6/1 | 0.015 | 0.002 | 0.00 | 1.9 |
+| plate_with_hole | 0.05 | 3/3 | 6/1 | 0.028 | 0.006 | 0.05 | 1.9 |
+| plate_with_hole | 0.2 | 3/3 | 6/1 | 0.078 | 0.023 | 0.18 | 1.9 |
+| boss | 0.0 | 3/3 | 7/1 | 0.027 | 0.006 | 0.01 | 1.9 |
+| boss | 0.05 | 3/3 | 7/1 | 0.031 | 0.007 | 0.04 | 2.5 |
+| boss | 0.2 | 3/3 | 7/1 | 0.090 | 0.016 | 0.07 | 3.5 |
+| ribbed_plate | 0.0 | 3/3 | 10/2 | 0.013 | 0.002 | 0.00 | 2.3 |
+| ribbed_plate | 0.05 | 3/3 | 10/2 | 0.028 | 0.005 | 0.03 | 2.4 |
+| ribbed_plate | 0.2 | 3/3 | 10/2 | 0.108 | 0.020 | 0.12 | 2.3 |
+| spherical_dome | 0.0 | 3/3 | 6/0/1 | 0.065 | 0.007 | 0.04 | 1.9 |
+| spherical_dome | 0.05 | 3/3 | 6/0/1 | 0.067 | 0.008 | 0.06 | 2.0 |
+| spherical_dome | 0.2 | 3/3 | 6/0/1 | 0.241 | 0.028 | 0.25 | 1.9 |
+| spherical_pocket | 0.0 | 3/3 | 6/0/1 | 0.046 | 0.004 | 0.02 | 2.2 |
+| spherical_pocket | 0.05 | 3/3 | 6/0/1 | 0.057 | 0.005 | 0.03 | 2.1 |
+| spherical_pocket | 0.2 | 3/3 | 6/0/1 | 0.075 | 0.015 | 0.13 | 2.1 |
+| countersunk_plate | 0.0 | 3/3 | 6/1/0/1 | 0.039 | 0.003 | 0.01 | 2.2 |
+| countersunk_plate | 0.05 | 3/3 | 6/1/0/1 | 0.055 | 0.009 | 0.07 | 2.1 |
+| countersunk_plate | 0.2 | 0/3 | 16322/0, 16324/0 | 0.745 | 0.130 | 0.17 | 41.6 |
+| tapered_boss | 0.0 | 3/3 | 7/0/0/1 | 0.086 | 0.011 | 0.05 | 1.8 |
+| tapered_boss | 0.05 | 3/3 | 7/0/0/1 | 0.095 | 0.012 | 0.05 | 1.7 |
+| tapered_boss | 0.2 | 3/3 | 7/0/0/1 | 0.123 | 0.025 | 0.22 | 1.7 |
+| chamfered_boss | 0.0 | 3/3 | 7/1/0/1 | 0.053 | 0.005 | 0.02 | 2.9 |
+| chamfered_boss | 0.05 | 1/3 | 16461/0, 16473/0, 7/1/0/1 | 0.160 | 0.035 | 0.08 | 31.6 |
+| chamfered_boss | 0.2 | 0/3 | 16565/0, 16569/0, 16570/0 | 0.703 | 0.133 | 0.26 | 34.8 |
+
+"exact" = valid B-rep with the true face count. All rows are the worst of 3 seeds from `python benchmarks/run.py` (v0.3.0). Rows with fewer than 3/3 exact fell back to a faceted solid (face count = triangles): noisy narrow countersink/chamfer bands don't constrain a cone yet, so CI checks those parts only where they are exact.
+
+## Limits (v0.3)
+Planes, cylinders, spheres, and cones. Cylinders must end on planes perpendicular to their axis (through-holes, bosses, pins). Spheres are exact when each boundary loop is a single plane ∩ sphere circle (domes, spherical bosses, spherical pockets, and zones between two parallel planes). A sphere that runs into an edge or a cylinder falls back to the hybrid solid — vertices are still snapped onto the fitted sphere, but the face stays faceted. Cones are exact when they run between two planes perpendicular to their axis: countersinks, chamfered holes and bosses, drafted/tapered bosses. Cones may also meet a coaxial cylinder (countersunk hole, chamfered pin). A cone that ends in its apex falls back to the hybrid solid, and countersinks are reliable up to ~0.05 mm noise, narrow chamfers on clean meshes only for now (wide drafted walls up to 0.2 mm). No fillets or free-form yet.
+
+## Development
+```
+pip install -e .[dev]
+pytest
+python docs/make_demo.py   # regenerate images (needs matplotlib, pillow)
+```
+
+License: MIT
+
+<!-- mcp-name: io.github.emirmrkaljevic/mesh2cad -->

mesh2cad-0.3.0/benchmarks/run.py

@@ -0,0 +1,44 @@
+"""Synthetic benchmark: cadquery parts -> mesh -> noise -> mesh2cad -> compare to truth.
+
+Usage: python benchmarks/run.py [--seeds 3]
+"""
+import argparse
+import pathlib
+import sys
+import time
+
+sys.path.insert(0, str(pathlib.Path(__file__).resolve().parents[1] / "tests"))
+from parts import CONES, CURVED, PLANAR, SPHERES, to_mesh  # noqa: E402
+
+from mesh2cad.convert import convert_mesh  # noqa: E402
+from mesh2cad.io import add_noise  # noqa: E402
+
+ap = argparse.ArgumentParser()
+ap.add_argument("--seeds", type=int, default=3)
+a = ap.parse_args()
+print("| part | noise σ (mm) | exact | faces (plane/cyl[/sph[/cone]]) | Hausdorff (mm) | RMS (mm) | vol err % | time s |")
+print("|---|---|---|---|---|---|---|---|")
+
+
+def _faces(f):
+    parts = [str(f.get("plane", 0)), str(f.get("cylinder", 0))]
+    if f.get("sphere") or f.get("cone") or f.get("other"):
+        parts.append(str(f.get("sphere", 0)))
+    if f.get("cone"):
+        parts.append(str(f.get("cone", 0)))
+    s = "/".join(parts)
+    return s + ("+" if f.get("other") else "")
+
+
+for name, fn in {**PLANAR, **CURVED, **SPHERES, **CONES}.items():
+    ref = to_mesh(fn())
+    for sigma in (0.0, 0.05, 0.2):
+        ok, H, R, V, T, F = 0, [], [], [], [], set()
+        for seed in range(a.seeds):
+            t = time.time()
+            r = convert_mesh(add_noise(ref, sigma, seed=seed), tol=max(sigma, 0.01), ref=ref).report
+            T.append(time.time() - t)
+            ok += r["method"] == "exact" and r["output"]["valid_solid"]
+            F.add(_faces(r["output"]["faces"]))
+            H.append(r["deviation"]["hausdorff"]); R.append(r["deviation"]["rms"]); V.append(r["volume_error_pct"])
+        print(f"| {name} | {sigma} | {ok}/{a.seeds} | {', '.join(sorted(F))} | {max(H):.3f} | {max(R):.3f} | {max(V):.2f} | {max(T):.1f} |")

mesh2cad-0.3.0/docs/make_demo.py

@@ -0,0 +1,84 @@
+"""Render README before/after images and demo GIF (needs matplotlib, pillow)."""
+import pathlib
+import sys
+
+import matplotlib
+matplotlib.use("Agg")
+import matplotlib.pyplot as plt  # noqa: E402
+import numpy as np  # noqa: E402
+from mpl_toolkits.mplot3d.art3d import Poly3DCollection  # noqa: E402
+from PIL import Image  # noqa: E402
+
+ROOT = pathlib.Path(__file__).resolve().parents[1]
+sys.path.insert(0, str(ROOT / "tests"))
+import cadquery as cq  # noqa: E402
+from parts import ribbed_plate as demo_part, to_mesh  # noqa: E402
+
+from mesh2cad.convert import convert_mesh  # noqa: E402
+from mesh2cad.io import add_noise  # noqa: E402
+from mesh2cad.report import shape_to_mesh  # noqa: E402
+from OCP.TopAbs import TopAbs_FACE  # noqa: E402
+from OCP.TopExp import TopExp_Explorer  # noqa: E402
+
+OUT = ROOT / "docs"
+
+
+def faces_of(shape):
+    exp, out = TopExp_Explorer(shape, TopAbs_FACE), []
+    while exp.More():
+        out.append(shape_to_mesh(exp.Current(), 0.05))
+        exp.Next()
+    return out
+
+
+def draw(ax, polys, colors, edges, title, azim):
+    ax.clear()
+    pc = Poly3DCollection(polys, facecolors=colors, edgecolors=colors if edges == "same" else edges,
+                          linewidths=0.15 if edges != "same" else 0.4)
+    ax.add_collection3d(pc)
+    ax.set_xlim(-20, 20); ax.set_ylim(-20, 20); ax.set_zlim(-10, 14)
+    ax.set_box_aspect((1, 1, 0.65))
+    ax.view_init(28, azim)
+    ax.set_axis_off()
+    ax.set_title(title, fontsize=11)
+
+
+def shade(normals, base):
+    light = np.array([0.4, -0.5, 0.75]); light /= np.linalg.norm(light)
+    k = 0.45 + 0.55 * np.clip(normals @ light, 0, 1)
+    return np.clip(np.asarray(base)[None, :3] * k[:, None], 0, 1)
+
+
+def main():
+    ref = to_mesh(demo_part())
+    noisy = add_noise(ref, 0.1, seed=0)
+    res = convert_mesh(noisy, tol=0.1, ref=ref)
+    r = res.report
+    print(r["method"], r["output"]["faces"], r["deviation"])
+    before = (noisy.triangles, shade(noisy.face_normals, (0.75, 0.75, 0.8)), (0.2, 0.2, 0.25, 0.6))
+    cmap = plt.get_cmap("tab20")
+    P, C = [], []
+    for i, fm in enumerate(faces_of(res.shape)):
+        P += list(fm.triangles)
+        C += list(shade(fm.face_normals, cmap(i % 20)))
+    after = (np.array(P), np.array(C), "same")
+    f = r["output"]["faces"]
+    t_before = f"input: {len(noisy.faces)} noisy triangles"
+    t_after = f"mesh2cad STEP: {f.get('plane', 0)} planes + {f.get('cylinder', 0)} cylinders"
+    fig = plt.figure(figsize=(10, 4.4), dpi=90)
+    a1, a2 = fig.add_subplot(121, projection="3d"), fig.add_subplot(122, projection="3d")
+    draw(a1, *before, t_before, -60); draw(a2, *after, t_after, -60)
+    fig.tight_layout(); fig.savefig(OUT / "before_after.png"); plt.close(fig)
+    frames = []
+    for az in range(-60, 300, 15):
+        fig = plt.figure(figsize=(8, 3.6), dpi=70)
+        a1, a2 = fig.add_subplot(121, projection="3d"), fig.add_subplot(122, projection="3d")
+        draw(a1, *before, t_before, az); draw(a2, *after, t_after, az)
+        fig.tight_layout(); fig.canvas.draw()
+        frames.append(Image.fromarray(np.asarray(fig.canvas.buffer_rgba())[..., :3]).convert("P", palette=Image.ADAPTIVE))
+        plt.close(fig)
+    frames[0].save(OUT / "demo.gif", save_all=True, append_images=frames[1:], duration=120, loop=0, optimize=True)
+
+
+if __name__ == "__main__":
+    main()

mesh2cad-0.3.0/glama.json

@@ -0,0 +1,4 @@
+{
+  "$schema": "https://glama.ai/mcp/schemas/server.json",
+  "maintainers": ["emirmrkaljevic"]
+}

mesh2cad-0.3.0/pyproject.toml

@@ -0,0 +1,19 @@
+[build-system]
+requires = ["hatchling"]
+build-backend = "hatchling.build"
+
+[project]
+name = "mesh2cad"
+version = "0.3.0"
+description = "Convert STL/OBJ/PLY meshes into B-rep CAD (STEP) with primitive fitting"
+readme = "README.md"
+license = "MIT"
+requires-python = ">=3.10"
+dependencies = ["trimesh", "rtree", "networkx", "numpy", "scipy", "cadquery-ocp"]
+
+[project.optional-dependencies]
+mcp = ["mcp"]
+dev = ["pytest", "cadquery", "mcp"]
+
+[project.scripts]
+mesh2cad = "mesh2cad.cli:main"

mesh2cad-0.3.0/server.json

@@ -0,0 +1,16 @@
+{
+  "$schema": "https://static.modelcontextprotocol.io/schemas/2025-09-29/server.schema.json",
+  "name": "io.github.emirmrkaljevic/mesh2cad",
+  "description": "Mesh (STL/OBJ/PLY) to editable STEP CAD: fitted planes, cylinders, spheres, cones + fidelity report",
+  "repository": {"url": "https://github.com/emirmrkaljevic/mesh2cad", "source": "github"},
+  "version": "0.3.0",
+  "packages": [
+    {
+      "registryType": "pypi",
+      "identifier": "mesh2cad",
+      "version": "0.3.0",
+      "transport": {"type": "stdio"},
+      "packageArguments": [{"type": "positional", "value": "mcp"}]
+    }
+  ]
+}

mesh2cad-0.3.0/src/mesh2cad/__init__.py

@@ -0,0 +1 @@
+__version__ = "0.3.0"