Rotated Slices of 3D volumes

Question

Ahmed Zankoor 2018-8-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/416163-rotated-slices-of-3d-volumes

评论： Ahmed Zankoor 2018-10-22

Hello, I have a 3D volume of voxels, this volume only contains integers 1,2 and 3. If I took slices parallel to principal planes ( x-y, x-z, y-z) the slice.CData only contains integers 1,2 and 3, which is correct. But if I used a rotated slice, slice.CData contains numbers inbetween and many NaNs. which may be because the slice does not cut exactly through the voxel. So my question is, is there a method to make slices that gives me the integer of the closest voxel which the slice cut through. To Clarify this , we can run the example given in Matlab documentation:

 [x,y,z] = meshgrid(-2:.2:2,-2:.25:2,-2:.16:2);
v = x.*exp(-x.^2-y.^2-z.^2);
figure
colormap hsv
for k = 0
   hsp = surf(linspace(-2,2,20),linspace(-2,2,20),...
      zeros(20) + k);
   rotate(hsp,[1,-1,1],30)
   xd = hsp.XData;
   yd = hsp.YData;
   zd = hsp.ZData;
   delete(hsp)
   slice(x,y,z,v,[-2,2],2,-2)
   hold on
   S = slice(x,y,z,v,xd,yd,zd)
   hold off
   view(-5,10)
   axis([-2.5 2.5 -2 2 -2 4])
   drawnow
end
A = S.CData;
unique(A)

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Rik 2018-8-25

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/416163-rotated-slices-of-3d-volumes#answer_334078

在 MATLAB Online 中打开

You can use round on the output data, although I don't see where your data would have the values you mention. The code below should help.

A = S.CData;
A(A<1)=1;
A(A>3)=3;
A=round(A);
S.CData=A;

3 个评论
显示 1更早的评论隐藏 1更早的评论

Rik 2018-8-26

Glad to be of help. I don't really understand if you still have a question or not, so if you have a follow-up question, please clarify what it is in a comment.

Ahmed Zankoor 2018-10-22

Hello rik, it has been a while since this discussion, but I just wanted to say that using 'nearest' option in slicing can be useful to avoid the previous interpolation. S = slice(V,Xp,Yp,Zp,'nearest');

请先登录，再进行评论。