function [LL LH HL HH]=haar_dwt2D(img)
[m n]=size(img);
for i=1:m %每一行进行分解
%[L H]=haar_dwt(img(i,:));
p=length(img(i,:));
p=p/2;
L=zeros(1,p); %低频分量
H=zeros(1,p); %高频分量
for i1=1:p
L(i1)=(img(i,2*i1-1)+img(i,2*i1))/sqrt(2);
H(i1)=(img(i,2*i1-1)-img(i,2*i1))/sqrt(2);
end
img(i,:)=[L H];
end
for j=1:n %每一列进行分解
%[L H]=haar_dwt(img(:,j));
q=length(img(:,j));
q=q/2;
L=zeros(1,q); %低频分量
H=zeros(1,q); %高频分量
for j1=1:p
L(j1)=(img(2*j1-1,j)+img(2*j1,j))/sqrt(2);
H(j1)=(img(2*j1-1,j)-img(2*j1,j))/sqrt(2);
end
img(:,j)=[L H];
end
%本来分解不应该加mat2gray的,不过为了有好的显示效果就加上了
LL=mat2gray(img(1:m/2,1:n/2)); %行列都是低频
LH=mat2gray(img(1:m/2,n/2+1:n)); %行低频列高频
HL=mat2gray(img(m/2+1:m,1:n/2)); %行高频列低频
HH=mat2gray(img(m/2+1:m,n/2+1:n)); %行列都是高频
end
没有合适的资源?快使用搜索试试~ 我知道了~
小波变换练习Code
共6个文件
m:3个
jpg:2个
asv:1个
需积分: 9 1 下载量 94 浏览量
2016-04-26
15:16:05
上传
评论
收藏 139KB RAR 举报
温馨提示
主要实现Haar小波变换,源码有一些小bug,已经做了修改,可以完美实现;资源中包括lena图片和Code。
资源推荐
资源详情
资源评论
收起资源包目录
小波变换练习Code.rar (6个子文件)
picture
lena1.jpg 86KB
lena.jpg 50KB
code
main.m 445B
haar_dwt.m 285B
haar_dwt2D.asv 963B
haar_dwt2D.m 1KB
共 6 条
- 1
资源评论
i胡说
- 粉丝: 185
- 资源: 31
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功