RI.rar_The Image

%load woman; I=imread('bubble.tif'); [cA1,cH1,cV1,cD1] = dwt2(I,'db1'); A1 = upcoef2('a',cA1,'db1',1); H1 = upcoef2('h',cH1,'db1',1); V1 = upcoef2('v',cV1,'db1',1); D1 = upcoef2('d',cD1,'db1',1); colormap(map); subplot(2,3,1); image(wcodemat(A1,192)); title('Approximation A1') subplot(2,3,2); image(wcodemat(H1,192)); title('Horizontal Detail H1') subplot(2,3,3); image(wcodemat(V1,192)); title('Vertical Detail V1') subplot(2,3,4); image(wcodemat(D1,192)); title('Diagonal Detail D1') subplot(2,3,5); image(I); title('Original Image') te=0;re=0; e=0; f=0; g=0; h=0; [m n]=size(cA1); for i=1:1:m for j=1:1:n e= e +int32( cA1(i,j)*cA1(i,j)); f= f +int32( cH1(i,j)*cH1(i,j)); g= g +int32( cV1(i,j)*cV1(i,j)); h= h +int32( cD1(i,j)*cD1(i,j)); end end e f g h te=e+f+g+h; te s=size(I); cD1 = zeros(s, s, 'uint32'); X = idwt2(cA1,cH1,cV1,cD1,'db1',s); [m n]=size(X); for i=1:1:m for j=1:1:n re= re +int32( X(i,j)*X(i,j)); end end re subplot(2,3,6); image(X); title('Reconstructed Image') %figure,imhist(X)