The input-image is a gray-scale image of soil that has some thin fractures. These fractures are detected by the application of the top-hat and area open operators.
The image of fractures in soil is read.
a = mmreadgray('soil.tif');
The fracture lines are enhanced by the close top-hat operator.
b = mmcloseth(a,mmsebox(2));
Small connected bright regions are removed by the gray-scale area open operator. Note the connectivity used (mmsebox : 8-connected).
c = mmareaopen(b,80,mmsebox());
The fracture lines are detected. This threshold is very robust.
d = mmthreshad(c,15);
Warning: upcasting image from scalar to uint8
Overlay of the fracture lines over the original image.
|Copyright (c) 1998-2008 by SDC Information Systems|