Шумоподавление для изображений

Содержание

Слайд 2

Виды и примеры шумов Шумы Стационарные Импульсные Смешанные Salt and pepper

Виды и примеры шумов

Шумы

Стационарные

Импульсные

Смешанные

Salt and pepper
Помехи в видео

Аддитивный белый Зерно пленки

Белый

шум – пиксели случайны, не коррелированны друг с другом
Гауссов/равномерный/треугольный шум – закон распределения амплитуд пикселей
Аддитивный шум – прибавляется к «чистому» изображению и не зависит от него

Шум + помехи в видео

Salt and pepper

AWGN

Слайд 3

Методы шумоподавления Шумы Стационарные Импульсные Смешанные Медианный фильтр Взвешенная медиана Ранговые

Методы шумоподавления

Шумы

Стационарные

Импульсные

Смешанные

Медианный фильтр
Взвешенная медиана
Ранговые фильтры

Bilateral filter
Non-Local Means
Wavelet thresholding
DCT, PCA, ICA
Анизотропная диффузия
Алгоритм

BM3D

Ранговые фильтры
Комбинированные методы

Salt and pepper

AWGN

Слайд 4

Импульсные шумы Медианная фильтрация Центральный элемент отсортированного массива яркостей Как быть с цветными изображениями? медиана 3х3

Импульсные шумы

Медианная фильтрация
Центральный элемент отсортированного массива яркостей
Как быть с цветными изображениями?

медиана

3х3
Слайд 5

Импульсные шумы Медианная фильтрация Центральный элемент отсортированного массива яркостей Как быть с цветными изображениями? медиана 7х7

Импульсные шумы

Медианная фильтрация
Центральный элемент отсортированного массива яркостей
Как быть с цветными изображениями?

медиана

7х7
Слайд 6

Импульсные шумы Медианная фильтрация Центральный элемент отсортированного массива яркостей Как быть с цветными изображениями? медиана 15х15

Импульсные шумы

Медианная фильтрация
Центральный элемент отсортированного массива яркостей
Как быть с цветными изображениями?

медиана

15х15
Слайд 7

Импульсные шумы Т.к. искажена лишь малая часть пикселей, то не надо

Импульсные шумы

Т.к. искажена лишь малая часть пикселей, то не надо фильтровать

все изображение!
Детектирование искаженных пикселей (простейший способ – анализ разности изображения с его медианной фильтрацией)
Фильтрация только искаженных пикселей

медиана 3х3

адаптивная медиана

Слайд 8

Стационарные шумы Простейшие методы Размытие изображения – вместе с шумом размывает

Стационарные шумы

Простейшие методы
Размытие изображения – вместе с шумом размывает детали
Размытие в

гладких областях – остается шум вблизи границ
Медианная фильтрация – хорошо подавляет импульсный шум, но удаляет мелкие детали
Слайд 9

Bilateral filter Адаптивные алгоритмы Bilateral filter усреднение окружающих пикселей с весами фотометрическая близость пространственная близость

Bilateral filter

Адаптивные алгоритмы
Bilateral filter
усреднение окружающих
пикселей
с весами

фотометрическая близость

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

Слайд 10

Bilateral filter Bilateral filter: художественное применение (при слишком сильном действии)

Bilateral filter

Bilateral filter: художественное применение

(при слишком сильном действии)

Слайд 11

Non-Local Means Адаптивные алгоритмы Non-local means (NL-means) – веса зависят от

Non-Local Means

Адаптивные алгоритмы
Non-local means (NL-means) – веса зависят от близости целых

блоков, а не отдельных пикселей

ν(xi,j) – блок вокруг
пикселя xi,j

Слайд 12

Non-Local Means Вычисление весов Способен сохранять текстуру изображения лучше, чем bilateral

Non-Local Means

Вычисление весов
Способен сохранять текстуру изображения лучше, чем bilateral filter

Иллюстрация из
Buades

et al 2005

Веса высоки для q1 и q2,
но не для q3

+

Слайд 13

Non-Local Means Достоинства и недостатки: Высокое качество результирующего изображения В исходном

Non-Local Means

Достоинства и недостатки:
Высокое качество результирующего изображения
В исходном варианте – очень

высокая вычислительная сложность
Ускоряющие расчет оптимизации:
Использование команд MMX/SSE для сравнения блоков
Разбиение изображения на блоки и усреднение целых блоков, а не отдельных пикселей
Промежуточный вариант между усреднением блоков и усреднением пикселей: усреднение маленьких блоков

+


Слайд 14

Non-Local Means Применение к видео Область поиска блоков можно расширить на

Non-Local Means

Применение к видео
Область поиска блоков можно расширить на соседние кадры

(сделать ее пространственно-временной)
Для ускорения просчета можно применять сравнение блоков только по Y-каналу в цветовой модели YCrCb (YUV)
Слайд 15

Вейвлетный метод Вейвлетное шумоподавление для изображений DWT Оценка уровня и спектра

Вейвлетный метод

Вейвлетное шумоподавление для изображений
DWT
Оценка уровня и спектра шума
Подавление вейвлет-коэффициентов (thresholding,

shrinkage)
Обратное DWT

+

Подавление шума различного масштаба


Отсутствие инвариантности к сдвигу
Плохая локализация энергии для наклонных границ

Слайд 16

Тригонометрические базисы Дискретные преобразования: Фурье, косинусное и синусное + Быстро вычисляются

Тригонометрические базисы

Дискретные преобразования: Фурье, косинусное и синусное

+

Быстро вычисляются


Зависимость от размера блока:

эффект Гиббса либо плохое подавление НЧ компонент шума

DFT

DCT

DST

Слайд 17

Метод главных компонент Идея: так выбрать базисные вектора, чтобы минимальным числом

Метод главных компонент

Идея: так выбрать базисные вектора, чтобы минимальным числом базисных

векторов можно было наилучшим образом приблизить данные из некоторого набора.
PCA – ортогональное линейное преобразование базиса, при котором первый вектор нового базиса соответствует направлению максимальной дисперсии данных, второй вектор – следующему направлению максимальной дисперсии и т.д.

Principal Component Analysis (PCA, KLT)

размерность данных = 2

Слайд 18

Метод главных компонент Решение: базисные вектора – собственные вектора ei матрицы

Метод главных компонент

Решение: базисные вектора – собственные вектора ei матрицы ковариации

Cx исходных данных x:
Альтернативное решение – с помощью сингулярного разложения матрицы исходных данных

Principal Component Analysis (PCA, KLT)

на рисунке размерность n = 2

Слайд 19

Шумоподавление Применение к блокам изображения 8x8: Метод главных компонент … … 64 базисных вектора

Шумоподавление

Применение к блокам изображения 8x8:

Метод главных компонент



64 базисных вектора

Слайд 20

Шумоподавление PCA-денойзинг изображений Блочное преобразование PCA Подавление коэффициентов в новом базисе

Шумоподавление

PCA-денойзинг изображений
Блочное преобразование PCA
Подавление коэффициентов в новом базисе
Обратное преобразование и наложение

блоков с перекрытием

Существующие подходы (Muresan/Parks, 2003)

+

Адаптация к деталям изображения (в т.ч. – к линиям любого наклона)


При больших блоках – эффект Гиббса,
при маленьких – не подавляется крупный шум

Слайд 21

Шумоподавление Результаты Зашумленное изображение Вейвлетное шумоподавление

Шумоподавление

Результаты

Зашумленное изображение

Вейвлетное шумоподавление