CC BY
289
УДК 004.42
0 ПОВЫШЕНИИ РЕЗКОСТИ ИЗОБРАЖЕНИЙ НА ОСНОВЕ СУБП0Л0СН0Г0 ДИФФЕРЕНЦИРОВАНИЯ
Е.Г. 9КИЛЯК0В Д.Д. ЧЕРНОМОРЕЦ Е.В. БОЛГОВА
В статье разработан метод повышении резкости изображений на основе субполосного дифференцирования, позволяющий учитывать частотные свойства изображений.
Белгородский государ■
Ключевые слова: изображение, контур, дифференцирование, субпо-лосный анализ, резкость изображений
ственный национальный
исследовательский
университет
e-mail:
zhilyakov@bsu.edu.ru
chernomorets@bsu.edu.ru
При решении многих задач анализа изображений, возникающих в процессе производственной деятельности или в научных исследованиях, зарегистрированные изображения, зачастую, не имеют необходимой резкости, четкости отображения различных объектов, вследствие этого проблеме улучшения резкости изображений и разработке соответствующих методов постоянно уделяется существенное внимание. Под резкостью изображений обычно понимают степень полноты воспроизведения мелких деталей, а также контуров объектов на изображениях. Основная цель повышения резкости изображений состоит в выявлении и выделении мелких деталей изображения или улучшении фрагментов и деталей изображений, которые не являются достаточно четкими для проведения анализа по причине технических характеристик регистрирующей аппаратуры и условий съемки. Как показывают исследования, во многих случаях наиболее информативными при выделении объектов на изображении являются параметры границ однородных областей - контуров данных объектов. Так, в биологических системах зрительного восприятия при выделении отдельных объектов используются, в основном, очертания контуров, а не значения яркости данных объектов [l]. Учитывая данный факт, в настоящей работе для повышения резкости изображений осуществляется поиск контуров объектов.
Контурам объектов (границы) изображения, как правило, соответствует существенный перепад яркости между двумя относительно однородными областями значений яркости. Поскольку математической операцией, позволяющей определить области изменения функции, является дифференцирование, то контуры изображения могут быть определены на основе вычисления производных.
В работе предложено операцию дифференцирования выполнять на основе субполосного анализа-синтеза, который в отличие от большинства известных методов дифференцирования изображений позволяет учитывать их частотные свойства.
Рассмотрим наиболее известные пространственные методы, использующие линейные фильтры, которые применяют для вычисления производных на изображении при подчеркивании перепадов яркости и контуров. Представим изображение в виде значений
двумерной функции яркости в отдельных равноудаленных точках
(хиУк) , i = 1,2,..., Nl ^ к = 1,2,..., N2 (при условии, что регистрация значений яркости изображения осуществлена в эквидистантных отсчетах). Выделение контуров (границ однородных областей) может быть достигнуто на основе операции дифференцирования в пространственной области изображения.
При выделении на изображении вертикальных перепадов осуществляется дифференцирование по строкам, горизонтальных перепадов - по столбцам. В пространственной области изображения шаг между отсчетами обычно принимают равным 1, следовательно,
для дифференцирования строк и столбцов изображения (//,*), i = 1,2,..., N 1?
к = 1,2,..., N2 можно применять следующие выражения:
г , _ _ f
у ду
1 - fЯ
G' = = /■ , , - 2 /, + /.
.XV Л .. 2 i, /с — 1 -'ус -'г,
9v:
к = 2,..., -1
(1)
При вьщелении на изображении контуров произвольной ориентации используют значение модуля градиента, вычисленного для двумерной функции, значение которой соответствуют значениям яркости пикселей изображения. Максимальное изменение двумерной функции происходит в направлении ее градиента. При вычислении первой производной изображения применяют различные приближения двумерного градиента и его модуля.
Известно, что для функции / О - У) градиент в точке (х - У) определяется как двумерный вектор-столбец
rdf
V/ =
G,
dx
#
dy
(2)
Соотношение для вычисления оценки модуля двумерного градиента изображения имеет следующий вид:
V / ( -т , V ) = -JI/ (х / (-Т - 1 , V ) } + [/ ( -Т , V ) - / ( -Т , V - 1) ]2 . (з)
При практической реализации процесса вычисления оценки модуля градиента применяют упрощенное выражение:
|v/.|= И|+ |g ; | s (4)
которое сохраняет относительные изменения в уровнях яркости.
С физической точки зрения рассмотренные выражения соответствуют применению фильтров высоких частот. На основе указанных выражений выделяют высокочастотные составляющие анализируемого изображения, определяющие перепады яркости и контуры, и подавляют «постоянную составляющую». Однако, их использование повышает уровень шума на изображении. Во многих случаях для повышения помехоустойчивости при выполнении операции выделения контуров применяют нелинейную обработку изображений. Особенность нелинейных методов обработки заключается в том, что они предполагают использование нелинейных операторов численного дифференцирования.
Методы численного дифференцирования, наиболее часто применяемые при проведении реальных расчетов по выделению контуров, используют дифференциальные операторы Робертса, Собела и Превитта.
В операторе Робертса для выделения границ используют диагонально смежные (перекрестные) пиксели
f ik | f i. к f i + 1. к + 1 | | f i. к + 1 f i + 1. к | , (S)
/ = 1,2,..., к = 1,2,..., N2- 1.
Превитт предложил для аппроксимации производных двумерной функции использовать следующие выражения, являющиеся некоторым обобщением соотношений (l):
Gx — + fi + \,k f г + 1,A- + 1 ) _ ( f / — 1,A" — 1 fi-\,k f г — 1,Ar +1 )
Gy — ( fj - 1, A- + 1 + fi,k+ 1 f / + 1, A" + 1 ) _ ( f / — 1, A" — 1 fi ,k — \ f / + 1, A" — 1 ) ,
(6)
/ = 2,..., J^-l, k = 2,..., N2- 1.
Для увеличения значимости элементов, используемых в соотношении (l), которое определяет значение производной, Собел ввел весовые коэффициенты в соотношениях
(6) и предложил для вычисления производных вдоль координатных осей применять следующие соотношения:
Gx = (// + 1Д--1 + 2 fi + l к + fi +1,А- +1 ) _ (//-1Д--1 + 2 /,_1Д- + /,-l,A-+l) 5
Gy = (//-1Д- + 1 + 2 fi к + \ + // + 1Д- + 1) _ (ft-l,k-l + 2 /,д--1 + /,+l,A--l) ,
(7)
/ = 2,..„ ^-1 £ = 2,..., 7V:
1
Весовые коэффициенты позволяют при применении метода (7) уменьшить эффект сглаживания за счет придания большего веса средним точкам окрестности.
Рассмотренные выше методы имеют существенный недостаток - при их реализации преобразование осуществляется на уровне отдельных пикселей в пространственной области, вследствие чего построенные контуры, зачастую, не являются гладкими линиями.
Предлагаемый в данной работе метод субполосного дифференцирования изображений разработан на основе результатов исследований, полученных при разработке метода субполосной интерполяции изображений [2]. Интерполяция в указанном методе осуществлена на основе реализации формулы Ньютона-Лейбница:
I
м(0 = М0 + | /(х)с& 5
о
примененной к двумерной функции. В [2] показано, что для интерполяции изображений на основе субполосного анализа-синтеза можно использовать приведенные ниже соотношения.
Исходное изображение представим в виде прямоугольной матрицы вещественных чисел ^ = (ит1т2), т\ = 1,2,..., М1 ^ т2 = 1,2,..., М2? соответствующих яркости отдельных пикселей интерполируемого изображения. Значения интерполирующего изображения и = (мИ1„2) 5 пх = 1,2,..., Nх 5 п2 = 1,2,..., Ы2 должны быть вычислены в ^ \ и В 2
точках между пикселями исходного изображения при условии выполнения интерполирующих равенств в узлах интерполяции:
ио1 (т1 -1) + 1,Д2 (т2-1)+1 — ит1т2 , (8)
т1 = 1,2,..., М1 ^т2 = 1,2,..., М 2 _
Размерности исходного и интерполирующего изображений связаны следующими соотношениями:
Л^ = Д(М1-1) + 1> (9)
Ы2 = В2(М2- 1) + 1 .
Для решения указанной задачи в [2] предложен метод субполосной интерполяции изображений, который позволяет избирательно задавать свойства результатов интерполяции в области пространственных частот (свойство гладкости). Достигается данная характеристика метода за счет применения собственных векторов субполосных матриц [3], соответствующих различным подобластям пространственных частот, а именно - для получения более гладких результатов интерполяции следует применять субполосные матрицы, соответствующие подобласти более низких пространственных частот. Так, в рамках теории субполосного анализа-синтеза [3] требование наибольшей гладкости результатов интерполяции изображения может быть сформулировано как требование размещения основной доли энергии результирующего изображения в области низких частот (в первой
подобласти пространственных частот ^ п [4]):
Оп ={{а,(3)\ -Ц <а< Ц, -02</?<02}5 (10)
где ^1 и ^ 2 - параметры, значение которых влияет на гладкость результата интерпо-
ляции:
П; = —- £22 = —
КХБХ и К2В2 ’
^1 и К 2 - коэффициенты, значения которых выбираются интерактивно.
Обозначим
ии={игк], г = 2,...,М1гк = 2,...,М2г ии = }, г = 2N1 , к = 2,..., N2 ^
Я 1 И б 2 - матрицы, состоящие ИЗ М 1 и М 2— 1 собственных векторов субполосных матриц и [3], соответствующих подобласти пространственных ча-
стот ^ п (ю);
^1 и Н 2 - квадратные нижние треугольные матрицы, состоящие из единиц и нулей, размерностей
(К -\)х(К и (^2-1)х(^2-1)
соответственно;
БелГУ
*
® 1 , ^ 2 - состоящие из единиц векторы, размерностей (^1 1)и(^2 1) соответственно;
У1 , У 2 - векторы, размерностей М 1 и М 2— 1 , состоящие из единиц;
В\, В 2 - матрицы размерностей 1_ 1)х(^ 1_ 1) и (М 2- 1)х(ТУ2-1) соответственно, состоящие из строк матриц ^ 1 и ^2 с номерами 1 +1 ,
2/) ! + 1 5...5 (М 1- 1)£> 1 + 1 и£>2 + 1;2£>2 + 1 5 ; (М 2- 1)/) 2 + 1 соответственно.
Тогда соотношение для вычисления интерполирующего изображения и и имеет вид [2]:
ии = щ&ё* + Ш(ШУ\ии -щЖхОгВ^у'о^в?. (11)
В выражении (п) обратные матрицы для матриц В х() х и (В202)т существуют, поскольку они являются произведением неособенных матриц соответственно.
Учитывая формулу Ньютона-Лейбница и тот факт, что умножение на матрицы В х
и В2 слева и справа соответствует операции интегрирования по строкам и столбцам, можно показать, следующий фрагмент выражения (п):
олАоу (ии - Штв2т у1 о2т,
которое определяет результат субполосного дифференцирования интерполирующего изображения.
Очевидно, что при выборе количества промежуточных точек
А = В2 = 1 , (12)
результат интерполяции (п) совпадает с исходным изображением. Следовательно, при
выполнении соотношений (12) матрица У
У = олтУ\ии-ипГгГ2тШтв2т)-1о2т (13)
определяет результат субполосного дифференцирования исходного изображения.
Для получения более резкого изображения и 5 чем исходное изображение и 5 можно использовать следующее выражение:
и = и +¥ , (14)
где Г - результат градационной коррекции изображения У , при которой диапазон значений элементов матрицы Г не выходит за пределы диапазона значений элементов матрицы и .
Выражения (13), (14) определяют новый метод увеличения резкости изображений на основе субполосного дифференцирования.
Для проверки работоспособности данного метода были проведены вычислительные эксперименты по улучшению резкости изображения на основе предложенного метода и операторов Робертса, Собела и Превитта. Так, при вычислительных экспериментах было использовано исходное изображение, приведенное на рис. 1.
Рис. 1. Исходное изображение
Результаты дифференцирования и улучшения резкости на основе анализируемых методов приведены на рис. 2 и рис. 3-
РЬш
Щф
'ІИІ
1 .у*® і II ,г 5' "V
- X • л ^ "А 4— Д ,
и . У
/ ,
в г
Рис. 2. Результаты дифференцирования исходного изображения: а - субполосное дифференцирование, б - оператор Робертса, в - оператор Собела, г - оператор Превитта
/
а
б
НАУЧНЫЕ ВЕДОМОСТИ В '• серия история, политология, экономика, информатика. -|gу
Рис. 3. Результаты улучшения резкости: субполосное дифференцирование, б - оператор Робертса, в - оператор Собела, г - оператор Превитта
а
Изображения, приведенные на рис. з, демонстрируют, что изображение на рис. за (получено при субполосном дифференцировании) является более резким по сравнению с изображениями на рис. 36, рис. зв и рис. 3г. Также показано, что при улучшении резкости на основе субполосного дифференцирования выделяется больше мелких объектов без существенного искажения изображения, в отличие от результатов применения других рассматриваемых методов.
Проведенный эксперимент, а также другие многочисленные вычислительные эксперименты продемонстрировали работоспособность предложенного в работе метода улучшения резкости и его преимущество по сравнению с методами, основанными на применении операторов Робертса, Собела и Превитта.
Работа выполнена при поддержке гранта РФФИ 12-07-00257-а.
Список литературы
1. Гонсалес Р., Вудс Р. Цифровая обработка изображений. - М.: Техносфера, 2006. - 1072 с.
2. Жиляков, Е.Г. Алгоритм масштабирования изображений [Текст] / Е.Г. Жиляков, А.А. Черноморец // Вопросы радиоэлектроники, Сер. ЭВТ. - 2013. - Вып. 1. - С. 155-1633. Жиляков, Е.Г. Вариационные методы анализа сигналов на основе частотных представлений [Текст] / Е.Г. Жиляков, С.П. Белов, А.А. Черноморец // Вопросы радиоэлектроники, Сер. ЭВТ. - 2010. - Вып. 1. - С. 10-25.
4. Черноморец, А.А. О формировании квазициклических компонент изображений с заданными частотными свойствами [Текст] / В.В. Красильников, А.А. Черноморец // Научные ведомости БелГУ. Сер. История. Политология. Экономика. Информатика. - 2010. - № 13 (84). -
Вып. 15/1. - С. 168-174.
ABOUT IMAGE SHARPNESS INCREASING ALGORITHM BASED ON SUBBANB DIFFERENTIATION
E.G.2HILYAK0V The image sharpness increasing algorithm based on subband dif-
A A PUERUMUMRETC ferentiation is given in the article. Given algorithm allows to take into
ftJt. linEnNUIHUIfE I a account the frequency properties of images.
E.V. BOLGOVA
Keywords: image, edge, differentiation, subband analysis, image
Belgorod National sharpness
Research University
Ц2Щ
e-mail:
zhilyakov@bsu.edu.ru chernomorets @bsu.edu.ru
