Produce simple iradon geometry
2 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
Brett Ruben
am 2 Jun. 2022
Bearbeitet: Matt J
am 17 Aug. 2022
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 Kommentare
Akzeptierte Antwort
Matt J
am 3 Jun. 2022
Bearbeitet: Matt J
am 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 Kommentare
Weitere Antworten (1)
Voss
am 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 Kommentare
Voss
am 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.
Siehe auch
Kategorien
Mehr zu Image Acquisition Support Packages for Hardware Adaptors (Generic Video Interface) finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!