我已经实现了 BackProjection 的代码并且它可以工作。现在我必须使用第一个 ramlak 过滤器和汉明窗口来实现过滤版本的代码。我的正弦图尺寸为 144x180。我的想法是
filter=zeros(1,144)';
filter(1:73)=[0:1/72:1];
filter(74:end)=[1-1/72:-1/72:1/72];
然后对于每个投影 s
filtered_proj=fft(s).*filter;
proj=real(ifft(filtered_proj));
然后执行反投影算法
你能告诉我它是否正确吗?
谢谢