zhijienitoux.m

来自「几个图象处理程序」· M 代码 · 共 27 行

%直接逆投的校正(逆投滤波）
I=phantom(256);
th=1:180;
G=radon(I,th);%G矩阵即投影函数,经过radon后，为一367×180的矩阵

F=zeros(256,256);
FF=zeros(256,256);
x=-1:2/255:1;
y=1:-2/255:-1;
clf
for m=1:256
    for n=1:256
        for th=1:180;
            R=(x(m)*cos(th*pi/180)+y(n)*sin(th*pi/180));
            r=round(R/1.414*183)+184;
            F(n,m)=F(n,m)+G(r,th);           
        end
    end
end
F1=fftshift(fft2(F));
for m=1:256
    for n=1:256
        FF(n,m)=F1(n,m)*sqrt((x(n)-128).^2+(y(m)+128).^2);
    end
end
ff=ifft2(fftshift(FF));
imagesc(real(ff))