matlab高斯1D模糊+差分（频域实现）

load Sinogram_Parallel
Sinogram=P.Projection;
[N_beam,N_angle]=size(Sinogram);


FWHM = 5; % in pixel
SIGMA = FWHM*Pixel_Size/2.0/sqrt(2.0*log(2.0));  % in physical dimension


L =2048;
Pixel_Size = 1e-3;
Fs = 1/Pixel_Size;
Kx = -Fs/2:Fs/(L*2):Fs/2-Fs/(L*2);   


F_diff = (1j*2*pi)*(Kx);
F_diff = fftshift(F_diff);


F_gauss = exp( - (2*pi^2)*SIGMA^2*(Kx.^2)) ;       
F_gauss = fftshift(F_gauss);
    
for i=1%1:360


    f = [Sinogram(:,i);zeros(L*2-N_beam,1)];        
    F_temp = fft(f);
    f_blur = ifft(F_temp.*F_gauss');
    f_blur_diff = ifft(F_temp.*F_gauss'.*F_diff');
        
end


%%


[temp,para]=phantom( );
head=