h_projection.m

平行束CT重建算法

function Q=h_projection()

close all;
clear all;


Y_Length=40;
L=64;
N=128;
d=Y_Length/N;

p=0;
q=0;





H=zeros(L,N,10);
h1=ones(L,1);

H1=sub_hfilter(N,d);
H2=h1*ifftshift(H1);

for j=1:10
H(:,:,j)=H2;
end

load P_projection;

S=fftn(P);
Q=ifftn(S.*H);

filename='h_projection';
save(filename,'Q');