mirdwt.m

在国外网站上找到的一个小波降噪的实例

function [x,L] = mirdwt(yl,yh,h,L);
%    function [x,L] = mirdwt(yl,yh,h,L);
% 
%    Function computes the inverse redundant discrete wavelet
%    transform x  for a 1D or 2D input signal. (Redundant means here
%    that the sub-sampling after each stage of the forward transform
%    has been omitted.) yl contains the lowpass and yl the highpass
%    components as computed, e.g., by mrdwt. In the case of a 2D
%    signal, the ordering in
%    yh is [lh hl hh lh hl ... ] (first letter refers to row, second
%    to column filtering).  
%
%    Input:
%       yl : lowpass component
%       yh : highpass components
%       h  : scaling filter
%       L  : number of levels. In the case of a 1D signal, 
%            length(yl) must  be divisible by 2^L;
%            in the case of a 2D signal, the row and
%            the column dimension must be divisible by 2^L.
%   
%    Output:
%	     x : finite length 1D or 2D signal
%	     L : number of levels
%
%  HERE'S AN EASY WAY TO RUN THE EXAMPLES:
%  Cut-and-paste the example you want to run to a new file 
%  called ex.m, for example. Delete out the % at the beginning 
%  of each line in ex.m (Can use search-and-replace in your editor
%  to replace it with a space). Type 'ex' in matlab and hit return.
%
%
%    Example 1:
%    xin = makesig('Leopold',8);
%    h = daubcqf(4,'min');
%    L = 1;
%    [yl,yh,L] = mrdwt(xin,h,L);
%    [x,L] = mirdwt(yl,yh,h,L)
%    x = 0.0000 1.0000 0.0000 -0.0000 0 0 0 -0.0000
%    L = 1
%  
%    Example 2:  
%    load lena;
%    h = daubcqf(4,'min');
%    L = 2;
%    [ll_lev2,yh,L] = mrdwt(lena,h,L); % lena is a 256x256 matrix
%    N = 256;
%    lh_lev1 = yh(:,1:N); 
%    hl_lev1 = yh(:,N+1:2*N); 
%    hh_lev1 = yh(:,2*N+1:3*N);
%    lh_lev2 = yh(:,3*N+1:4*N); 
%    hl_lev2 = yh(:,4*N+1:5*N); 
%    hh_lev2 = yh(:,5*N+1:6*N);
%    figure; colormap(gray); imagesc(lena); title('Original Image');
%    figure; colormap(gray); imagesc(ll_lev2); title('LL Level 2');
%    figure; colormap(gray); imagesc(hh_lev2); title('HH Level 2');
%    figure; colormap(gray); imagesc(hl_lev2); title('HL Level 2');
%    figure; colormap(gray); imagesc(lh_lev2); title('LH Level 2');
%    figure; colormap(gray); imagesc(hh_lev1); title('HH Level 1');
%    figure; colormap(gray); imagesc(hl_lev2); title('HL Level 1');
%    figure; colormap(gray); imagesc(lh_lev2); title('LH Level 1');
%    [lena_Hat,L] = mirdwt(ll_lev2,yh,h,L);
%    figure; colormap(gray); imagesc(lena_Hat); 
%                            title('Reconstructed Image');
%
%    See also: mdwt, midwt, mrdwt
%
%    Warning! min(size(yl))/2^L should be greater than length(h)
%

%File Name: mirdwt.m
%Last Modification Date: 08/07/95	15:14:21
%Current Version: mirdwt.m	2.4
%File Creation Date: Wed Oct 19 10:51:58 1994
%Author: Markus Lang  <[e-mail]>
%
%Copyright (c) 2000 RICE UNIVERSITY. All rights reserved.
%Created by Markus Lang, Department of ECE, Rice University. 
%
%This software is distributed and licensed to you on a non-exclusive 
%basis, free-of-charge. Redistribution and use in source and binary forms, 
%with or without modification, are permitted provided that the following 
%conditions are met:
%
%1. Redistribution of source code must retain the above copyright notice, 
%   this list of conditions and the following disclaimer.
%2. Redistribution in binary form must reproduce the above copyright notice, 
%   this list of conditions and the following disclaimer in the 
%   documentation and/or other materials provided with the distribution.
%3. All advertising materials mentioning features or use of this software 
%   must display the following acknowledgment: This product includes 
%   software developed by Rice University, Houston, Texas and its contributors.
%4. Neither the name of the University nor the names of its contributors 
%   may be used to endorse or promote products derived from this software 
%   without specific prior written permission.
%
%THIS SOFTWARE IS PROVIDED BY WILLIAM MARSH RICE UNIVERSITY, HOUSTON, TEXAS, 
%AND CONTRIBUTORS AS IS AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, 
%BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS 
%FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL RICE UNIVERSITY 
%OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, 
%EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, 
%PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; 
%OR BUSINESS INTERRUPTIONS) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, 
%WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR 
%OTHERWISE), PRODUCT LIABILITY, OR OTHERWISE ARISING IN ANY WAY OUT OF THE 
%USE OF THIS SOFTWARE,  EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
%
%For information on commercial licenses, contact Rice University's Office of 
%Technology Transfer at [e-mail] or [phone]
%
%Change History:
% 
%Modification #1
%Mon Aug  7 15:09:51 CDT 1995
%Rebecca Hindman <[e-mail]>
%Added L to function line so that it can be displayed as an output
% 
%Modification #2  
%Thursday Mar 2 2000
% Added Example 2
% Felix Fernandes <[e-mail]>