--- a/src/Translate.jl Wed Dec 25 17:42:32 2019 +0200 +++ b/src/Translate.jl Fri Dec 27 20:59:44 2019 +0200 @@ -132,18 +132,18 @@ vyʹ = (imy-by)/2 - v[2] # Target image indices within source image - px = max(ceil(Int, max(1, vxʹ + 1) - vxʹ), 1) - py = max(ceil(Int, max(1, vyʹ + 1) - vyʹ), 1) - qx = min(floor(Int, min(imx, vxʹ + bx) - vxʹ), bx) - qy = min(floor(Int, min(imy, vyʹ + by) - vyʹ), by) + px = min(max(ceil(Int, max(1, vxʹ + 1) - vxʹ), 1), bx) + py = min(max(ceil(Int, max(1, vyʹ + 1) - vyʹ), 1), by) + qx = max(min(floor(Int, min(imx, vxʹ + bx) - vxʹ), bx), 1) + qy = max(min(floor(Int, min(imy, vyʹ + by) - vyʹ), by), 1) @inbounds begin b[1:px-1, :] .= 0 - b[qx+1:end, :] .= 0 + b[qx+1:bx, :] .= 0 Threads.@threads for i=px:qx b[i, 1:py-1] .= 0 - b[i, qy+1:end] .= 0 + b[i, qy+1:by] .= 0 @simd for j=py:qy b[i, j] = interpolate2d_quadrants(im, (i+vxʹ, j+vyʹ)) end