2017-02-16 10 views
2

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

Line artifact

Я в настоящее время его обнаружение пути сравнения среднего значения каждого столбца. Если среднее значение меньше половины среднего значения как левой, так и правой сторон столбцов, то столбец считается артефактом линии. Затем он интерполируется как максимум левого и правого соседних пикселей.

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

Можете ли вы придумать лучший способ обнаружить эти «артефакты линий»? Обратите внимание на сильные границы объектов на изображениях и большой шум, включенный также в артефакт.

+1

Ах, эти извещатели деформируются ... Так часто в CT! –

+1

Поскольку эти изображения являются КТ, я считаю, что всегда будут вертикальные линии. Пробовали ли вы фильтр вертикальной фильтрации? Что-то вроде 'conv2 (img, [- 1 0 1])'? Это создаст две вертикальные белые линии рядом с вашими черными линиями. С ними, вероятно, легче справиться. Кроме того, и жаль, что я говорю это снова: Итерационные алгоритмы, вероятно, помогут уменьшить эффект этих линий в recosntruction –

+1

Конечно, я согласен! Я также наблюдаю за тем, как я люблю поведение, хотя у меня есть спектральная модель, и у меня есть некоторые проблемы с NLPV. Поэтому я также рассматриваю различные методы реконструкции. Но одна проблема за раз. Вертикальный фильтр - хороший старт, я редактирую вопрос. – Tapio

ответ

3
  1. Мы хотим найти вертикальные линии на изображении, чтобы сначала сверлить изображение с фильтром [1 -2 1]. Это даст большие значения для пикселей, которые ниже их вертикальных соседей.
  2. Суммировать все столбцы изображения.
  3. Найдите индекс столбца с максимальным значением. Эта колонка является проблематичной.
+0

«Суммировать все столбцы», вы имеете в виду сумму всех строк? – Florian

+0

Нет. Я имею в виду столбцы. После шага 1 пиксели на вертикальных линиях будут иметь высокие значения, поэтому суммирование всех этих значений по вертикали (colum) даст высокое значение. –

+0

Хорошо, я неправильно понял. Сумма всех столбцов звучит как форма суммы по всем столбцам (т. Е. По индексу столбца). Вы имеете в виду: для каждого столбца суммируйте все записи столбца (т. Е. Суммируйте по строкам). Немного неоднозначный для меня, но, возможно, это только я. – Florian

 Смежные вопросы

  • Нет связанных вопросов^_^