dynamic display of radon transform

Hi everyone I have a problem: I'm trying to do a dynamic display of a ct-reconstruction. Basically what I'd like to do is to display on the same image the iradon for different angles. My problem is the following:

im = imread('image');

p = radon(im,[0:180]);

imshow(iradon(p,[0:180],'none')

this displays all the projection in one figure and this projection are at different rotation. But I'd like to display them first the one at 0 then on the same image, keeping the ones at 0 and adding also the one at 1 degree and so one. I tried with the following code but I got a single figure for each angle:

figure,

for theta = 0:180

p = radon(im,theta);

imshow(iradon(p,theta,'non'),[]);

end

[Rik]

So you want to animate the backprojection? I believe iradon has optimizations making it impossible to see intermediate steps (similar to how conv works). You can implement something yourself.

[LORENZO D'AMICO]

I solved the problem something like this:

theta = []
figure,
for i = 0:20:180
    theta = [theta,i]
    p = radon (im,theta)
    imshow(iradon(p,theta,'none'))
end