Produce simple iradon geometry
7 次查看(过去 30 天)
显示 更早的评论
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 个评论
采纳的回答
Matt J
2022-6-3
编辑:Matt J
2022-6-3
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 个评论
更多回答(1 个)
Voss
2022-6-2
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 个评论
Voss
2022-6-3
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.
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!