Исследование методов обработки изображений с помощью пакета IMAGE PROCESSING TOOLBOX: Методические указания к лабораторным работам, страница 15

Все многообразие разностных критериальных функций задается общей формулой

, где                   (23)

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

Наибольшее распространение получили критериальные разностные функции вида

   ;                               (24)

   .                                (25)

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

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

Экспериментальные исследования показали, что такой способ адаптации к возможным изменениям оказался неработоспособным из-за накопления погрешностей измерения координат, которые вызывают "случайные блуждания" объекта в эталоне и даже его уход из эталона. Это объясняется тем, что каждый акт определения координат путем вычисления экстремума той или иной критериальной функции сопровождается некоторой ошибкой, вызываемой эффектами пространственной дискретизации изображения, шумами, влиянием неравномерного фона. При многократной смене эталона очевидно неизбежное накопление ошибки и срыв слежения. Для борьбы с этим явлением было предложено несколько подходов. Один из них состоит в проверке целесообразности смены эталонного изображения путем сравнения глобального экстремума критериальной функции с некоторой пороговой величиной, например со значением критериальной функции, вычисленной в момент последней смены эталона. Это позволяет сократить число смен эталонного изображения, но принципиально не решает проблему накопления ошибки.

Эффективным методом обновления эталонного изображения объекта является использование в качестве эталона усредненного изображения объекта. Усредненное изображение может быть получено с использованием алгоритма экспоненциального сглаживания

 ,                 (26)

где ,  - координаты центра объекта, найденные в n-м кадре; b - параметр экспоненциального сглаживания, выбираемый в диапазоне от 0 до 1. При b=1 эталонное изображение не обновляется, при b=0 осуществляется покадровая смена эталона, а при 0<b<1 эталонное изображение формируется путем межкадрового усреднения изображения объекта. Достоинством такого подхода является малая чувствительность к ошибкам дискретизации, в отличие от алгоритмов с покадровой сменой эталона. Кроме этого, использование межкадровой фильтрации позволяет повысить точность за счет подавления шумов, присутствующих на изображении.

С точки зрения требуемого объема вычислений более предпочтительным выглядит использование критерия (24), так как он не требует возведения в квадрат. Но при больших размерах зоны поиска необходимость отыскания глобального минимума критерия (24) может привести к неприемлемо большим вычислительным затратам. Необходимость использования КЭА в реальном масштабе времени заставляет искать возможность снизить требуемое количество вычислений без заметного ухудшения качества решения задачи слежения.

Одним из возможных методов достижения этой цели является использование метода последовательных испытаний. Идея метода заключается в уменьшении количества вычислений в тех областях изображения, где искомый объект отсутствует. Эти области характеризуются быстрым ростом значения функции отличия с увеличением числа анализируемых точек. Считается, что объект отсутствует в анализируемой области изображения, если функция отличия растет быстрее, чем так называемая пороговая функция [2]. В данной работе реализован упрощённый метод последовательных испытаний. Это связано с особенностями реализации алгоритмов измерения координат в системе MATLAB. Суть этого метода состоит в  следующем: вычисление критериальной разностной функции происходит несколько раз, причём сначала по участку, составляющему небольшую часть от эталона. Для этого участка выполняется проверка следующего условия:

   ,                            (27)

где  – минимум критериальной функции, вычисленный в предыдущем (n – 1) кадре; – критериальная функция, вычисленная по частичному участку эталона на i-й ступени алгоритма в n–ом кадре;  k – весовой коэффициент;  – площадь (количество точек)  всего эталона; - площадь   частичного участка эталона.

Если условие (27) выполняется, то алгоритм фиксирует отсутствие в данной области текущего кадра искомого объекта  и переходит к рассмотрению следующей области. В противном случае алгоритм переходит на следующую ступень, увеличивая вчетверо площадь частичного участка эталона. Таким образом, изменяя коэффициент k и количество ступеней, можно добиться  отсева большей части точек без потери устойчивости работы алгоритма.

3. Описание программы лабораторной работы

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