Считайте 2D полутоновое изображение в рабочую область. Отобразите изображение. Предметы интереса являются темными потоками против светлого фона.

I = imread('threads.png');
imshow(I)

Скелетирование требует двухуровневого изображения, в котором приоритетными пикселями является (белый) 1, и фоном является (черный) 0. Чтобы сделать оригинальное изображение подходящим для скелетирования, возьмите дополнение изображения так, чтобы объекты были легки, и фон является темным. Затем binarize результат.

Icomplement = imcomplement(I);
BW = imbinarize(Icomplement);
imshow(BW)

Выполните скелетирование двухуровневого изображения с помощью bwskel.

out = bwskel(BW);

Отобразите скелет по оригинальному изображению при помощи функции labeloverlay. Скелет появляется как голубая строка 1 пиксель шириной по темным потокам.

imshow(labeloverlay(I,out,'Transparency',0))

Сократите маленькие шпоры, которые появляются на скелете и просматривают результат. Одно короткое ответвление сокращено от потока около центра изображения.

out2 = bwskel(BW,'MinBranchLength',15);
imshow(labeloverlay(I,out2,'Transparency',0))

Считайте двухуровневое изображение в рабочую область.

BW1 = imread('circbw.tif');

Скелетонайз возражает в изображении с помощью функции bwskel.

BW2 = bwskel(BW1);

Просмотрите оригинальное изображение и изображение skeletonized рядом.

montage({BW1,BW2},'BackgroundColor','blue','BorderSize',5)

Загрузите объемный набор данных в рабочую область. Именем набора данных является spiralVol. Отобразите объем с помощью volshow.

load spiralVol.mat;
volshow(spiralVol);

Преобразуйте набор данных spiralVol в двоичный формат, который требуется функцией bwskel.

spiralVolLogical = imbinarize(spiralVol);

Skeletonize спиральная форма в наборе данных. Отобразите skeletonized объем с volshow.

spiralVolSkel = bwskel(spiralVolLogical);

.