Сегментация изображений в градациях серого с использованием GMM

Я ищу функции для выполнения сегментации зашумленных медицинских изображений (оттенки серого) с помощью GMM (модели смесей Гаусса).

Я нашел в MATLAB:

gm = gmdistribution(mu,sigma)
idx = cluster(gm,X)

учитывая X, мое изображение в градациях серого.

Как бы вы определили mu и sigma? Какого размера они должны быть? И как бы вы их инициализировали?

Я пробовал следующее (учитывая размер изображения (576x720)):

mu = rand(5,size(X,2));
sigma = ones(720,720);
gm = gmdistribution(mu,sigma);
idx = cluster(gm,X);

но я получаю сообщение об ошибке:

Ошибка при использовании wdensity (строка 29)
Создана неправильная ковариация.

Ошибка в gmdistribution/cluster (строка 59)
log_lh=wdensity(X,obj.mu, obj.Sigma, obj.PComponents, obj.SharedCov, CovType);

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


matlab image-segmentation gmm
person SheppLogan    schedule 11.01.2019    source источник
comment
нет ответа? любой?   -  person SheppLogan    schedule 11.01.2019

Ответы (1)


arrow_upward
0
arrow_downward

Неправильная функция. Вы ищете fitgmdist(X,k) и вводите свою оценку количества объектов, которые нужно сегментировать, как k. Затем программа попытается вычислить мю и сигму, используя EM-алгоритм.

Эта плохо обусловленная ковариация, созданная предупреждением, является типичной, вы увидите это часто, если ваши данные зашумлены. Я предлагаю регуляризацию, изменив параметр "RegularizationValue", возможно, установив ограничения на ковариационные структуры и/или фильтрация зашумленного изображения. У меня всегда были хорошие результаты при использовании BM3D (для 2D изображений) и BM4D -фильтры (для 3D-изображений).

Я буду рад помочь, если у вас есть какие-либо конкретные вопросы, но вам также придется выполнить домашнюю работу по этому вопросу. Обработка изображений сложна, вам нужно обработать кучу движущихся частей, прежде чем даже основные вещи начнут надежно работать.

person Tapio    schedule 15.01.2019
comment
Спасибо за Ваш ответ! На самом деле я уже не студент и когда-то видел эти вещи (недавно закончил магистратуру), но иногда я чувствую себя немного обескураженным из-за огромной сложности и количества вещей, хотя я нахожу их очень интересными. Вот почему я иногда прихожу просить помощи для кикстарта, и обычно я могу продолжить сам... Я знаю, что обработка изображений сложна, но это то, что меня больше всего интересует. - person SheppLogan; 16.01.2019
comment
Без проблем! И я не хотел быть грубым, я могу подчеркнуть с чувством подавленности. Как я уже сказал, я рад помочь, если вам нужна дополнительная помощь. - person Tapio; 17.01.2019