%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Image decomposition by nonsubsampled contourlet transform (NSSC).% This is the iterated filter bank that computes the nonsubsampled% contourlet transform.  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Parameteters:nlevels = 4 ;        % Decomposition level%pfilter = 'pkva' ;              % Pyramidal filterdfilter = 'dmaxflat7'; %'cd' ;              % Directional filter% Nonsubsampled Contourlet decompositioncoeffs = nssdfbdec( double(im), dfilter, nlevels );disp( nlevels); disp(dfilter);% Display the coefficients%disp('Displaying the contourlet coefficients...') ;%shownssc( coeffs ) ;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Nonsubsampled Contourlet transform (NSSC) reconstruction.% This is the inverse of nsscdec, i.e.% imrec = nsscrec(coeffs, dfilter, pfilter);% would reconstruct imrec = im%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Reconstruct imageimrec = nssdfbrec( coeffs, dfilter ) ;% Show the reconstruction image and the original imageif 0figure;subplot(1,2,1), imagesc( im, [0, 255] ); title('Original image' ) ;colormap(gray);axis image off;subplot(1,2,2), imagesc( imrec, [0, 255] );title('Reconstructed image' ) ;colormap(gray);axis image off;endmse = sum( sum( (imrec - double(im)).^2 ) );mse = mse / prod(size(im));disp( sprintf('The mean square error is: %f', mse ) );disp(' ');