From 30817f2626172d8a8b38c01ad5d721d47c19b1c2 Mon Sep 17 00:00:00 2001
From: Alex Matulich <anachronist@gmail.com>
Date: Wed, 21 May 2025 13:48:26 -0700
Subject: [PATCH 1/4] Add size value to top of img2scad.html output file

---
 scripts/img2scad.html | 8 +++++---
 1 file changed, 5 insertions(+), 3 deletions(-)

diff --git a/scripts/img2scad.html b/scripts/img2scad.html
index 4830684c..6e3c72a0 100644
--- a/scripts/img2scad.html
+++ b/scripts/img2scad.html
@@ -8,8 +8,9 @@ Version 6: 23 April 2025 - added cropping UI
 Version 7: 25 April 2025 - added contrast and threshold sliders
 Version 8: 26 April 2025 - added file size estimate to output section
 Version 9: 20 May 2025 - improved appearance UI, added Sobel edge detection
+Verskl  10: 21 May 2025 - Added array_name_size value at top of output file
 -->
-  <title>Image to OpenSCAD array, v9</title><!-- REMEMBER TO CHANGE VERSION -->
+  <title>Image to OpenSCAD array, v10</title><!-- REMEMBER TO CHANGE VERSION -->
   <meta charset="UTF-8">
   <style>
     body { font-family: sans-serif; padding-left:1em; padding-right:1em;}
@@ -772,9 +773,10 @@ Alpha channel is ignored. After processing the image as desired, you may save it
         return " [" + row.map(val => useUnit ? parseFloat((val/255.0).toFixed(3)) : val).join(",") + "]";
       }).join(",\n");
       const introcomment = " = [ // " + cropDim.width + "×" + cropDim.height + "\n";
-      const dimSuffix = "_"+cropDim.width + "x" + cropDim.height
-      const openscadArray = (arrayName.value.length>0 ? arrayName.value : 'image_array') + introcomment + arrayContent + "\n];";
+	  const sizevar = (arrayName.value.length>0 ? arrayName.value : 'image_array')+"_size = [" + cropDim.width + "," + cropDim.height + "];\n";
+      const openscadArray = sizevar + (arrayName.value.length>0 ? arrayName.value : 'image_array') + introcomment + arrayContent + "\n];";
       const blob = new Blob([openscadArray], { type: "text/plain" });
+      const dimSuffix = "_"+cropDim.width + "x" + cropDim.height;
       let filename = (arrayName.value.length>0 ? arrayName.value : "image_array") + dimSuffix + '.scad';
       if (window.showSaveFilePicker) {
         saveWithFilePicker(blob, filename);

From f5bc61841b1f1b5a6ba2d9079c9e07acaa90e153 Mon Sep 17 00:00:00 2001
From: Alex Matulich <anachronist@gmail.com>
Date: Thu, 22 May 2025 20:02:50 -0700
Subject: [PATCH 2/4] Added sharpening filter, fixed edge artifact, fixed
 contour topic in isosurface.scad

---
 isosurface.scad       |  2 +-
 scripts/img2scad.html | 82 +++++++++++++++++++++++++++----------------
 2 files changed, 52 insertions(+), 32 deletions(-)

diff --git a/isosurface.scad b/isosurface.scad
index 599df2b9..081805fe 100644
--- a/isosurface.scad
+++ b/isosurface.scad
@@ -3523,7 +3523,7 @@ function _showstats_isosurface(voxsize, bbox, isoval, cubes, triangles, faces) =
 // Function&Module: contour()
 // Synopsis: Creates a 2D contour from a function or array of values.
 // SynTags: Geom,Path,Region
-// Topics: Isosurfaces, Path Generators (2D), Regions
+// Topics: Contours, Path Generators (2D), Regions
 // Usage: As a module
 //   contour(f, isovalue, bounding_box, pixel_size, [pixel_count=], [use_centers=], [smoothing=], [exact_bounds=], [show_stats=], [show_box=], ...) [ATTACHMENTS];
 // Usage: As a function
diff --git a/scripts/img2scad.html b/scripts/img2scad.html
index 6e3c72a0..2487be5e 100644
--- a/scripts/img2scad.html
+++ b/scripts/img2scad.html
@@ -8,9 +8,10 @@ Version 6: 23 April 2025 - added cropping UI
 Version 7: 25 April 2025 - added contrast and threshold sliders
 Version 8: 26 April 2025 - added file size estimate to output section
 Version 9: 20 May 2025 - improved appearance UI, added Sobel edge detection
-Verskl  10: 21 May 2025 - Added array_name_size value at top of output file
+Version 10: 21 May 2025 - Added array_name_size value at top of output file
+Version 11: 22 May 2025 - Fixed filter artifacts at image edges, added sharpening filter
 -->
-  <title>Image to OpenSCAD array, v10</title><!-- REMEMBER TO CHANGE VERSION -->
+  <title>Image to OpenSCAD array, v11</title><!-- REMEMBER TO CHANGE VERSION -->
   <meta charset="UTF-8">
   <style>
     body { font-family: sans-serif; padding-left:1em; padding-right:1em;}
@@ -223,13 +224,15 @@ Alpha channel is ignored. After processing the image as desired, you may save it
       <div style="margin-top:8px;">
         <label><input type="checkbox" id="invertBrightness"> Invert brightness</label>
       </div>
-      <div style="margin:8px 0;">
+      <fieldset style="margin:8px 0;">
+	  <legend style="font-size:medium;">Filter cascade</legend>
         <label for="blurRadius">Gaussian blur radius (pixels):</label>
         <input type="number" id="blurRadius" size="5" min="0" max="20" value="0"><br>
-        <label for="sobelRadius" class="tooltip">Edge detect radius (pixels):
-        <span class="tooltiptext">Sobel filter uses own radius if Gaussian blur=0</span></label>
+        <label for="sharpenRadius">Sharpen radius (pixels):
+        <input type="number" id="sharpenRadius" size="5" min="0" max="20" value="0"><br>
+        <label for="sobelRadius">Edge detect radius (pixels):
         <input type="number" id="sobelRadius" size="5" min="0" max="20" value="0">
-      </div>
+      </fieldset>
 
       <div class="slider-row">
         <label for="contrast" class="slider-label tooltip">Contrast
@@ -293,6 +296,7 @@ Alpha channel is ignored. After processing the image as desired, you may save it
     const cropBottom = document.getElementById('cropBottom');
     const blurRadiusInput = document.getElementById('blurRadius');
     const sobelRadiusInput = document.getElementById('sobelRadius');
+    const sharpenRadiusInput = document.getElementById('sharpenRadius');
     const contrastInput = document.getElementById('contrast');
     const contrastValue = document.getElementById('contrastValue');
     const thresholdInput = document.getElementById('threshold');
@@ -355,47 +359,47 @@ Alpha channel is ignored. After processing the image as desired, you may save it
       return kernel.map(v => v / norm);
     }
 
-    function convolve1DHorizontal(matrix, kernel, normalize=true) {
+    function convolve1DHorizontal(matrix, kernel) {
       const width = matrix[0].length;
       const height = matrix.length;
       const r = Math.floor(kernel.length / 2);
       const result = [];
+	  let indx, nx;
       for (let y = 0; y < height; y++) {
         result[y] = [];
         for (let x = 0; x < width; x++) {
           let sum = 0;
-          let weightSum = 0;
           for (let k = -r; k <= r; k++) {
-            const nx = x + k;
+			indx = x+k;
+            nx = indx<0 ? -indx : (indx>=width ? 2*(width-1)-indx : indx); // reflect edges
             if (nx >= 0 && nx < width) {
                 sum += matrix[y][nx] * kernel[k+r];
-                weightSum += kernel[k+r];
             }
           }
-          result[y][x] = normalize ? (weightSum !== 0 ? sum / weightSum : 0) : sum;
+          result[y][x] = sum;
         }
       }
       return result;
     }
 
-    function convolve1DVertical(matrix, kernel, normalize=true) {
+    function convolve1DVertical(matrix, kernel) {
       const width = matrix[0].length;
       const height = matrix.length;
       const r = Math.floor(kernel.length / 2);
       const result = [];
+	  let indx, ny;
       for (let y = 0; y < height; y++) {
         result[y] = [];
         for (let x = 0; x < width; x++) {
           let sum = 0;
-          let weightSum = 0;
           for (let k = -r; k <= r; k++) {
-            const ny = y + k;
+		    indx = y+k;
+            ny = indx<0 ? -indx : (indx >= height ? 2*(height-1)-indx : indx); // reflect edges
             if (ny >= 0 && ny < height) {
                 sum += matrix[ny][x] * kernel[k+r];
-                weightSum += kernel[k+r];
             }
           }
-          result[y][x] = normalize ? (weightSum !== 0 ? sum / weightSum : 0) : sum;
+          result[y][x] = sum;
         }
       }
       return result;
@@ -416,18 +420,34 @@ Alpha channel is ignored. After processing the image as desired, you may save it
     }
 
     function applyGaussianBlur(matrix, blurRadius) {
+	  if (blurRadius <= 0) return matrix;
       const gKernel = gaussianKernel1D(blurRadius)
       g1 = convolve1DVertical(matrix, gKernel);
       return convolve1DHorizontal(g1, gKernel);
     }
 
-    function applySobel(matrix, sobelRadius, blurRadius) {
+    function applySharpen(original, blurMatrix, radius, blurRadius, k) {
+	  if (radius <= 0) return blurMatrix;
+      const height = original.length;
+      const width = original[0].length;
+	  blurred = blurRadius === 0 ? applyGaussianBlur(original, radius) : blurMatrix;
+      const result = [];
+      for (let y = 0; y < height; y++) {
+        result[y] = [];
+        for (let x = 0; x < width; x++) {
+          result[y][x] = original[y][x] + k * (original[y][x] - blurred[y][x]);
+        }
+      }
+      return result;
+	}
+
+    function applySobel(matrix, sobelRadius, sharpenRadius, blurRadius) {
       if (sobelRadius <= 0) return matrix; // No edge detection
       const sobelSize = 2 * sobelRadius + 1;
       const dKernel = sobelDerivativeKernel(sobelSize);
-      let gblur = blurRadius === 0 ? applyGaussianBlur(matrix, sobelRadius) : matrix;
-      gx = convolve1DHorizontal(gblur, dKernel, false);
-      gy = convolve1DVertical(gblur, dKernel, false);
+      let gblur = blurRadius === 0 && sharpenRadius === 0 ? applyGaussianBlur(matrix, sobelRadius) : matrix;
+      gx = convolve1DHorizontal(gblur, dKernel);
+      gy = convolve1DVertical(gblur, dKernel);
       return computeEdgeMagnitude(gx, gy);
     }
 
@@ -492,13 +512,13 @@ Alpha channel is ignored. After processing the image as desired, you may save it
         brightnessMatrix.push(row);
       }
 
-      // apply blurring to the grayscale image
+      // apply filter cascade
       const blurRadius = parseInt(blurRadiusInput.value) || 0;
-      const blurredMatrix = applyGaussianBlur(brightnessMatrix, blurRadius);
-
-      // apply Sobel edge detection
+	  const sharpenRadius = parseInt(sharpenRadiusInput.value) || 0;
       const sobelRadius = parseInt(sobelRadiusInput.value) || 0;
-      const sobelMatrix = applySobel(blurredMatrix, sobelRadius, blurRadius);
+      const blurMatrix = applyGaussianBlur(brightnessMatrix, blurRadius);
+	  let filteredMatrix = applySharpen(brightnessMatrix, blurMatrix, sharpenRadius, blurRadius, 1.0);
+      filteredMatrix = applySobel(filteredMatrix, sobelRadius, sharpenRadius, blurRadius);
 
       // crop the matrix, gather min and max values in crop area
       const cropMatrix = [];
@@ -506,14 +526,14 @@ Alpha channel is ignored. After processing the image as desired, you may save it
       let cropx2 = parseInt(cropID[edgeID[0]].value) || 0;
       let cropy1 = parseInt(cropID[edgeID[1]].value) || 0;
       let cropy2 = parseInt(cropID[edgeID[3]].value) || 0;
-      let min = 255;
-      let max = 0;
+      let min = 32000;
+      let max = -32000;
       for (let y=cropy1; y<uncropDim.height-cropy2; y++) {
         const row = [];
         for(let x=cropx1; x<uncropDim.width-cropx2; x++) {
-          row.push(sobelMatrix[y][x]);
-          min = Math.min(min, sobelMatrix[y][x]);
-          max = Math.max(max, sobelMatrix[y][x]);
+          row.push(filteredMatrix[y][x]);
+          min = Math.min(min, filteredMatrix[y][x]);
+          max = Math.max(max, filteredMatrix[y][x]);
         }
         cropMatrix.push(row);
       }
@@ -626,7 +646,7 @@ Alpha channel is ignored. After processing the image as desired, you may save it
 
     // set up event listeners for all the input gadgets
 
-    [blurRadiusInput, sobelRadiusInput, contrastInput, thresholdInput,
+    [blurRadiusInput, sobelRadiusInput, sharpenRadiusInput, contrastInput, thresholdInput,
         ...document.querySelectorAll('input[name="grayModel"]')].forEach(el => el.addEventListener('input', processImage));
 
     resizeWidthInput.addEventListener('input', function () {

From 4afc72cb22afbb64a9f9ee80c2d488dd553b39ac Mon Sep 17 00:00:00 2001
From: Alex Matulich <anachronist@gmail.com>
Date: Fri, 30 May 2025 16:38:01 -0700
Subject: [PATCH 3/4] Made filters mutually exclusive

---
 scripts/img2scad.html | 52 +++++++++++++++++++++++++++++++------------
 1 file changed, 38 insertions(+), 14 deletions(-)

diff --git a/scripts/img2scad.html b/scripts/img2scad.html
index 2487be5e..843baca0 100644
--- a/scripts/img2scad.html
+++ b/scripts/img2scad.html
@@ -10,8 +10,9 @@ Version 8: 26 April 2025 - added file size estimate to output section
 Version 9: 20 May 2025 - improved appearance UI, added Sobel edge detection
 Version 10: 21 May 2025 - Added array_name_size value at top of output file
 Version 11: 22 May 2025 - Fixed filter artifacts at image edges, added sharpening filter
+Version 12: 30 May 2025 - Made filters mutually exclusive
 -->
-  <title>Image to OpenSCAD array, v11</title><!-- REMEMBER TO CHANGE VERSION -->
+  <title>Image to OpenSCAD array, v12</title><!-- REMEMBER TO CHANGE VERSION -->
   <meta charset="UTF-8">
   <style>
     body { font-family: sans-serif; padding-left:1em; padding-right:1em;}
@@ -225,11 +226,14 @@ Alpha channel is ignored. After processing the image as desired, you may save it
         <label><input type="checkbox" id="invertBrightness"> Invert brightness</label>
       </div>
       <fieldset style="margin:8px 0;">
-	  <legend style="font-size:medium;">Filter cascade</legend>
+	  <legend style="font-size:medium;">Filter</legend>
+	    <input type="radio" name="filterSelect" value="blur" checked>
         <label for="blurRadius">Gaussian blur radius (pixels):</label>
         <input type="number" id="blurRadius" size="5" min="0" max="20" value="0"><br>
+		<input type="radio" name="filterSelect" value="sharp">
         <label for="sharpenRadius">Sharpen radius (pixels):
         <input type="number" id="sharpenRadius" size="5" min="0" max="20" value="0"><br>
+		<input type="radio" name="filterSelect" value="edge">
         <label for="sobelRadius">Edge detect radius (pixels):
         <input type="number" id="sobelRadius" size="5" min="0" max="20" value="0">
       </fieldset>
@@ -294,6 +298,7 @@ Alpha channel is ignored. After processing the image as desired, you may save it
     const cropLeft = document.getElementById('cropLeft');
     const cropRight = document.getElementById('cropRight');
     const cropBottom = document.getElementById('cropBottom');
+	const filterSelect = document.getElementById('filterSelect');
     const blurRadiusInput = document.getElementById('blurRadius');
     const sobelRadiusInput = document.getElementById('sobelRadius');
     const sharpenRadiusInput = document.getElementById('sharpenRadius');
@@ -426,11 +431,11 @@ Alpha channel is ignored. After processing the image as desired, you may save it
       return convolve1DHorizontal(g1, gKernel);
     }
 
-    function applySharpen(original, blurMatrix, radius, blurRadius, k) {
-	  if (radius <= 0) return blurMatrix;
+    function applySharpen(original, radius, k=1.0) {
+	  if (radius <= 0) return original;
       const height = original.length;
       const width = original[0].length;
-	  blurred = blurRadius === 0 ? applyGaussianBlur(original, radius) : blurMatrix;
+	  blurred = applyGaussianBlur(original, radius);
       const result = [];
       for (let y = 0; y < height; y++) {
         result[y] = [];
@@ -441,11 +446,11 @@ Alpha channel is ignored. After processing the image as desired, you may save it
       return result;
 	}
 
-    function applySobel(matrix, sobelRadius, sharpenRadius, blurRadius) {
-      if (sobelRadius <= 0) return matrix; // No edge detection
-      const sobelSize = 2 * sobelRadius + 1;
+    function applySobel(matrix, radius) {
+      if (radius <= 0) return matrix; // No edge detection
+      const sobelSize = 2 * radius + 1;
       const dKernel = sobelDerivativeKernel(sobelSize);
-      let gblur = blurRadius === 0 && sharpenRadius === 0 ? applyGaussianBlur(matrix, sobelRadius) : matrix;
+      let gblur = applyGaussianBlur(matrix, radius);
       gx = convolve1DHorizontal(gblur, dKernel);
       gy = convolve1DVertical(gblur, dKernel);
       return computeEdgeMagnitude(gx, gy);
@@ -512,13 +517,29 @@ Alpha channel is ignored. After processing the image as desired, you may save it
         brightnessMatrix.push(row);
       }
 
-      // apply filter cascade
+      // apply filter
       const blurRadius = parseInt(blurRadiusInput.value) || 0;
 	  const sharpenRadius = parseInt(sharpenRadiusInput.value) || 0;
       const sobelRadius = parseInt(sobelRadiusInput.value) || 0;
-      const blurMatrix = applyGaussianBlur(brightnessMatrix, blurRadius);
-	  let filteredMatrix = applySharpen(brightnessMatrix, blurMatrix, sharpenRadius, blurRadius, 1.0);
-      filteredMatrix = applySobel(filteredMatrix, sobelRadius, sharpenRadius, blurRadius);
+	  let filteredMatrix = [];
+      switch(document.querySelector('input[name="filterSelect"]:checked').value) {
+	  // any of the filters return the original if the radius=0
+	    case "blur":
+		console.log("blur");
+		  filteredMatrix = applyGaussianBlur(brightnessMatrix, blurRadius);
+		  break;
+		case "sharp":
+		console.log("sharp");
+	      filteredMatrix = applySharpen(brightnessMatrix, sharpenRadius);
+		  break;
+		case "edge":
+		console.log("edge");
+		  filteredMatrix = applySobel(brightnessMatrix, sobelRadius);
+		  break;
+		default:
+		console.log("none");
+		  filteredMatrix = brightnessMatrix;
+      }
 
       // crop the matrix, gather min and max values in crop area
       const cropMatrix = [];
@@ -647,7 +668,10 @@ Alpha channel is ignored. After processing the image as desired, you may save it
     // set up event listeners for all the input gadgets
 
     [blurRadiusInput, sobelRadiusInput, sharpenRadiusInput, contrastInput, thresholdInput,
-        ...document.querySelectorAll('input[name="grayModel"]')].forEach(el => el.addEventListener('input', processImage));
+        ...document.querySelectorAll('input[name="grayModel"]'),
+        ...document.querySelectorAll('input[name="filterSelect"]')
+		].forEach(el => el.addEventListener('input', processImage)
+		);
 
     resizeWidthInput.addEventListener('input', function () {
         let min = parseInt(this.min);

From cb9c933e8ccc165f547e736977401b12ded4c004 Mon Sep 17 00:00:00 2001
From: Adrian Mariano <avm4@cornell.edu>
Date: Fri, 13 Jun 2025 17:43:36 -0400
Subject: [PATCH 4/4] fix workflows

---
 .github/workflows/gen_docs.yml      | 2 +-
 .github/workflows/gen_tutorials.yml | 2 +-
 .github/workflows/main.yml          | 2 +-
 3 files changed, 3 insertions(+), 3 deletions(-)

diff --git a/.github/workflows/gen_docs.yml b/.github/workflows/gen_docs.yml
index 13324589..f46d6fa1 100644
--- a/.github/workflows/gen_docs.yml
+++ b/.github/workflows/gen_docs.yml
@@ -21,7 +21,7 @@ jobs:
       run: sudo apt-get install python3-pip python3-dev python3-setuptools python3-pil gifsicle libfuse2
 
     - name: Install OpenSCAD-DocsGen package.
-      run: sudo pip3 install openscad-docsgen
+      run: sudo pip3 install openscad-docsgen imageio
 
     - name: Install OpenSCAD
       run: |
diff --git a/.github/workflows/gen_tutorials.yml b/.github/workflows/gen_tutorials.yml
index 0bc8e61e..f67684e0 100644
--- a/.github/workflows/gen_tutorials.yml
+++ b/.github/workflows/gen_tutorials.yml
@@ -21,7 +21,7 @@ jobs:
       run: sudo apt-get install python3-pip python3-dev python3-setuptools python3-pil gifsicle libfuse2
 
     - name: Install OpenSCAD-DocsGen package.
-      run: sudo pip3 install openscad-docsgen
+      run: sudo pip3 install openscad-docsgen imageio
 
     - name: Install OpenSCAD
       run: |
diff --git a/.github/workflows/main.yml b/.github/workflows/main.yml
index 8a5cab9f..803adcaa 100644
--- a/.github/workflows/main.yml
+++ b/.github/workflows/main.yml
@@ -43,7 +43,7 @@ jobs:
       run: sudo apt-get install python3-pip python3-dev python3-setuptools python3-pil libfuse2
 
     - name: Install OpenSCAD-DocsGen package.
-      run: sudo pip3 install openscad-docsgen
+      run: sudo pip3 install openscad-docsgen imageio
 
     - name: Install OpenSCAD
       run: |