[转帖]二代小波harr变换，整数小波变换

happy

1. %文件名：liftwavedec2.m
   
2. %编 写：郭林庚
   
3. %时 间：2005/11/01
   
4. %函数功能：二代小波harr变换，整数小波变换
   
5. %输入格式举例：imgwave=liftwavedec2(image,256,3)
   
6. %参数说明：
   
7. % image--输入的图像矩阵，要为方阵
   
8. % m --输入的图像矩阵大小
   
9. % n --小波变换次数
   
10. %测试用例：
    
11. % img=imread('lena.jpg');
    
12. % [m,m]=size(img);
    
13. % imgwave=liftwavedec2(img,m,3);
    
14. % imshow(imgwave);
    
15. 
    
16. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
17. function imgwave=liftwavedec2(image,m,n)
    
18. img=image;
    
19. M=m;
    
20. for i=1:n
    
21. imgwave1=lwavedec2(img,M);
    
22. imgwave(1:M,1:M)=imgwave1;
    
23. M=M/2;
    
24. img=imgwave1(1:M,1:M);
    
25. end
    
26. 
    
27. %
    
28. % 二代小波harr变换，整数小波变换
    
29. %
    
30. function f_row=lwavedec2(image,N)
    
31. f=image;
    
32. T=N/2; % 子图像维数
    
33. 
    
34. 
    
35. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
36. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
37. %%%% 1.列变换
    
38. 
    
39. % A.分裂（奇偶分开）
    
40. 
    
41. f1=f([1:2:N-1],:); % 奇数
    
42. f2=f([2:2:N],:); % 偶数
    
43. 
    
44. % f1(:,T+1)=f1(:,1); % 补列
    
45. % f2(T+1,:)=f2(1,:); % 补行
    
46. 
    
47. % B.预测
    
48. 
    
49. for i_hc=1:T;
    
50. high_frequency_column(i_hc,:)=f1(i_hc,:)-f2(i_hc,:);
    
51. end;
    
52. 
    
53. % high_frequency_column(T+1,:)=high_frequency_column(1,:); % 补行
    
54. 
    
55. % C.更新
    
56. 
    
57. for i_lc=1:T;
    
58. low_frequency_column(i_lc,:)=f2(i_lc,:)+1/2*high_frequency_column(i_lc,:);
    
59. end;
    
60. 
    
61. % D.合并
    
62. f_column([1:1:T],:)=low_frequency_column([1:T],:);
    
63. f_column([T+1:1:N],:)=high_frequency_column([1:T],:);
    
64. 
    
65. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
66. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
67. %%%% 2.行变换
    
68. 
    
69. % A.分裂（奇偶分开）
    
70. 
    
71. f1=f_column(:,[1:2:N-1]); % 奇数
    
72. f2=f_column(:,[2:2:N]); % 偶数
    
73. 
    
74. 
    
75. % f2(:,T+1)=f2(:,1); % 补行
    
76. 
    
77. % B.预测
    
78. 
    
79. for i_hr=1:T;
    
80. high_frequency_row(:,i_hr)=f1(:,i_hr)-f2(:,i_hr);
    
81. end;
    
82. 
    
83. % high_frequency_row(:,T+1)=high_frequency_row(:,1); % 补行
    
84. 
    
85. % C.更新
    
86. 
    
87. for i_lr=1:T;
    
88. low_frequency_row(:,i_lr)=f2(:,i_lr)+1/2*high_frequency_row(:,i_lr);
    
89. end;
    
90. 
    
91. % D.合并
    
92. f_row(:,[1:1:T])=low_frequency_row(:,[1:T]);
    
93. f_row(:,[T+1:1:N])=high_frequency_row(:,[1:T]);

复制代码

[ 本帖最后由 suffer 于 2006-10-9 20:04 编辑 ]

happy

1. function imgwave=liftwaverec2(image,m,n)
   
2. M=m/(2^n);
   
3. imgwave=image;
   
4. for i=n:-1:1
   
5. M=M*2;
   
6. img=imgwave(1:M,1:M);
   
7. imgwave1=lwaverec2(img,M);
   
8. imgwave(1:M,1:M)=imgwave1;
   
9. end
   
10. 
    
11. 
    
12. function f_column=lwaverec2(f_row,N);
    
13. T=N/2;
    
14. %figure(2);
    
15. %subplot(221),imshow(f_row),title('变换后图像');
    
16. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
17. %%%%%%%%%%%%%%%%%%%反变换%%%%%%%%%%%%%%%%%%%%%%%%%
    
18. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
19. 
    
20. 
    
21. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
22. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
23. %%%% 1.行变换
    
24. 
    
25. % A.提取（低频高频分开）
    
26. 
    
27. f1=f_row(:,[T+1:1:N]); % 奇数
    
28. f2=f_row(:,[1:1:T]); % 偶数
    
29. 
    
30. 
    
31. % f2(:,T+1)=f2(:,1); % 补行
    
32. 
    
33. % B.更新
    
34. 
    
35. for i_lr=1:T;
    
36. low_frequency_row(:,i_lr)=f2(:,i_lr)-1/2*f1(:,i_lr);
    
37. end;
    
38. 
    
39. % C.预测
    
40. 
    
41. for i_hr=1:T;
    
42. high_frequency_row(:,i_hr)=f1(:,i_hr)+low_frequency_row(:,i_hr);
    
43. end;
    
44. 
    
45. % high_frequency_row(:,T+1)=high_frequency_row(:,1); % 补行
    
46. 
    
47. 
    
48. % D.合并(奇偶分开合并)
    
49. f_row(:,[2:2:N])=low_frequency_row(:,[1:T]);
    
50. f_row(:,[1:2:N-1])=high_frequency_row(:,[1:T]);
    
51. 
    
52. %subplot(223),imshow(f_row),title('行变换图像');
    
53. 
    
54. 
    
55. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
56. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
57. %%%% 2.列变换
    
58. 
    
59. % A.提取（低频高频分开）
    
60. 
    
61. f1=f_row([T+1:1:N],:); % 奇数
    
62. f2=f_row([1:1:T],:); % 偶数
    
63. 
    
64. % f1(:,T+1)=f1(:,1); % 补列
    
65. % f2(T+1,:)=f2(1,:); % 补行
    
66. 
    
67. % B.更新
    
68. 
    
69. for i_lc=1:T;
    
70. low_frequency_column(i_lc,:)=f2(i_lc,:)-1/2*f1(i_lc,:);
    
71. end;
    
72. 
    
73. % C.预测
    
74. 
    
75. for i_hc=1:T;
    
76. high_frequency_column(i_hc,:)=f1(i_hc,:)+low_frequency_column(i_hc,:);
    
77. end;
    
78. 
    
79. % high_frequency_column(T+1,:)=high_frequency_column(1,:); % 补行
    
80. 
    
81. % D.合并(奇偶分开合并）
    
82. f_column([2:2:N],:)=low_frequency_column([1:T],:);
    
83. f_column([1:2:N-1],:)=high_frequency_column([1:T],:);
    
84. 
    
85. %subplot(224),imshow(f_column),title('行变换图像');

复制代码

[ 本帖最后由 suffer 于 2006-10-9 20:05 编辑 ]

graduate

一楼代码有问题吗 ？？

范科举

我在厌学论坛上看到过这个代码