Deriche edge detector

Differences between Canny and Deriche edge detector

Deriche edge detector, like Canny edge detector, consists of the following 4 steps:

1. Smoothing
2. Calculation of magnitude and gradient direction
3. Non-maximum suppression
4. Hysteresis thresholding (using two thresholds)

The essential difference is in the implementation of the first two steps of the algorithm. Unlike the Canny edge detector, Deriche edge detector uses the IIR filter in the form:

: f(x)=\frac{S}{\omega}e^{-\alpha|x|}\sin\omega x

The filter optimizes the Canny criteria. As is evident from the preceding formula, the most effective filter is obtained when the value of \omega approaches 0. Such filter then uses the formula:

: f(x)=Sxe^{-\alpha|x|}

The advantage of such a filter is that it can be adapted to the characteristics of the processed image using only one parameter. If the value of α is small (usually between 0.25 and 0.5), it results in better detection. On the other hand, better localization is achieved when the parameter has a higher value (around 2 or 3). For most of the normal cases parameter value of around 1 is recommended.

Using the IIR filter makes sense especially in cases where the processed image is noisy or a large amount of smoothing is required (which leads to large convolution kernel for FIR filter). In these cases, the Deriche detector has considerable advantage over the Canny detector, because it is able to process images in a short constant time independent of the desired amount of smoothing.

Deriche detector implementation

It is possible to separate the process of obtaining the value of a two-dimensional Deriche filter into two parts. In first part, image array is passed in the horizontal direction from left to right according to the following formula:

: y_{ij}^1=a_1x_{ij}+a_2x_{ij-1}+b_1y_{ij-1}^1+b_2y_{ij-2}^1

and from right to left according to the formula:

: y_{ij}^2=a_3x_{ij+1}+a_4x_{ij+2}+b_1y_{ij+1}^2+b_2y_{ij+2}^2

The result of the computation is then stored into temporary two-dimensional array:

: \theta_{ij}=c_1(y_{ij}^1+y_{ij}^2)

The second step of the algorithm is very similar to the first one. The two-dimensional array from the previous step is used as the input. It is then passed in the vertical direction from top to bottom and bottom-up according to the following formulas:

: y_{ij}^1=a_5\theta_{ij}+a_6\theta_{i-1j}+b_1y_{i-1j}^1+b_2y_{i-2j}^1

: y_{ij}^2=a_7\theta_{i+1j}+a_8\theta_{i+2j}+b_1y_{i+1j}^2+b_2y_{i+2j}^2

: \Theta_{ij}=c_2(y_{ij}^1+y_{ij}^2)

The description of the algorithm implies that the processed rows and columns are independent of each other. As a result, the solution based on the IIR filter is often employed in embedded systems and architectures which support a high level of parallelization.

The mathematical properties of the algorithm are often used in practical implementation of the Deriche detector. It is sufficient to implement only one part of the algorithm, which is then called twice, while performing a transposition of the resulting matrix.