How to display second level dwt?
5 views (last 30 days)
Show older comments
Ömer Faruk GÜNAYDIN
on 23 Dec 2021
Commented: Imtiyaz Ahmad
on 7 Aug 2023
Hi. I want to display second level dwt. If I display seperatly to frequency components, it works fine. But when I try to merge and display, It appears just black color. Here is my code
close all;
clear;
cd 'C:\Users\LENOVO\OneDrive\Masaüstü' %% Uploading picture to Matlab
X=imread('Coleen.jpg','jpg');%Reading picture by Matlab
[LoD,HiD] = wfilters('haar','d');
[cA,cH,cV,cD] = dwt2(X,LoD,HiD,'mode','symh');
[cA2,cH2,cV2,cD2]=dwt2(cA,LoD,HiD,'mode','symh');
subplot(2,2,1)
imagesc(cA)
colormap gray
title('Approximation')
subplot(2,2,2)
imagesc(cH)
colormap gray
title('Horizontal')
subplot(2,2,3)
imagesc(cV)
colormap gray
title('Vertical')
subplot(2,2,4)
imagesc(cD)
colormap gray
title('Diagonal')
figure
Level2=[cA2,cH2; cV2,cD2];
imagesc([Level2,cH; cV,cD]);
colormap gray

It appears like this. But I want to display it like below.

Please. Could you help me?
0 Comments
Accepted Answer
yanqi liu
on 24 Dec 2021
Edited: yanqi liu
on 24 Dec 2021
yes,sir,may be use the dimension information,such as
close all;
clear all; clc;
%Reading picture by Matlab
X=imread('cameraman.tif');
[LoD,HiD] = wfilters('haar','d');
[cA,cH,cV,cD] = dwt2(X,LoD,HiD,'mode','symh');
[cA2,cH2,cV2,cD2]=dwt2(cA,LoD,HiD,'mode','symh');
figure
Level2=[mat2gray(cA2),mat2gray(cH2); mat2gray(cV2),mat2gray(cD2)];
d = [Level2,mat2gray(cH); mat2gray(cV),mat2gray(cD)];
% get dimension
rects1 = [1 1 size(cA2)
1 2 size(cA2)
2 1 size(cV2)
2 2 size(cD2)];
rects2 = [1 1 size(Level2)
1 2 size(cH)
2 1 size(cV)
2 2 size(cD)];
rects = {rects1, rects2};
imagesc(d);
colormap gray
axis equal
axis off
hold on;
rectangle('position', [1 1 size(d,2)-1 size(d,1)-1], 'EdgeColor', 'y', 'LineWidth', 2)
for k = 1 : length(rects)
rectsk = rects{k};
for i = 1 : size(rectsk, 1)
recti = [zeros(1, 2) rectsk(i,3:4)];
% get location
rowi = rectsk(i,1);
coli = rectsk(i,2);
for j = 1 : coli-1
recti(1) = recti(1) + rectsk(j, 4);
end
for j = 1 : rowi-1
recti(2) = recti(2) + rectsk(j, 3);
end
rectangle('position', recti, 'EdgeColor', 'y', 'LineWidth', 2)
end
end
1 Comment
Imtiyaz Ahmad
on 7 Aug 2023
Hi,
Dimensions of arrays being concatenated are not consistent.... but its working fine while using haar wavelet, can you provide the solution for that also ..
More Answers (0)
See Also
Categories
Find more on Image Analysis in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!