Сегментация двухуровневого изображения с помощью быстро идущий метод
[ возвращается нормированное геодезическое расстояние сопоставляют BW,D] =
imsegfmm(___)D вычисленное использование быстрого метода похода. BW пороговая версия D, где все пиксели, которые нормировали геодезические значения расстояния меньше, чем thresh рассматриваются приоритетными пикселями и устанавливаются в trueD может быть порог на разных уровнях, чтобы получить различные результаты сегментации.
В этом примере показано, как сегментировать объект в изображении, использующем Быстро идущий Метод на основе различий в полутоновой интенсивности по сравнению с местоположениями seed.
Readimage.
I = imread('cameraman.tif'); imshow(I) title('Original Image')
Создайте маску и задайте местоположение seed. Можно также использовать roipoly создать маску в интерактивном режиме.
mask = false(size(I)); mask(170,70) = true;
Вычислите массив весов на основе полутоновых различий в интенсивности.
W = graydiffweight(I, mask, 'GrayDifferenceCutoff', 25);Сегментируйте изображение с помощью весов.
thresh = 0.01;
[BW, D] = imsegfmm(W, mask, thresh);
figure
imshow(BW)
title('Segmented Image')
Вы можете порог геодезическая матрица расстояния D использование различных порогов, чтобы получить различные результаты сегментации.
figure
imshow(D)
title('Geodesic Distances')
Этот пример сегментирует мозг из данных MRI человеческой головы.
Загрузите данные MRI.
load mri
V = squeeze(D);Визуализируйте данные.
sizeO = size(V); figure; slice(double(V),sizeO(2)/2,sizeO(1)/2,sizeO(3)/2); shading interp, colormap gray; title('Original');
Установите местоположения seed.
seedR = 75; seedC = 60; seedP = 10;
Вычислите веса на основе полутоновых различий в интенсивности.
W = graydiffweight(V, seedC, seedR, seedP , 'GrayDifferenceCutoff', 25);Сегментируйте изображение с помощью весов.
thresh = 0.002; BW = imsegfmm(W, seedC, seedR, seedP, thresh);
Визуализируйте сегментированное изображение с помощью поверхности ISO.
figure; p = patch(isosurface(double(BW))); p.FaceColor = 'red'; p.EdgeColor = 'none'; daspect([1 1 27/64]); camlight; lighting phong;
imsegfmm операции плавающей точки двойной точности использования для внутренних расчетов для всех классов кроме класса single. Если W имеет класс single, imsegfmm использует операции с плавающей точкой с одинарной точностью внутренне.
imsegfmm пиксели наборов с 0 или NaN значения веса к Inf в геодезическом расстоянии отображают D. Эти пиксели являются частью фона (логическая ложь) в сегментированном изображении BW.
[1] Сифианин, Дж. А. Методы установки уровня и быстро идущие методы: развивая интерфейсы в вычислительной геометрии, гидроаэромеханике, компьютерном зрении, и материаловедении, издательстве Кембриджского университета, 2-м выпуске, 1999.
Image Segmenter | activecontour | gradientweight | graydiffweight | graydist