Считайте 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);

.