mirdwt.m
上传用户:speoil
上传日期:2022-06-23
资源大小:224k
文件大小:5k
源码类别:

波变换

开发平台:

Matlab

  1. function [x,L] = mirdwt(yl,yh,h,L);
  2. %    function [x,L] = mirdwt(yl,yh,h,L);
  3. %    Function computes the inverse redundant discrete wavelet
  4. %    transform x  for a 1D or 2D input signal. (Redundant means here
  5. %    that the sub-sampling after each stage of the forward transform
  6. %    has been omitted.) yl contains the lowpass and yl the highpass
  7. %    components as computed, e.g., by mrdwt. In the case of a 2D
  8. %    signal, the ordering in
  9. %    yh is [lh hl hh lh hl ... ] (first letter refers to row, second
  10. %    to column filtering).  
  11. %
  12. %    Input:
  13. %       yl : lowpass component
  14. %       yh : highpass components
  15. %       h  : scaling filter
  16. %       L  : number of levels. In the case of a 1D signal, 
  17. %            length(yl) must  be divisible by 2^L;
  18. %            in the case of a 2D signal, the row and
  19. %            the column dimension must be divisible by 2^L.
  20. %   
  21. %    Output:
  22. %      x : finite length 1D or 2D signal
  23. %      L : number of levels
  24. %
  25. %  HERE'S AN EASY WAY TO RUN THE EXAMPLES:
  26. %  Cut-and-paste the example you want to run to a new file 
  27. %  called ex.m, for example. Delete out the % at the beginning 
  28. %  of each line in ex.m (Can use search-and-replace in your editor
  29. %  to replace it with a space). Type 'ex' in matlab and hit return.
  30. %
  31. %
  32. %    Example 1:
  33. %    xin = makesig('Leopold',8);
  34. %    h = daubcqf(4,'min');
  35. %    L = 1;
  36. %    [yl,yh,L] = mrdwt(xin,h,L);
  37. %    [x,L] = mirdwt(yl,yh,h,L)
  38. %    x = 0.0000 1.0000 0.0000 -0.0000 0 0 0 -0.0000
  39. %    L = 1
  40. %  
  41. %    Example 2:  
  42. %    load lena;
  43. %    h = daubcqf(4,'min');
  44. %    L = 2;
  45. %    [ll_lev2,yh,L] = mrdwt(lena,h,L); % lena is a 256x256 matrix
  46. %    N = 256;
  47. %    lh_lev1 = yh(:,1:N); 
  48. %    hl_lev1 = yh(:,N+1:2*N); 
  49. %    hh_lev1 = yh(:,2*N+1:3*N);
  50. %    lh_lev2 = yh(:,3*N+1:4*N); 
  51. %    hl_lev2 = yh(:,4*N+1:5*N); 
  52. %    hh_lev2 = yh(:,5*N+1:6*N);
  53. %    figure; colormap(gray); imagesc(lena); title('Original Image');
  54. %    figure; colormap(gray); imagesc(ll_lev2); title('LL Level 2');
  55. %    figure; colormap(gray); imagesc(hh_lev2); title('HH Level 2');
  56. %    figure; colormap(gray); imagesc(hl_lev2); title('HL Level 2');
  57. %    figure; colormap(gray); imagesc(lh_lev2); title('LH Level 2');
  58. %    figure; colormap(gray); imagesc(hh_lev1); title('HH Level 1');
  59. %    figure; colormap(gray); imagesc(hl_lev2); title('HL Level 1');
  60. %    figure; colormap(gray); imagesc(lh_lev2); title('LH Level 1');
  61. %    [lena_Hat,L] = mirdwt(ll_lev2,yh,h,L);
  62. %    figure; colormap(gray); imagesc(lena_Hat); 
  63. %                            title('Reconstructed Image');
  64. %
  65. %    See also: mdwt, midwt, mrdwt
  66. %
  67. %    Warning! min(size(yl))/2^L should be greater than length(h)
  68. %
  69. %File Name: mirdwt.m
  70. %Last Modification Date: 08/07/95 15:14:21
  71. %Current Version: mirdwt.m 2.4
  72. %File Creation Date: Wed Oct 19 10:51:58 1994
  73. %Author: Markus Lang  <lang@jazz.rice.edu>
  74. %
  75. %Copyright (c) 2000 RICE UNIVERSITY. All rights reserved.
  76. %Created by Markus Lang, Department of ECE, Rice University. 
  77. %
  78. %This software is distributed and licensed to you on a non-exclusive 
  79. %basis, free-of-charge. Redistribution and use in source and binary forms, 
  80. %with or without modification, are permitted provided that the following 
  81. %conditions are met:
  82. %
  83. %1. Redistribution of source code must retain the above copyright notice, 
  84. %   this list of conditions and the following disclaimer.
  85. %2. Redistribution in binary form must reproduce the above copyright notice, 
  86. %   this list of conditions and the following disclaimer in the 
  87. %   documentation and/or other materials provided with the distribution.
  88. %3. All advertising materials mentioning features or use of this software 
  89. %   must display the following acknowledgment: This product includes 
  90. %   software developed by Rice University, Houston, Texas and its contributors.
  91. %4. Neither the name of the University nor the names of its contributors 
  92. %   may be used to endorse or promote products derived from this software 
  93. %   without specific prior written permission.
  94. %
  95. %THIS SOFTWARE IS PROVIDED BY WILLIAM MARSH RICE UNIVERSITY, HOUSTON, TEXAS, 
  96. %AND CONTRIBUTORS AS IS AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, 
  97. %BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS 
  98. %FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL RICE UNIVERSITY 
  99. %OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, 
  100. %EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, 
  101. %PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; 
  102. %OR BUSINESS INTERRUPTIONS) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, 
  103. %WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR 
  104. %OTHERWISE), PRODUCT LIABILITY, OR OTHERWISE ARISING IN ANY WAY OUT OF THE 
  105. %USE OF THIS SOFTWARE,  EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
  106. %
  107. %For information on commercial licenses, contact Rice University's Office of 
  108. %Technology Transfer at techtran@rice.edu or (713) 348-6173
  109. %
  110. %Change History:
  111. %Modification #1
  112. %Mon Aug  7 15:09:51 CDT 1995
  113. %Rebecca Hindman <hindman@ece.rice.edu>
  114. %Added L to function line so that it can be displayed as an output
  115. %Modification #2  
  116. %Thursday Mar 2 2000
  117. % Added Example 2
  118. % Felix Fernandes <felixf@rice.edu>