Produce simple iradon geometry
3 ビュー (過去 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!