%%
clear all
figure(1); subplot(2,1,1)

d   = double(imread('obj_observed.gif'));
psf = double(imread('psf_recovered_2.gif'));

figure(1); subplot(2,1,1)
image(d)
subplot(2,1,2)
image(psf)
psf = psf./(sum(sum(psf))); % Normalisierung

psf  = fftshift(psf);  % Was macht fftshift?