Radon Transformation - Digital Signal Processing for Medical Imaging Using Matlab - page 9

Fig. 1.7 Illustration of filtered parallel beam backprojected images obtained for θ = 0 , 20 . 0557 ,

4457 ◦ (Note that parallel-beam

40 . 1114 , 60 . 1617 , 80 . 2228 , 100 . 2786 , 120 . 3343 , 140 . 3900 , 160 .

ranging from 0 to 180 ◦ with the resolution of 0 5014 ◦ )

reconstructed images are computed for all θ .

temp=imrotate(temp,90);
D=imrotate(temp,-theta,'nearest','crop');
RECONSTRUCTEDIMAGE=RECONSTRUCTEDIMAGE+D;
if(Z(z)==k)
SNAPSHOT_RECONSTRUCTED{z}=RECONSTRUCTEDIMAGE;
SNAPSHOP_DATA{z}=D;
z=z+1;
end
S=S+1;
end
figure
colormap(gray(256))
imagesc(RECONSTRUCTEDIMAGE)