Считывание изображения с объектом в сегмент. Преобразуйте изображение в градации серого и отобразите результат.

A = imread('hands1.jpg');
I = im2gray(A);
figure
imshow(I)
title('Original Image')

Используйте активные контуры (змеи) для сегментации руки.

mask = false(size(I));
mask(25:end-25,25:end-25) = true;
BW = activecontour(I, mask, 300);

Прочитайте в наземной сегментации правды.

BW_groundTruth = imread('hands1-mask.png');

Вычислите индекс Dice для сегментации активных контуров в зависимости от истинности грунта.

similarity = dice(BW, BW_groundTruth);

Отображение масок друг на друге. Цвета указывают на различия в масках.

figure
imshowpair(BW, BW_groundTruth)
title(['Dice Index = ' num2str(similarity)])

В этом примере показано, как сегментировать изображение на несколько областей. Затем в примере вычисляется коэффициент подобия Диса для каждой области.

Считывание изображения с несколькими областями для сегментации.

RGB = imread('yellowlily.jpg');

Создайте писанки для трех областей, отличающих их типичные цветовые характеристики. Первая область классифицирует жёлтый цветок. Вторая область классифицирует зелёный стебель и листья. Последний регион классифицирует коричневую грязь по двум отдельным пятнам изображения. Области задаются 4-элементным вектором, элементы которого обозначают координату x и y верхнего левого угла ROI, ширину ROI и высоту ROI.

region1 = [350 700 425 120]; % [x y w h] format
BW1 = false(size(RGB,1),size(RGB,2));
BW1(region1(2):region1(2)+region1(4),region1(1):region1(1)+region1(3)) = true;

region2 = [800 1124 120 230];
BW2 = false(size(RGB,1),size(RGB,2));
BW2(region2(2):region2(2)+region2(4),region2(1):region2(1)+region2(3)) = true;

region3 = [20 1320 480 200; 1010 290 180 240]; 
BW3 = false(size(RGB,1),size(RGB,2));    
BW3(region3(1,2):region3(1,2)+region3(1,4),region3(1,1):region3(1,1)+region3(1,3)) = true;
BW3(region3(2,2):region3(2,2)+region3(2,4),region3(2,1):region3(2,1)+region3(2,3)) = true;

Отображение начальных областей поверх изображения.

imshow(RGB)
hold on
visboundaries(BW1,'Color','r');
visboundaries(BW2,'Color','g');
visboundaries(BW3,'Color','b');
title('Seed Regions')

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

L = imseggeodesic(RGB,BW1,BW2,BW3,'AdaptiveChannelWeighting',true);

Загрузите сегментацию истинности на земле изображения.

L_groundTruth = double(imread('yellowlily-segmented.png'));

Визуально сравните результаты сегментации с основной истинностью.

figure
montage({label2rgb(L),label2rgb(L_groundTruth)})
title('Comparison of Segmentation Results (Left) and Ground Truth (Right)')

Вычислите индекс подобия Dice для каждой сегментированной области. Индекс подобия Диса заметно меньше для второй области. Этот результат согласуется с визуальным сравнением результатов сегментации, которая ошибочно классифицирует грязь в правом нижнем углу изображения как листья.

similarity = dice(L, L_groundTruth)

similarity = 3×1

    0.9396
    0.7247
    0.9139