Удаление концентрических кругов с изображения в Matlab

У меня есть данные изображения ОКТ, которыми я пытаюсь манипулировать:

введите здесь описание изображения

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

Первоначально я думал, что могу использовать imfindcircles() для идентификации этих кругов, но алгоритм не может обнаружить большинство из них независимо от указанных радиусов. Методы, основанные на связанных компонентах (bwareaopen, regionprops), также не работают, поскольку концентрические круги кажутся связанными с более крупной структурой.

У меня была одна мысль, что я могу преобразовать в полярные координаты, так как центр концентрических кругов всегда будет в центре изображения, создавая это изображение:

width = 512
greyImg = ImToPolar(rgb2gray(rgb_img), 0, 1, width, width);

введите здесь описание изображения

Это должно быть значительно проще в обработке. У кого-нибудь есть идеи?


comment
Вы вставили неправильные изображения.   -  person Yves Daoust    schedule 12.11.2015
comment
Ну стыдно, исправили   -  person    schedule 12.11.2015
comment
Декартово преобразование в полярное не является хорошей идеей. Это должно быть выполнено с шестом, помещенным в центр круга, который вы не знаете.   -  person Yves Daoust    schedule 12.11.2015
comment
Центр круга всегда будет в центре изображений, созданных с использованием этого модальности изображения. Извините, я должен был указать это в своем исходном сообщении.   -  person    schedule 12.11.2015


Ответы (1)


Накопить все столбцы изображения по горизонтали. Тогда вертикальный профиль позволит вам легче находить радиусы окружностей (они отображаются в виде пиков). Зная внешний радиус, вы можете извлечь круглую область интереса из исходного изображения.

введите здесь описание изображения

person Yves Daoust    schedule 12.11.2015