CT backprojection of a sinogram without using iradon

20 次查看(过去 30 天)
Hi,
Can somebody please give me a reference code on how to perform simple backprojection of a sinogram without using IRADON ???

回答(1 个)

Matt J
Matt J 2013-11-20
  2 个评论
ram
ram 2013-11-20
编辑:Matt J 2013-11-20
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')

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Biomedical Imaging 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by