I want to use compressive sensing to reconstruct an image from fewer samples.
My problem is with Psi matrix which I want to be Biorthogonal wavelet coefficients but I don't know how to define it.
I have used Fourier basis and Haar wavelet and it worked well. Here is my code with Haar wavelet.
Can anyone tell me how to define Psi matrix as Biorthogonal wavelet transform?
Thanks in advance.
clc
close all
clear all
[fn,fp]=uigetfile({'*.*'});
tic
A=im2double(rgb2gray(imread([fp,fn])));
figure(1),imshow(A)
xlabel('original')
x=A(:);
n=length(x);
m=0.3*n;
Phi=randn(m,n); %Measurment Matrix
Psi=generate_haar(n); %sensing Matrix which is made by generate_haar function
y=Phi*x; %compressed signal
Theta=Phi*Psi;
%Initial Guess: y=Theta*s => s=Theta\y
s2=Theta\y;
%Solution
s1=OMP( Theta, y, 1e3);
%Reconstruction
x1=Psi*s1;
figure,imshow(reshape(x1,size(A))),xlabel('OMP')
toc
