Pruning (morphology) explained

The pruning algorithm is a technique used in digital image processing based on mathematical morphology.[1] 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.[2]

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

B1=\begin{bmatrix}x&0&0\ 1&1&0\x&0&0\end{bmatrix} B2=\begin{bmatrix}x&1&x\ 0&1&0\ 0&0&0\end{bmatrix} B3=\begin{bmatrix}0&0&x\ 0&1&1\ 0&0&x\end{bmatrix} B4=\begin{bmatrix}0&0&0\ 0&1&0\x&1&x\end{bmatrix}

B5=\begin{bmatrix}1&0&0\ 0&1&0\ 0&0&0\end{bmatrix} B6=\begin{bmatrix}0&0&1\ 0&1&0\ 0&0&0\end{bmatrix} B7=\begin{bmatrix}0&0&0\ 0&1&0\ 0&0&1\end{bmatrix} B8=\begin{bmatrix}0&0&0\ 0&1&0\ 1&0&0\end{bmatrix}

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.

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

X1=A\{B\}

Step 2: Find End Points

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

X2=

8
cup
k=1

(X1\circledastBk)

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.

X3=(X2H)\capA

Step 4: Union of X1 & X3

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

X4=X1\cupX3

MATLAB Code

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

% Image read inimg = 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.

width=250px Original Image width=250px Image Skeletonwidth=250px Skeleton Image Pruned
align=center

See also

External links

Notes and References

  1. Book: Russ. John C.. The image processing handbook. 2011. CRC Press. Boca Raton. 978-1-4398-4045-0. 6th.
  2. Book: Digital image processing. Gonzalez. Rafael C.. Woods. Richard E.. 2008. Prentice Hall. 978-0131687288. 3rd. Upper Saddle River, N.J..