CC BY
141
УДК 535.416.3
КОНСТРУИРОВАНЕ БИОРТОГОНАЛЬНЫХ И КОМПЛЕКСНЫХ ВЕЙВЛЕТ-БАЗИСОВ ДЛЯ ОБРАБОТКИ ОПТИЧЕСКИХ ИЗОБРАЖЕНИЙ
Ю.Н. Исаев
Томский политехнический университет; Институт оптического мониторинга СО РАН E-mail: Isaev_Yusup@mail.ru
Описывается алгоритм синтеза биортогональных и комплексных вейвлетов. Приводятся примеры восстановления и сжатия сигналов на основе сконструированных вейвлетов. Описывается их сходство с базисом Карунена-Лоева. Представлена картина пространно-временной развёртки случайного сигнала по масштабам неоднородности. Приводятся примеры разложения и сжатия двумерных сигналов с фильтрацией по направлениям, примеры выделения контуров и усиления контраста изображений.
Ранее [1] мы рассматривали синтез биортогональных вейвлетов, в данной работе мы продолжим рассмотрение биортогональны вейвлетов и приведём алгоритм построения комплексных вейвлетов, которые полезны при анализе и синтезе амплитуды и фазы оптических комплексных сигналов, а также весьма эффективны для выделения контуров изображения. В заключение приведем примеры сжатия и повышения контрастности двумерных изображений.
Критерии синтеза вейвлет-базиса
Для полученных вейвлетов, приведённых в [1], автором проверялся критерий локализации, который выражается соотношением
АхАк> 2, (1)
где
[|^(х)|2(х - х0)<±с №>(к)\ (к-к0)с1к
Ах =}' '-2--, Дк = ^—-I— ------------. (2)
||Ч/(.г)|" с1х Цчч*-)!" йх-
При увеличении числа нулевых моментов вейвлетов неравенство (1) быстро стремилось к равенству. То есть вейвлеты модулируются гауссовой функцией, значит они приближаются к собственным функциями оператора Фурье [2] и становятся оптимальными по критерию (1). Это означает, что сигнал ограниченной протяжённости при фиксированном числе членов разложения будет иметь наибольшую энергию, если он представлен вейвлет-базисом, удовлетворяющем равенству (1), это -сходство с базисом Карунена-Лоева. С помощью таких вейвлетов удобно определять скрытую периодичность исследуемых сигналов и определять их мультифрактальную структуру. Приведём вейвлет-разложение случайного сигнала (скорости ветра), приведенного на рис. 1 (вверху). В качестве базиса выберем вейвлет, изображённый на рис. 2, в из работы [1]. На рис. 1 (второй и третий сверху) приведены абсолютные значения коэффициентов вейвлет-разложения случайного сигнала. Это - картина, показывающая относительный вклад неоднородностей разного масштаба в формирование сиг-
нала. Из картины видно, что крупномасштабные неоднородности рассыпаются на более мелкие, в свою очередь мелкие неоднородности рассыпаются в ещё более мелкие, и это каскадный процесс продолжается до самой мелкой фракции неоднородности о чём свидетельствуют характерные разветвления - "вилки".
Рис. 2. Трапвциидальная скейлинговая функция
Этот пример демонстрирует сходство случайного сигнала с мултифрактальной структурой неоднородного канторовского ряда, а, значит, среду
распространения сигнала описывать степенным спектром типа К можно, если только считать, что а является функцией к. Такою физическую картину описания неоднородностей удаётся получить благодаря фрактальным свойствам скейлинговых функций вейвлетов, то есть ур. (4) в работе [1], смысл которого графически представлен на рис. 2, где для простоты в качестве скейлинговой функции представлена трапециидальная функция
ф(х)=0,5р(2х)+0,5ф(2х-1)+0,5ф(2х-2)+0,5ф(2х-3). (3)
Каждая трапеция состоит из четырёх себе подобных трапеций, этот каскадный процесс, (напоминающий каскадный процесс взаимодействия вихрей в турбулентном процессе) можно продолжать до бесконечности в обоих направлениях, как в сторону увеличения, так и в сторону уменьшения размеров трапеций. Можно без труда заметить, что подобный базис более точно соответствует структуре турбулентного потока и, следовательно, и его адекватному описанию.
р0 = 0,003536
- ) 0,049875 Р1 = 0,703571 +] 0,049875
Р2 = 0,703571 +] 0,049875 Р3= 0,003536
- } 0,049875
Действительная часть
В
Мнимая часть
Гя11п д ЛтсИо:
г
а)
р0= 0,25 - ] 0,25 р1=0,25 + ] 0,25 р2= 0,25 + } 0,25 Р3= 0,25 - ) 0,25
Действительная часть
Итсио:
Мнимая часть
tt4.oth.er \vavelet
б)
р0 =-,0785674 - } 0,1111111 Р1= 0,1571348 - ) 0,1111111 Р2 = 0,6285394 +] 0,222222 Р3 = 0,6285394 +] 0,222222 Р4= 0,1571348 - ) 0,1111111 Р5=-0,0785674 - } 0,1111111
Л
/ \
/
р*/ V—^
Действительная часть
Мнимая часть
в)
Рис. 3. Комплексные вейвлеты. Слева - скейлинговые функции, справа - вейвлеты
При компрессии сигнала необходимо в его разложении иметь наименьшее число членов разложения сохраняя при этом высокой степени аппроксимации. В этом случае при выборе базиса необходимо руководствовать энтропийным критерием.
N
Е= -Х ак ъ&(ак) ^ тп-1
Чем меньше энтропия коэффициентов разложения, тем меньше коэффициентов остаются в разложении. Этот критерий означает, что удалось сконструировать вейвлет, наилучшим образом отражающий профиль анализируемого сигнала. Пример использования этого критерия будет продемонстрирован на двумерных сигналах.
Комплексные вейвлеты
В ряде случаев анализ амплитуды и фазы оптического сигнала удобно производить в его комплексном представлении. В этом случае удобно для разложения использовать комплексные вейвлеты. Скейлинговая функция ф(х) в этом случае будет имеет комплексные коэффициенты, а двойственная ей функция будет сопряжённой ей функцией ф{х)=ф{х). (4)
Приведём алгоритм синтеза комплексных вейвлетов.
Задаёмся вещественной частью коэффициентов в разложении (4), сумма которых равна единице, например
Лф„) = ЯФО =Лф2)= Яв(р3) =1/4, (5)
затем находим мнимую часть коэффициентов 1т(рК) из уравнения ортогональности (8) из работы [1] и с учётом условия (4)
Р„ = Рз =1/4+/1/4, Р1= Р2 =1/4+/'1/4. (6)
После определения коэффициентов находим скейлинговые и вейвлет-функции с помощью выражений (4) и (5), приведённых в [1], соответственно. Комплексные вейвлеты можно получить также, используя комплексное решение уравнения Безу. Например, при N=3 решение можно записать так
В{х) = (х - zj)(х - ^ X
У=1
(7)
тогда в выражения для фильтров т(к) и т(к) будут входить различное сочетание комплексных множителей из (7). На рис. 3 приведены несколько комплексных вейвлетов, полученных по выше описанному алгоритму.
2D-вейвлет-базисы
Приведём пример использования вейвлет-базиса для обработки двумерных изображений. Запишем тензорное произведение одномерных вейвлетов и перейдем к определению трёх двумерных вейвлетов
Уг{х,у) =ф(х)х¥(у), х¥,(х,у)=х¥(х)ф(у), Х¥,(х,у)=ф(х) <р(у). (8)
Такое представление удобно тем, что осуществляется фильтрация изображения по горизонтальным вертикальны и диагональном направлениях, соответствующие вейвлеты снабжены индексами. Путем сдвигов и сжатия образуется базис
^Мй(*,У)= 2-'Т(2-х-и1,2-у-и2). (9)
Все соотношения (1-5), приведённые в [1] для одномерного случая, остаются верными для двумерного базиса.
В качестве примера работы вейвлетов приведём пример сжатия 2D-изображения (рис. 4) с помощью вейвлета, изображённого на рис. 2, б в работе [1]. Исходное количество коэффициентов было 256x256 с последующем уменьшением в 4 раза. Для выделения краевых эффектов изображения с учётом направления выберем вейвлет, изображённый на рис. 6, б в работе [1]. Результаты восстановления опробованы на двух различных структурах и приведены на рис. 5. Как рис. 4, так и рис. 5 отчётливо демонстрируют полезное свойство вейвлетов фильтровать по горизонтальным, вертикальным и диагональным направлениях.
Приведём пример выбора наиболее оптимального вейвлета. На рис. 6 представлено изображение разложенное по трём различным вейвлетам. Изображения расположены снизу вверх по увеличению степени сжатия, а, следовательно, оптимальности вейвлета. Оптимальность разложения определя-
N
лась по энтропийному критерию Е = -£ ак ).
1
Оптимальность разложения можно увидеть "нево-
Рис. 4. Сжатие и фильтрация изображения по направлениям
Рассмотрим применение комплексного вейвлета для выделения контуров изображения. При анализе функций с сингулярными особенностями хорошо известен метод её регуляризации. Сингулярную функцию аналитически продолжают в комплексную область. Тогда функция разделяется на обобщённую (гладкая версия исходной) и сингулярный вклад который является мнимым [3, 4]. Если рассматривать изображение, как функцию с
mother Wavelet
f-rk *
HI
Рис. 6. Выбор оптимального вейвлета для сжатия изображения
Рис. 5. Фильтрация текстового изображения по направлениям
оружённым глазом". Заметим, что первое изображение (снизу вверх) после разложения второго уровня (т.е. после сжатия изображения в 4 раза) потеряло значительно меньше информации, чем второе, а второе изображение потеряло меньше, чем третье. За приделами не оптимально сжатого изображения остаётся узнаваемый рельеф изображения. Особенно это хорошо видно после фильтрации (рис. 7).
Рис. 7. Выбор оптимального вейвлета для сжатия изображения с фильтрацией высокочастотного шума
Re(mother wavelet)
Л)
V
Im(mother wavelet)
Im(Y1)
Im(Y2)
H
Y1 Y2
Рис. 8. Использование комплексных вэйлетов для выделения контуров изображения
особенностями, и, используя комплексное вейвлет-разложение изображения, можно ожидать, что сильные изменения градиента интенсивности изображения окажутся в мнимой части разложения. Следует только правильно сконструировать комплексный вейвлет. Необходимо, чтобы мнимая часть вейвлета имела как можно больше нулевых моментов, тогда его мнимая часть будет быстро осциллирующей [5, 6]. Автором сконструирован комплексный вейвлет, изображённый на рис. 8 на основе решения уравнений (7), (8) из работы [1] и с дополнительным условием
Яе[£(-1)кктрк | = 0, т = 0...N/2.
V к=0 )
На том же рисунке продемонстрировано действие вейвлета на изображение.
Первое изображение сверху разложено в действительную часть вейвлета, второе изображение сверху разложено в мнимую часть вейвлета. Третий изображение сверху получено из второго после фильтрации шума.
Рассмотрим использование вейвлетов для повышения контрастности изображения. Одной из важных проблем видения изображений, наблюдаемых в рассеивающих средах и описываемых в терминах уравнения сверки, является задача восстановления замутненных изображений. Класс задач, связанных с этой проблемой, например, в спутниковом мониторинге подстилающей поверхности Земли, необозримо велик. Сложность решения этой задачи определяется тем, что функция размытия точки на момент регистрации изображений, как правило, неизвестна, кроме того, необходимо
—1
1
□ 50 1 □□ 150
Рис. 9. Вейвлет-разложение сигнала по уровням разрешения
Рис. 10. Сужение функции рассеяния точки
Рис. 11. Усиление контрастности изображения
учитывать некорректность операции обращения свертки. Среди многообразия реальных атмосферно-оптических ситуаций нередки случаи, когда мы наблюдаем достаточно четкое, хотя и слабоконтрастное изображение на фоне общего "размытия". Это соответствует такому виду функции размытия точки, которая содержит четко выраженную дельта-составляющую и "размытое основание". Если
/- м м
построить преобразование, "усиливающее" роль высокочастотных составляющих наблюдаемого размытого изображения и подавить фоновую составляющую, то получим восстановленное высококонтрастное изображение. Технология использования вейвлет-преобразований для решения этой задачи заключается в следующем. Выбирается объект, имеющий резкие границы, например объект с прямоугольными краями (рис. 9). Производится его вейвлет-разложение по уровням разрешения. Далее, для усиления роли высоких частот вейвлет-коэффициенты высокого уровня разрешения помножаются на множитель выше единицы, а низкочастотные - на множитель меньше единицы. Для примера на рис. 9 было выбрано 7 уровней разложения. Самый низкий уровень разложения с коэффициентом £7=„,9 и дальше по уровням возрастания разрешения, соответственно, £6=„,8, £5=1,„, £4=1,4, £3=1,5, £2=1,5, &1=1,5. Затем производится обратное вейвлет-преобразование и получатся кар-
тина, приведённая на рис. 1„, где заметно усилена роль краёв изображения. Таким образом, перераспределяя относительный вклад различных уровней разложения можно сузить функцию рассеяния точки. После того, как коэффициенты уровней определены, можно переходить к обработке 2D-изобра-жения. В качестве примера рассмотрено спутниковое изображение Новосибирского водохранилища (рис. 11). Множители по уровням разрешения были выбраны те же, что и для рис. 1„. Видно, что в обработанном изображении повышенная контрастность, чётче узнаётся рельеф местности.
В настоящей работе представлены группа новых биортогональных и комплексных вейвлетов, приведён алгоритм их синтеза. Приведены примеры синтеза и анализа модельных сигналов. Приведена возможность визуализации фрактальной структуры сигнала, варьирования масштабов неоднородности среды распространения сигнала с помощью коэффициентов вейвлет-разложения. Продемонстрированы полезные свойства полученных вейвлетов фильтровать по направлениям и выделять тонкую структуру изображения. Полученные комплексные сигналы позволяют выделять контуры изображения и оптимально их представлять. Полученные вейвлеты позволяют улучшить контрастность изображения.
СПИСОК ЛИТЕРАТУРЫ
Исаев Ю.Н. Конструирование биортогональных вейвлет-базисов для оптимального представления сигналов // Известия Томского политехнического университета. - 2„„4. - Т. 3„7. -№ 1. -С. 37-42.
Хургин Я.И., Яковлев В.П. Финитные функции в физике и технике. -М.: Наука, 1971. — 4„6 с.
Гельфанд И.М., Шилов ГЕ. Обобщенные функции и действия над ними. - М: Физматгиз, 1959. - 47„ с.
Владимиров В.С. Обобщённые функции в математической физике. - М.: Наука, 1976. - 347 с.
Исаев Ю.Н. Синтез вейвлет-базиса для анализа оптических сигналов. Ч. I. Ортогональный вейвлет-базис // Оптика атмосферы и океана. — 2„„2. -Т. 15. - № 11. - С. 974-981.
Исаев Ю.Н. Синтез вейвлет-базиса для анализа оптических сигналов. Ч. III. Представление дифференциальных и обратных им операторов в вейвлет-базисах. Вейвлет пакеты // Оптика атмосферы и океана. - 2„„3. - Т. 16. - № 8. - С. 688-694.
