%Test Radontransformation und Berechnung über FFT (Fourier Slice Theorem)
%
%Alexander Liebhold 2009
%
%Input file enthält ein Testbild der Größe 50x50x3 Pixel (Grauwert)


%Bild laden
I=imread('C:\test.png')
%Neue Ausgabeform
figure(1)
%Bild plotten
imagesc(I)
title("Input")

%Intensitätsbild aus Grauwertbild entnehmen (Alle 3 Farbkanäle sind bei Grauwert gleich)
im(1:50,1:50,1)=I(1:50,1:50,1)

%2-D FFT berechnen
FT=fft2(im)
%Neue Ausgabeform
figure(2)
%Betragsspektrum
A=abs(FT)
%Phasenspektrum
phase=angle(FT)
%Ausgabe
subplot(1,2,1)
imagesc(A)
title("Betragsspektrum")
subplot(1,2,2)
imagesc(phase)
title("Phasenspektrum")

%Ausgabe in Polarkoordinaten
%Neue Ausgabeform
figure(3)
polar(phase,A)


%Berechnung Radontransformation
%Neue Ausgabeform
figure(4)
R=radon(im)
imagesc(R)
title("Radontransformation")
%Farbwertetabelle hinzufügen
colorbar