Fix polygonize with LinRange

JuliaGeo · Jul 30, 2024 · 45c6780 · 45c6780
1 parent ca428c6
commit 45c6780
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/methods/polygonize.jl b/src/methods/polygonize.jl
@@ -121,8 +121,8 @@ function _polygonize(f::Base.Callable, xs::AbstractRange, ys::AbstractRange, A::
     ybounds = first(ys) - yhalf : step(ys) : last(ys) + yhalf
     Tx = eltype(xbounds)
     Ty = eltype(ybounds)
-    xvec = similar(Vector{Tuple{Tx,Tx}}, xs)
-    yvec = similar(Vector{Tuple{Ty,Ty}}, ys)
+    xvec = similar(Vector{Tuple{Tx,Tx}}, axes(xs))
+    yvec = similar(Vector{Tuple{Ty,Ty}}, axes(ys))
     for (xind, i) in enumerate(eachindex(xvec))
         xvec[i] = xbounds[xind], xbounds[xind+1]
     end