Скачать презентацию
Идет загрузка презентации. Пожалуйста, подождите
Презентация была опубликована 10 лет назад пользователемАртем Чуров
1 Оптический поток
2 PDE = уравнения в частных производных Видео = …
3 Компьютерное зрение Одна из основных задач: Получить информацию о движении в кадре
4 Компьютерное зрение Одна из основных задач: Получить информацию о движении в кадре Идентифицировать движение Определить его направление Определить скорость
5 Компьютерное зрение Одна из основных задач: Получить информацию о движении в кадре Положение и перемещения камеры обычно неизвестны Идентифицировать движение Определить его направление Определить скорость
6 Компьютерное зрение Одна из основных задач: Получить информацию о движении в кадре Положение и перемещения камеры обычно неизвестны Идентифицировать движение Определить его направление Определить скорость Ограничиваем задачу поиском относительного движения объектов в системе отсчета, связанной с камерой
7 Компьютерное зрение Относительное движение между двумя последовательными кадрами можно представить как векторное поле – оптический поток
11 Оптический поток Применение Визуальные эффекты Отслеживание объектов Автономная навигация роботов Системы видеонаблюдения
12 Оптический поток Методы нахождения Фазовая корреляция преобразование Фурье только прямолинейное движение все точки перемещаются одинаково
13 Взаимная корреляция сигналов
14 Оптический поток Методы нахождения Фазовая корреляция преобразование Фурье только прямолинейное движение все точки перемещаются одинаково Блочные методы поиск похожих блоков все точки блока перемещаются одинаково
15 Оптический поток Методы нахождения Фазовая корреляция преобразование Фурье только прямолинейное движение все точки перемещаются одинаково Блочные методы поиск похожих блоков все точки блока перемещаются одинаково Вариационные методы минимизация некоторого функционала
16 Оптический поток Вариационные методы формулировка в виде задачи оптимизации плотное поле скоростей для каждого пикселя высокая точность смещение не обязательно кратно размеру пикселя
17 Оптический поток Построение модели Дано: Непрерывная последовательность изображений Надо найти: Поле смещений в точкев момент временияркость
18 Оптический поток
19 Построение модели Предположение 1 Постоянство яркости пикселей Предположение 2 – малы, – гладкая функция В этом случае предыдущее равенство можно линеаризовать в точке
20 Построение модели Одно уравнение для двух неизвестных Некорректная задача с бесконечным числом решений Известно как проблема апертуры(aperture problem) наблюдая лишь за небольшой частью кадра невозможно однозначно определить направление движения
21 Построение модели Как движется шаблон? Одно уравнение для двух неизвестных Некорректная задача с бесконечным числом решений Известно как проблема апертуры(aperture problem) наблюдая лишь за небольшой частью кадра невозможно однозначно определить направление движения
22 Построение модели Как движется шаблон? По диагонали вправо? Одно уравнение для двух неизвестных Некорректная задача с бесконечным числом решений Известно как проблема апертуры(aperture problem) наблюдая лишь за небольшой частью кадра невозможно однозначно определить направление движения
23 Построение модели Как движется шаблон? По диагонали вправо? Только вправо? Одно уравнение для двух неизвестных Некорректная задача с бесконечным числом решений Известно как проблема апертуры(aperture problem) наблюдая лишь за небольшой частью кадра невозможно однозначно определить направление движения
24 Построение модели Как движется шаблон? По диагонали вправо? Только вправо? Или только вниз? Одно уравнение для двух неизвестных Некорректная задача с бесконечным числом решений Известно как проблема апертуры(aperture problem) наблюдая лишь за небольшой частью кадра невозможно однозначно определить направление движения
25 Как визуализировать векторное поле При помощи стрелок, предварительно понизив разрешение При помощи цвета: Направление – цвет Абсолютное значение – яркость
27 Как оценить качество рассчитанного потока Настоящий поток (ground truth) Вычисленный поток (estimated) Размер изображения Средняя угловая ошибка (AAE – Average Angular Error) Средняя ошибка по конечной точке (AEE – Average Endpoint Error)
28 Как оценить качество рассчитанного потока AAE = AEE = 7.22 AAE = 2.76 AEE = 0.37 AAE = 3.14 AEE = 1.53 Ground truth
29 Вариационные методы Основная идея поле смещений минимизирует некоторый энергетический функционал (data term) – штраф за отклонения от предположений о постоянстве какой-либо величины (например, яркости) (smooth term) – штраф за отклонения от предположений гладкости векторного поля параметр регуляризации – определяет гладкость получаемого векторного поля
30 МЕТОД HORN-SCHUNK B.K.P. Horn and B.G. Schunck, "Determining optical flow." Artificial Intelligence, vol 17, pp , 1981
31 Метод Horn-Schunck Предположение Поле смещений – гладкая функция в области Поле смещений минимизирует функционал data term – штраф за отклонения от предположения постоянства яркости smooth term – штраф за отклонение от предположения гладкости поля
32 Минимизация функционала Уравнения Эйлера-Лагранжа С граничными условиями
33 Метод Horn-Schunck Уравнения Эйлера-Лагранжа Граничные условия
34 Метод Horn-Schunck Аппроксимация оператора Лапласа
35 Метод Horn-Schunck Пространственные производные удобно вычислить заранее свертка столбцов(строк) с ядром 5x1 (1x5) Ядро
36 Метод Horn-Schunck Дискретизация уравнений Граничные условия уже учтены в дискретизованном уравнении Шаг сетки h обычно принимают равным 1 Система с 2xNxM неизвестными – метод Гаусса неприменим из-за высокой сложности и плохой устойчивости
37 Метод Якоби для решения СЛАУ Рассмотрим СЛАУ вида Итерационный метод
38 Метод Гаусса-Зейделя для решения СЛАУ Рассмотрим СЛАУ вида Итерационный метод
39 Метод релаксации (SOR) для решения СЛАУ Рассмотрим СЛАУ вида Итерационный метод
40 Метод Horn-Schunck Метод Якоби Просто распараллеливается Приближение на следующей итерации полностью определяется по приближению на текущей итерации – нет зависимости между разными пикселями
41 Метод Horn-Schunck Метод Гаусса-Зейделя i,j+1 i-1,ji,ji+1,j i,j-1 Нельзя напрямую использовать для параллельных вычислений
42 Красно-черная схема Гаусса- Зейделя Новое значение в красном узле зависит только от текущих значений в самом узле и в его «черных соседях» Новое значение в черном узле зависит только от текущих значений в самом узле и в его «красных соседях» Новые значения в узлах одного цвета можно вычислять параллельно Схема одной итерации Обновить значения в красных узлах Обновить значения в черных узлах Аналогично можно сделать параллельную версию SOR
43 Метод Horn-Schunck Схема вычислений Вычислить производные изображения (пространственные и временные) Установить начальное приближение - нулевой поток Выполнить некоторое количество итераций одного из итерационных методов Новые значения сохраняются в массив, отличный от массива со старыми значениями. В противном случае возникает конфликт чтения-записи. Ядро выполняет одну итерацию итерационного метода критерий останова
44 МЕТОД BROX ET AL T. Brox, A. Bruhn, N. Papenberg, J. Weickert High accuracy optical flow estimation based on a theory for warping, T. Pajdla and J. Matas (Eds.), European Conference on Computer Vision (ECCV) Prague, Czech Republic, Springer, LNCS, Vol. 3024, , May 2004
45 Метод Brox et al Рассмотренные методы чувствительны к изменению освещения
46 Метод Brox et al Рассмотренные методы чувствительны к изменению освещения Идея Рассмотреть другой data term Предположение Сохраняется градиент изображения Плюс: градиент не чувствителен к аддитивному изменению яркости Минус: чувствительность к шуму
47 Метод Brox et al Рассмотренные методы чувствительны к изменению освещения Идея Рассмотреть другой data term Предположение Сохраняется градиент изображения С самого начала ищутся только небольшие векторы смещения
48 Метод Brox et al Рассмотренные методы чувствительны к изменению освещения Идея Рассмотреть другой data term Предположение Сохраняется градиент изображения С самого начала ищутся только небольшие векторы смещения Идея Использовать нелинеаризованное уравнение сохранения яркости
49 Метод Brox et al Тогда штраф за отклонение от предположений о сохранении яркости и градиента яркости примет вид
50 Метод Brox et al Проблема: квадратичный штраф чувствителен к выбросам
51 Метод Brox et al Проблема: квадратичный штраф чувствителен к выбросам Идея Заменить квадратичный штраф выпуклой возрастающей функцией малый положительный параметр
52 Метод Brox et al Получим новый функционал
53 Метод Brox et al Штраф за отклонения от предположения о гладкости потока Предположение Поток кусочно-гладкий
54 Метод Brox et al Задача Найти u и v, минимизирующие функционал
55 Метод Brox et al Уравнения Эйлера-Лагранжа Введем обозначения: Использование z вместо t говорит о том, что соответствующее выражение НЕ производная, а разность, которую будем минимизировать
56 Метод Brox et al Уравнения Эйлера-Лагранжа Граничные условия Неймана
57 Метод Brox et al Численный метод Уравнения Эйлера-Лагранжа нелинейные Используем метод неподвижной точки для w Если есть смещения на расстояние, большее размера пикселя, то функционал может иметь множество локальных минимумов Идея Использовать пирамиду изображений разного разрешения
58 Метод неподвижной точки
59 Метод Brox et al Численный метод Совместим метод неподвижной точки с масштабированием Масштабировать будем с коэффициентом Для более плавного перехода от разрешения к разрешению коэффициент выберем близким к 1 Используем полную пирамиду изображений, начиная с наименьшего возможного разрешения (например, с 10x10)
60 Метод Brox et al Численный метод Итерации по разрешению назовем внешними Пусть k – номер изображения в пирамиде (0 – самое низкое разрешение) – решение на k-ой внешней итерации Осталась нелинейность в
61 Метод Brox et al Численный метод Устранение нелинейности в производных изображения
62 Метод Brox et al Численный метод Введем обозначения
63 Метод Brox et al Численный метод Аппроксимация слагаемого с дивергенцией
64 Метод Brox et al Численный метод Аппроксимация smooth term
65 Метод Brox et al Численный метод Аппроксимация smooth term
66 Метод Brox et al Численный метод Аппроксимация smooth term
67 Метод Brox et al Численный метод Аппроксимация smooth term
68 Метод Brox et al Численный метод Аппроксимация smooth term
69 Метод Brox et al Численный метод Оставшаяся нелинейность в Ѱ устраняется повторным применением метода неподвижной точки Начальные значения Полученная система уравнений линейная может быть решена одним из рассмотренных методов
70 Метод Brox et al Решение линейной системы Используем красно-черный метод релаксации (Red-Black SOR) Граничные условия
71 Метод Brox et al Решение линейной системы Используем красно-черный метод релаксации (Red-Black SOR)
72 Метод Brox et al Билинейная интерполяция для I(x+w)
73 Метод Brox et al Реализация Кадры храним в виде текстур Для продолжения решения с разрешения k на разрешение k+1 используем текстуры Общая схема метода Подготовить текстуры с изображениями Задать начальное значения для u и v Для каждого разрешения, начиная с наименьшего Установить начальное значение du, dv Подготовить данные для итераций SOR Выполнить некоторое число итераций SOR Обновить u, v Продолжить решение на большее разрешение Сохранить результат
74 Метод Brox et al Реализация Кадры храним в виде текстур Для продолжения решения с разрешения k на разрешение k+1 используем текстуры Общая схема метода Подготовить текстуры с изображениями Задать начальное значения для u и v Для каждого разрешения, начиная с наименьшего Установить начальное значение du, dv Подготовить данные для итераций SOR Выполнить некоторое число итераций SOR Обновить u, v Продолжить решение на большее разрешение Сохранить результат warping FP iteration lagged nonlinearity FP iteration solver FP iteration
75 Метод Brox et al Сложные шаблоны доступа к памяти Большая часть данных не меняется во время итераций Используем текстуры Текстуры + shared memory работает медленнее чем просто текстуры Для объединения запросов при записи индексы начала строк округляем до значений, кратных 16
77 Дополнительные материалы:
Еще похожие презентации в нашем архиве:
© 2024 MyShared Inc.
All rights reserved.