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

最新推荐文章于 2022-12-13 23:09:56 发布

wisdompeak

最新推荐文章于 2022-12-13 23:09:56 发布

阅读量949

点赞数

CC 4.0 BY-SA版权

分类专栏： Matlab 文章标签： fft matlab

本文链接：https://blog.youkuaiyun.com/wisdompeak/article/details/43231585

这篇博客介绍了如何在MATLAB中使用频域方法进行1D高斯模糊和差分操作。通过设置高斯滤波器的标准差，计算频域中的高斯函数和差分矩阵，然后对投影数据应用快速傅里叶变换(FFT)进行模糊和差分。实验证明，这种方法与直接差分的结果一致。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

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=phantom(para,2048);
imagesc(head);