mmdsoil
Detect fractures in soil.

The image of fractures in soil is read.

a = mmreadgray('soil.tif');

mmshow(a);

a

The fracture lines are enhanced by the close top-hat operator.

b = mmcloseth(a,mmsebox(2));

mmshow(b);

b

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);

mmshow(c);

c

The fracture lines are detected. This threshold is very robust.

d = mmthreshad(c,15);

Warning: converting image from scalar to uint8

mmshow(d);

d

Overlay of the fracture lines over the original image.

mmshow(a,d);

a,d