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);