# Radon Transform reducing the number of parallel beams (MATLAB)

I am trying to understand different filters' and other things' effects on Radon Transform by using MATLAB. First I upload an image, then take its Radon transform with radon function and then reconstruct it with iradon function.

Now I wanted to see the difference between undersampling the projection data by reducing the number of parallel beams by a factor of 4 and reducing the number of rotation angles by a factor of 4. However according to what I understood from radon and iradon functions' descriptions, I can only change the number of rotation angles. Is not it true? If not, what is the way to do it?

Here is what I do to decrease the number of rotation angles by a factor of 4:

I = imread ('asdasd.tif');


Unfortunately you can't do this with the radon function as it comes. However, the function is written in MATLAB form (as opposed to being a compiled function), so you could extend the function yourself. This would also be a good way to learn more about the process.

To reconstruct properly, you would also need to perform edits on the iradon function.

An alternative way to get what you're looking for would be to get projections as they come from radon, then reduce the array size by averaging the projection rays in groups of 4. You would probably want the reconstruction to be the same size, so you'd need to repeat the averaged values 4 times in the argument to iradon. So maybe something like this:

p = phantom ();