clc;
clear;
close all;
warning off;
addpath(genpath(pwd));
I=imread('peppers_gray.bmp');
I=uint8(I)
subplot(131)
imshow(I)
title('Original Image')
dim=size(I)
N=dim(1)
n=log2(N)
%Compute the Slant Transform Matrix
S=1/sqrt(2)*[1 1; 1 -1];
a=1
for i=2:n
b=1/sqrt(1+4*a^2)
a=2*b*a
if i==2
%%
%
% $$e^{\pi i} + 1 = 0$$
%
% $$e^{\pi i} + 1 = 0$$
%
%
S=(1/sqrt(2))*[[1 0; a b] [1 0;-a b];
[0 1;-b a] [0 -1;b,a]]*[S zeros(size(S));zeros(size(S)),S]
else
S=(1/sqrt(2))*[[1 0;a b] [zeros(2,(2^i-4)/2)] [1 0;-a b] [zeros(2,(2^i-4)/2)];
[zeros((2^i-4)/2,2)] [eye((2^i-4)/2)] [zeros((2^i-4)/2,2)] [eye((2^i-4)/2)];
[0 1;-b a] [zeros(2,(2^i-4)/2)] [0 -1;b,a] [zeros(2,(2^i-4)/2)];
[zeros((2^i-4)/2,2)] [eye((2^i-4)/2)] [zeros((2^i-4)/2,2)] [-eye((2^i-4)/2)] ]*[S zeros(size(S));zeros(size(S)),S]
end
end
%Compute Slant transform
S_t=S*double(I)*S';
subplot(132)
imshow(uint8(S_t),[])
title('Slant Transform ')
%Image Reconstruction
r_i=S'*S_t*S
subplot(133)
imshow(r_i,[])
title('Reconstructed Image')
matlab-图像的Haar变换,hadamard变换,KL变换,slant变换matlab对比仿真-源码
版权申诉
174 浏览量
2021-09-15
21:02:31
上传
评论
收藏 205KB RAR 举报
mYlEaVeiSmVp
- 粉丝: 1889
- 资源: 19万+
最新资源
- 随机森林RF葡萄酒分类,随机森林二分类(代码完整,数据齐全)
- 香蕉新鲜度等级识别香蕉成熟度识别分类数据集1186张3类别.7z
- 香蕉成熟度识别分类数据集13478张4类别.7z
- 安卓Android图书馆座位占座app设计可导入Studio源码案例设计.zip
- 安卓Android汽车票售票订票app设计可导入Studio源码案例设计.zip
- 安卓android高校开放式实验室管理app可导入Studio源码案例设计.zip
- 安卓Android车辆运输信息管理App可升级Studio源码案例设计.zip
- Android Studio移动应用开发.zip
- android studio-GIT.zip
- 优先队列式分支限界法解决旅行商问题.zip
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈