Produce simple iradon geometry
3 views (last 30 days)
Show older comments
Below I have a very basic script that creates a very small circular geometry when using the radon transform of my data. I have my "data =" defined as a few values that would be the values from one projection, I just want to take that projection and rotate it around 360 degrees such that it would produce a circle when I use the iradon() funciton. When I do this, it produces what I think I'd like, but it is so very small and pixelated. Is there a way to make the image larger? It looks like I took a CT scan of a tiny, dense, pea-sized object, but I would like it to appear on the screen as a dense torso-sized image. Any help would be great.
Thank you.
data = [0,1,2,3,4,3,2,1,0]
R1 = radon(data)
figure(1)
imagesc(R1)
colormap(hot)
colorbar
figure(2)
I1 = iradon(R1,0:179,'linear','none')
imshow(I1,[])
0 Comments
Accepted Answer
Matt J
on 3 Jun 2022
Edited: Matt J
on 3 Jun 2022
data = [0,1,2,3,4,3,2,1,0];
theta=linspace(0,180,1000); theta(end)=[];
R1 = radon(data,theta);
figure(2)
N=numel(theta);
R2=repmat(R1(:,1),1,N);
for skip=[1000, 500, 250,125,60,1];
nexttile
I = iradon(R2(:,1:skip:end),theta(1:skip:end),'spline','none',1,30);
imshow(I,[])
end
2 Comments
More Answers (1)
Voss
on 2 Jun 2022
data = [0,1,2,3,4,3,2,1,0];
R1 = radon(data);
figure(1)
imagesc(R1)
colormap(hot)
colorbar
figure(2)
I1 = iradon(R1,0:179,'linear','none',20)
imshow(I1,[])
2 Comments
Voss
on 3 Jun 2022
All I can suggest is to increase the output_size in iradon further, to reduce the appearance of pixelation.
data = [0,1,2,3,4,3,2,1,0];
R1 = radon(data);
figure(1)
imagesc(R1)
colormap(hot)
colorbar
figure(2)
I1 = iradon(R1,0:179,'linear','none',64);
imshow(I1,[])
I'm not sure what those 8 different images all are nor what the numbers (1,2,4,8,16,32,64,180) mean. If you want to reproduce those 8 images, it seems like you'll have to use something besides (or in addition to) iradon.
See Also
Categories
Find more on 2-D and 3-D Plots in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!