CT backprojection of a sinogram without using iradon

Hi,

Can somebody please give me a reference code on how to perform simple backprojection of a sinogram without using IRADON ???

Answer 1

[Matt J]

http://sourceforge.net/apps/mediawiki/stir/index.php?title=Main_Page

[ram]

Hi,

Thank you for your reply.I have the following code which takes only a single backprojection and i wish to take multiple projections.

I would be really grateful if you could please figure out the problem.

I have a projection data and the following code which performs the sinogram and takes only a single backprojected view.

the code is as follows:--

    clc

    load proja(1)

        
    colormap(bone);

    imagesc(proj)

    axis('image')

    Num_of_projs = size(proj,1)

    Image_size = size(proj,2)

    % backproject projection data

    bpimage = zeros(Image_size,Image_size);  % Result image initialization

    delta = pi / Num_of_projs;

        
    view=1          %%Backprojection of first view

    for iy=1:Image_size      

     for ix=1:Image_size  %% project many views.
      val = proj(view,ix);

      bpimage(ix,iy) = bpimage(ix,iy) + val;
     end
    end

    % normalize and display result

    bpimage = bpimage * delta;

    figure(2)

    colormap(bone)

    imagesc(bpimage)

    axis('image')

[Rifqi Naufal]

What is (proja1) ??