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 input image from int32 to uint8.

>>> mmshow(d);

d

Overlay of the fracture lines over the original image.

>>> mmshow(a,d);

a,d