### strel class

The

`strel`

class was something that I wanted to implement a long time ago but had never found the time. The name strel comes from *str*ucturing

*el*ement (SE), the shapes used in morphological operations such as dilation and erosion. I have only seen it as a standard way to create SEs, but is actually much more. Specially, SE decomposition can have a really nice increase in performance.

Roberto Metere submitted his own implementation of the class last month and we have been working on it, slowly adding it to the other functions of the package. It started as a single .m function, no OOP at all, but he managed to implement

`@strel`

with the old `@class`

style while keeping matlab compatibility. All the basic methods have been implemented: the object constructor, `getnhood`

, `getheight`

and `reflect`

.The idea behind SE decomposition is that morphology operations take as much space as the number of pixels in a SE. The bigger the SE, the slower it will be. However, some SE can be replaced by a sequence of smaller ones so it's in our best interest to use them. For an example on performance, see how the use of a square compares to use of 1 row and 1 column of the same size:

octave>im1 = im2 = randp (5, 2000) > 15;octave>t = cputime ();octave>im1 = conv2 (im1, true (20), "same") > 0; # dilation by 1 squareoctave>cputime () - t ans = 2.6402octave>t = cputime ();octave>im2 = conv2 (im2, true (20, 1), "same") > 0; # dilation by 1 columnoctave>im2 = conv2 (im2, true (1, 20), "same") > 0; # dilation by 1 rowoctave>cputime () - t ans = 0.52803octave>isequal (im1, im2) ans = 1

At the moment, decomposition is only being done for rectangular, square and cube shapes but other will come with time. Functions that can gain from SE decomposition, are written so that it does not matter if it has been implemented for a specific shape. This means that when it is implemented for another shape, its effect will be immediate across the whole image package. The only file where this is done is

`inst/@strel/getsequence.m`

so send us patches.Going through

`imdilate()`

and `imerode()`

to make them use `strel`

, brought up a bunch of other matlab incompatibilities that I hope are now fixed, as well as other improvements. I'm a bit interested in morphology of volumes so one of the changes made was making them work for N-dimensional images (think MRI scans).On top of the matlab shapes for

`strel`

, I implemented the cube as an optional shape. I also wanted to implement ball as a volume but unfortunately, matlab has already done it wrong as a non-flat ellipsoid. Note that non-flat is unrelated to volumes.### Non-flat morphology

This has confused me for a very long time. Because of the name (non-flat), and because 3D images are the norm for me, I have always assumed that a non-flat SE was one used for volumes. The fact that the only non-flat standard shape in matlat is named ball, which immediately brings up the idea of volume, also helped to the confusion.

Actually, non-flat morphology is something that only makes sense for grayscale operations. A non-flat SE is defined by two different matrices, one defining the neighbourhood (same as a flat SE) and another defining the height of each neighbour. These heights are added to the image pixels values before the erosion and dilation, in the same way as the variable

`S`in

`ordfiltn`

.Basically, not useful for volumes at all but the image package can do this now. To create a non-flat SE, use the arbitrary shape of

`strel`

.### Spatial transformations

Pantxo Diribarne has also submitted a set of functions for spatial transformations:

`maketform`

, `cp2tform`

, `tformfwd`

, and `tforminv`

. These are still not completely implemented and generally restricted to 2D transforms. Adding the missing options should now be much easier.
Hello,

ReplyDeleteI have installed this 2.1.0 version, but strel doesn'T exist and have got error with function "medfilt2":

error: padarray: invalid string in padval parameter.

error: called from:

error: /Users/p/octave/image-2.1.0/padarray.m at line 190, column 13

error: /Users/p/octave/image-2.1.0/private/pad_for_spatial_filter.m at line 28, column 8

error: /Users/p/octave/image-2.1.0/ordfiltn.m at line 78, column 5

error: /Users/p/octave/image-2.1.0/ordfilt2.m at line 36, column 5

error: /Users/p/octave/image-2.1.0/medfilt2.m at line 56, column 10

Regards,

Pai

Hi Pai

Deletestrel should be there. What do you mean does not exist? Does the code

strel ("square", 5)

works? Can you see the directory inst/@strel on the tarball?

Also, could you test medfilt2 again with the current revision (11849)? It should be already fixed.

Carnë

hi....can anyone give me the code for image dilation and erosion in octave?

ReplyDeleteDownload the latest version of the image package and you will find the source code for both those operations on the file named imerode.cc

DeleteInteresting functionality can be put to use in signal processing of 1-D signals too.

ReplyDelete