# Haskell image processing - easy and automatically parallel

@ 2011-06-20 by João Paulo Pizani Flor

One of the topics we learn in the Pattern Recognition course here at UFSC is pattern recognition devired from image processing. Our final course project is to implement a selective intrusion detector that uses only a webcam… This is pretty complicated, and involves several processing steps, much more dificult tham it might seem at a first glance :)

One of the most common processing steps performed to identify objects in an image is edge detection. Our first course assignment is to implement the Roberts Cross edge detection algorithm. Because of my wish to get up-to-speed and review my Haskell knowledge, and also because this sort of algorithm fits perfectly the functional paradigm, I decided to give Haskell a shot! I’ll then show you the code, the pretty images and report my (almost painless) experience implementing it :P

One of the strongest selling points of functional programming is that code written according to the paradigm can be easily parallelized. Roberts edge detection consists, for example, of nothing more than scanning an image, pixel by pixel, and applying a small matrix of coefficients that calculates the weighted average of brightness (it’s called Convolution by the guys at Wikipedia). Well, this sort of algorithm is INHERENTLY parallel, and we can divide the task among as many threads as we want, without changing the semantics of the code. This is exactly what the Repa library does: it allows you to easily express several algorithms over multidimensional arrays and have them parallelized automatically.

Without any more explanation, here goes an example of executing the “Roberts” program on 1 thread, for a 1024x768 image:

Original image:

Edges:

Execution time (some of it is IO-bound):

real	0m7.112s
user	0m7.044s
sys	0m0.020s

Now executing it on 2 and 4 threads, the only difference is that you need to run the program with some arguments on the command line:

Roberts wow.png wow_bordas.bmp +RTS -N{2,4}

real	0m4.000s
user	0m7.176s
sys	0m0.076s

real	0m3.096s
sys	0m0.456s