CC BY
63
Лазарев ИЮ.
Аспирант, Самарский государственный технический университет РАЗНОСТНО - КОРРЕЛЯЦИОННЫЙ МЕТОД ИЗМЕРЕНИЯ ДВИЖЕНИЯ ИЗОБРАЖЕНИЯ
Аннотация
В статье представлены результаты компьютерного моделирования автокорреляционных и взаимно корреляционных функций разностных изображений видеоряда движущейся сцены и оценки их использования для измерения вектора скорости движения. Ключевые слова: разностное изображение, вектор скорости, корреляционная функция.
Lazarev N.Y.
Postgraduate student, Samara State Technical Univetsity
RAZNOSTNO - THE CORRELATION METHOD OF MEASUREMENT OF MOVEMENT OF THE IMAGE
Abstract
Results of computer modeling of autocorrelated and mutually correlation functions of differential images of a video series of a moving scene and assessment of their use are presented in article for measurement of a vector of speed of movement.
Keywords: differential image, speed vector, correlation function.
Движение является одним из главных факторов анализа динамических сцен. Не случайно биологические зрительные сенсоры используют движение для выделения интересующих их объектов из фона. Например, лягушки реагируют только на движущихся насекомых. Некоторые устройства технического зрения работают по такому же принципу. Известны видеорегистраторы, включающие автоматическую запись видеосигнала только при наличии динамической составляющей на контролируемом объекте. При этом важен факт самого движения, а не изменения освещенности или других параметров изображения.
Появление новых стандартов видеокодирования, таких как H.264 и MPEG-4, основано на широком использовании векторов движения динамических блоков изображения, остающихся практически неизменными на протяжении ряда кадров [1]. Закодированная информация о таком блоке на начальном кадре используется далее путем переноса на следующих кадрах. Не менее важна задача автоматической стабилизации изображений на основе контроля смещений [2]. Другая важная задача - восстановление смазанных движением изображений методом деконволюции [5]. Для ее эффективного решения необходимо знать функцию рассеяния точки (ФРТ размытия), вычисляемую на основе параметров движения.
Большой класс задач, связанных с определением движения изображения, относится к управлению движением автономных роботов, автономной навигации беспилотных летательных и космических аппаратов [7-9].
В общем виде оценка движения локальных областей изображения и отдельных точек решается с помощью уравнения оптического потока [3]. Для условия неизменности освещенности сцены за время анализа это уравнение имеет вид:
VE-V + — = 0,
dt
где: VE = [—, —J - пространственный градиент яркости в точке (x, y);
V - вектор скорости оптического потока.
Если предположить, что в окрестности каждой точки изображения значение оптического потока одинаково, то можно записать основное уравнение оптического потока для всех пикселей окрестности и решить полученную систему для всех пикселей окрестности методом наименьших квадратов, как это сделано в алгоритме Лукаса - Канаде [4].
Использование локальных областей (блоков, эталонов) лежит в основе большого класса корреляционно-экстремальных систем (КЭС). Работа таких систем подробно рассмотрена в литературе [2, 8, 9] и здесь не анализируется. Отметим только, что КЭС не позволяют измерять скорость при движении, как в прямом, так и в обратном направлении. Для реализации такой функции требуется два коррелятора, один из которых работает в прямом, а другой - в обратном направлении [2].
Из всего перечисленного множества задач рассмотрим те из них, которые требуют оценки глобального движения, когда вектор движения элементов изображения одинаков для всех точек. Такая ситуация имеет место, когда достаточно удаленная сцена движется без вращения относительно неподвижной камеры или наоборот, камера движется относительно сцены.
Пусть видеокамера формирует последовательность кадров (фреймов) полутоновых черно-белых изображений движущейся сцены с межкадровым периодом Тк и временем экспозиции г << Тк. Примем в первом приближении, что освещенность сцены достаточна, т мало, и смазом изображения за время экспозиции можно пренебречь. Тогда все элементы кадров, кроме краевых, сдвигаются на величину Д = V-N Тк, а скорость движения изображения равна
V= Д /(N- Тк), (1)
где N - число кадров, между которыми определяется сдвиг.
Отметим, что все параметры привязаны к изображениям: сдвиг измеряется в пикселях (пик), скорость - в пик/сек.
Один из наиболее простых и распространенных способов оценки межкадрового сдвига между фреймами /;(х,у, t;) и lj(x,y, ГД взятыми в моменты времени t; и tj, основан на вычислении поэлементной разности соответствующих фреймов [7]. Такая разность вычисляется по формуле
Ri,j(x,y) =
(2)
1, если \li(x,y,ti) — lj{x,y,tj)\ > 0 0, в противном случае где 0 - значение порогового уровня.
Формула (2) позволяет формировать так называемые разностные изображения. Очевидно, если исходное изображение 1(х,у) содержит значительные однородные области, то разность яркостей элементов, принадлежащих этим областям, будет минимальной, меньше порога 0. На разностном изображении Rij(x,y) они принимают значение 0. Для областей с высоким градиентом (кромки, 33
границы объектов) элементы в результате движения существенно отличаются по яркости, и они принимают значение 1. Таким образом, при формировании разностных изображений осуществляется своего рода двухуровневая сегментация (бинаризация) движущихся областей на изображении.
Результаты, представленные в данной публикации, получены путем компьютерного моделирования в среде MATLAB. На рис. 1 показана последовательность трех полутоновых изображений f1, f2, f3 размером MxM (M=256) движущегося ландшафта.
33
----------------------------------------------------------------Уиз
Рис. 1. Исходные изображения. Показана ось времени t и вектор скорости изображения.
Предполагается, что наблюдатель с камерой находится на достаточно большом расстоянии от объектов ландшафта и все точки изображения движутся с одинаковой скоростью Уиз. В результате движения элементы соседних изображений сдвигаются на А пикселей (на рисунке А = 20) в направлении вектора движения.
На рисунке 2а представлены разностные изображения Й21 и й32 , полученные при 0 =50. Как видно из примера, разностные изображения содержат много единичных пикселей и мелких областей. Хотя для соседних фреймов существенная разница в шумовых компонентах маловероятна, тем не менее, мелкие градиентные детали проявляются существенно. Для их устранения применим морфологические операции [6]. Так, для устранения изолированных точек и мелких областей хороший эффект дает операция эрозии Re = R © В = {г: Br с R] со структурным элементом B. После эрозии целесообразно нарастить оставшиеся объекты с помощью операции дилатации: Rd = R ф В = {г: Вг П R ^0}.
Последовательное применение этих операций называется размыканием бинарного изображения RoB = (R© В) фВ. Его результат показан на рис. 26. В качестве структурного элемента использован квадрат 5x5.
Рис. 2. а) Разностные изображения Й21 и R32 б) размыкание R21 и R32 Исследуем корреляционные функции разностных изображений. Сформируем суммарное изображение Rs = Й21 + й32 и двойное разностное Rs = R21 — R32 и определим их автокорреляционные функции (АКФ). Поскольку разностные изображения R2i и R32 сильно коррелированны, справедливы следующие соотношения:
АКФ (йг) = АКФ (й21) + АКФ (й32) + ВКФ (й21,й32) + ВКФ(Й32,Й21),
АКФ(Й5) = АКФ (й21) + АКФ (й32) — ВКФ (й21,й32) — ВКФ(Й32,Й21), где через ВКФ обозначена взаимно корреляционная функция. (3)
Рис.3. Сечения АКФ соответственно Rz и Rs в плоскости, совпадающей с вектором движения.
Первые два слагаемых в формулах (3) определяют основную форму и абсолютный максимум АКФ в центральной точке, вторые два слагаемых дают локальные максимумы для Rs и минимумы для Rs, смещенные соответственно влево и вправо относительно центра АКФ на величину Д. Очевидно, минимумы на АКФ(Й5) выделить проще на фоне абсолютного максимума, расстояние между ними равно 2Д и его измерение может быть использовано в формуле (1) для определения скорости движения изображения. Тем не менее, в силу симметрии АКФ, возникает проблема неоднозначности при определении направления вектора скорости (его начала и конца). Решить ее можно используя взаимно корреляционные функции. Форма ВКФ (йг, й32) имеет два ярко выраженных пика, соответствующих началу и концу вектора движения изображения. Один из пиков, соответствующий начальной точке вектора V , расположен всегда в центре корреляционной плоскости. Второй пик соответствует конечной точке вектора V и перемещается по корреляционной плоскости при изменении модуля и направления вектора скорости, что позволяет однозначно измерить вектор движения изображения в диапазоне 0^360гр. 34
34
Рис. 4. Корреляционная плоскость ВКФ (Дг, Д32). Яркие точки соответствуют корреляционным пикам.
Очевидно, для расчета вектора скорости достаточно найти положение второго пика, т.к. первый пик всегда фиксирован в центре корреляционной плоскости, как показано на рис. 4. В этом случае отпадает необходимость суммирования разностных изображений, что упрощает задачу. Достаточно найти ВКФ (Д21, Д32) и определить положение максимума относительно центральной точки корреляционной плоскости.
Проведенное моделирование показало эффективность использования описанной методики для измерения вектора скорости движения изображения. Выбор порога 0 не столь критичен, как при яркостной бинаризации исходных изображений. В модели порог менялся в диапазоне от 0,05 до 0,8 максимальной яркости.
Для повышения быстродействия вместо вычисления ВКФ по классической формуле в модели опробована более простая метрика сравнения разностных изображений логического типа
SXOR = Ef=iSf=i-(fi2i © Я32),
что позволило сократить временные затраты на поиск максимума приблизительно в 6 раз. Диапазон измеряемых сдвигов Д составляет Н40% от M при сохранении погрешности в пределах одного пикселя. При больших сдвигах существенное влияние оказывают краевые эффекты.
Литература
1. Ян Ричардсон. Видеокодирование. Н.264 и MPEG-4 — стандарты нового поколения. - М.: Техносфера, 2005. - 368 с.
2. Еськов Д.Н., Ларионов Ю.П. и др. Автоматическая стабилизация оптического изображения / Под общ. ред. Д.Н. Еськова, В.А. Новикова. - Л.: Машиностроение, 1988. - 240 с.
3. Шапиро Л., Стокман Дж. Компьютерное зрение / Пер. с англ. - М.: БИНОМ. Лаборатория знаний, 2006. - 752с.
4. Lucas B.D. An Iterative Image Registration Technique with an Application to Stereo Vision / B.D. Lucas, T. Kanade //
Proceedings of the 7th international joint conference on Artifical intelligence. - 1981. - Vol. 2. - P. 674-679.
5. Гонсалес Р., Вудс Р., Эддинс С. цифровая обработка изображений в среде MATLAB. - М.: Техносфера, 2006. -616 с.
6. Фу К., Гонсалес Р., Ли К. Робототехника: Пер. с англ. - М.: Мир 1989. - 624 с.
7. Грузман И.С., Киричук В.С. и др. Цифровая обработка изображений в информационных системах: Учеб. Пособие. -Новосибирск: Изд-во НГТУ, 2002. - 352 с.
8. Красовский А.А., Белоглазов И.Н., Чигин Г.П. Теория корреляционно- экстремальных навигационных систем. - М.: Наука, 1979. - 408 с.
9. Козубовский С.Ф. Корреляционные экстремальные системы. - Киев: Наукова думка 1973. - 223 с.
Кадыров А.С.1, Максутова З.А.2, Бахриденов А.Б.3
'Доктор технических наук, профессор, Карагандинский государственный технический университет 2Магистрант, Карагандинский государственный технический университет 3Магистрант, Карагандинский государственный технический университет НОВАЯ КОНСТРУКЦИЯ ДЛЯ ОЧИСТКИ ВЫХЛОПНЫХ ГАЗОВ
Аннотация
В статье акцентируется внимание на проблеме загрязнения окружающей среды выхлопными газами автомобиля, предлагается автомобильный глушитель с ультразвуковыми устройствами. Авторами приведены предварительные расчеты основных параметров ультразвукового глушителя, являющиеся базой для планирования эксперимента.
Ключевые слова: выхлопные газы, глушитель, ультразвук.
Kadyrov A.S.1, Maxutova Z.A.2, Bakhridenov A.B.3
1 Doctor of technical science, professor, Karaganda State Technical University; 2 undergraduate, Karaganda State Technical University;
3 undergraduate, Karaganda State Technical University A NEW CONSTRUCTION OF EXHAUST GAS PURIFICATION
Abstract
The article focuses on the problem of environmental pollution by exhaust gases of the car, the car is offered with a muffler ultrasonic devices. The authors presented preliminary calculations of the basic parameters of ultrasonic muffler which are the basis for the planning of the experiment.
Keywords: exhaust gase, muffler, ultrasound.
Эффективность и удобства при использовании автомобиля во всех хозяйственных и деловых областях, а также личной жизни человека привели к исключительно широкому его применению.
Первые автомобили появились в конце XIX столетия. В 1910 г. их насчитывалось около 300 тыс., а перед началом второй мировой войны уже 46 млн. Засорок с лишним послевоенных лет автомобильный парк вырос более чем в десять раз и в 1987 г. превысил полумиллиардный рубеж. В 1998г. число автомобилей достигло 700 млн.
Автомобильный парк растет быстрее, чем народонаселение. В настоящее время с конвейеров автозаводов всего мира ежегодно сходит около 50 млн автомобилей, то есть в среднем при двухсменной работе 170 автомобилей каждую минуту [1].
Ожидается, что к концу первого десятилетия XXI века парк автомобилей достигнет миллиардной отметки.
Практически все современные автомобили снабжены двигателями внутреннего сгорания, качествам которых автомобиль и обязан своим повсеместным распространением. При сравнительно небольшой массе этот двигатель развивает значительную
35
