GPU calculation : step 2

GPU calculation becomes interesting when you have to run the same calculation a large number of times. Our aim is to split the iradon transform of one sinogram in many identical tasks.

If we manage to do that, we will be able to run these task in parallel on the GPU and, we hope, gain calculation time.

In one sinogram, each row corresponds to a given angle in the imaging plane. We will create a sinogram of a phantom for 18 projections angles and compute a reconstructed image using the inverse radon transform.

Fig .2 : Original phantom (left), reconstructed image (right)

Now we will try to do the same, but we will perform one inverse radon transform per projection angle.

Here are the results we get :

Fig .3 : Reconstructed images for different theta values

And after calculating a sum of the images :

Fig .4 : Sum of the images from Fig3

The next step will be to run in parallel the reconstruction for different theta values and see if this as an impact on our calculation time.