CC BY
65
БЫСТРЫЙ СПОСОБ ПОСТРОЕНИЯ МОМЕНТНЫХ ИНВАРИАНТОВ НА ИЗОБРАЖЕНИИ С ИСПОЛЬЗОВАНИЕМ АППРОКСИМАЦИИ
Сергеев1 В.В., Титова2 О.А. 1 Институт систем обработки изображений РАН 2 Самарский государственный аэрокосмический университет
Аннотация
В данной работе описан быстрый способ построения характеристик изображения. В качестве признаков берутся моментные инварианты Ху. Предлагается использовать параллельно-рекурсивный алгоритм для расчета обобщенных моментов в "скользящем окне" по изображению, а затем аппроксимировать по ним значения инвариантов. Данный метод вычислительно более эффективен, чем прямой расчет.
1. Введение
Одним из хорошо известных подходов в распознавании образов на цифровом изображении является использование моментов различных порядков, а также моментных инвариантов. Для обработки цифровых изображений используются дискретные аналоги моментных характеристик. Поэтому формула момента порядка (к, s) записывается в следующем виде:
Ф1 = «20 + «02 , Ф2 = («20 -«02 )2 + «Ь
фз =(«зо - 3«12 )2 + (3«21 -«03 ) Ф4 = («30 + «12 )2 + («21 + «03 )
Ф5 =(«30 - 3«12 Х/^30 + «12 )[(«30 + «12 )2 - 3(/"21 + ^03 ) Ф6 = («20 — «02 )[(«30 +«12 )2 -(«21 +«03 )2 ]++4«11 («3 Ф7 = («21 - «03 )(«30 +«12 )[(«^30 +«12 )2 - 3(«21 +«03)
Основным отличием соотношений (2) является их инвариантность к повороту 2 — Б объекта. С тех пор множество работ было посвящено обобщению и улучшению инвариантных признаков, предложенных Ху, а также попыткам применить их в приложениях из разных областей знаний. Например, в [2] и [1] было предложено использовать моментные инварианты для распознавания силуэтов самолетов на изображении, в [9] моментные инварианты применяются при сравнении с эталоном и совмещении снимков, полученных со спутника. Ряд авторов используют инварианты для распознавания символов. В [7] и [6] предложено использовать инварианты для изменений контраста, а в [8] для классификации текстур, в [3] описаны инварианты к линейным преобразованиям.
Много усилий исследователи прилагают для нахождения эффективных алгоритмов вычисления моментов. Возможность параллельно-рекурсивного вычисления локальных моментных характеристик для обработки изображений при использовании информационной технологии "скользящего окна" показана в монографии [10] (глава 8). Целью данной работы стала разработка быстрого способа вычисления моментных инвариантов, построенных на основе локальных характеристик, рассчитанных в "скользящем окне".
M-1N-1
«ks = ZZ
mknsx(m, n), k, s = 0,1,...
(1)
Здесь х(т, п) - матрица, отсчетами которой являются яркости пикселей изображения; М х N -размеры изображения.
За последние сорок лет моментные инварианты, предложенные Ху, стали классическим инструментом распознавания образов. Приведенный им в [5] набор характеристик включает семь инвариантов:
(2)
+ (3«21 -«03 Х«21 +«03 )(«30 +«12 )2 "(«21 + «03 f 1 + «12 Х«21 +«03 ^
" (30 " 3«12 Х«21 +«03 Х^30 +«12 )2 — («21 +«03 )2 I
2. Построение моментных инвариантов в режиме "скользящего окна" Во многих случаях использование моментов, вычисленных по всему изображению, является необязательным. Можно ограничиться "локальными" моментами, рассчитанными только на ограниченной области - по так называемому «скользящему окну» обработки (например, при обнаружении и распознавании объектов ограниченность области вытекает из конечности размеров объектов). Обычно окно задается прямоугольным и симметричным
D :— M < m < M ;— N < n < N , где M ,N
па-
раметры, задающие границы окна по координатам. Тогда для каждой точки изображения с координатами (т, п) локальные моментные характеристики порядка (к, 5) определяются выражением:
«к*(т,п)= X Ет1п1х(т -т1,п -п), (3)
(т1,п1 )еБ
где х(т, п) - матрица, отсчетами которой являются яркости пикселей изображения; Б - область, задающая окно обработки. Размеры окна выбираются, исходя из условий решаемой задачи.
m=0n=0
Для такого окна можно применить параллельно-рекурсивный алгоритм вычисления моментных характеристик [4], в соответствии с которым организуется покоординатная рекурсия с использованием моментов младших порядков, определенных для предыдущего положения окна: к
як^ЬЕ^и^ чпгМ-м^
1=0
-(( +]))( -М
= ЕС^Я3(1,п2 -1)+(-М2)?Я(п1,П2 +М2)
}=0
+) ( -М -1),
0 < к < К,0 < 5 < К - к. (4)
В таком случае алгоритм построения моментных инвариантов будет выглядеть следующим образом:
• вычисление "скользящим окном" по изображению степенных моментов с применением параллельно-рекурсивного алгоритма,
• пересчет инвариантов по известным формулам.
Общую сложность и такого алгоритма можно оценить как сумму двух слагаемых:
и = ии + и1,
(5)
где
и и = иЯ + иЯ,
и я = (К +1)( К + 3)(К + 8) а 6 ; ия = К (К + 2)(К + 7)
6
Здесь К - порядок моментов; и - оценка сложности расчета треугольной матрицы степенных моментов, необходимых для вычисления инвариантов;
иЯ - количество операций сложения на каждый отсчет обрабатываемого изображения [4]; и Я -количество операций умножения на отсчет [4]; и1 -
сложность собственно пересчета инварианта, она зависит от того, какой именно инвариант выбран.
3. Вычисление моментных инвариантов через обобщенные моменты
Как было показано в работе [10], наиболее эффективным рекурсивным способом вычисляются обобщенные моменты с некоторым классом полиномиальных базисов:
Як5 п2 )= Е Е Чк (т1)Ч, (т2) *
(1,т2 (6)
* х( - п2 - т2 ),
где Чк(т) = Е<
одномерный полином к -го
г=0
порядка с коэффициентами {ак }}=0 (акк ф 0); В -
область, задающая окно обработки. Однако, несмотря на вычислительную эффективность, напрямую использовать обобщенные моменты для нахождения моментных инвариантов (2) нельзя. Возможный способ использования этих моментов состоит в том, чтобы сначала пересчитать по ним степенные моменты, а затем уже искать инварианты. Однако этот способ практически обесценивает выигрыш в вычислительных затратах, который имелся за счет использования обобщенных моментов.
В данной работе предлагается следующая процедура нахождения моментных инвариантов:
• расчет обобщенных моментов,
• вычисление инвариантов через обобщенные моменты.
Таким образом, сначала по изображению рассчитываются обобщенные моменты для полиномиального базиса (6), коэффициенты разложения которого по базису подбираются таким образом, чтобы вычислительная сложность алгоритма была минимальной [10]. При этом локальные обобщенные моменты порядка (к, 5) для каждой точки изображения с координатами (т, п) рассчитываются по формулам:
Я0 (т,п) = я0 (т -1,п)+ х(т + М*,п) - х(т -М* -1,п),
Як (т, п) = як (т -1, п)+ ¡ик-\ (т -1, п) + чк (-М *)х(т + М *, п), 0 < к < К;
Як 0 (т, п)= Як 0 (т, п - 1)+Яэ (, п + N')-я (, п - N * -1),
Як5(m, п) = Як5(m, п -1) + Як,5-1(m, п -1) + (-N* )Як ( п + N*) 0 < 5 < Б,
(7)
где х(т, п) - матрица, отсчетами которой являются яркости пикселей изображения; К, Б - наибольшие порядки моментов соответственно по двум координатам; М * х N * - размеры «скользящего окна».
Далее, чтобы получить моментные инварианты, предлагается аппроксимировать значение каждого инварианта с помощью вычисленных характеристик.
Будем использовать древовидную аппроксимацию, которая предполагает следующий порядок действий. В процессе предварительной обработки область значений используемых обобщенных моментов, представляющая собой многомерный ги-
перкуб, последовательно разбивается по осям и порождает в памяти ЭВМ древовидную структуру. В каждой из областей, полученных в результате разбиения, осуществляется аппроксимация функции решений. В простейших вариантах древовидная аппроксимация может быть кусочно-постоянной и кусочно-линейной [11]. Области с малой ошибкой аппроксимации принимаются за терминальные вершины дерева. Те области, в которых ошибка велика, подвергаются дальнейшему разбиению.
Общая сложность и снова может быть оценена как сумма двух слагаемых:
U = Vм + Ua
(8)
где
U "= U»+ U^iï,
U^= 2K + 2, Uiïiï= K,
Ua = Uakp or Uki,
Ukp = L +1, Ua = 2M + L.
Здесь K - порядок моментов; L - максимальная глубина построенного дерева (мы строим оценку «сверху»); M - количество обобщенных моментов,
используемых при аппроксимации; Uм - оценка сложности вычисления обобщенных моментов с полиномиальным базисом специального вида; -количество операций сложения на каждый отсчет обрабатываемого изображения [10]; - количество операций умножения на отсчет [10]; U^p - сложность кусочно-постоянной древовидной аппроксимации; - сложность кусочно-линейной древовидной аппроксимации.
4. Экспериментальное исследование Были проведены экспериментальные исследования погрешности полученных инвариантов в зависимости от сложности метода их аппроксимацион-ного вычисления. В качестве моментных инвариантов использовались семь инвариантов Ху. Эксперименты проводились на тестовом изображении "Лена" размерами 256x256 пикселов. На рисунке 1 приведен график относительной погрешности кусочно-постоянной аппроксимации третьего моментного инварианта (2) в зависимости от сложности расчета инварианта (количества операций на отсчет). £ 12 | 10 S 8
ж ¡6
о 4 s о
10 20 30 40 50 60 70 80 90 100 110
Сложность, операций/пиксел
Рис.1. График зависимости относительной погрешности расчета инварианта от сложности вычисления.
Отдельной точкой показана сложность метода точного вычисления моментных инвариантов в режиме "скользящего окна" (5). Можно заметить, что с увеличением сложности точность аппроксимации инварианта приближается к этому значению. Однако зачастую изображения содержат некоторый уровень шумов, который делает бессмысленным вы-
числение точных инвариантов. Поэтому можно воспользоваться предложенным методом, допустив наличие небольшой погрешности вычисления.
5. Заключение
В работе рассмотрен аппроксимационный метод вычисления моментных инвариантов в "скользящем окне" по изображению. Указанный метод позволяет снизить затраты времени при вычислении инвариантов. Заметим, что необязательно аппроксимировать именно инварианты Ху. В общем случае можно брать достаточно общие интегральные характеристики, рассчитанные по изображению.
6. Благодарность Работа выполнена при поддержке Министерства Образования РФ, Администрации Самарской области и Американского фонда гражданских исследований и развития (CRDF Project SA-014-02) в рамках российско-американской программы "Фундаментальные исследования и высшее образование" (BRHE), а также при поддержке Российского фонда фундаментальных исследований (РФФИ), проект № 04-01-96507.
Литература
1. Belkasim S.O., Shridhar M., Ahmadi M. Pattern recognition with moment invariants: a comparative study and new results// Pattern Recognition, 1991. Vol. 24. P. 1117-1138.
2. Dudani S.A., Breeding K.J., McGhee R.B. Aircraft identification by moment invariants// IEEE Trans. Comput., 1977. Vol. 26. P. 39-45.
3. Flusser J., Suk T. Degraded image analysis: an invariant approach// IEEE Trans. Pattern Anal. Mach. Intell.,1998. Vol. 20. P. 590-603.
4. Glumov N.I., Krainukov N.I., Sergeev V.V., Khramov A.G. The Fast Algorithm of Image Approximation in a Sliding Window // Pattern Recognition and Image Analysis, 1991. Vol. 4. P. 424-426.
5. Hu M.K. Visual pattern recognition by moment invariants// IRE Trans. Inform. Theory, 1962. Vol. 8. P. 179187.
6. Hupkens T.M., de Clippeleir J. Noise and intensity invariant moments// Pattern Recognition, 1995. Vol. 16. P. 371376.
7. Maitra S. Moment invariants// Proc. IEEE, 1979. Vol. 67. P. 697-699.
8. Wang L., Healey G. Using Zernike moments for the illumination and geometry invariant classifcation of multis-pectral texture// IEEE Trans. Image Process, 1998. Vol. 7. P. 196-203.
9. Wong R.Y., Hall E.L. Scene matching with invariant moments// Comput. Graphics Image Process., 1978. Vol. 8. P. 16-24.
10. Методы компьютерной обработки изображений // Под редакцией Сойфера В.А. М.: Физматлит, 2001. 784 c.
11. Сергеев В.В., Чернов А.В. Методы восстановления изображений, основанные на принципах теории распознавания образов // Тезисы докладов третьей Всероссийской конференции "Распознавание образов и анализ изображений: новые информационные технологии". Н.Новгород, 1997. Часть 2. С. 3-4.
