Серия История. Политология. Экономика. Информатика. 2010. № 19 (90). Выпуск 16/1
УДК 621.397
МЕТОД АНАЛИЗА РАСПРЕДЕЛЕНИЯ ЭНЕРГИЙ ИЗОБРАЖЕНИЙ ПО ЗАДАННЫМ ЧАСТОТНЫМ ИНТЕРВАЛАМ
Белгородский
государственный
университет
А.А. ЧЕРНОМОРЕЦ О.Н. ИВАНОВ
В работе изложен метод анализа распределения энергий изображений по заданным частотным интервалам, который основывается на вычислениях точных значений соответствующих долей энергий.
Ключевые слова: изображения, частотные интервалы, доли энергии, субполосные матрицы.
е-таїЬ chernomorets @bsu.edu.ru
Использование частотных представлений в задачах цифровой обработки изображений позволяет естественным образом сформулировать задачу анализа свойств изображений на основе соотнесения их с выделенными, исходя из тех или иных соображений, частотными интервалами специального вида. Одной из наиболее важных и часто определяемых характеристик является доля энергии изображения, попадающей в заданный частотный интервал. В частотной области традиционно обработка изображений предполагает применение дискретного преобразования Фурье (ДПФ) или быстрого преобразования Фурье (БПФ). Следует отметить, что получаемые оценки долей энергий затем используются для принятия тех или иных решений. Достижение максимально возможной обоснованности возможно только тогда, когда минимальны погрешности вычислений. Вместе с тем, применение алгоритма БПФ приводит к большим погрешностям при оценивании попадающей в выбранный частотный интервал доли энергии анализируемых изображений, особенно в случае интервала малой ширины. Поэтому целесообразно вычислять точные значения долей энергий, при этом время получения результата можно уменьшить при использовании распараллеливания вычислений при применении соответствующих архитектур вычислителей, что без особых проблем позволяют осуществить современные технические средства.
1. Частотные представления в задачах анализа изображений.
Возможность проведения анализа изображений на основе частотных представлений определяется тем, что в визуальных данных зачастую наблюдается периодичность или квазипериодичность отображаемых процессов. На изображении могут присутствовать повторяющиеся объекты, которые задают некоторую периодичность изменению яркости изображения.
В процессе преобразований изображение можно определить как двумерную функцию /(х, у) на плоскости. Известно, что любая функция, периодически воспроизводящая свои значения и удовлетворяющая условиям Дирихле, может быть представлена в виде суммы синусов и/или косинусов различных частот, умноженных на некоторые коэффициенты. В этом случае рассматривается преобразование Фурье [1].
Пусть /(х, у) - некоторая функция с конечной или неограниченной областью определения. Тогда при выполнении указанных условий [2, 3] справедливо следующее представление:
Равенства
еІіХ = сов(их ) + і 8Іп(и.х) и е1УУ = соБ(уу) + і БІп(уу)
определяют набор базисных функций, с помощью которых представляется функция (функция f (х, у) представляется в базисе Фурье).
Таким образом, функция f (х, у) , описывающая некоторое изображение, определяется через множество базисных функций, являющихся функциями синус и косинус различных аргументов. Каждая синусоидальная базисная функция характеризуется своей частотой, что позволяет говорить о частотных представлениях функций, задающих различные изображения.
Значение трансформанты Фурье F(и, V) функции f (х, у) определяется выражением [1, 4]:
ад ад
F(и, V) = 11 f (х, у)е~ 1ихе-тйхйу. (2)
-ад-ад
В содержательном отношении указанным соотношениям (1) и (2) можно придать разный смысл, который и определяет роль частотных представлений в задачах анализа и обработки изображений. В частности, представление (1) является выражением принципа суперпозиции в общем случае бесконечного количества периодических компонент, на которые может быть, согласно (2), разложена исходная функция. При этом справедливо равенство Парсеваля [2, 4, 5]:
Ь d 1 ад ад
Я f 2(х, у^у=732 Ц | F(и, V) |2 dudv,
~ ~ ад ад
которое можно преобразовать к виду
Ь d 1 ад ад
Я f2(ху^ф=222 ц 1F(и,v)|2
2
а с Г1 =0 г2 =0 (иУ)ЕПпп
где интервалы О определяют разделение частотной области на непересекающиеся интервалы.
Таким образом, оказывается возможным осуществить анализ энергетических характеристик исследуемой функции на основе частотных представлений, так как интегралы
Srír2 = № (и, V)2 dudv (3)
4п (и^еО^
определяют части энергии, попадающие в выбранные частотные интервалы. В частности, можно выделить частотные интервалы, в которых сосредоточена подавляющая доля энергии, или выделить квазипериодические компоненты исходной функции, энергии которых сосредоточены в разных интервалах.
В настоящее время для нахождения энергетических характеристик в большинстве случаев используется традиционный подход, основанный на использовании преобразования Фурье, а именно: пусть Ф = (fik ) - некоторое изображение, заданное матрицей яркости точек изображения (пикселей) размерностью Мх^ i = 1,2,...,М, к = 1,2,...,N, F (и, V) - его Фурье-образ.
Тогда дискретное преобразование Фурье определяется следующим выражением [1, 4]:
МД - 0'-1)( и 1) - ]—(к—1 )(v—1)
Р(и,V) = 22f,^e М е » ,
i=1 к=1
и = 1,2,...,М, V = 1,2,...,N .
Энергетический спектр [1] преобразования Фурье для дискретных значений частот и,V , и = 1,2,...,М, V = 1,2,...,N определяется как
Р(и, V) = ^ (и, V))2 = R2 (и, V) +12 (и, V), (4)
где величины R(u, V) и I(и, V) обозначают действительную и мнимую части величины F (и, V) соответственно.
Подход, основанный на использовании выражения (4) для нахождения энергетического спектра преобразования Фурье, обладает существенной вычислительной сложностью (трудоемкостью). В настоящее время для построения энергетического спектра традиционно используют быстрое преобразование Фурье (БПФ). Исследования [6, 7] показали, что алгоритмы, использующие преобразование Фурье и БПФ, не позволяют находить точные значения энергетического спектра одномерного сигнала в заданных частотных диапазонах. Одна из проблем ДПФ - приближенное нахождение части энергии в некотором частотном интервале. В данном случае оно строится как сумма дискретных значений, попадающих в интервал, что приводит к неточным результатам.
2. Метод анализа распределения энергий изображений по заданным частотным интервалам.
Для нахождения точных значений доли энергии РО изображения Ф в заданном частотном интервале О предлагается использовать выражение
Е
Ро= тО, (5)
Ео
где ЕО - часть энергии изображения Ф, соответствующая частотному интервалу О, Е0 -
энергия изображения Ф.
Частотный интервал, в котором предлагается определять точные значения долей энергий, является симметричной двумерной частотной областью О (субинтервал), задаваемой следующим выражением:
О :{О(и, V) | (и е[а,а2 ], V е[Д,А ]) и и (и е[а,а2 ] V е[-в2,-Д ]) и
и (и е[-а-а 1 V е[-вг,-А ]) и и (и е[-а-а ], V е[д,А ])}, (6)
где 0 < а,а2, Д,в2 < п.
Изображение будем рассматривать в виде непрерывной двумерной функции А(х, у) . Обозначим F(и, V) - Фурье-преобразование функции А(х, у) , то есть имеем
ад ад
F(и, V) = 11 А(х, у)е-!ихе-lvydxdy.
-ад-ад
Из выражения (3) следует, что точное значение части энергии ЕО изображения А (х, у) в частотной области О при известном Фурье-образе F (и, V) определяется выражением
ЕО = —П Ц|F(и, V)|2dudv. (7)
—П (и,у )еО
Значение энергии функции А(х, у) на всей области определения задается соотношениями
адад
Ео = 11АI2 = 11 /2(х у)^Ф.
-ад -ад
Можно показать, что соотношение для вычисления точного значения доли энергии функции А (х, у) в выбранном частотном интервале О (6) имеет следующий вид
ад ад ад ад
ни /(*1, У1)/(Х2 , У2 ) Аа (Х1, Х2 ) Ав (У1 , У2 )dXldУldX2dУ:
Л) ___ —ад—ад—ад—ад
го =
I I 12(Х, У)dxdy
(8)
где
Аа(*1, Х2) = <
2Cos
(а2 +а1)(*1 — *2) (а2 —а1)(Х1 — *2)
п( *1 — *2)
п
х1 = Х2
Ав(У1, У2) = <
2£05 (в2 + Д)(у1 — у2) £Цп (в — Д)(у1 — у2)
п( У1 — У2)
в2 — в1
п
, У1 = У 2
Следует отметить, что функции Аа( Хl, *2) и Ав( у1, у2) обладают
1
Аа (Х1, Х2 ) = Аа ( Х1,— Х2 ) ,
Ав( У^ У2) = Ав(—У1,—У2) и при конкретных значениях а1, а2, Д, в2 представляют собой волну (рис. 1)
свойством
(9)
(10)
(11)
ад —ад
2
2
а2 — а1
2
2
Рис. 1. Графическое представление функции Аа(Х1, Х2) : а - функция в виде изображения, б - значения функции при конкретном значении Х1
Рассмотрим далее дискретный случай. Дискретную конечную двумерную функцию будем задавать в виде матрицы, что соответствует форме описания изображений в цифровом виде.
Рассмотрим некоторое изображение Ф в виде матрицы Ф = (/1к), г = 1,2,...,М, к = 1,2,...,N, элементы которой представляют собой значения яркости в равноотстоящих точках пространственной области изображения.
В дискретном случае непрерывным переменным Х1, у1 и Х2, у2 в выражениях (8, 9, 10) соответствуют порядковые номера г1,к1 и г2,к2 элементов матрицы Ф. Функциям Аа(Х1, Х2) и Ав(у1, у2) (9, 10) соответствуют матрицы А = (аг 2) и Ав = (а\к2), размерно-
сти МхМ и NxN (называемые субполосными матрицами [2]), элементы которых задаются следующими соотношениями:
, (а + а )(/1 - г2 к (а - а Х!- ь)
а“,
2Соя 2 1/41 2/ Sin -
1: * ¿2,
п(11 - 12) (12)
¿1 ¿2,
а^ k:k2 —<
п
2Со^и2' "2; Sin
(в2 + Д )(k1 - ^ ) О, ™ (в2 - Д )(k1 - Ю
П№ - k2) (13)
в2 - Д
п
Данное представление элементов субполосных матриц позволяет разработать быструю процедуру приближенного вычисления значений элементов матриц.
Таким образом, для определения точного значения доли энергии изображения Ф — (fiк) в частотном субинтервале О (6), используя выражения (12, 13), можем указать следующее соотношение:
Гг(Аат Ф- Ав Фт)
Ро — ---- , (14)
О Гг(ФФ т )
где Ф - матрица исходного изображения, Аа и Ар - субполосные матрицы, значения
элементов которых определяются на основании выражения (12, 13), функция «£г» - след матрицы.
Представление (14) определяет метод анализа распределения энергий изображений по заданным частотным интервалам на основе частотных представлений и позволяет разработать вычислительный алгоритм нахождения доли энергии Ро дискретного
изображения Ф в частотном двумерном субинтервале О, не вычисляя при этом трансформанту Фурье.
Выводы
Возможность вычисления точных значений указанных долей энергий без перехода в пространство трансформант Фурье определяет основное преимущество предложенного метода перед используемыми в настоящее время методами, основанными на вычислении ДПФ. Значимость такой возможности определяется тем, что эти доли являются одними из важнейших физических характеристик, и интерес к ним проявляется почти всегда, когда для анализа изображений используются частотные представления.
С целью проверки работоспособности разработанного метода был проведен ряд вычислительных экспериментов, результаты которых показали, что данный метод анализа распределений энергий изображений по заданным частотным интервалам является инвариантным к повороту изображений на 90° и к изменениям яркости; данный факт позволяет применять разработанный метод в качестве эффективного инструмента выявления частотных свойств различных изображений.
Метод анализа распределения энергий изображений по заданным частотным интервалам может быть использован при создании алгоритмов и программ, позволяющих решать проблему выделения квазициклических компонент, соответствующих частотным субинтервалам, в которых сосредоточена подавляющая доля энергии изображения, в задачах сжатия при хранении, передаче и повышении скрытности визуальной информации.
2
2
а2 - а
2
2
Литература
1. Гонсалес Р. Цифровая обработка изображений / Р. Гонсалес, Р. Вудс. - М.: Техносфера, 2006. - 1072 с.
2. Жиляков Е. Г. Методы анализа и построения функций по эмпирическим данным на основе частотных представлений / Е. Г. Жиляков. - Белгород: Изд-во БелГУ, 2007. - 160 с.
3. Хургин Я. И. Финитные функции в физике и технике / Я. И. Хургин, В. П. Яковлев. -М.: Либроком, 2010. - 414 с.
4. Сойфер В. А. Методы компьютерной обработки изображений / В. А. Сойфер. - М.: ФИЗМАТЛИТ, 2003. - 784 с.
5. Рабинер Л. Теория и применение цифровой обработки сигналов / Л. Рабинер, Б. Гоулд. -М.: Мир, 1978. - 327 с.
6. Жиляков Е. Г. Вариационные алгоритмы анализа и обработки изображений на основе частотных представлений / Е. Г. Жиляков, А.А. Черноморец. - Белгород: Изд-во «ООО ГиК», 2009. - 146 с.
7. Жиляков Е.Г. Об эффективности метода оценивания значений долей энергии изображений на основе частотных представлений / Е.Г. Жиляков, А.А. Черноморец, А.Н. Заливин // Известия ОрелГТУ. Информационные системы и технологии. - 2009. - № 2/52 (563). С. 12-22.
METHOD OF ANALYSIS OF IMAGE ENERGY DISTRIBUTION IN SPECIFIED FREQUENCY INTERVALS
A.A. CHERNOMORETS O.N. IVANOV
In this work we propose a method of analysis of image energy distribution in specified frequency intervals based on calculation of exact values of corresponding energy parts.
Belgorod State University
Key words: image, frequency interval, energy parts, subband matrices.
e-mail :
chernomorets@bsu.edu.ru