MATLAB2 维小波变换经典程序
% FWT_DB.M;
% 此示意程序用 DWT 实现二维小波变换
%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;clc;
T=256; % 图像维数
SUB_T=T/2; % 子图维数
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 1.调原始图像矩阵
load wbarb; % 下载图像
f=X; % 原始图像
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 2.进行二维小波分解
l=wfilters('db10','l'); % db10（消失矩为 10)低通分解滤波器冲击响应（长度为 20）
L=T-length(l);
l_zeros=[l,zeros(1,L)]; % 矩阵行数与输入图像一致，为 2 的整数幂
h=wfilters('db10','h'); % db10（消失矩为 10)高通分解滤波器冲击响应（长度为 20）
h_zeros=[h,zeros(1,L)]; % 矩阵行数与输入图像一致，为 2 的整数幂
for i=1:T; % 列变换
row(1:SUB_T,i)=dyaddown( ifft( fft(l_zeros).*fft(f(:,i)') ) ).'; % 圆周卷积<->FFT
row(SUB_T+1:T,i)=dyaddown( ifft( fft(h_zeros).*fft(f(:,i)') ) ).'; % 圆周卷积<->FFT
end;
for j=1:T; % 行变换
line(j,1:SUB_T)=dyaddown( ifft( fft(l_zeros).*fft(row(j,:)) ) ); % 圆周卷积<->FFT
line(j,SUB_T+1:T)=dyaddown( ifft( fft(h_zeros).*fft(row(j,:)) ) ); % 圆周卷积<->FFT
end;
decompose_pic=line; % 分解矩阵
% 图像分为四块
lt_pic=decompose_pic(1:SUB_T,1:SUB_T); % 在矩阵左上方为低频分量--fi(x)*fi(y)
rt_pic=decompose_pic(1:SUB_T,SUB_T+1:T); % 矩阵右上为--fi(x)*psi(y)
lb_pic=decompose_pic(SUB_T+1:T,1:SUB_T); % 矩阵左下为--psi(x)*fi(y)
rb_pic=decompose_pic(SUB_T+1:T,SUB_T+1:T); % 右下方为高频分量--psi(x)*psi(y)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 3.分解结果显示
figure(1);
colormap(map);
subplot(2,1,1);
image(f); % 原始图像
资源评论
hhappy0123456789
- 粉丝: 72
- 资源: 5万+
最新资源
- 伯克利大学机器学习-12Time series&sequential hypothesis testing&anomaly de
- 上市公司企业网络舆论关注数据集
- 伯克利大学机器学习-11Bootstrap&cross-validation&ROC plots Michael Jordan
- java项目工时统计成本核算管理系统源码数据库 MySQL源码类型 WebForm
- arnabdhar YOLOv8-Face-Detection Dateset Source
- CSnet-tls1.3-packet数据集
- 创建vue+electron项目流程
- Python-基于Pygame的贪吃蛇
- C#ASP.NET高校移动考勤webapp源码数据库 SQL2008源码类型 WebForm
- （2000-2023年）中国各、省、市、县、乡镇基尼系数数据（全新整理）
资源上传下载、课程学习等过程中有任何疑问或建议，欢迎提出宝贵意见哦~我们会及时处理！
点击此处反馈
