10448

Использование фильтров и медианной фильтрации для подавления шумов различных видов

Лабораторная работа

Информатика, кибернетика и программирование

Использование фильтров и медианной фильтрации для подавления шумов различных видов. Подавление шумов одна из наиболее часто встречающихся задач в обработке изображений. Как правило шум является дельта-коррелированным. Исключением может являться лишь шум связанный ...

Русский

2013-03-26

46 KB

33 чел.

Использование фильтров и медианной фильтрации для подавления шумов различных видов.

Подавление шумов – одна из наиболее часто встречающихся задач в обработке изображений. Как правило, шум является дельта-коррелированным. Исключением может являться лишь шум, связанный с макронеоднородностью чувствительность фотоприемника, различными дефектами и подобными явлениями. В силу дельта-коррелированности шума в его спектре содержится гораздо больше высоких пространственных частот, по сравнению с незашумленным изображением. Поэтому любой фильтр высоких частот будет с той или иной степенью эффективности подавлять шум. Понятно, что этот фильтр не должен быть слишком "сильным", чтобы его влияние на полезный сигнал было минимальным. Чаще всего при подавлении шумов с использованием фильтров используется свертка исходного изображения со следующими матрицами:

,

,

.

Из этих матриц наиболее эффективно подавляет шумы матрица Н1, но она же и вносит в изображение наибольшие искажения. Матрица Н3 вносит наименьшие искажения, но и с шумами она борется хуже, чем остальные матрицы. Матрица Н2 является промежуточной между ними.

Рассматривая борьбу с шумами необходимо учесть, что нередко в изображениях присутствует шум вида "соль и перец". Это означает, что в изображении присутствуют пикселы, сигнал в которых резко больше или резко меньше относительно незашумленного изображения. Появление таких шумов связано с точечными дефектами ПЗС, сбоями, радиационными повреждениями. Уже визуальный анализ таких изображений подсказывает, что можно предложить более эффективный алгоритм с такими шумами, чем частотная фильтрация. Чаще всего в таких случаях применяется медианная фильтрация.

Сначала рассмотрим одномерный случай. Медианный фильтр представляет собой скользящее окно, содержащее нечетное количество элементов. Медианой последовательность из 2N+1 элементов называется элемент, значение которого меньше значений N элементов в последовательность и больше N элементов в последовательности. По смыслу медиана близка к среднему значению, однако в ряде случаев выгоднее использовать именно медиану. Рассмотрим это на примерах, в которых сравнивается медианный фильтр и одномерный усредняющий фильтр, который является одномерным аналогом фильтра H1.

Входной сигнал

Сигнал на выходе усредняющего фильтра

Сигнал на выходе медианного фильтра

1 1 1 1 3 3 3 3

1 1 1,33 1,66 3 3

1 1 1 1 3 3 3 3

1 1 2 3 4 5 6 6

1 1,33 2 3 4 5 5,66 6

1 1 2 3 4 5 6 6

1 1 1 2 2 1 1 1

1 1 1,33 1,66 1,66 1,33 1 1

1 1 1 2 2 1 1 1

1 1 1 2 1 1 1 1

1 1 1,33 1,33 1,33 1 1 1

1 1 1 1 1 1 1 1

Как видим, в 3 примерах из 4 медианный фильтр обеспечивает лучшие результаты по сравнению с усредняющим фильтром.

В двумерном случае окно превращается в прямоугольное окно размером 3 на 3, 5 на 5, 7 на 7 и т. д. элементов. В более сложных случаях иногда оказывается целесообразно использовать крестообразное окно, содержащее 9, 13, 17 … элементов.

0 0 1 0 0

0 0 1 0 0

1 1 1 1 1

0 0 1 0 0

0 0 1 0 0

Рассмотрим пример.

Исходное изображение

1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1

1 1 3 3 3 3 1 1

1 1 3 3 3 3 1 1

1 1 3 3 3 3 1 1

1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1

Изображение после медианного фильтра с окном размером 3 на 3 элемента

1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1

1 1 1 3 3 1 1 1

1 1 3 3 3 3 1 1

1 1 1 3 3 1 1 1

1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1

Изображение после медианного фильтра с крестообразным окном.

1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1

1 1 3 3 3 3 1 1

1 1 3 3 3 3 1 1

1 1 3 3 3 3 1 1

1 1 1 1 1 1 1 1

1 1 1 1 1 1 1 1

Список полезных функций MATLAB.

imread – чтение изображения из графического фвойла

rand – возвращает случайную величину из диапазона от 0 до 1 с равномерным распределением

randn – возвращает случайную величину из диапазона от -∞ до ∞ с нулевым математическим ожиданием и гауссовским распределением с дисперсией 1.

conv2 – двумерная свертка,

image, imagesc, imshow – вывод изображения на экран,

double – преобразование числе к типу double


Задание для лабораторной работы.

  1.  Получить у преподавателя файл с изображением, дисперсию для гауссовского шума, вероятность сбоя для шума "соль и перец"
  2.  Внести в файл следующие шумы (по отдельности):
    1.  гауссовский,
    2.  соль и перец.
  3.  Применить к каждому зашумленному изображению фильтры H1, H2, H3, медианный с окном 3х3, медианный с окном 5х5, медианный с крестообразным окном 5х5.
  4.  Сделать выводы об эффективности каждого фильтра для подавления тог или иного вида шумов.


 

А также другие работы, которые могут Вас заинтересовать

3604. Відношення і пропорції. 149.18 KB
  Тема. Відношення і пропорції. Мета: узагальнити і систематизувати знання учнів з теми; показати застосування математичних пропорцій у мистецтві та архітектурі;розвивати вміння застосовувати математику в проблемних ситуаціях; виховувати ерудова...
3605. КАЧЕСТВО САХАРА И ПУТИ ЕГО ПОВЫШЕНИЯ 4.64 MB
  Введение Увеличение производства сахара в мире, в том числе в Республике Беларусь, связано одновременно с возрастающими требованиями к его качеству. Имея отличные вкусовые качества и высокую калорийность, сахар является одним из самых важных пр...
3606. Водні ресурси України 8.51 MB
  Водні ресурси – це природне багатство, яке вимагає збереження і охорони особливо в теперішній час коли суспільство змінює свої погляди, розширює можливості. У сільськогосподарському виробництві використувується 72,2% території країни, а 57,5%...
3607. ВАНТАЖОПІДЙОМНА, ТРАНСПОРТУЮЧА ТА ТРАНСПОРТНА ТЕХНІКА 9.59 MB
  В посібнику описані сучасні конструкції вантажопідйомної, транспортуючої та транспортної техніки, яка використовується при переміщенні великої кількості вантажів в процесі виробництва в різних галузях народного господарства, в будівництві промислово...
3609. Банкротство предприятий и антикризисный менеджмент в современных российских условиях 887.5 KB
  Объективным процессом рыночной экономики, основанной на конкуренции, является постоянный переток капиталов в наиболее доходные сферы, перераспределение собственности от неэффективных хозяйствующих субъектов к эффективным. Осуществляется дан...
3610. Инженерная графика 9.34 MB
  Учебно-методическое пособие представляет базовый курс инженерной графики. Приводится необходимая информация для освоения курса инженерной графики и выполнения расчетно-графических работ. Содержатся основные положения нормативно-технической документа..
3611. Исследование некоторых эксплуатационных показателей трелевочных тракторов ОТЗ различной энергонасыщенности 8.82 MB
  Введение Основной задачей технического прогресса в лесозаготовительной промышленности на перспективный период является увеличение производительности труда за счет интенсификации общественного производства, т.е., за счет роста его энерговооруженности...
3612. Характеристика суждения 59.5 KB
  Характеристика суждения Придать логическую форму следующим высказываниям (выразить через суждения вида А, Е, I, О): а) Рыба дышит жабрами. А: Все S есть P б) Лень никогда не приводит к добру; E: все S не есть P в) Иногда металлы жидкие. I: Некоторые ...