В этом примере показано, как извлечь и отобразить изображения деталей уровня разложения вейвлетов.

Загрузка изображения. Выполните вейвлет изображения уровня 2 с помощью haar вейвлет.

load woman
[c,s]=wavedec2(X,2,'haar');

Извлеките приближение уровня 1 и коэффициенты детализации.

[H1,V1,D1] = detcoef2('all',c,s,1);
A1 = appcoef2(c,s,'haar',1);

Использование wcodemat для пересмотра коэффициентов на основе их абсолютных значений. Отображение пересмотренных коэффициентов.

V1img = wcodemat(V1,255,'mat',1);
H1img = wcodemat(H1,255,'mat',1);
D1img = wcodemat(D1,255,'mat',1);
A1img = wcodemat(A1,255,'mat',1);

subplot(2,2,1)
imagesc(A1img)
colormap pink(255)
title('Approximation Coef. of Level 1')

subplot(2,2,2)
imagesc(H1img)
title('Horizontal Detail Coef. of Level 1')

subplot(2,2,3)
imagesc(V1img)
title('Vertical Detail Coef. of Level 1')

subplot(2,2,4)
imagesc(D1img)
title('Diagonal Detail Coef. of Level 1')

Извлеките приближение уровня 2 и коэффициенты детализации.

[H2,V2,D2] = detcoef2('all',c,s,2);
A2 = appcoef2(c,s,'haar',2);

Использование wcodemat для пересмотра коэффициентов на основе их абсолютных значений. Отображение пересмотренных коэффициентов.

V2img = wcodemat(V2,255,'mat',1);
H2img = wcodemat(H2,255,'mat',1);
D2img = wcodemat(D2,255,'mat',1);
A2img = wcodemat(A2,255,'mat',1);

figure
subplot(2,2,1)
imagesc(A2img)
colormap pink(255)
title('Approximation Coef. of Level 2')

subplot(2,2,2)
imagesc(H2img)
title('Horizontal Detail Coef. of Level 2')

subplot(2,2,3)
imagesc(V2img)
title('Vertical Detail Coef. of Level 2')

subplot(2,2,4)
imagesc(D2img)
title('Diagonal Detail Coef. of Level 2')

Этот пример показывает структуру wavedec2 выходные матрицы.

Загрузка и отображение изображения.

load woman
imagesc(X)
colormap(map)

Сохраните текущий режим расширения дискретного вейвлет.

origMode = dwtmode('status','nodisplay');

Изменение периодической обработки границ. The dwtmode функция отображает сообщение об изменении режима расширения DWT.

dwtmode('per')

 
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!  WARNING: Change DWT Extension Mode  !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
                                         
*****************************************
**  DWT Extension Mode: Periodization  **
*****************************************
                                         

Выполните разложение изображения уровня 3 с помощью db1 (Haar) вейвлет.

[c,s] = wavedec2(X,3,'db1');

Возвращает количество элементов в изображении X и вектор коэффициентов c. Подтвердите, что количество элементов в каждом равно.

numel(X)

ans = 65536

numel(c)

ans = 65536

Отобразите матрицу бухгалтерии s. В первой строке отображаются размерности грубой шкалы приближения изображения. В последней строке отображаются размерности оригинального изображения. Промежуточные строки отображают размерности коэффициентов детализации на трех уровнях разложения, исходя из грубой и мелкой шкалы.

s

s = 5×2

    32    32
    32    32
    64    64
   128   128
   256   256

Установите режим расширения дискретного вейвлет в исходный режим.

dwtmode(origMode,'nodisplay')

Обратитесь к разделу Поддержка GPU по версии (Parallel Computing Toolbox), чтобы узнать, какие графические процессоры поддерживаются.

Загрузка изображения. Поместите изображение на графический процессор с помощью gpuArray. Сохраните текущий режим расширения.

load mask
imgg = gpuArray(X);
origMode = dwtmode('status','nodisp');

Использование dwtmode для изменения режима внутреннего абонента на нулевое заполнение. Получите трехуровневый DWT изображения на графическом процессоре с помощью db4 вейвлет.

dwtmode('zpd','nodisp')
[c,s] = wavedec2(imgg,3,'db4');

Текущий режим расширения zpd не поддерживается для gpuArray вход. Поэтому DWT выполняется вместо этого с помощью sym режим расширения. Чтобы подтвердить это, установите режим внутреннего абонента на sym и возьмите DWT из noisdoppg, затем сравните с предыдущим результатом.

dwtmode('sym','nodisp')
[csym,ssym] = wavedec2(imgg,3,'db4');
[max(abs(c-csym)) max(abs(s-ssym))]

ans =

     0     0     0

Установите текущий режим расширения на per и получите трехуровневый DWT imgg. Режим расширения per поддерживается для gpuArray вход. Подтвердите, что результат отличается от sym результаты.

dwtmode('per','nodisp')
[cper,sper] = wavedec2(imgg,3,'db4');
[length(csym) ; length(cper)]

ans = 2×1

       71542
       65536

ssym

ssym = 5×2

    38    38
    38    38
    69    69
   131   131
   256   256

sper

sper = 5×2

    32    32
    32    32
    64    64
   128   128
   256   256

Восстановите режим расширения к исходной настройке.

dwtmode(origMode,'nodisp')