Исследование сходимости итерационной процедуры определения параметров движения изображений методом функционализации

Размещено на http://www.allbest.ru/

Самарский государственный технический университет

ИССЛЕДОВАНИЕ СХОДИМОСТИ ИТЕРАЦИОННОЙ ПРОЦЕДУРЫ ОПРЕДЕЛЕНИЯ ПАРАМЕТРОВ ДВИЖЕНИЯ ИЗОБРАЖЕНИЙ МЕТОДОМ ФУНКЦИОНАЛИЗАЦИИ

П.К. Кузнецов, Е.Ю. Чекотило, Б.В. Мартемьянов

Аннотация

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

Ключевые слова: метод функционализации, итерационная процедура, параметры движения, изображение.

Annotation

RESEARCH ITERATIVE PROCEDURE CONVERGENCE FOR COMPUTING IMAGE MOTION PARAMETERS BY FUNCTIONALI ZATION METHOD

P.K. Kuznesov, E.Yu. Chekotilo, B.V. Martemyanov Samara State Technical University

In this paper, we present an iterative procedure for computing image motion parameters. The proposed procedure is based on a functionalization method for computing velocity of moving objects by image analysis.

Key words: functionalization method, iterative procedure, motion parameters, image.

Введение

Одной из актуальных задач технического зрения автономных мобильных роботов является задача восстановления изображений окружающей обстановки на основе данных, получаемых от видео-сенсоров. Ее можно решить путем перехода от анализа изображения окружающей обстановки к анализу динамики векторного поля скоростей движения изображения подвижных яркостных объектов окружающей среды методом функционализации. Этот метод является обобщением известных градиентных методов [1-3] и основан на использовании соотношения, функционально связывающего параметры движения с изменениями измеримых характеристик изображения (ФС-уравнения).

Метод функционализации состоит в формировании с помощью изобразительной системы на чувствительной поверхности ПИ изображения движущегося яркостного объекта, выделении на изображении (по параметрам пространства или времени) окон анализа, измерении в выделенных окнах обобщенных производных по времени и пространственным координатам от функционала специального вида [4], определенного в окне анализа на функциях распределения освещенности изображения, вычислении для каждого окна анализа коэффициентов ФС-уравнений с использованием этих производных, составлении системы ФС-уравнений на множестве окон анализа в количестве, достаточном для получения единственного решения этой системы уравнений.

Практика использования систем технического зрения автономных мобильных роботов, построенных на алгоритмах, реализующих этот метод, показала высокое быстродействие и хорошие метрологические характеристики. Однако при всей универсальности метода функционализации не удается избежать погрешности измерения скорости движения изображения, вызванной неточностью аппроксимации производной по времени первой разностью. Эта погрешность может быть уменьшена до допустимого значения за счет применения компенсационного варианта метода функционализации с использованием физического или «электронного» (режим ВЗН ПЗС матриц) подслеживания движения изображения [4-5].

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

изображение итерационный функционализация движение

Размещено на http://www.allbest.ru/

В компенсационном методе используется итерационная процедура совмещения изображений (рис. 1). Задача совмещения формулируется следующим образом: на анализируемом изображении, заданном функцией освещенности ( - вектор смещения в системе координат плоскости изображения), выделяются две пересекающиеся прямоугольные области анализа и с изображениями и соответственно, т.е. , причем область является результатом преобразования области :

, (1)

где - матрица преобразования вращения.

Согласно методу функционализации, на каждом шаге итерационной процедуры формируется система алгебраических уравнений, по которой определяются оценки составляющих вектора смещения изображения [6]. Для этого в каждой области анализа выделяются идентичные системы N окон анализа D_ij размером (i = 1, 2 - номер области; j = 1, 2, …, N - номер окна), образующих покрытие предполагаемой зоны пересечения рассматриваемых областей анализа (см. рис. 1).

Для каждой пары окон составляются уравнения функциональной связи [6] относительно компонентов вектора смещения:

(2)

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

Для N пар окон анализа формируется переопределенная система ФС-уравнений:

, (3)

решение которой находится методом квазиобращений:

. (4)

В итерационной процедуре на каждом шаге определяются оценки составляющих вектора смещения изображения:

(5)

где = 0, 1, 2,…, n - номер итерации; ;

- матрица «коэффициентов обратной связи», значения которых подлежат определению из условий сходимости процесса (5);

- предполагаемая оценка смещения изображения области анализа относительно области .

После вычисления оценки смещения изображения на i-том шаге производится сдвиг области анализа на величину получаемой оценки.

Описанные процедуры повторяются до тех пор, пока компоненты вектора смещения на n-ном шаге не станут меньше заданных величин , и (в идеальном случае ).

В реализации итерационного метода удобно использовать «пирамидальную» функцию веса [5]:

(6)

Параметрами функции веса являются размеры окон анализа .

Получим аналитические оценки сходимости итерационного метода с «пирамидальной» функцией веса без учета аддитивного шума для тестового моночастотного изображения, которое зададим функцией освещенности вида:

, (7)

где - первая и вторая строки матрицы .

В силу сложности получения аналитических оценок сходимости итерационного метода для общего случая (7) рассмотрим отдельно плоско-параллельное движение и вращение в плоскости движения изображения. Для случая плоско-параллельного движения () система уравнений (5) принимает вид:

(8)

На основе известного критерия сходимости итерационных процессов [7] получим достаточный признак сходимости итерационной процедуры (8):

(9)

Соотношение (9) дает оценку области сходимости в пространстве начальных данных для тестового изображения (7).

Аналогичным образом получим оценку сходимости итерационного метода для случая чистого вращения в плоскости движения (угол крена и угол тангажа равны нулю). Итерационная процедура имеет вид:

. (10)

Достаточным признаком сходимости итерационной процедуры (10) является условие

. (11)

Условие сходимости (11) будет выполняться для всех таких, что

(12)

где - функция нескольких переменных: . Таким образом, сходимость итерационного процесса (11) при функции веса принятого вида зависит от размеров окон анализа и пространственной частоты тестового сигнала (7). Обозначим через максимальное (критическое) значение, отвечающее достаточному признаку сходимости (12). Считаем, что окна анализа имеют форму квадрата, т.е. .

В табл. 1 приведена аналитическая оценка зависимости количества итераций итерационного процесса (10) от частоты тестового сигнала при критическом значении коэффициента обратной связи . Условие (12) является лишь достаточным.

Нетрудно показать, что значение может быть увеличено, при этом скорость сходимости увеличивается.

Результаты исследования влияния на процесс сходимости итерационной процедуры (10) для тестового изображения вида (7) представлены в табл. 2.

Таблица 1

При увеличении сначала увеличивается количество итераций, а затем процесс (10) теряет устойчивость. Таким образом, можно сделать вывод, что наилучшая сходимость с точки зрения максимального быстродействия итерационного процесса достигается при . В табл. 2 приведены результаты, соответствующие допустимой погрешности . При других значениях значение остается в тех же пределах.

Таблица 2

Размещено на http://www.allbest.ru/

Для подтверждения достоверности высказанных предположений были проведены компьютерные эксперименты по оценке области сходимости итерационного метода по начальным данным и скорости сходимости для случая сложного движения (плоско-параллельное движение, вращение в плоскости движения). Эксперименты проводились в области начального смещения 50х50 пикселей для изображения типа «Горы» (рис. 2).

Результаты экспериментов сравнивались с результатами, полученными корреляционным методом с использованием для поиска экстремума взаимно-корреляционной функции совмещаемых изображений.

Результаты оценки сравнительного анализа областей сходимости представлены на рис. 3, скорости сходимости - на рис. 4.

Размещено на http://www.allbest.ru/

Проведенные эксперименты показали работоспособность метода и его существенное преимущество по точностным характеристикам и диапазону измеряемых величин по сравнению с корреляционным методом.

Заключение

Размещено на http://www.allbest.ru/

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

Библиографический список

1. Horn B., Schunck B. Determining Optical Flow// Artificial Intelligence 17. 1981. P. 185-203.

2. Black M.J., Anandan P. A framework for the robust estimation of optical flow // ICCV'93, May. 1993. P. 231-236.

3. Schalkov R.J. Image Motion Analysis Using Concept of Weak Solution to Distributed Parametr Systems // Proc. 1983 IEEE Computer Vision and Pattern Recognition Conf., Washington, DC. 1983. P. 232-239.

4. Кузнецов П.К., Семавин В.И. Метод определения параметров движения яркостного поля // Известия вузов. Приборостроение. 1990. №6. С. 26-30.

5. Кузнецов П.К., Семавин В.И., Мишин В.Ю., Владимиров М.В. Метод функциональных преобразований в задаче определения скорости движения яркостных полей/ / Вестник Самарского технического университета. 1994. №1. С. 66-76.

6. Кузнецов П.К., Мартемьянов Б.В., Семавин В.И., Чекотило Е.Ю. Метод определения вектора скорости движения подстилающей поверхности // Вестник Самарского технического университета. Серия Технические науки. 2008. № 2(22). С. 96-110.

7. Волков Е.А. Численные методы. М.: Наука, 1987. 248 с.

Размещено на Allbest.ru