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

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

如果您想在傅里叶域中进行无需过滤的反Radon变换,则列出的公式是中间结果.另一种方法是使用空间域中的卷积来完成整个滤波反投影算法,这可能在40年前更快;你最终会重新发布你发布的公式.但是,我现在不建议,特别是不适合你的第一次重建;你应该先了解希尔伯特变换.

无论如何,这里有一些Matlab代码,它执行强制性的Shepp-Logan模拟滤波反投影重建.我将展示如何使用Ram-Lak过滤器进行自己的过滤.如果我真的很有动力,我会用一些interp2命令和总结来代替radon / iradon.

phantomData =幻象();

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 do inverse FBP without a filter