Этот раздел берет вас через функции 2D дискретного анализа вейвлета с помощью программного обеспечения Wavelet Toolbox™. Тулбокс обеспечивает эти функции для анализа изображения. Для получения дополнительной информации смотрите страницы ссылки на функцию.
В этом разделе представление и примеры используют 2D массивы, соответствующие индексируемым представлениям изображений. Однако описанные функции также доступны при использовании изображений истинного цвета, которые представлены m-by-n-by-3 массивы uint8. Для получения дополнительной информации о форматах изображения смотрите Вейвлеты: Работа с Изображениями.
FunctionName | Цель |
|---|---|
Обеспечьте значения по умолчанию для шумоподавления и сжатия | |
Оштрафованный порог для вейвлета 1D или 2D шумоподавление | |
Пороги для вейвлета 2D использование стратегия Birgé-Massart | |
Шумоподавление вейвлета и сжатие | |
Пороговый менеджер по настройкам |
В этом разделе вы будете учиться
Как загрузить изображение
Как анализировать изображение
Как выполнить одноуровневые и многоуровневые разложения изображений и реконструкции (только командная строка)
Как использовать функции режима Square и Tree (только графический интерфейс пользователя)
Как увеличить масштаб детали (только графический интерфейс пользователя)
Как сжать изображение
В этом примере мы покажем, как можно использовать 2D анализ вейвлета, чтобы сжать изображение эффективно, не жертвуя его ясностью.
Вместо того, чтобы непосредственно использовать image(I), чтобы визуализировать изображение I, мы use image(wcodemat(I)), который отображает перемасштабированную версию продвижения I к более ясному представлению деталей и приближениям (см. страницу с описанием wcodemat).
Загрузите изображение.
От подсказки MATLAB® ввести
load wbarb; whos
| Имя | Размер | Байты | Класс |
|---|---|---|---|
X | 256x256 | 524288 | double array |
map | 192x3 | 4608 | double array |
Отобразите изображение. Ввод
image(X); colormap(map); colorbar;
Преобразуйте индексируемое изображение в полутоновое изображение.
Если палитра сглаженна, вейвлет преобразовывают, может быть непосредственно применен к индексируемому изображению; в противном случае индексируемое изображение должно быть преобразовано в полутоновый формат. Для получения дополнительной информации смотрите Вейвлеты: Работа с Изображениями.
Поскольку палитра сглаженна в этом изображении, можно теперь выполнить разложение.
Выполните одноуровневое разложение вейвлета.
Чтобы выполнить одноуровневое разложение изображения с помощью вейвлета bior3.7, ввести
[cA1,cH1,cV1,cD1] = dwt2(X,'bior3.7');
Это генерирует содействующие матрицы уровня одно приближение (cA1) и горизонталь, вертикальные и диагональные детали (cH1, cV1, cD1, соответственно).
Создайте и отобразите приближения и детали от коэффициентов.
Чтобы создать уровень одно приближение и детали (A1, H1, V1 и D1) от коэффициентов cA1, cH1, cV1 и cD1, вводят
A1 = upcoef2('a',cA1,'bior3.7',1);
H1 = upcoef2('h',cH1,'bior3.7',1);
V1 = upcoef2('v',cV1,'bior3.7',1);
D1 = upcoef2('d',cD1,'bior3.7',1);
или
sx = size(X); A1 = idwt2(cA1,[],[],[],'bior3.7',sx); H1 = idwt2([],cH1,[],[],'bior3.7',sx); V1 = idwt2([],[],cV1,[],'bior3.7',sx); D1 = idwt2([],[],[],cD1,'bior3.7',sx);
Чтобы отобразить результаты разложения уровня 1, ввести
colormap(map);
subplot(2,2,1); image(wcodemat(A1,192));
title('Approximation A1')
subplot(2,2,2); image(wcodemat(H1,192));
title('Horizontal Detail H1')
subplot(2,2,3); image(wcodemat(V1,192));
title('Vertical Detail V1')
subplot(2,2,4); image(wcodemat(D1,192));
title('Diagonal Detail D1')
Регенерируйте изображение одноуровневым Обратным Преобразованием Вейвлета.
Чтобы найти обратное преобразование, ввести
Xsyn = idwt2(cA1,cH1,cV1,cD1,'bior3.7');
Это восстанавливает или синтезирует оригинальное изображение от коэффициентов приближения уровня 1 и деталей.
Выполните многоуровневое разложение вейвлета.
Чтобы выполнить разложение уровня 2 изображения (снова использующий вейвлет bior3.7), ввести
[C,S] = wavedec2(X,2,'bior3.7');
где X является матрицей оригинального изображения, и 2 уровень разложения.
Коэффициенты всех компонентов разложения второго уровня (то есть, приближения второго уровня и первых двух уровней детализации) возвращены конкатенированные в один вектор, C. Аргумент S является бухгалтерской матрицей, которая отслеживает размеры каждого компонента.
Извлеките приближение и детализируйте коэффициенты.
Чтобы извлечь коэффициенты приближения уровня 2 от C, ввести
cA2 = appcoef2(C,S,'bior3.7',2);
Чтобы извлечь первое - и коэффициенты детали второго уровня от C, ввести
cH2 = detcoef2('h',C,S,2);
cV2 = detcoef2('v',C,S,2);
cD2 = detcoef2('d',C,S,2);
cH1 = detcoef2('h',C,S,1);
cV1 = detcoef2('v',C,S,1);
cD1 = detcoef2('d',C,S,1);
или
[cH2,cV2,cD2] = detcoef2('all',C,S,2);
[cH1,cV1,cD1] = detcoef2('all',C,S,1);
где первый аргумент (' h ', ' v ' или ' d ') определяет тип детали (горизонталь, вертикальная, диагональная) извлеченный, и последний аргумент определяет уровень.
Восстановите приближение Уровня 2 и детали Уровня 1 и 2.
Чтобы восстановить приближение уровня 2 от C, ввести
A2 = wrcoef2('a',C,S,'bior3.7',2);
Чтобы восстановить детали уровня 1 и 2 от C, ввести
H1 = wrcoef2('h',C,S,'bior3.7',1);
V1 = wrcoef2('v',C,S,'bior3.7',1);
D1 = wrcoef2('d',C,S,'bior3.7',1);
H2 = wrcoef2('h',C,S,'bior3.7',2);
V2 = wrcoef2('v',C,S,'bior3.7',2);
D2 = wrcoef2('d',C,S,'bior3.7',2);
Отобразите результаты многоуровневого разложения.
Со всеми подробностями, вовлеченными в многоуровневое разложение изображений, это целесообразно импортировать разложение в Вейвлет 2D графический инструмент по порядку, чтобы более легко отобразить его. Для получения информации о том, как сделать это, смотрите Разложения Загрузки.
Чтобы отобразить результаты разложения уровня 2, ввести
colormap(map);
subplot(2,4,1);image(wcodemat(A1,192));
title('Approximation A1')
subplot(2,4,2);image(wcodemat(H1,192));
title('Horizontal Detail H1')
subplot(2,4,3);image(wcodemat(V1,192));
title('Vertical Detail V1')
subplot(2,4,4);image(wcodemat(D1,192));
title('Diagonal Detail D1')
subplot(2,4,5);image(wcodemat(A2,192));
title('Approximation A2')
subplot(2,4,6);image(wcodemat(H2,192));
title('Horizontal Detail H2')
subplot(2,4,7);image(wcodemat(V2,192));
title('Vertical Detail V2')
subplot(2,4,8);image(wcodemat(D2,192));
title('Diagonal Detail D2')
Восстановите оригинальное изображение от многоуровневого разложения.
Чтобы восстановить оригинальное изображение от структуры разложения вейвлета, ввести
X0 = waverec2(C,S,'bior3.7');
Это восстанавливает или синтезирует оригинальное изображение от коэффициентов C многоуровневого разложения.
Сожмите изображение и отобразите его.
Чтобы сжать оригинальное изображение X, используйте команду ddencmp, чтобы вычислить параметры по умолчанию и команду wdencmp, чтобы выполнить фактическое сжатие. Ввод
[thr,sorh,keepapp]= ddencmp('cmp','wv',X);
[Xcomp,CXC,LXC,PERF0,PERFL2] = ...
wdencmp('gbl',C,S,'bior3.7',2,thr,sorh,keepapp);
Обратите внимание на то, что мы передаем в wdencmp результаты разложения (C и S), мы вычислили на шаге 7. Мы также задаем вейвлеты bior3.7, потому что мы использовали этот вейвлет, чтобы выполнить исходный анализ. Наконец, мы задаем глобальную опцию пороговой обработки ' gbl '. Смотрите ddencmp и страницы с описанием wdencmp для получения дополнительной информации об использовании этих команд.
Чтобы просмотреть сжатое изображение бок о бок с оригиналом, ввести
colormap(map);
subplot(121); image(X); title('Original Image');
axis square
subplot(122); image(Xcomp); title('Compressed Image');
axis square
PERF0
PERF0 =
49.8076
PERFL2
PERFL2 =
99.9817
Эти возвращенные значения говорят, соответственно, какой процент коэффициентов вейвлета был обнулен и какой процент энергии изображения был сохранен в процессе сжатия.
Обратите внимание на то, что, даже при том, что сжатое изображение создается только из приблизительно вдвое меньше ненулевых коэффициентов вейвлета, чем оригинал, в качестве изображения нет почти никакого обнаруживаемого ухудшения.
В этом разделе мы исследуем то же изображение как в предыдущем разделе, но мы используем приложение Wavelet Analyzer, чтобы анализировать изображение.
Запустите 2D аналитический инструмент вейвлета.
От посдказки MATLAB введите waveletAnalyzer.
Главное меню Инструмента Вейвлета появляется.
Кликните по пункту меню Wavelet 2-D. Дискретный аналитический инструмент вейвлета для 2D данных изображения появляется.
Загрузите изображение.
В подсказке команды MATLAB ввести
load wbarb
X. Нажмите ОК, чтобы импортировать изображение.Изображение загружается в инструмент Wavelet 2-D.
Анализируйте изображение.
Используя меню Wavelet и Level, расположенные к верхнему правому углу, определите семейство вейвлетов, тип вейвлета и количество уровней, которые будут использоваться для анализа.
Для этого анализа выберите вейвлет bior3.7 на уровне 2.
Нажмите кнопку Analyze. После паузы для вычисления инструмент Wavelet 2-D отображает свой анализ.
Использование квадратных функций режима
По умолчанию анализ появляется в “Квадратном Режиме”. Этот режим включает четыре различных отображения. В верхнем левом углу оригинальное изображение. Ниже этого изображение, восстановленное от различных приближений и деталей. К нижнему правому углу разложение, показывающее самые грубые коэффициенты приближения и всю горизонталь, диагональ и вертикальные коэффициенты детали. Наконец, пробел визуализации в правых верхних отображениях любой компонент анализа, на который вы хотите посмотреть более тесно.
Нажмите на любой компонент разложения в нижнем правом окне.
Синяя граница подсвечивает выбранный компонент. В нижнем правом углу окна Wavelet 2-D существует набор трех кнопок, маркировал “Operations on selected image”. Обратите внимание на то, что, если вы щелкнете еще раз по тому же компоненту, вы отмените выбор его, и синяя граница исчезает.
Нажмите кнопку Visualize.
Выбранное изображение отображено в области визуализации. Вы видите сырые данные, неисправимые 2D коэффициенты вейвлета. Используя другие кнопки, можно отобразить восстановленную версию выбранного компонента изображения, или можно просмотреть выбранный компонент в полноэкранном разрешении.
Использование древовидных функций режима
Выберите Tree из меню View Mode.
Ваше отображение изменяется, чтобы показать следующее.
Это - та же информация, показанная в квадратном режиме, с, кроме того, всеми коэффициентами приближения, но расположенный, чтобы подчеркнуть древовидную структуру разложения. Различные кнопки и меню работают все равно, как они делают в квадратном режиме.
Увеличивание масштаб детали
Перетащите поле круглой резинки (удержанием левой кнопки мыши) по фрагменту изображения, которое вы хотите увеличить.
Нажмите кнопку XY + (расположенный внизу экрана), чтобы масштабировать горизонтально и вертикально.
Инструмент Wavelet 2-D увеличивает отображенные изображения.
Чтобы масштабировать назад к исходному увеличению, нажмите History <<-кнопка.
Нажмите кнопку Compress, расположенную к верхнему правому углу окна Wavelet 2-D. Окно Wavelet 2-D Compression появляется.
Инструмент автоматически выбирает уровни пороговой обработки, чтобы обеспечить, хороший начальный баланс между сохранением энергии изображения при минимизации количества коэффициентов должен был представлять изображение.
Однако можно также настроить пороги вручную с помощью опцией пороговой обработки Уровня, и затем ползунками или редактированиями, соответствующими каждому уровню.
В данном примере выберите опцией пороговой обработки Уровня и выберите Remove около 0 методов из Избранного меню метода задания порога.
Следующее окно отображено.
Выберите из меню направления, хотите ли вы настроить пороги для горизонтали, диагональных или вертикальных деталей. Чтобы внести фактические корректировки для каждого уровня, используйте ползунки или используйте левую кнопку мыши, чтобы непосредственно перетащить пунктирные линии.
Чтобы сжать оригинальное изображение, нажмите кнопку Compress. После паузы для вычисления сжатое изображение отображено около оригинала. Заметьте, что сжатие устраняет почти половину коэффициентов, все же никакое обнаруживаемое ухудшение изображения не появляется.
Покажите невязки.
От инструмента Wavelet 2-D Compression нажмите кнопку Residuals. Окно More on Residuals for Wavelet 2-D Compression появляется.
Отображенные статистические данные включают меры тенденции (среднее значение, режим, медиана) и дисперсия (область значений, стандартное отклонение). Кроме того, инструмент предоставляет схемы плотности распределения (гистограммы и совокупные гистограммы). Тот же инструмент существует для инструмента Wavelet 2-D Denoising.
Статистические данные, отображенные в вышеупомянутой фигуре, связаны с отображенным изображением, но не с исходным. Обычно этой информацией является то же самое, но в некоторых случаях, краевые эффекты могут заставить оригинальное изображение быть обрезанным немного. Чтобы видеть точную статистику, используйте функции командной строки, чтобы получить желаемое изображение и затем применить желаемый MATLAB статистическая функция (функции).
2D графический инструмент Вейвлета позволяет вам импортировать информацию из и экспортировать информацию в диск, если вы придерживаетесь соответствующих форматов файлов.
Можно сохранить синтезируемые образы, коэффициенты и разложения от инструмента Wavelet 2-D до диска, где информацией можно управлять и позже повторно импортировать в графический инструмент.
Сохранение Синтезируемых Образов. Можно обработать изображение в инструменте Wavelet 2-D, и затем сохранить обработанный образ к MAT-файлу (с дополнительным mat или другим).
Например, загрузите анализ в качестве примера:
Файл> Анализ В качестве примера> Индексируемые Изображения> на уровне 3, с Детайлем Дюрером sym4
и выполните сжатие на оригинальном изображении. Когда вы закрываете окно Wavelet 2-D Compression, обновляете синтезируемое изображение путем нажимания кнопку Да в диалоговом окне, которое появляется.
Затем от инструмента Wavelet 2-D выберите пункт меню File> Save> Synthesized Image. Диалоговое окно кажется позволяющим вам, чтобы выбрать папку и имя файла для MAT-файла (с дополнительным mat или другим). В данном примере выберите имя symage.
Чтобы загрузить изображение в вашу рабочую область, ввести
load symage whos
| Имя | Размер | Байты | Класс |
|---|---|---|---|
X | 359x371 | 1065512 | double array |
map | 64x3 | 1536 | double array |
valTHR | 1x1 | 8 | double array |
wname | 1x4 | 8 | char array |
Синтезируемое изображение дано X, и map содержит палитру. Кроме того, параметры процесса шумоподавления или сжатия даны именем вейвлета (wname) и глобальный порог (valTHR).
Сохранение Дискретных Коэффициентов Преобразования Вейвлета. Инструмент Wavelet 2-D позволяет вам сохранить коэффициенты дискретного вейвлета преобразовывает (DWT) к диску. Тулбокс создает MAT-файл в текущей папке с именем, которое вы выбираете.
Чтобы избавить коэффициенты DWT от существующего анализа, используйте пункт меню File> Save> Coefficients.
Диалоговое окно появляется, который позволяет вам задать папку и имя файла для хранения коэффициентов.
Рассмотрите анализ в качестве примера:
Файл> Анализ В качестве примера> Индексируемые Изображения> на уровне 3, с Детайлем Дюрером sym4
После сохранения дискретных коэффициентов вейвлета к файлу cfsdurer.mat загрузите переменные в свою рабочую область:
load cfsdurer whos
| Имя | Размер | Байты | Класс |
|---|---|---|---|
coefs | 1x142299 | 1138392 | double array |
map | 64x3 | 1536 | double array |
sizes | 5x2 | 80 | double array |
valTHR | 0x0 | 0 | double array |
wname | 1x4 | 8 | char array |
Переменный map содержит палитру. Переменный wname содержит имя вейвлета, и valTHR пуст, поскольку синтезируемое изображение совпадает с исходным.
Переменные coefs и sizes содержат дискретные коэффициенты вейвлета и связанные матричные размеры. Более точно, в вышеупомянутом примере, coefs 1 142299 вектор конкатенированных коэффициентов, и sizes дает длину каждого компонента.
Сохранение Разложений. Инструмент Wavelet 2-D позволяет вам сохранить целый набор данных от дискретного анализа вейвлета до диска. Тулбокс создает MAT-файл в текущей папке с именем, которое вы выбираете, сопровождаемый дополнительным wa2 (2D анализ вейвлета).
Откройте инструмент Wavelet 2-D и загрузите анализ в качестве примера:
Файл> Анализ В качестве примера> Индексируемые Изображения> на уровне 3, с Детайлем Дюрером sym4 .
Чтобы сохранить данные от этого анализа, используйте пункт меню File> Save> Decomposition.
Диалоговое окно появляется, который позволяет вам задать папку и имя файла для того, чтобы хранить данные о разложении. Введите имя decdurer.
После сохранения данных о разложении к файлу decdurer.wa2 загрузите переменные в свою рабочую область:
load decdurer.wa2 -mat whos
| Имя | Размер | Байты | Класс |
|---|---|---|---|
coefs | 1x142299 | 1138392 | double array |
data_name | 1x6 | 12 | char array |
map | 64x3 | 1536 | double array |
sizes | 5x2 | 80 | double array |
valTHR | 0x0 | 0 | double array |
wave_name | 1x4 | 8 | char array |
Переменные coefs и sizes содержат структуру разложения вейвлета. Другие переменные содержат имя вейвлета, палитру и имя файла, содержащее данные. Переменный valTHR пуст, поскольку синтезируемое изображение совпадает с исходным.
Сохраните опции, также доступны при выполнении шумоподавления или сжатия в инструменте Wavelet 2-D. В окне Wavelet 2-D Denoising можно сохранить denoised образ и разложение. То же самое сохраняется для окна Wavelet 2-D Compression. Таким образом, можно сохранить много различных испытаний из окон Denoising и Compression, не возвращаясь к основному окну Wavelet 2-D во время процесса подстройки. При сохранении синтезируемого сигнала, разложения или коэффициентов к MAT-файлу, расширение файла mat не необходимо. Можно сохранить приближения индивидуально для каждого уровня или сохранить их целиком.
Можно загрузить изображения, коэффициенты или разложения в приложение Wavelet Analyzer. Информация, которую вы загружаете, возможно, ранее экспортировалась из приложения Wavelet Analyzer, и затем управлялась в рабочей области; или это, возможно, была информация, которую вы сгенерировали первоначально из командной строки.
Или в случае, необходимо наблюдать строгие форматы файлов и структуры данных, используемые инструментом Wavelet 2-D, или иначе в ошибки закончатся, когда вы попытаетесь загрузить информацию.
Загрузка Изображений. Этот тулбокс поддерживает только индексированные изображения. Индексируемое изображение является матрицей, содержащей только целые числа от 1 до n, где n является количеством, раскрашивает изображение.
Это изображение может опционально сопровождаться n-by-3 матрица под названием map. Это - палитра, сопоставленная с изображением. Когда MATLAB отображает такое изображение, он использует значения матрицы, чтобы искать требуемый цвет в этой палитре. Если палитра не дана, инструмент Wavelet 2-D использует монотонную палитру с цветами max(max(X))min(min(X))+1.
Чтобы загрузить изображение, вы создали в своем рабочем пространстве MATLAB в инструмент Wavelet 2-D, сохраните образ (и опционально, переменная map) в MAT-файле (с дополнительным mat или другим).
Например, предположите, что вы создали названный мозг изображения и хотите анализировать его в инструменте Wavelet 2-D. Ввод
X = brain; map = pink(256); save myfile X map
Чтобы загрузить это изображение в инструмент Wavelet 2-D, используйте пункт меню File> Load> Image.
Диалоговое окно появляется, который позволяет вам выбрать соответствующий MAT-файл, чтобы загрузиться.
Графические инструменты позволяют вам загружать изображение, которое не содержит целые числа от 1 до n. Вычисления правильны, потому что они действуют непосредственно на матрицу, но отображение изображения является странным. Значения меньше чем 1 оценен как 1, значения, больше, чем n, оценены как n и действительное значение в интервале [1, n] оценен как самое близкое целое число.
Коэффициенты, приближения и детали, произведенные разложением вейвлета, не индексируются матрицы изображений.
Чтобы отобразить эти изображения подходящим способом, инструмент Wavelet 2-D следует этим правилам:
Восстановленные приближения отображены с помощью палитры map.
Коэффициенты и восстановленные детали отображены с помощью палитры, map применился к перемасштабированной версии матриц.
Первая 2D переменная, с которой сталкиваются в файле (кроме переменной map, которая резервируется для палитры), рассматривается изображением. Переменные осмотрены в алфавитном порядке.
Загрузка Дискретных Коэффициентов Преобразования Вейвлета. Чтобы загрузить коэффициенты дискретного вейвлета преобразовывает (DWT) в инструмент Wavelet 2-D, сначала сохраните соответствующие данные в MAT-файле, который должен содержать, по крайней мере, эти две переменные:
coefs, содействующий вектор
sizes, бухгалтерская матрица
Для индексируемого изображения матричным sizes является n+2-by-2 массив:
Для изображения истинного цвета матричным sizes является n+2-by-3:
Переменный coefs должен быть вектором конкатенированных коэффициентов DWT. Вектор coefs для n - разложение уровня содержит разделы 3n+1, состоя из коэффициентов приближения уровня-n, сопровождаемых горизонталью, вертикальными, и диагональными коэффициентами детали, в том порядке, для каждого уровня. Переменный sizes является матрицей, строки которой задают размер cAn, размер cHn (или cVn или cDn)..., размер cH1 (или cV1 или cD1), и размер оригинального изображения X. Размеры вертикальных и диагональных деталей совпадают с горизонтальной деталью.
После построения или редактирования соответствующих данных в вашей рабочей области, ввести
save myfile coefs sizes
Используйте пункт меню File> Load> Coefficients от инструмента Wavelet 2-D, чтобы загрузить данные в графический инструмент.
Диалоговое окно появляется, позволяя вам выбрать папку и файл, в котором находятся ваши данные.
Загрузка Разложений. Чтобы загрузить дискретный вейвлет преобразовывают данные о разложении в инструмент Wavelet 2-D, необходимо сначала сохранить соответствующие данные в MAT-файле (с дополнительным wa2 или другим).
MAT-файл содержит эти переменные.
| Переменная | Состояние | Описание |
|---|---|---|
coefs | Необходимый | Вектор конкатенированных коэффициентов DWT |
sizes | Необходимый | Матричные размеры определения компонентов |
wave_name | Необходимый | Имя определения вектора символов вейвлета используется для разложения (например, |
map | Дополнительный |
|
data_name | Дополнительный | Имя определения вектора символов разложения |
После построения или редактирования соответствующих данных в вашей рабочей области, ввести
save myfile.wa2 coefs sizes wave_name
Используйте пункт меню File> Load> Decomposition от инструмента Wavelet 2-D, чтобы загрузить данные о разложении изображений.
Диалоговое окно появляется, позволяя вам выбрать папку и файл, в котором находятся ваши данные.
При загрузке изображения, разложения или коэффициентов из MAT-файла, расширение этого файла свободно. Расширение mat не необходимо.