CC BY
25RESEARCH RESULT
Zhilyakov E.G., Chernomorets A.A., Bolgova E.V., Goloschapova V.A. On the method of optimal subband filtration // Сетевой журнал «Научный результат». Серия «Информационные технологии». - Т.1, №1,2016.
УДК 621.397
DOI 10.18413/2518-1092-2016-1-1-58-64
Жиляков Е.Г.
О МЕТОДЕ ОПТИМАЛЬНОЙ СУБПОЛОСНОЙ ФИЛЬТРАЦИИ
Черноморец А. А. Болгова Е.В. Голощапова В.А.
1) заведующий кафедрой информационно-телекоммуникационных систем и технологий, доктор технических наук, профессор. Белгородский государственный национальный исследовательский университет, ул. Победы д.85,
г. Белгород, 308015, Россия. e-mail: Zhilyakov@bsu.edu.ru 2) профессор кафедры прикладной информатики и информационных технологий, кандидат технических наук, доцент Белгородский государственный национальный исследовательский университет, ул. Победы д.85, г. Белгород, 308015,
Россия. e-mail: Chernomorets@bsu.edu.ru 3) аспирант кафедры прикладной информатики и информационных технологий. Белгородский государственный национальный исследовательский университет, ул. Победы д.85, г. Белгород, 308015, Россия.
e-mail: Bolgova_e@bsu.edu.ru 4) старший преподаватель кафедры информационно-телекоммуникационных систем и технологий Белгородский государственный национальный исследовательский университет, ул. Победы д.85, г. Белгород,
308015, Россия. e-mail: vgoloschapova@bsu.edu.ru
Аннотация
В данной работе изложены результаты экспериментальных исследований эффективности применения метода оптимальной субполосной фильтрации, обеспечивающего минимальный уровень просачивания энергии за пределы выбранной для фильтрации подобласти пространственных частот.
В ходе экспериментов показано, что энергия входного изображения, сосредоточенная вне выбранной подобласти пространственных частот, влияет на выходную частотных фильтров, основанных на применении ДПФ, в значительно большей степени, чем на результат, полученный на основании метода оптимальной фильтрации. Данное влияние проявляется как на графиках значений соответствующих трансформант Фурье, так и на изображениях, полученных в результате фильтрации.
Ключевые слова: изображение; оптимальная фильтрация; трансформанта Фурье; подобласть пространственных частот.
Chernomorets A.A. Bolgova E.V. Goloschapova V.A.
1) Doctor of Technical Sciences, Professor, Head of Department of Information and Telecommunication Systems and
Technologies, Belgorod State National Research University, 85 Pobedy St., Belgorod, 308015, Russia
e-mail: Zhilyakov@bsu.edu.ru
2) Candidate of Technical Sciences, Associate Professor. Department of Applied Informatics and Information Technologies,
Belgorod State National Research University, 85 Pobedy St., Belgorod, 308015, Russia e-mail: Chernomorets@bsu.edu.ru 3) Postgraduate Student. Department of Applied Informatics and Information Technologies. Belgorod State National Research University, 85 Pobedy St., Belgorod, 308015, Russia. e-mail: Bolgova_e@bsu.edu.ru 4) Senior Lecturer. Department of Information and Telecommunication Systems and Technologies Belgorod State National Research University, 85 Pobeda St., Belgorod, 308015, Russia. e-mail: vgoloschapova@bsu.edu.ru
Abstract
The paper presents the results of experimental studies of the effectiveness of the optimal subband method. This method allows to achieve the minimal level of energy leakage out of the subdomain of spatial frequencies selected for filtration.
Zhilyakov E.G.
ON THE METHOD OF OPTIMAL SUBBAND FILTRATION
Zhilyakov E.G., Chernomorets A.A., Bolgova E.V., Goloschapova V.A. On the method of optimal subband filtration // Сетевой журнал «Научный результат». Серия «Информационные технологии». - Т.1, №1,2016.
The experiments show that the energy of the input image, being out of the selected subdomain of spatial frequencies, affects in a greater degree the output sequence of frequency filters based on DFT than the results based on the optimal filtration method. This influence is evident both on the graphs of corresponding Fourier transforms, and on the images obtained as a result of filtering. Keywords: image; optimal filtration; Furie transform; subdomain of spatial frequencies.
Введение
Одно из направлений интенсивного развития современных информационно-
телекоммуникационных систем определяется тенденцией использования наиболее
естественных для человека форм информационного обмена, и, прежде всего, в виде визуальных данных. Обработка изображений является важной задачей. Существующие в настоящее время методы выделения (фильтрации) значимых характеристик отдельных компонент изображения, некоторых периодических структур изображений не являются оптимальными с точки зрения аппроксимации трасформант Фурье в заданной подобласти пространственных частот (11114), в которых осуществляется фильтрации. Поэтому актуальной проблемой является создание математических моделей и методов фильтрации, позволяющих адекватно учитывать энергетические характеристики изображений в выбранной подобласти пространственных частот. В работе [1] разработан и теоретически обоснован новый метод оптимальной линейной фильтрации изображений на основе частотных представлений, который является оптимальным в том смысле, что спектр получаемого в результате фильтрации изображения имеет наименьшее
среднеквадратическое отклонение от спектра фильтруемого изображения в заданной двумерной подобласти пространственных частот, а вне этой подобласти пространственных частот имеет наименьшее отклонение от нуля. В работе представлены результаты экспериментальных исследований эффективности метода
оптимальной фильтрации изображений на основе частотных представлений в заданной подобласти пространственных частот на примере модельных изображений.
Теоретические основы метода
оптимальной фильтрации изображений на основе частотных представлений
О частотных представлениях при обработке изображений [2] можно говорить, поскольку изображение, как функция с конечной или
неограниченной областью определения, может быть представлено в виде суммы синусов и косинусов различных частот, умноженных на некоторые весовые коэффициенты
1 М N
/к =~Т2II Р(и,1е™ 1,I = 1,2,..М,
и= 1 Г=1
к= 1,2,...,N (1)
где Р(и, V) - значение трансформанты Фурье функции
F (u, v) = 11 fü
-ju(i-1) -jv(k-1) T1 _
e
, u = 1,2,...,M,
=1 k=1
v= 1,2,., N. (2)
Метод оптимальной фильтрации
изображений на основе частотных представлений [1,3] позволяет для нахождения результатов УО фильтрации изображения Ф=(/к), г=1,2,...,М, к=1,2,...,^ в частотной двумерной области О построить вычислительную процедуру, не вычисляя при этом трансформанту Фурье. Для фильтрации изображений предложено использовать следующий вариационный принцип: трансформанта Фурье Z(u,v) результата фильтрации УО, является оптимальной в смысле евклидовой нормы ее отклонения в заданной подобласти пространственных частот О от трансформанты Фурье F(u,v) исходного изображения Фи от нуля - вне данной 11114, т.е.
jj \F(u, v) - Z(u, v)| dudv +
(m,V)gQ
+ jj |Z (u, v)|2 dudv ^ min. (3)
(u,v
Указанному вариационному принципу соответствует следующее изображение Yq
Yn= AT-Ф-B, (4)
где матрицы А=(аШ2), ii,i2=l,2,...,M, и B=(bkm), ki,k2=l,2,...,N (субполосные матрицы [2]), определяются соотношениями
Zhilyakov E.G., Chernomorets A.A., Bolgova E.V., Goloschapova V.A. On the method of optimal subband filtration // Сетевой журнал «Научный результат». Серия «Информационные технологии». - Т.1, №1,2016.
аи =
¥2
Sin(а2(ij -i2)) -Sin«(ij -i2)) .
, ij
x(ij - ¿2)
kjk-2
Sin(ß2(kj -^2))- Sin(ßj(kj -k2)) k , ,
, kj ^ k2,
ß2 -ßj
7r(kl - k2) , kj = k2.
(5)
Значения
0<а1,а2,Д,р2 <ж (6)
задают границы подобласти пространственных частот О.
Важным свойством определенного таким образом результата фильтрации является то, что указанный метод не допускает растекания энергии двумерного сигнала за пределы заданной частотной области.
Концептуальные основы
экспериментальных исследований
В ходе вычислительных экспериментов проводился сравнительный анализ оценки доли энергии изображения У0р1, полученного в результате оптимальной фильтрации,
SOOpt
Opt\\
- Pa (Yopt )
l|2
= j -
Pa (Yop, )
(7)
Opt\\
Optn
где ||YOpJ -энергия изображения YOvt в области
Opt
определения трансформанты Фурье,
Ро (У0р,) - доля энергии изображения У0р{ в
подобласти пространственных частот О, определяемая на основании соотношения
Ро (Уор,) = *ес(К?ор,РХр;)- (8)
с соответствующими характеристиками изображений, полученных при реализации известных частотных фильтров [4]:
- идеальный фильтр (результат фильтрации
Ideal
il, u е[-ж,4 v e[-ß2,-ßj]Ußß],
HIdeal(u,v) = Г — (9)
[О, в протиен ом случае.
- фильтр Баттерворта порядка n=2 (результат
фильтрации YButterw X
1
HBu„rrw(U V) = ■
( Г> С,, т,\Л2п ( n с,, „W
(j0)
j +
Dj(u, v)
D,
+
0 J
D2 (u,v)
D
0 J
где О0 - частота среза, Д(ы,V), 02(ы,V) -расстояние между точкой (ы, V) и осями симметрии I, 12 фильтра,
- Гауссов фильтр (результат фильтрации
Gaussian
^GausiaXU v) = ®Ф
Dj2(u, v) ^
2D,
+ exp
0 J
D22(u, v) ^
2D
, (Л)
0J
где Б0 - частота среза, Д(ы,V), 02(ы,V) -расстояния между точкой (ы, V) и осями симметрии ^, 12 фильтра.
Построение фильтров осуществлено на основании параметров подобласти
пространственных частот, в которой выполнялась оптимальная фильтрация [5,6].
Соответствующие оценки доли энергии результатов фильтрации (изображения УЫеа1,
УЕи„ег„, YGaыSSian), которая сосредоточена вне выбранной подобласти пространственных частот О, имеют следующий вид
S = ! -
aideal j
Pa (Yideal)
Y
idealU
S = j -
OButterw
Pa (YButterw )
Bu tterw
S = j -
OGuassian
P (Y )
О v Guassian
(j2) (j3) (j4)
У II
|| Gыassiaц\
Отметим, что вычисления как относительной доли «просачивания» энергии за пределы подобласти пространственных частот оказались доступными, благодаря разработанному в [1] методу вычисления долей энергии (8).
Для получения большей наглядности при оценивании различий оптимальной фильтрации и фильтрации с помощью известных фильтров представляется естественным построить графики соответствующих характеристик выходных последовательностей фильтров.
Сравнительный анализ в работе выполнен, исследуя параметры фильтров и результаты их применения в различных подобластях двумерной области (подобласти получены при фиксированных значениях одной из независимых переменных частотного пространства). Однако, данный факт не снижает важности полученных результатов в виду известного свойства преобразования Фурье (преобразование Фурье
2
а^-а
j2
к
л
2
2
2
2
2
2
Zhilyakov E.G., Chernomorets A.A., Bolgova E.V., Goloschapova V.A On the method of optimal subband filtration // Сетевой журнал «Научный результат». Серия «Информационные технологии». - Т.1, №1,2016.
можно записать так, что переменные в нем оказываются разделенными).
Результаты вычислительных экспериментов
Экспериментальные исследования
целесообразно осуществлять на основе обработки
модельных сигналов. Генерирование значений модельного изображения (рис. 1) Ф = (/к), г = 1,2,...,512, k = 1,2,...,512, размерностью 512x512 пикселей, осуществляется на основе соотношения
где
fk = a1 sin(2nf1k) + о> sin(2nf2k) + a3 sin(2nf3k) + a4 sin(2^f4k) + a5 sin(2^f5k), (15) i = 1,2,...,512, k = 1,2,...,512,
f1 = 0.0193; f2 = 0.0208; f3 = 0.0279, f4 = 0.0295, f5 = 0.0311, a = 3.2; a2 = 2.9; a3 = 3.5, a4 = 2.1, a5 = 2.2.
Такие значения частот позволяют осуществить исследование различных аспектов фильтрации и, прежде всего, влияний энергии исходного сигнала за пределами подобласти пространственных частот на ее результаты.
Фильтрация модельного сигнала (15) на основании метода оптимальной фильтрации изображений на основе частотных представлений
осуществлялась в следующей подобласти пространственных частот:
О,: ^ = 0, а2 = л, & = &+Д&(г-1), & = &+¥(»-1)}, (16)
г=1,2,.,15,
где Д& - величина смещения границ очередной
ппч.
Рис. 1. Модельное изображение: а - в виде изображения, б - профиль Fig. 1. Model image: a - in the form of images, b - profile
В ходе экспериментов были использованы следующие значения границ подобласти пространственных частот
&01 = 0.035 л, &02 = 0.04я; Д& = 0.0025^. (17)
В таблице представлены значения долей энергии выходных последовательностей оптимального и других анализируемых фильтров, находящихся вне заданной подобласти пространственных частот, вычисленные согласно (7), (12)-(14).
Как показано в ходе экспериментов выходная последовательность оптимального фильтра
определяется только трансформантой Фурье в выбранной подобласти пространственных частот. В отличие от этого, на выходную последовательность идеального, Гауссова фильтров и фильтра Баттерворта влияет энергия входного изображения, сосредоточенная вне выбранной подобласти пространственных частот. Это различие легко обнаружить на рис. 2, на котором представлены абсолютные значения трансформант Фурье выходных последовательностей оптимального и других фильтров.
Zhilyakov E.G., Chernomorets A.A., Bolgova E.V., Goloschapova V.A On the method of optimal subband filtration // Сетевой журнал «Научный результат». Серия «Информационные технологии». - Т.1, №1,2016.
Оценка долей энергии результатов фильтрации вне подобласти пространственных частот Q Estimation of the energy share of filtered results outside the subdomain of spatial frequencies Q
Таблица 1
Table 1
№ Границы ППЧ Q àQOpt àQldeal С QButterw с- QGaussian
A A
1 0.035л 0.04л 0.13211 0.28323 0.27313 0.27189
2 0.0375л 0.0425л 0.3492 0.30822 0.31119 0.31036
3 0.04л 0.045л 0.13101 0.55406 0.55503 0.55472
4 0.0425л 0.0475л 0.34099 0.72638 0.55996 0.60419
5 0.045 л 0.05л 0.41356 0.86981 0.87022 0.87001
6 0.0475л 0.0525л 0.23196 0.56816 0.57263 0.57082
7 0.05л 0.055л 0.34246 0.58755 0.6305 0.62484
8 0.0525л 0.0575л 0.12423 0.22294 0.22436 0.22396
9 0.055л 0.06л 0.33214 0.63152 0.6312 0.63129
10 0.0575л 0.0625л 0.40426 0.85964 0.7762 0.8095
11 0.06л 0.065л 0.1968 0.41986 0.4204 0.42016
12 0.0625л 0.0675л 0.24863 0.57319 0.57095 0.57158
13 0.065 л 0.07л 0.4328 0.92163 0.92157 0.92155
14 0.0675л 0.0725л 0.47163 0.69156 0.69313 0.69208
15 0.07л 0.075л 0.20767 0.7969 0.68474 0.7815
Далее представлены результаты фильтрации модельного изображения (рис. 1) в подобласти пространственных частот
О :{а1 = 0, а2 =л, & = 0.04л, & = 0.045л}.
Рис. 2. Трансформанты Фурье (абсолютные значения) в подобласти пространственных частот а,1=0, а2=я, & = 0.04л, /3, = 0.045л при и=0.22п: - исходного изображения; выходных последовательностей фильтров: —Ф-Ф— „ ___*_ „ + ■ ^ ■■+......ф-
---------------- » _ лдеадьныл
~ - идеальный, ' ' ~ ' - Батгерворта, т Гауссов
Fig. 2. Fourier transforms (absolute values) in a subdomain of spatial frequencies a1=0, а2=л, Д = 0.04л, ¡32 = 0.045л at
-Ф-Ф-
- оптимальный.
и=0.22л:
- original image; output sequences of filters: Butterworth, + +
- optimal,
- ideal,
- Gaussian
Zhilyakov E.G., Chernomorets A.A., Bolgova E.V., Goloschapova V.A On the method of optimal subband filtration // Сетевой журнал «Научный результат». Серия «Информационные технологии». - Т.1, №1,2016.
На рис. 2 приведены абсолютные значения трансформант Фурье при и=0.22л, ve [0.03 л, 0.075 л] исходного изображения и результатов фильтрации, полученных при применении оптимального и других фильтров в заданной подобласти пространственных частот О.
На рис. 3 влияние дополнительной энергии исходного изображения, сосредоточенной вне
выбранной подобласти пространственных частот у ее левой границы (рис. 2), проявляется в наличии существенных изменений яркости пикселей изображений (рис. 3в, 3д, 3ж), являющихся результатами фильтрации с помощью фильтров, отличных от оптимального фильтра, и в наличии биений на графиках соответствующих профилей изображений.
SO 100 130
Рис. 3. Результаты фильтрации (изображение и профиль) на основании фильтров: а, б - оптимальный; в, г -
идеальный; д, е - Батгерворта, ж, з - Гауссов
Fig. 3. The filter results (image and profile) based on the filters: a, b - optimal; c, d - ideal; e, f Butterworth; g, h - Gaussian
Выводы
Результаты вычислительных экспериментов как с модельными изображениями, так и с реальными фотоснимками в цифровом виде
показывают, что относительные доли «просачивания» энергии за пределы указанной подобласти пространственных частот у выходных последовательностей частотных фильтров,
Zhilyakov E.G., Chernomorets A.A., Bolgova E.V., Goloschapova V.A. On the method of optimal subband filtration // Сетевой журнал «Научный результат». Серия «Информационные технологии». - Т.1, №1,2016.
основанных на применении ДПФ, существенно больше, чем у разработанного метода.
В ходе экспериментов показано, что энергия входного изображения, сосредоточенная вне выбранной подобласти пространственных частот, влияет на выходную частотных фильтров, основанных на применении ДПФ, в значительно большей степени, чем на результат, полученный на основании метода оптимальной фильтрации. Данное влияние проявляется как на графиках значений соответствующих трансформант Фурье, так и на изображениях, полученных в результате фильтрации.
Работа выполнена при поддержке гранта РФФИ 15-07-01570 (Субполосная скрытная интеграция/извлечение дополнительной
информации в аудио или видео контентах).
Список литературы
1. Жиляков Е.Г., Черноморец A.A. Вариационные алгоритмы анализа и обработки изображений на основе частотных представлений. Белгород: Издательство БелГУ, 2008. 146 с.
2. Жиляков Е.Г. Методы анализа и построения функций по эмпирическим данным на основе частотных представлений. Белгород: Издательство БелГУ, 2007. 160 с.
3. Жиляков Е.Г., Черноморец A.A. Optimal filtering images based on frequency representations // Вопросы радиоэлектроники. Сер. ЭВТ, 2008. Вып. 1. С. 118-131.
4. Гонсалес Р., Вудс Р. Цифровая обработка изображений. М.: Техносфера, 2006. 1072 с.
5. Жиляков Е.Г., Черноморец A.A., Болгова Е.В. О субполосной фильтрации изображений с применением компенсирующих коэффициентов // Техника и технологии - 2014: сборник статей Международной научно-технической конференции. 2014. С. 30-34.
6. Черноморец A.A., Лысенко ИВ., Болгова Е.В. Компьютерная реализация алгоритма взвешенной оптимальной фильтрации изображений. Вопросы радиоэлектроники. Сер. ЭВТ, 2012. Вып. 1. С. 103-111.
References
1. Zhilyakov E.G., Chernomorets A.A. Variation Analysis Algorithms and Image Processing Based on the Frequency Representations. Belgorod: BSU publishing house, 2008. 146 p.
2. Zhilyakov E.G. Methods of Analysis and Construction of Functions from the Empirical data Based on Frequency Representations. Belgorod: BSU Publishing House, 2007. 160 p.
3. Zhilyakov E.G., Chernomorets A.A. Optimal Filtering Images Based on Frequency Representations. Electronic Engineering. 2008. Part 1. Pp. 118-131.
4. Gonzalez R., Woods R. Digital Image Processing. Moscow: Tekhnosfera, 2006. 1072 p.
5. Zhilyakov E.G., Chernomorets A.A., Bolgova E.V. About Subband Filtering Images Based on the Use of Compensating Coefficients. Engineering and Technologies - 2014: collected works of the International Scientific and Technical Conference. 2014. Pp. 30-34.
6. Chernomorets A.A., Lysenko I.V., Bolgova E.V. Computer Implementation of Algorithm of the Weighed Optimal Image Filtration. Electronic Engineering. 2012. Part 1. Pp. 103-111.
