0
$\begingroup$

I am using matlab image processing , I have a 3DBodyphantom image of a CT scan which displays the axial,coronal and saggital slices, I am to apply to the image the radon function to get a sinogram and then the iradon function, I keep getting an error message while trying to apply the radon function within the loop.

This is the code for viewing the slices, how to apply the radon and iradon is what I am trying to figure out, is there a way to apply it or how am I suppose to apply it? My theta for applying the radon is that projections are taken every 2 degrees of rotation. This is the code I am using;

dim = [256 256 131]; %[256 256 116];
fid = fopen('3DBodyPhantom.bin'); %fopen('output42_atn_1.bin');
x= fread(fid, 'float32');
y=reshape(x, dim);

%to view sagittal slices
for i=1:dim(1)   
    imagesc(squeeze(y(i,:,:))); 
    pause(0.05); 
end

%to view coronal slices
for i=1:dim(2)
  img1=  imagesc(squeeze(y(:,i,:))); 
    pause(0.05); 
end

%to view axial slices
for i=1:dim(3)
    imagesc(squeeze(y(:,:,i))); 
    pause(0.05); 
end

fclose(fid);
$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.