fix in gaussian, now consistent with python version

jonocarroll · Feb 3, 2024 · f51c97f · f51c97f
1 parent d65549b
commit f51c97f
Show file tree

Hide file tree

Showing 8 changed files with 5,185 additions and 6,557 deletions.
diff --git a/examples/tri_bird_fast.svg b/examples/tri_bird_fast.svg
diff --git a/examples/tri_bird_small.svg b/examples/tri_bird_small.svg
diff --git a/examples/tri_bird_small_100_refine.svg b/examples/tri_bird_small_100_refine.svg
diff --git a/examples/tri_bird_small_500.svg b/examples/tri_bird_small_500.svg
diff --git a/examples/tri_bird_small_debug.svg b/examples/tri_bird_small_debug.svg
diff --git a/examples/tri_bird_small_debug_noimage.svg b/examples/tri_bird_small_debug_noimage.svg
diff --git a/examples/tri_everest.svg b/examples/tri_everest.svg
diff --git a/src/TriangulArt.jl b/src/TriangulArt.jl
@@ -60,20 +60,20 @@ function generate_weighted_points(image::Matrix{RGB{N0f8}}, n_points::Integer=10
     radius = floor(Int, length_scale)
 
     amp = 3
-    entropy_width = 20
+    entropy_width = 0.2 * length_scale
     im2 = skimage.filters.rank.entropy(im1, skimage.morphology.disk(entropy_width))
 
     epts = Array{Int}(undef, (n_points, 2))
     for n in 1:n_points
         maxpos = last(findmax(im2))
         x, y = getindex(maxpos, 1), getindex(maxpos, 2)
         epts[n, :] = [y x]
-        for i in (x-radius):(x+radius)
-            for j in (y-radius):(y+radius)
+        for i in (x-2radius):(x+2radius)
+            for j in (y-2radius):(y+2radius)
                 if !checkbounds(Bool, im2, i, j)
                     continue
                 end
-                im2[i, j] = im2[i, j] - amp * exp(-sqrt((i - x)^2 + (j - y)^2) / 2sigma^2)
+                im2[i, j] = im2[i, j] - amp * exp(-((i - x)^2 + (j - y)^2) / 2sigma^2)
             end
         end
     end