Radon Transform of image

I want radon transform in this steps: 1. Use CT scan image. 2. apply X an Y projections on image 3. Rotate image by 1 degree. 4. Apply Projection by X axis and Y axis on image. 5. Repeat II to IV steps up-to 360 degree.

Using built in radon transform function, I only rotate image with degree but how can I apply all this steps.

I use these code: %1. Create a Shepp-Logan head phantom image. P = phantom(512); figure,imshow(P),title('Original image');

%2. Compute the Radon transform of the phantom brain for three different sets of theta values. R1 has %18 projections, R2 has 36 projections, and R3 has 90 projections. theta1 = 0:5:360; [R1,xp] = radon(P,theta1); theta2 = 0:2:358; [R2,xp] = radon(P,theta2); theta3 = 0:1:359; [R3,xp] = radon(P,theta3);

%3.Display a plot of one of the Radon transforms of the Shepp-Logan head phantom. The following %figure shows R3, the transform with 180 projections. figure, imagesc(theta3,xp,R3); colormap(hot); colorbar, title('Radon transfomr of 360 degree projections'); xlabel('\theta'); ylabel('x\prime');

%4. Reconstruct the head phantom image from the projection data created in step 2 and display the results. I1 = iradon(R1,5); I2 = iradon(R2,2); I3 = iradon(R3,1);

figure, imshow(I1),title('reconstruct of image with 36 degree projections'); figure, imshow(I2),title('reconstruct of image with 180 degree projections'); figure, imshow(I3), title('reconstruct of image with 360 degree projections');

Kindly give me answer, how can I provide projections?

Thanks