% simple example for FFST
% computes the shearlet transform of some geometric image
%
%--------------------------------------------------------------------------
% 2012-01-20, v1.0, (c) Sören Häuser
addpath ./create
addpath ./helper
addpath ./shearlets
% create image
%A = myPicture2();
A = double(imread('../lena.tif'));

% shearlet transform
[ST, Psi] = shearletTransformSpect(A);

% inverse shearlet transform
C = inverseShearletTransformSpect(ST,Psi);

num_shear = 43;
% plot results
subplot(2,2,1)
imagesc(A)
axis image off
colormap(gray)
title('original image')

subplot(2,2,2)
imagesc(abs(ST(:,:,num_shear))) %18
axis image off
colormap(gray)
title('shearlet coefficients')

X = Psi(:,:,num_shear);
subplot(2,2,3)
imagesc(X); %Psi(:,:,10))
axis image off
colormap(gray)
title('shearlet')

%modyfikacje

X = ifft2(fftshift(fftshift(Psi,1),2)); %ifft2(fftshift(X));
theta = 0:179;%theta = 0:179;
[R,xp] = radon(X(:,:,num_shear),theta);
subplot(2,2,4)
imagesc(R)
axis image off
colormap(gray)
title('radon-shearlet image')
% koniec

% subplot(2,2,4)
% imagesc(C)
% axis image off
% colormap(gray)
% title('reconstructed image')