Пример 10.11. Сегментация полутонового изображения с помощью градиента и преобразования по водоразделу.

f = imread('C: \Users\Alex\Desktop\small-blobs.tif');

h = fspecial('sobel');

fd = double(f);

g = sqrt(imfilter(fd, h, 'replicate').^2 +...

imfilter(fd, h', 'replicate').^2);

L = watershed(g);

wr = L == 0;

g2 = imclose(imopen(g, ones(3, 3)), ones(3, 3));

L2 = watershed(g2);

wr2 = L2 == 0;

f2 = f;

f2(wr2) = 255;

imshow(f2)

Пример 10.12. Иллюстрация использования маркеров при сегментации по водоразделам.

f = imread('C: \Users\Alex\Desktop\gel.tif');

h = fspecial('sobel');

fd = double(f);

g = sqrt(imfilter(fd, h, 'replicate').^2 +...

imfilter(fd, h', 'replicate').^2);

L = watershed(g);

wr = L == 0;

rm = imregionalmin(g);

im = imextendedmin(f, 2);

fim = f;

fim(im) = 175;

Lim = watershed(bwdist(im));

em = Lim == 0;

g2 = imimposemin(g, im | em);

L2 = watershed(g2);

f2 = f;

f2(L2 == 0) = 255;

imshow(f2)

Пример 11.1. Простые иллюстрации работы со смешанным массивом.

f = imread('C: \Users\Alex\Desktop\gel.tif');

G = image_stats(f);

H = image_stats2(f);

G

G =

[1x2 double] [129.9775] [256x1 double] [1x256 double]

H

[1x2 double] [129.9775]

[256x1 double] [1x256 double]

v = G{1}

v =

256 256

v = H{1, 1}

v =

256 256