КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ
УДК 004.932
ПОВЫШЕНИЕ ЧЕТКОСТИ МАСШТАБИРОВАННЫХ ИЗОБРАЖЕНИЙ НА ОСНОВЕ ВАРИАЦИОННОГО МЕТОДА ЧИСЛЕННОГО ДИФФЕРЕНЦИРОВАНИЯ
Одной из задач повышения визуального качества изображений является увеличение его размеров. Такие вопросы обработки изображений достаточно часто возникают при решении различных задач в науке и технике. В качестве примера можно упомянуть увеличение размеров космоснимков при формировании карт земной поверхности, при наблюдении за различными объектами в тех случаях, когда технические средства не позволяют получить изображения достаточно большого размера.
В таких случаях для увеличения размеров изображения используются различные методы интерполяции, наиболее распространенным среди которых является метод интерполяции кубическими сплайнами.
Одной из проблем при увеличении размеров цифровых изображений является недостаточная четкость получаемых масштабированных снимков, что проявляется в смазанных контурах объектов изображения, размытии мелких деталей изображения и т.д. Повысить визуальное качество изображений, в том числе и масштабированных, можно, используя градиентную обработку, которая заключается в комлексировании изображения с его оценками производных. Для оценки производных, как правило, используется дифференцирование с использованием тех или иных операторов над разностями значений функций. Недостатком таких операторов является чувствительность к воздействиям так называемых шумов измерений, что приводит к неустойчиво-стям получаемых оценок производных. Особенно ярко данный эффект наблюдается при повышении четкости увеличенных в размерах изображений. Таким образом, возникает необходимость разработки иных методов интерполяции и численного дифференцирования дискретных двумерных сигналов, к которым можно отнести цифровые изображения.
Т.Н. БАЛАБАНОВА О.Н. ИВАНОВ
В работе изложен новый метод вычисления производных сигнала по его дискретным значениям, основанный на частотных представлениях, применение данного метода для градиентной обработки цифровых изображений с целью повышения их четкости.
Белгородский
государственный университет
e-mail:
sozonova@bsu.edu.ru
Ключевые слова: дифференцирование, четкость изображения, частотные представления, вариационный принцип.
В данной работе предлагается иной подход к интерполяции и численному дифференцированию цифровых изображений. В основе таких методов предлагается использовать принцип минимизации евклидовых норм оценок первых производных из класса функций с финитными областями трансформант Фурье, при дополнительных условиях совпадения соответствующих определённых интегралов (формула Ньютона- Лейбница) с разностями зарегистрированных значений исходной функции. Использование таких принципов позволяет получить устойчивые оценки производных, осуществлять интерполяцию сигналов и использовать разработанные методы для масштабирования и увеличения четкости изображений.
Предлагаемый метод основан на использовании известной из математического анализа формулы, позволяющей выразить дифференцируемую функцию через производную (обозначения очевидны):
I
) = и^0) + | /(x)dx,I > 1{).
и(1) = „ ^
ч
Понятно, что при известном начальном значении и известном способе вычисления производной искомая функция также может быть вычислена с любой заранее оговоренной точностью.
Пусть задан вектор и = (и0,и1з...,)т отсчётов дискретного сигнала, где и{ = u(iAt), г = 1,..., N, At - интервал дискретизации. Обозначим V = (у1з..., )т, где
V = иг -иг-1>г = 1,...>К- (!)
Введём частотные интервалы:
0 = (-О 2,-00 и [П1,О 2), (2)
о = [-02,-о1) и [Ц, о2). 01 = At * 01 = 02 = At * 02 = д2*п. (з)
В основе дальнейших построений используется представление интерполирующей функций через производную
t
й{1) = и-1 + | I (т^т (4)
(г-1) At
для At (г -1) < t < iAt.
Тогда для первых разностей исходных данных должно выполняться равенство:
iAt
V. = и - и
1I (т^т, (5)
4-1
I(т)- первая производная интерполирующей функции, которая является оценкой первой производной неизвестной функции и(1), выборка из которой обрабатывается.
Для достижения устойчивости оценки первой производной используем представление
I(т) = -1- |F(aycoтda, (6)
где F(®) - трансформанта Фурье
F(®) = 11(т)е-Jютdт .
-ОТ
В качестве области определения трансформанты Фурье предлагается использовать частотный интервал, в котором сосредоточена основная доля энергии сигнала.
Соотношение для интерполирующей функции на основе трансформанты Фурье производной можно получить путем подстановки представления (6) в правую часть (4).
и(г) = иг-1 + — |F(а)(ехр^Ш) - ехр(]ок1 (г - 1))^о / jо, (7)
так что интерполирующие равенства представимы в виде
. (оАг Л
Sinl
1 I о ] /
— Г F(о)—^^ ехр(7'оАг(г - 0,5)^о = V. (8)
1 -л- тА1./ / Аг.
2п Jn (Аг
ю&и /2
Ясно, что такие интерполирующие функции тоже относятся к классу целых. Вместе с тем имеется возможность использовать дополнительные ограничения.
Можно привести достаточно много аргументов использования вариационного принципа
ад \
Г / 2(т^т = — || F (о)2 dю = тт
2п п
-ад ю&и
Один из аргументов заключается в целесообразности построения функции с наименьшей в смысле евклидовой нормы производной скорости изменения значений.
Другим важным соображением может служить необходимость повышения устойчивости вычислений к воздействиям случайных ошибок измерений (регуляризация).
Искомое решение вариационной задачи (8),(7) представимо в виде
. ( оАг Л
N
F (о) = (0^/ ехР(- jоАt (г - 0,5)), (9)
г=\ /2 когда ое^ F (о) = 0 в противном случае.
Общую формулу для вычисления оценки производной получаем при подстановке по следнего представления в соотношение (6)
\ п.2 )
/ (т) = ^вк * - х(АТ - к + 0.5)^. (10)
к=\ пй\ (х2) Аг
Коэффициенты р должны удовлетворять системе уравнений (8), на основании чего получаем
Ар = V,
где А = {ак } - матрица учета исходных данных (УИД), элементы которой определяются из соотношения
А. п.^т2(х2)
агк = — I-у-^— cos(х(г - к))dx; г, к = 1,..., N. (и)
(х/2)2
В общем случае матрица УИД может быть особенной, так что для нахождения коэффициентов р необходимо использовать псевдообращение
в = А++V, (12)
А++ = ^Ц^т2 , (13)
где G - матрица собственных векторов.
AG = GL ; G = (д 1,..., д N ),
L = diag (Л,...Л);
4 = diag (Л1,...,Лр ), (14)
если
Лр+1 = Лр+2 = ... = Лм = 0, (15)
где Р - оценка ранга матрицы УИД.
G = (Ч,..., Чр ). (16)
Если заранее выбрать точки в виде
Т = (/ - 0.5)ДГ,i = 1,...,N (17)
области определения, где необходимо вычислять оценку производной то из (ю) получим
N
1 % )
/ = /(Т) = - Г-cos(- к. (18)
Или для вектора / = (/,..., /ы )Т, / = /(т,),
/ = В1А++ V, (19)
гдеВ! = [Ь)к },
1 1 ) Ьк = — I-со!(х(, - кШх. (20)
Ч (х/2)
Старшие производные в тех же точках вычисляются на основе дифференцирования выражения (ю)
^М = ц€-2)(т) = -Увк— Г^^/2.хsin(х(—- к + 0.5))dx. (21)
^ ^ пМ £ (^) М
В тех же точках области определения получим
1 °.2 !!П( Х2 )
в2 = Ь}: Ь2к =-—1-х !!п(х(, - к)^. (22)
А (*/2)
Вектор оценок вторых производных вычисляется на основе соотношения
/(1) = С/1(1),.., = В2 А++V = В2р. (23)
В рамках данной работы предлагается использовать новый вариационный метод оценки производных для увеличения четкости масштабированных изображений.
В качестве экспериментальных данных использовались изображения небольшого размера ЫхМ (К, М от 100 до 200 пикселей), содержащие оттенки серого, полученные при помощи цифровой техники. Выбор изображений, содержащих оттенки серого, обусловлен простотой реализации алгоритма для данного вида изображений. При повышении четкости цветных изображений обработке подвергается каждая из отдельных изображений-компонент (КОВ).
На первом этапе эксперимента проводилось увеличение размера изображения при помощи вариационного алгоритма интерполяции в К=2, 3, 5 раз [1]. В результате чего изображение получало большие размеры, однако контуры объектов на масштабированном изображении оказывались нечеткими.
На втором этапе эксперимента осуществлялось повышение четкости масштабированного изображения:
Осуществлялось вычисление матриц Вх = [Ьш }; к = 1,...,N; I = 1,...,N и
Ву = {ь. }; к = 1,...,М; I = 1,...,М с элементами вида (20) и осуществлялось вычисление второй смешанной производной по выражению
дхду
= ВхА - • / • В А -
где Ау = \ау } - матрица с элементами вида (и)/- исходное изображение.
Затем к исходному изображению добавлялось значение второй смешанной производной, то есть
€=у
дхду
что позволило получить более четкие, в смысле субъективного восприятия, изображения. Результаты эксперимента представлены на рисунках 1-2.
Рис. 1. Исходное изображение
Рис. 2. Масштабированное изображение (слева) и изображение после градиентной обработки (справа)
Рис. 4. Масштабированное изображение
Рис. 5. Масштабированное изображение с повышенной четкостью
По результатам эксперимента видно, что предлагаемый метод повышения визуального качества масштабированных изображений, в частности повышения четкости цифровых изображений, позволяет получить изображения, в которых более четко наблюдаются границы перехода от одного объекта к другому. Так же следует отметить,
что повышение четкости масштабированных изображений позволяет увидеть объекты небольшого размера.
Работа выполнена в рамках государственного контракта П2613 Федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» на 2009-2013 гг.
Литература
1. Жиляков, Е.Г. Вариационный метод оценивания производных и интерполяции сигналов по эмпирическим данным / Т.Н. Созонова, И.Ю. Мисливец// «Вестник ВГУ», сер. Системный анализ и информационные технологии, 2006, № 2. - С. 70-73.
2. К Де Бор Практическое руководство по сплайнам [Текст]: - М.: Радио и связь, 1985.- 304 с.
3. Ланцош, К. Практические методы прикладного анализа [Текст] : справ, рук. / К Ланцош ; пер. с англ. М. 3. Кайнера. - М.: Физматгиз, 1961. - 524 с.
4. Созонова, Т.Н. Применение вариационных алгоритмов интерполяции и оценки первой производной для некоторых аспектов обработки изображений [Текст]/ Т.Н. Созонова, Н.С. Титова, Н.В. Щербинина // Научные ведомости БелГУ. Сер. История. Политология. Экономика. Информатика. - 2008. - № ю (50). Вып. 8/1. - С. 4 - 12.
5. Хургин, Я. И. Финитные функции в физике и технике [Текст] / Я.И. Хургин, В.П. Яковлев. - М.: Наука, 1971. - 408 с.: ил.
SCALED IMAGES SHARPNESS ENHANCEMENT WITH VARIATIONAL APPROACH
TO NUMERICAL DIFFERENTIATION
The paper describes a new method for signal derivatives calculation from its discrete values, based on frequency representations. An application of the method to gradient processing of digital images aimed at increasing sharpness is covered.
Key words: differentiation, image sharpness, frequency representations, variational principle.
T.N. BALABANOVA O.N. IVANOV
Belgorod State University e-mail:
sozonova@bsu.edu.ru