МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ И ЧИСЛЕННЫЕ МЕТОДЫ
УДК 336.153
ПРИМЕНЕНИЕ ВАРИАЦИОННЫХ АЛГОРИТМОВ ИНТЕРПОЛЯЦИИ И ОЦЕНКИ ПЕРВОЙ ПРОИЗВОДНОЙ ДЛЯ НЕКОТОРЫХ АСПЕКТОВ ОБРАБОТКИ ИЗОБРАЖЕНИЙ
Т.Н. СОЗОНОВА Н.С. ТИТОВА21 Н.В. ЩЕРБИНИНА31
1) Белгородский государственный университет
e-mail: Sozonova@bsu.edu.ru
2) Белгородский государственный университет
e-mail: NTitova@bsu.edu.ru
2) Белгородский государственный университет
e-mail:Shcherbinina @bsu. edu. ru
В настоящее время в различных областях науки и техники широкое развитие получили системы, ориентированные на получение, обработку и хранение информации, представленной в виде изображений. В данной работе рассмотрены некоторые аспекты обработки изображений, в частности, передискретизация и увеличение четкости, на основе вариационных алгоритмов интерполяции и оценивание первой производной на основе частотных представлений.
Ключевые слова: интерполяция, оценка производной, частотное представление, передискретизация, четкость изображения.
Многие отрасли науки и техники, имеющие отношение к получению, обработке, хранению и передаче информации, в значительной степени ориентируются в настоящее время на развитие систем, в которых информация имеет характер изображений. Изображение, которое рассматривается как двумерный сигнал, является значительно более емким носителем информации, чем обычный одномерный (временной) сигнал. Вместе с тем, при решении многих научных и инженерных задач при работе с визуальными данными требуется в качестве промежуточного этапа проводить пересчет цифровых изображений на другую растровую сетку. В качестве примера таких задач можно привести коррекцию погрешностей аппаратуры (апертурные искажения камеры), перевод аэрокосмических изображений в географические проекции, приведение нескольких изображений к единому размеру, устранение перспективных искажений.
Одной из задач обработки изображений является передискретизация. К задаче передискретизации предъявляются высокие требования по точности, так как каждое преобразование ведет к накоплению ошибки. Таким образом, наилучшим является алгоритм, который производит набор необходимых преобразований за один этап и использует наиболее точный для данной задачи метод аппроксимации. В настоящее время во многих редакторах обработки изображений используются в основном билинейная и бикубическая интерполяция. Отметим, что интерполяция изображения может привести к сильному размытию последнего.
Еще одной из задач обработки изображений является увеличение четкости. Недостаточная четкость изображения может иметь множество объяснений: недостаточное разрешение используемой фототехники, увеличение или уменьшение размера изображения и др.
Таким образом, возникает необходимость при помощи обработки придать изображению такие качества, благодаря которым его восприятие человеком было бы по возможности комфортным. Часто бывает полезным подчеркнуть, усилить какие-то черты, особенности, нюансы наблюдаемой картины с целью улучшения ее субъективного восприятия.
Последнее - субъективность восприятия - сильно усложняет применение формализованного подхода в достижении данных целей. Поэтому, при обработке изображений для визуализации получили распространение методы, в которых часто отсутствуют строгие математические критерии оптимальности. Их заменяют качественные представления о целесообразности той или иной обработки, опирающиеся на субъективные оценки результатов.
В настоящее время известно много методов увеличения четкости изображений, например, такие как: шарпинг, линейное контрастирование, соляризация, препарирование, эквализация и т.д.
В данной работе предлагается использование новых вариационных алгоритмов интерполяции (для увеличения размера изображения) и оценивания производных (для увеличения четкости изображения), основанных на использовании класса целых аппроксимирующих функций ограниченной степени роста. К ним, в частности, относятся функции с финитной областью определения трансформант Фурье, то есть с финитными спектрами.
Учитывая, что изображение можно рассматривать как двумерный сигнал, интерполяция изображения может производиться в два этапа: вдоль оси ОХ и вдоль оси ОУ.
------------------------>
У:
Рис. 1. Исходное изображение
Рассмотрим математические основы предлагаемых алгоритмов интерполяции и оценивания производных по одному направлению (вдоль оси ОХ).
Пусть в эквидистантных точках ti = г = 0,...,N дан вектор значений сигнала
и = ^..^UN)Т,
и1 = u(iAt), г = 0,..., N,
где символ Т означает транспонирование.
Представляется целесообразным областью определения аппроксимирующей функции считать всю числовую ось, т.е.
— ж < t < ж, (2)
При этом потребуем, чтобы в любой точке области определения существовали и были непрерывными производные любого порядка, причем выполнялись неравенства
| \и('к)^)| Л <<х>, к = 0,1,...
—ж
В основе дальнейших построений используется представление
t
) = и 0 +| / (т)Лт,
(3)
(4)
которое позволяет по производной вычислить интерполирующую функцию. Очевидно, что при этом должны выполняться интерполяционные равенства
и = и(/А^) = и, / = 0,1,..., N. (5)
Для оценки производной будем использовать класса моделей целых функций, которые представимы в виде
I(О = | F(ю)Лю, (6)
2ж „
оєО
где
0 = [-02,-01 )^[О1,О2); 01 < ж; О2 <ж
^ (ю) - трансформанта Фурье, связанная с модельной функцией дуальным (7) соотношением
Е(о) = |/(т)е }отёт (7)
Выбор области определения О трансформанты Фурье может быть продиктован априорными сведениями о свойствах сигнала.
Подстановка представления (7) в правую часть (6) позволяет получить соотношение для вычисления интерполирующей функции на основе трансформанты Фурье производной
' Ю Л
2П оЄо °2
і 81ПI 2 і о
и(ї) = и0 Н-----------------Г Е(о)------------------т— е 2 йо, (8)
9 7Т 1 ГІЇ /
так что условиям интерполяции (5) нетрудно придать вид
. ( (оАИ Л
1 г 81П V 2 і —і у /
^ IЕ(о)4лИе 2 <*» = %, (9)
2ж ^ оАї
ОЄО /2
где
V = (иі -и0),і = !,•••,N. (10)
Для отбора интерполирующей функции из класса целых предлагается использовать вариационный принцип минимизации евклидовой нормы производной. На основе представления (7) и равенства Планшереля (Парсеваля), этому принципу нетрудно придать вид
+ж і
I/2(т)йт =— I Е (о)2 йо = ШІП (11)
-ж оєО
Таким образом, задача сводится к поиску решения вариационной изоперимет-рической задачи, определяемой условием (11) и ограничениями вида (9). В работе показано, что искомое решение представимо в виде
. ( оАїі ^
^ 81П і 2 і -—і Е(о)-2е 2 • (12) г /2
когда о є О и Е(о) = 0нулю в противном случае.
Для вычисления вектора множителей Лагранжа ¡3 = 3,--,3Ы)т следует воспользоваться подстановкой представления (12) в левые части равенств (9). В результате нетрудно получить систему линейных алгебраических уравнений (СЛАУ), которым должен удовлетворять набор коэффиц иентов
А0 = V ^у^-,уы)т , (13)
где
А = К};
а,„ ак = —
-
-1
7Г А
хк
2
О281П| ----- I 81П|
XI
-008
2(к -О
ёх; О = А/О , г = 1,2. (14)
В строгом смысле симметричная матрица с элементами вида (14) является положительно определенной. Отметим, что если интервал интегрирования в (14) удовлетворяет условию
О2 - О1 = NАt(О 2 - О1) > 2-,
(15)
то определитель матрицы А будет заметно больше нуля, следовательно решение СЛАУ вида (13) будет устойчивым. Отметим, что неравенство вида (15) всегда можно выполнить за счёт соответствующего определения размеров области определения трансформанты Фурье моделирующей (аппроксимирующей производную) целой функции. Иными словами, в указанном смысле решение сформулированных задач оценивания производной и интерполяции на основе соотношения (8) всегда существует.
Если же интервал интегрирования удовлетворяет условию
N (О 2 - О1) < -,
(16)
то значения элементов матрицы А будут близки и ее определитель будет близок к нулю, следовательно решение СЛАУ вида (13) будет неустойчивым. В этом случае при вычисления вектора ¡3 предлагается использовать псевдообратную матрицу А+, вычисляемую через ненулевые собственные числа и соответствующие им собственные векторы
А+ = бА1^, (17)
где
¿1 = (А-.АX Л> А2 >... > А > 0 01 = Чз).
На основе соотношений (12) и (14) нетрудно получить вычислительные формулы для интерполирующей функции (18) и оценки производной (19).
. (ют
„ О 81ПI
N °2 I 2
{от/ |юА^ )
008
О{т -)
(18)
О 81ПI
1 N О2
I (т)=1 X з, Н
- ,■=1 К 0
юАt
о А//
008
ю
Т - гАУ2 )]ю
(19)
/2
Отметим еще одну возможность организации вычислений. Если заранее известен набор тк, к = 1,...,Ь интервала [0,NАt], в которых предполагается в дальнейшем
вычисление интерполирующих функций. Тогда, используя (18) и (19) с учетом СЛАУ (13) можно получить следующие вычислительные формулы:
Для интерполирующей функции
где
и = (и(т1),...,и(ть))т = и0еь + СА+V ,
= (1,...,1)Т
с = К-};
(20)
81П
^от,г Л . (оАt
81П| -------/
I 2
-
О1
от,
{°А/2)
008
00 {тк-/А/)
(21)
2
О
2
Т
к
к ~
Для оценки производной
где
I = (I(т1),...,I(ть))Т = ВА+У , (22)
В = {¿к,-}; к = 1,..., Ь; / = 1,..., N.
юА/ .Л О 81ПI -----/ I
^ = - 1 юд2/ С08° Т - /А/21° . (23)
°1 /2
Удобство применения этих формул заключается в том, что матрицы А, В и С могут быть вычислены заранее, храниться в памяти ПК и многократно использоваться.
Рассмотрим применение предложенных алгоритмов интерполяции и оценки производной для обработки изображений. Пусть в качестве исходных данных имеем изображение размером МхМ (y=N, х=М) пикселей.
1. Увеличение изображений в К раз.
Для увеличения изображения будем использовать предложенный вариационный алгоритм интерполяции.
На первом этапе осуществляем вычисление матриц Ах = {а/},
/ = (1: N -1), ] = (1: N-1) , с элементами вида (14) и Сх ={е/] },
/ = (1: (N -1)К), ] = (1: N -1) с элементами вида (21).
Затем осуществляем интерполяцию исходного изображения предложенным вариационным методом вдоль оси ОХ по выражению
I = IСхАхх+;
где I - изображение, интерполированное вдоль оси ОХ, f — исходное изображение.
В результате получаем изображение, интерполированное вдоль оси ОХ, приведенное на рисунке 2.
Рис. 2. Интерполяция изображения вдоль оси ОХ
На втором этапе эксперимента осуществляем интерполяцию изображения, приведенного на рисунке 2 вдоль оси ОУ, для чего предварительно производим вычисление матриц Ау = {а, }, / = (1: М -1),у = (1: М -1) с элементами вида (14) и
Су = {су }, / = (1: (М -1)К), у = (1: М -1) с элементами вида (21). Интерполяцию осуществляем по выражению
I = СуАу+/х,
где I - интерполированное изображение.
В результате получим изображение размером 5ЫхбМ пикселей. Результаты эксперимента представлены на рисунке 3.
Рис. 3. Интерполированное изображение
Из приведенных данных видно, что предложенный алгоритм интерполяции может быть применен для увеличения размера изображения с приемлемой точностью.
2. Увеличение четкости изображений.
Для увеличения четкости изображений предлагается использовать предложенный метод оценки производных.
На первом этапе эксперимента производилось вычисление матриц Вх = {¿и }; к = 1,..., N; / = 1,..., N и Ву ={ЬЫ }; к = 1,..., М; / = 1,...,М с элементами вида (23) и осуществлялось вычисление второй смешанной производной по выражению
5 2I
дхду
= Вх А-1 • I • Ву А-1,
где Ау = {агу} - матрица с элементами вида (14), f — исходное изображение.
Затем к исходному изображению добавлялось значение второй смешанной производной, то есть
д21
I = I +
дхду
что позволило получить более четкие, в смысле субъективного восприятия, изображения. Результаты эксперимента приведены на рисунках 4 — 6.
а б
Рис. 4. Исходное изображение (а), изображение после обработки (б)
После обработки изображения наблюдается увеличение четкости границ перехода от светлого тона к более темному.
а б
Рис. 5. Исходное изображение (а), изображение после обработки (б)
На данном рисунке наблюдается уменьшение «туманности» изображения, что дает возможность наблюдать мелкие детали изображения и увеличить субъективное восприятие изображения.
«ДО
а
б
Рис. 6. Исходное изображение (а), изображение после обработки (б)
По результатам эксперимента видно, что предлагаемый алгоритм увеличения четкости изображений позволяет повысить детальность наблюдаемой картины, дает возможность наблюдения мелких деталей на всех участках изображения. А так же позволяет получить более четкие границы отдельных объектов, то есть производить оконтуривание, что может быть полезно в геологии, при обработке аэроснимков земной поверхности.
Очевидным достоинством предлагаемого алгоритма увеличения четкости изображений является его предельная простота. Следует отметить быстродействие предлагаемого алгоритма, которое достигается за счет предварительного вычисления используемых матриц и хранения их в памяти компьютера. Вместе с тем, предлагаемый алгоритм обработки изображений приводит к очевидному субъективному улучшению визуального качества.
Литература
1. Ланцош, К. Практические методы прикладного анализа [Текст] : справ. рук. / К. Лан-цош; пер. с англ. М. З. Кайнера. - М. : Физматгиз, 1961. - 524 с.
2. Хургин, Я. И. Финитные функции в физике и технике [Текст] / Я.И. Хургин, В.П. Яковлев. - М. : Наука, 1971. - 408 с. : ил.
3. Жиляков, Е.Г. Вариационный метод оценивания производных и интерполяции сигналов по эмпирическим данным [Текст] / Е.Г. Жиляков, Т.Н. Созонова, И.Ю. Мисливец // Вестник Воронежского государственного университета, Серия: Системный анализ и информационные технологии. - Воронеж, 2006. - выпуск 2. - С.70-73.
VARIATIONAL ALGORITHMS OF THE INTERPOLATION AND THE ESTIMATION OF THE FIRSR DERIVED FOR IMAGE PROCESSING
T.N. SOZONOVA1 N.S. TITOVA2)
N.V. Shcherbinina3)
Nowadays in different areas of the science and technology broad development have got the systems, oriented on reception, processing and keeping the information, presented as the images. In the article some aspects of the image processing are considered, in particular, resampling and the sharpness increase, on the variational algorithms of interpolation and estimation of the first derived on the frequency presentations base.
1) Belgorod state university e-mail: Sozonova@bsu.edu.ru
2) Belgorod state university e-mail: NTitova@bsu.edu.ru
3) Belgorod state university
Key words: Interpolation, estimation of the first derived, frequency pres entation, resampling, image sharpness
e-mail: Shcherbinina@bsu. edu.ru