Обработка изображений. Общая информация Страница курса Этот курс подготовлен и читается при поддержке. - презентация

1 Обработка изображений

2 Общая информация Страница курса Этот курс подготовлен и читается при поддержке

3 На предыдущей лекции S ML Wavelength Power Камера-обскура Фотоаппарат, глаз Цвет Психологическое свойство человека Свет описывается спектром Сетчатка глаза Колбочки 3х видов Трихроматическая теория 3 канала для пиксель Адаптация зрения Цветовой баланс, «баланс белого» Серые карточки, «серый мир»

4 Обработка изображений Семейство методов и задач, где входной и выходной информацией являются изображения. Примеры : Устранение шума в изображениях Улучшение качества изображения Усиления полезной и подавления нежелательной (в контексте конкретной задачи) информации

5 Обработка изображений Зачем обрабатывать? 1.Улучшение изображения для восприятия человеком цель – чтобы стало «лучше» с субъективной точки зрения человека 2.Улучшение изображения для восприятия компьютером цель – упрощение последующего распознавания 3.Развлечение (спецэффекты) цель – получить эстетическое удовольствие от красивого эффекта

6 Цифровое изображение Вспоминаем процесс получения цифрового изображения…

7 Почему оно может получиться плохо? Ограниченный диапазона чувствительности датчика Плохой функции передачи датчика

8 Что такое гистограмма? Гистограмма – это график распределения яркостей на изображении. На горизонтальной оси - шкала яркостей тонов от белого до черного, на вертикальной оси - число пикселей заданной яркости

9 Изменение контраста изображения Что может не устраивать в полученном изображении: Узкий или смещенный диапазон яркостей пикселей (тусклое или «пересвеченое» изображение) Концентрация яркостей вокруг определенных значений, неравномерное заполнение диапазона яркостей (узкий диапазон - тусклое изображение) Коррекция - к изображению применяется преобразование яркостей, компенсирующий нежелательный эффект: y – яркость пикселя на исходном изображении, x – яркость пикселя после коррекции.

10 Линейная коррекция Компенсация узкого диапазона яркостей – линейное растяжение: График функции f -1 (y)

11 Линейная коррекция Компенсация узкого диапазона яркостей – линейное растяжение:

12 Линейная коррекция Линейное растяжение – «как AutoContrast в Photoshop»

13 Линейная коррекция Линейная коррекция помогает не всегда!

14 Нелинейная коррекция График функции f -1 (y) y x

15 Нелинейная коррекция Нелинейная компенсация недостаточной контрастности Часто применяемые функции: Гамма-коррекция Изначальная цель – коррекция для правильного отображения на мониторе. Логарифмическая Цель – сжатие динамического диапазона при визуализации данных

16 Гамма-коррекция Графики функции f -1 (y)

17 Нелинейная коррекция График функции f -1 (y)

18 Цветовая коррекция Изменение цветового баланса Компенсация: –Неверного цветовосприятия камеры –Цветного освещения Ряд алгоритмов рассмотрели на предыдущей лекции

19 Цветовая коррекция изображений Растяжение контрастности (autolevels) Идея – растянуть интенсивности по каждому из каналов на весь диапазон; Метод: Найти минимум, максимум по каждому из каналов: Преобразовать интенсивности:

20 Растяжение контрастности

21

22 Шумоподавление Причины возникновения шума: –Несовершенство измерительных приборов –Хранение и передача изображений с потерей данных Шум фотоаппарата Сильное сжатие JPEG

23 Цель: подавление шума Пусть дана камера и статичная сцена, требуется подавить шум. Простейший вариант: усреднить несколько кадров Source: S. Seitz

24 Заменим каждый пиксель взвешенным средним по окрестности Веса обозначаются как ядро фильтра Веса для усреднения задаются так: Усреднение box filter Source: D. Lowe

25 Пусть f – изображение, g -ядро. Свертка изображения f с помощью g обозначается как f * g. Определение свертки f Source: F. Durand Соглашение: ядро перевернуто MATLAB: conv2 vs. filter2 (also imfilter)

26 Основные свойства Линейность: filter(f 1 + f 2 ) = filter(f 1 ) + filter(f 2 ) Инвариантность к сдвигу: не зависит от сдвига пиксела: filter(shift(f)) = shift(filter(f)) Теория: любой линейный оператор, инвариантный к сдвигу, может быть записан в виде свертки Slide by S. Lazebnik

27 Свойства Коммутативность: a * b = b * a Нет никакой разницы между изображением и ядром фильтра Ассоциативность: a * (b * c) = (a * b) * c Последовательное применение фильтров: (((a * b 1 ) * b 2 ) * b 3 ) Эквивалентно применению такого фильтра: a * (b 1 * b 2 * b 3 ) Дистрибутивность по сложению: a * (b + c) = (a * b) + (a * c) Домножение на скаляр можно вынести за скобки: ka * b = a * kb = k (a * b) Единица: e = […, 0, 0, 1, 0, 0, …], a * e = a Slide by S. Lazebnik

28 Детали реализации Размер результирующего изображения? MATLAB: filter2(g, f, shape) shape = full: output size is sum of sizes of f and g shape = same: output size is same as f shape = valid: output size is difference of sizes of f and g f gg gg f gg g g f gg gg full samevalid Slide by S. Lazebnik

29 Детали реализации Как происходит фильтрация по краям? Окно фильтра выходит за границы изображения Необходимо экстраполировать изображение Варианты: –clip filter (black) –wrap around –copy edge –reflect across edge Source: S. Marschner

30 Простейшие фильтры Original ? Source: D. Lowe

31 Простейшие фильтры Original Filtered (no change) Source: D. Lowe

32 Простейшие фильтры Original ? Source: D. Lowe

33 Простейшие фильтры Original Shifted left By 1 pixel Source: D. Lowe

34 Простейшие фильтры Original ? Source: D. Lowe

35 Простейшие фильтры Original Blur (with a box filter) Source: D. Lowe

36 Сглаживание с box-фильтром Результат сглаживание с помощью усреднения отличается от разфокусированного изображения Точка света, наблюдаемая с расфокусированного объектива, выглядит как кружок света, а усреднение дает квадратик Source: D. Forsyth

37 Сглаживание Точка света, наблюдаемая с расфокусированного объектива, выглядит как кружок света, а усреднение дает квадратик Другой способ: взвешивает вклад пикселей по окрестности с учетом близости к центру: fuzzy blob Slide by S. Lazebnik

38 Point Spread Function (PSF) Point spread function (PSF) – отклик оптической системы на точечный источник света (объект)

39 Ядро фильтра гаусса x 5, = 1 Source: C. Rasmussen

40 Выбор размера ядра Размер ядра дискретного фильтра ограничен Source: K. Grauman

41 Выбор размера ядра Эмпирика: полуразмер фильтра равен 3σ Slide by S. Lazebnik

42 Сглаживание фильтром гаусса

43 Сравнение Slide by S. Lazebnik

44 Свойства фильтра Гаусса Свертка с сами собой дает тоже фильтр гаусса Сглаживание несколько раз фильтром с маленьким ядром дает результат, аналогичный свертке с большим ядром Свертка 2 раза с фильтром радиуса σ дает тот же результат, что с фильтром радиуса σ2 Source: K. Grauman

45 Маленькая экскурсия к Фурье + Низкие частоты Высокие частоты

46 Фильтр Гаусса (gaussian blurring) Результат свертки фильтром гаусса и усреднения Исходное изображениеФильтр Гаусса с Sigma = 4 Усреднение по 49 пикселям (7x7) Важное свойство фильтра Гаусса – он по сути является фильтром низких частот.

47 Сепарабельность Source: D. Lowe Сепарабельное ядро Раскладывается в произведение двух одномерных фильтром гаусса

48 Пример * * = = 2D свертка Source: K. Grauman Фильтр раскладывается в произведение двух 1D фильтров: Свертка по строкам: Затем свертка по столбцу:

49 Сепарабельность Почему сепарабельность полезна на практике? Slide by S. Lazebnik

50 Виды шума Соль и перец: случайные черные и белые пиксели Импульсный: случайные белые пиксели Гауссов: колебания яркости, распределенные по нормальному закону Source: S. Seitz

51 Гауссов шум Мат.модель: сумма множества независимых факторов Подходит при маленьких дисперсиях Предположения: независимость, нулевое матожидание Source: M. Hebert

52 Сглаживание фильтрами большого радиуса подавляет шум, но размывает изображение Подавление гауссова шума Slide by S. Lazebnik

53 Подавление шума «соль и перец» Чем результат плох? 3x35x57x7 Slide by S. Lazebnik

54 Медианный фильтр Выбор медианы из выборки пикселей по окрестности данного Является ли фильтр линейным? Source: K. Grauman

55 Медианный фильтр В чем преимущество медианного фильтра перед фильтром гаусса? Устойчивость к выбросам (outliers) Source: K. Grauman

56 Медианный фильтр Salt-and-pepper noise Median filtered Source: M. Hebert MATLAB: medfilt2(image, [h w])

57 Медианный фильтр Результат применения медианного фильтра с радиусом в 7 пикселей к изображению с шумом и артефактами в виде тонких светлых окружностей.

58 Сравнение фильтров 3x35x57x7 Гауссов Медианный Slide by S. Lazebnik

59 Повышение резкости Что теряется при сглаживании? original smoothed (5x5) – detail = sharpened = Добавим дополнительно высокие частоты: originaldetail + α Slide by S. Lazebnik

60 Фильтр Unsharp Гауссов Единичный фильтр Лапласиан гауссиана изображение сглаженное изображение Единичный Фильтр Slide by S. Lazebnik

61 Пример Ядро свертки

62 Компенсация разности освещения Пример

63 Компенсация разности освещения Идея: Формирование изображения: Плавные изменения яркости относятся к освещению, резкие - к объектам. объектосвещение Изображение освещенного объекта

64 Выравнивание освещения Алгоритм Single scale retinex (SSR) Получить приближенное изображение освещения путем низочастотной фильтрации Восстановить изображение по формуле

65 Обрезание по порогу

66 Выравнивание освещения Пример

67 Компенсация разности освещения Пример /= Gauss 14.7 пикселей

68 Многомасштабный вариант Чаще всего выбирают 3 масштаба Веса одинаковые (1/3)

69 Multi-scale retinex Source by Z.Rahman et.al.

70 Метрики качества Как измерить похожесть двух изображений? Для оценки качества подавления шума, например исходное изображение искаженное изображение

71 Метрики качества Среднеквадратичная ошибка (MSE) Пиковое отношение сигнал/шум (PSNR) N – число пикселей M – максимальное значение пикселя

72 Метрики качества PSNR и MSE не учитывают особенности человеческого восприятия! Оригинал Далее будут использованы рисунки из статьи Wang, Bovik, Lu WHY IS IMAGE QUALITY ASSESMENT SO DIFFICULT?

73 Метрики качества У этих изображений одинаковые PSNR с оригиналом (примерно 25 dB) Повышена контрастностьДобавлен белый гауссов шум

74 Метрики качества И у этих – тоже примерно 25 dB! Добавлен импульсный шумРазмытие

75 Метрики качества И у этого – тоже! Артефакт блочности после JPEG

76 Метрики качества Вывод: PSNR не всегда отражает реальный видимый уровень искажений. Как улучшить? –Использовать функцию чувствительности глаза к различным частотам (CSF) –Использовать свойство маскировки –Использовать равномерные к восприятию цветовые пространства (CIE Lab, CIEDE2000) HVS models (human visual system)

77 Спецэффекты Рассмотрим Тиснение Негатив «Светящиеся» края Геометрические эффекты –Перенос/поворот –Искажение «Эффект стекла»

78 Тиснение Фильтр + сдвиг яркости, нормировка…

79 Цифровой негатив

80 Светящиеся края Медианный фильтра + выделение краев + фильтр «максимума»

81 Перенос/поворот Перенос: x(k; l) = k + 50; y(k; l) = l; Поворот: x(k; l) = (k. x0)cos(µ) + (l. y0)sin(µ) + x0; y(k; l) =.(k. x0)sin(µ) + (l. y0)cos(µ) + y0; x0 = y0 = (центр поворота), µ = /6

82 «Волны» Волны 1: x(k; l) = k + 20sin(2 l / 128); y(k; l) = l; Волны 2: x(k; l) = k + 20sin(2 k / 30); y(k; l) = l;

83 «Эффект стекла» x(k; l) = k + (rand(1, 1) – 0.5) * 10; y(k; l) = l + (rand(1, 1) – 0.5) * 10;

84 На следующей лекции Старые-добрые методы распознавания объектов Сопоставление шаблонов Выделение краёв Выделение контрастных объектов Геометрические и фотометрические инварианты