I am starting to study 3D reconstruction from multiple CT slices. I used the following coding to create a 3D array :
X = zeros(512,512,n); for i = 1:n fileName = ['D:\Image Folder\', srcFiles(i).name]; I=imread(fileName); X(:,:,i) = I; end
I would like to know just stacking 2D images is enough to create a 3D image or not. Do I need to use some algorithms e.g – interpolation (I am not sure)? And for visualization, what kind of method is used for the following coding?
load X.mat; %%3D array
map = hsv(90);XR =Y;Ds = smooth3(XR);hiso = patch(isosurface(Ds,5),'FaceColor','blue','EdgeColor','none');hcap = patch(isocaps(XR,5),'FaceColor','interp','EdgeColor','none');colormap(map)daspect(gca,[1,1,.4])lightangle(305,30);fig = gcf;fig.Renderer = 'zbuffer';lighting phongisonormals(Ds,hiso)hcap.AmbientStrength = .6;hiso.SpecularColorReflectance = 0;hiso.SpecularExponent = 50;ax = gca;ax.View = [215,30];ax.Box = 'On';axis tighttitle('Original Data');
Please kindly help me to answer it. Thank you very much.
Best Answer