Edge detection using Sujoy Filter

using Images
using Statistics
using TestImages
using ImageBinarization

Sujoy's Edge Detection Algorithm is a better & more generic approach (first derivative) for edge detection than the other commonly used first-derivative methods (like Robert’s operator, Prewitt operator, Sobel operator etc.).

"""
edges = sujoy(img; four_connectivity=true)

Compute edges of an image using the Sujoy algorithm.

# Parameters

* img (Required): any gray image
* four_connectivity=true: if true, kernel is based on 4-neighborhood, else, kernel is based on
8-neighborhood,

# Returns

* edges : gray image
"""
function sujoy(img; four_connectivity=true)
img_channel = Gray.(img)

min_val = minimum(img_channel)
img_channel = img_channel .- min_val
max_val = maximum(img_channel)

if max_val == 0
return img
end

img_channel = img_channel./max_val

if four_connectivity
krnl_h = centered(Gray{Float32}[0 -1 -1 -1 0; 0 -1 -1 -1 0; 0 0 0 0 0; 0 1 1 1 0; 0 1 1 1 0]./12)
krnl_v = centered(Gray{Float32}[0 0 0 0 0; -1 -1 0 1 1;-1 -1 0 1 1;-1 -1 0 1 1;0 0 0 0 0 ]./12)
else
krnl_h = centered(Gray{Float32}[0 0 -1 0 0; 0 -1 -1 -1 0; 0 0 0 0 0; 0 1 1 1 0; 0 0 1 0 0]./8)
krnl_v = centered(Gray{Float32}[0 0 0 0 0;  0 -1 0 1 0; -1 -1 0 1 1;0 -1 0 1 0; 0 0 0 0 0 ]./8)
end

end

img = testimage("cameraman")
img_edge = sujoy(img, four_connectivity=true)
img_edge₀₁ = binarize(img_edge, Otsu()) # or use other binarization methods provided in ImageBinarization

mosaicview(img, img_edge, img_edge₀₁; nrow = 1)