Pruning (morphology)

The pruning algorithm is a technique used in digital image processing based on mathematical morphology. It is used as a complement to the skeleton and thinning algorithms to remove unwanted parasitic components (spurs). In this case 'parasitic' components refer to branches of a line which are not key to the overall shape of the line and should be removed. These components can often be created by edge detection algorithms or digitization. Common uses for pruning include automatic recognition of hand-printed characters. Often inconsistency in letter writing creates unwanted spurs that need to be eliminated for better characterization.

Mathematical Definition

The standard pruning algorithm will remove all branches shorter than a given number of points. If a parasitic branch is shorter than four points and we run the algorithm with ''n = 4'' the branch will be removed. The second step ensures that the main trunks of each line are not shortened by the procedure.

Structuring Elements

$B_1 = \begin x & 0 & 0 \\ 1 & 1 & 0 \\ x & 0 & 0 \end B_2 = \begin x & 1 & x \\ 0 & 1 & 0 \\ 0 & 0 & 0 \end B_3 = \begin 0 & 0 & x \\ 0 & 1 & 1 \\ 0 & 0 & x \end B_4 = \begin 0 & 0 & 0 \\ 0 & 1 & 0 \\ x & 1 & x \end$

$B_5 = \begin 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 0 \end B_6 = \begin 0 & 0 & 1 \\ 0 & 1 & 0 \\ 0 & 0 & 0 \end B_7 = \begin 0 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end B_8 = \begin 0 & 0 & 0 \\ 0 & 1 & 0 \\ 1 & 0 & 0 \end$

The x in the arrays indicates a “don’t care” condition i.e. the image could have either a 1 or a 0 in the spot.

Step 1:

Thinning

Apply this step a given (n) times to eliminate any branch with (n) or less pixels.

$X_1 = A \otimes \$

Step 2: Find End Points

Wherever the structuring elements are satisfied, the center of the 3x3 matrix is considered an endpoint.

$X_2 = \bigcup_^8 (X_1\circledast B^k)$

Step 3:

Dilate End Points

Perform dilation using a 3x3 matrix (H) consisting of all 1's and only insert 1's where the original image (A) also had a 1. Perform this for each endpoint in all direction (n) times.

$X_3 = (X_2 \oplus H) \cap A$

Step 4:

Union of X₁ & X₃

Take the result from step 1 and union it with step 3 to achieve the final results.

$X_4 = X_1 \cup X_3$

MATLAB Code

%% ---------------
% Pruning
% ---------------
clear; clc;

% Image read in
img = imread('Pruning.tif');

b_img_skel = bwmorph (img, 'skel', 40);
b_img_spur = bwmorph(b_img_skel, 'spur', Inf);

figure('Name', 'Pruning');
subplot(1,2,1);
imshow(b_img_skel);
title(sprintf('Image Skeleton'));
subplot(1,2,2);
imshow(b_img_spur);
title(sprintf('Skeleton Image Pruned'));

MATLAB Example

In the MATLAB example below, it takes the original image (below left) and skeletonize it 40 times then prunes the image to remove the spurs as per the MATLAB code above. As shown (below right) this removed the majority of all spurs resulting in a cleaner image.

See also

* Morphological image processing

External links

Morphological Pruning function
in Mathematica

Morphological Operations
in MATLAB

References

{{DEFAULTSORT:Pruning (Morphology)
Digital geometry
Mathematical morphology