Уменьшите все объекты до линий в 2-м бинарном изображении или 3D двойном томе
B = bwskel(A)A до криволинейных линий шириной 1 пиксел без изменения существенной структуры изображения. Этот процесс, называемый скелетонизацией, извлекает осевую линию с сохранением топологии и числа Эйлера (также известного как характеристика Эйлера) объектов.
Считывание 2-D изображения в градациях серого в рабочую область. Отображение изображения. Представляющими интерес объектами являются темные нити на светлом фоне.
I = imread('threads.png');
imshow(I)
Для скелетизации требуется двоичное изображение, в котором пикселы переднего плана 1 (белый) и фон 0 (черный). Чтобы сделать исходное изображение пригодным для скелетизации, возьмите дополнение изображения так, чтобы объекты были светлыми, а фон темным. Затем выполните бинаризацию результата.
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);
Просмотрите исходное изображение и скелетированное изображение бок о бок.
montage({BW1,BW2},'BackgroundColor','blue','BorderSize',5)
Загрузка набора объемных данных в рабочую область. Имя набора данных: spiralVol. Отображение тома с помощью volshow.
load spiralVol.mat;
volshow(spiralVol);
Преобразовать spiralVol в двоичном формате, который требуется для bwskel функция.
spiralVolLogical = imbinarize(spiralVol);
Каркасировать форму спирали в наборе данных. Отображение скелетированного тома с помощью volshow.
spiralVolSkel = bwskel(spiralVolLogical);
.
В то время как оба bwskel и bwmorph может скелетировать 2-D изображения, вы можете получить различные результаты, используя bwmorph чем при использовании bwskel. Потому что они используют разные алгоритмы, bwskel функция использует 4-связность с 2-D изображениями; bwmorph использует 8-связность.
bwskel предполагает, что объекты переднего плана в двоичном изображении являются белыми (логический true). Если изображение имеет белый фон и черные объекты, используйте дополнение изображения в качестве входных данных для bwskel. Вы можете вычислить дополнение с помощью imcomplement.
bwskel функция использует преобразование медиальной оси.
[1] Та-Чих Ли, Рангасами Л. Кашьяп и Чонг-Нам Чу. Строительство скелетных моделей с помощью 3D средних алгоритмов утончения поверхности/оси. Компьютерное зрение, графика и обработка изображений, 56 (6): 462-478, 1994.
[2] Kerschnitzki, M, Kollmannsberger, P, Burghammer, M. et al. Архитектура сети остеоцитов коррелирует с качеством костного материала. Журнал исследований костей и минералов, 28 (8): 1837-1845, 2013.