matlab怎样滤波器,MATLAB如何在频域实现Ram

matlab怎样滤波器,MATLAB如何在频域实现Ram

你列出的公式是：如果你想要做逆Radon变换而不会在傅立叶域中滤波的中间结果。另一种方法是在空间域中使用卷积来完成整个滤波反投影算法，这在40年前可能更快;你最终会重新发布你发布的公式。但是，我现在不会推荐它，特别是不适合你的第一次重建。你应该首先理解希尔伯特变换。

总之，这里的一些Matlab代码它执行强制性Shepp - 洛根幻像滤波反投影重建。我展示了如何使用Ram-Lak过滤器进行自己的过滤。如果我真的有动力，我会用一些interp2命令和总结替换氡/ iradon。

phantomData=phantom();

N=size(phantomData,1);

theta = 0:179;

N_theta = length(theta);

[R,xp] = radon(phantomData,theta);

% make a Ram-Lak filter. it's just abs(f).

N1 = length(xp);

freqs=linspace(-1, 1, N1).';

myFilter = abs(freqs);

myFilter = repmat(myFilter, [1 N_theta]);

% do my own FT domain filtering

ft_R = fftshift(fft(R,[],1),1);

filteredProj = ft_R .* myFilter;

filteredProj = ifftshift(filteredProj,1);

ift_R = real(ifft(filteredProj,[],1));

% tell matlab to