% odstraneni rozmazani obrazku pohybem svislem smeru

i=5;
im=double(imread('test.tif'));

%velikost obrazku
[I,J]=size(im);

%matice pro blur (rozmazani)
bl_matrix=zeros(I,J);
bl_matrix(1:i,1)=1;
bl_matrix(end-i+1:end,1)=1;

% fourierova transformace rozmazavaci matice
fftblur=fft2(bl_matrix); 
cnst=1e-5;

% vypocet filru, ktery odstrani rozmazani
invfilter=conj(fftblur)./(abs(fftblur).^2 + cnst);

% vypocet nerozmazaneho obrazku
im_deblur=abs(ifft2(fft2(im).*invfilter));
im_spect=abs(fftshift(fft2(im))).^2;
im_deblur_spect=abs(fftshift(fft2(im_deblur))).^2;

%zobrazeni vysledku
figure(1); imshow(im,[]); title('Puvodni obrazek');
figure(2); imshowl(im_spect);title(['Spektrum puvodniho obrazku']);
figure(3); imshow(im_deblur,[]); title('Obrazek s odstranenym rozmazanim');
figure(4); imshowl(im_deblur_spect); title(['Spektrum vysledneho obrazku']);s