以下代码构建2个CT图像,然后进行重建过程,
%reconstruction of an phantom
img = phantom(256) % generate the phantom
maxvalimg=max(max(img)) %obtain maximum value of 2D array
minvalimg=min(min(img)) %obtain minimum value of 2d array
imshow(img,[minvalimg,maxvalimg]); %display the generated image
theta = 0:359 % generate angle vector
R = radon(img, theta) % generate the sin construction
maxvalR=max(max(R))
minvalR=min(min(R))
imshow(R,[minvalR,maxvalR]); % display the image
当代码运行时,第一个图像覆盖第二个图像,我如何分别获取这两个图像,(当然我可以单独运行代码块),但我需要在单次运行中运行它们,然后需要获取构建的图像)