mirdwt.m
上传用户:speoil
上传日期:2022-06-23
资源大小:224k
文件大小:5k
- 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 <lang@jazz.rice.edu>
- %
- %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 techtran@rice.edu or (713) 348-6173
- %
- %Change History:
- %
- %Modification #1
- %Mon Aug 7 15:09:51 CDT 1995
- %Rebecca Hindman <hindman@ece.rice.edu>
- %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 <felixf@rice.edu>