Erosion (morphology)

From HandWiki
Short description: Basic operation in mathematical morphology


The erosion of the dark-blue square by a disk, resulting in the light-blue square.

Erosion (usually represented by ) is one of two fundamental operations (the other being dilation) in morphological image processing from which all other morphological operations are based. It was originally defined for binary images, later being extended to grayscale images, and subsequently to complete lattices. The erosion operation usually uses a structuring element for probing and reducing the shapes contained in the input image.

Binary erosion

In binary morphology, an image is viewed as a subset of a Euclidean space [math]\displaystyle{ \mathbb{R}^d }[/math] or the integer grid [math]\displaystyle{ \mathbb{Z}^d }[/math], for some dimension d.

The basic idea in binary morphology is to probe an image with a simple, pre-defined shape, drawing conclusions on how this shape fits or misses the shapes in the image. This simple "probe" is called structuring element, and is itself a binary image (i.e., a subset of the space or grid).

Let E be a Euclidean space or an integer grid, and A a binary image in E. The erosion of the binary image A by the structuring element B is defined by:

[math]\displaystyle{ A \ominus B = \{z\in E \mid B_{z} \subseteq A\} }[/math],

where Bz is the translation of B by the vector z, i.e., [math]\displaystyle{ B_z = \{b+z \mid b\in B\} }[/math], [math]\displaystyle{ \forall z\in E }[/math].

When the structuring element B has a center (e.g., a disk or a square), and this center is located on the origin of E, then the erosion of A by B can be understood as the locus of points reached by the center of B when B moves inside A. For example, the erosion of a square of side 10, centered at the origin, by a disc of radius 2, also centered at the origin, is a square of side 6 centered at the origin.

The erosion of A by B is also given by the expression: [math]\displaystyle{ A \ominus B = \bigcap_{b\in B} A_{-b} }[/math], where A−b denotes the translation of A by -b.

This is more generally also known as a Minkowski difference.

Example

Suppose A is a 13 x 13 matrix and B is a 3 x 3 matrix:

    1 1 1 1 1 1 1 1 1 1 1 1 1        
    1 1 1 1 1 1 0 1 1 1 1 1 1    
    1 1 1 1 1 1 1 1 1 1 1 1 1    
    1 1 1 1 1 1 1 1 1 1 1 1 1   
    1 1 1 1 1 1 1 1 1 1 1 1 1                
    1 1 1 1 1 1 1 1 1 1 1 1 1               1 1 1
    1 1 1 1 1 1 1 1 1 1 1 1 1               1 1 1
    1 1 1 1 1 1 1 1 1 1 1 1 1               1 1 1
    1 1 1 1 1 1 1 1 1 1 1 1 1        
    1 1 1 1 1 1 1 1 1 1 1 1 1    
    1 1 1 1 1 1 1 1 1 1 1 1 1    
    1 1 1 1 1 1 1 1 1 1 1 1 1   
    1 1 1 1 1 1 1 1 1 1 1 1 1

Assuming that the origin B is at its center, for each pixel in A superimpose the origin of B, if B is completely contained by A the pixel is retained, else deleted.

Therefore the Erosion of A by B is given by this 13 x 13 matrix.

    0 0 0 0 0 0 0 0 0 0 0 0 0
    0 1 1 1 1 0 0 0 1 1 1 1 0
    0 1 1 1 1 0 0 0 1 1 1 1 0
    0 1 1 1 1 1 1 1 1 1 1 1 0
    0 1 1 1 1 1 1 1 1 1 1 1 0
    0 1 1 1 1 1 1 1 1 1 1 1 0 
    0 1 1 1 1 1 1 1 1 1 1 1 0
    0 1 1 1 1 1 1 1 1 1 1 1 0 
    0 1 1 1 1 1 1 1 1 1 1 1 0 
    0 1 1 1 1 1 1 1 1 1 1 1 0
    0 1 1 1 1 1 1 1 1 1 1 1 0
    0 1 1 1 1 1 1 1 1 1 1 1 0
    0 0 0 0 0 0 0 0 0 0 0 0 0

This means that only when B is completely contained inside A that the pixels values are retained, otherwise it gets deleted or eroded.

Properties

  • The erosion is translation invariant.
  • It is increasing, that is, if [math]\displaystyle{ A\subseteq C }[/math], then [math]\displaystyle{ A\ominus B \subseteq C\ominus B }[/math].
  • If the origin of E belongs to the structuring element B, then the erosion is anti-extensive, i.e., [math]\displaystyle{ A\ominus B\subseteq A }[/math].
  • The erosion satisfies [math]\displaystyle{ (A\ominus B)\ominus C = A\ominus (B\oplus C) }[/math], where [math]\displaystyle{ \oplus }[/math] denotes the morphological dilation.
  • The erosion is distributive over set intersection

Grayscale erosion

Example of erosion on a grayscale image using a 5x5 flat structuring element. The top figure demonstrates the application of the structuring element window to the individual pixels of the original image. The bottom figure shows the resulting eroded image.

In grayscale morphology, images are functions mapping a Euclidean space or grid E into [math]\displaystyle{ \mathbb{R}\cup\{\infty,-\infty\} }[/math], where [math]\displaystyle{ \mathbb{R} }[/math] is the set of reals, [math]\displaystyle{ \infty }[/math] is an element larger than any real number, and [math]\displaystyle{ -\infty }[/math] is an element smaller than any real number.

Denoting an image by f(x) and the grayscale structuring element by b(x), where B is the space that b(x) is defined, the grayscale erosion of f by b is given by

[math]\displaystyle{ (f\ominus b)(x)=\inf_{y\in B}[f(x+y)-b(y)] }[/math],

where "inf" denotes the infimum.

In other words the erosion of a point is the minimum of the points in its neighborhood, with that neighborhood defined by the structuring element. In this way it is similar to many other kinds of image filters like the median filter and the gaussian filter.

Erosions on complete lattices

Complete lattices are partially ordered sets, where every subset has an infimum and a supremum. In particular, it contains a least element and a greatest element (also denoted "universe").

Let [math]\displaystyle{ (L,\leq) }[/math] be a complete lattice, with infimum and supremum symbolized by [math]\displaystyle{ \wedge }[/math] and [math]\displaystyle{ \vee }[/math], respectively. Its universe and least element are symbolized by U and [math]\displaystyle{ \emptyset }[/math], respectively. Moreover, let [math]\displaystyle{ \{ X_{i} \} }[/math] be a collection of elements from L.

An erosion in [math]\displaystyle{ (L,\leq) }[/math] is any operator [math]\displaystyle{ \varepsilon: L\rightarrow L }[/math] that distributes over the infimum, and preserves the universe. I.e.:

  • [math]\displaystyle{ \bigwedge_{i}\varepsilon(X_i)=\varepsilon\left(\bigwedge_{i} X_i\right) }[/math],
  • [math]\displaystyle{ \varepsilon(U)=U }[/math].

See also

References

  • Image Analysis and Mathematical Morphology by Jean Serra, ISBN:0-12-637240-3 (1982)
  • Image Analysis and Mathematical Morphology, Volume 2: Theoretical Advances by Jean Serra, ISBN:0-12-637241-1 (1988)
  • An Introduction to Morphological Image Processing by Edward R. Dougherty, ISBN:0-8194-0845-X (1992)
  • Morphological Image Analysis; Principles and Applications by Pierre Soille, ISBN:3-540-65671-5 (1999)
  • R. C. Gonzalez and R. E. Woods, Digital image processing, 2nd ed. Upper Saddle River, N.J.: Prentice Hall, 2002.