Neighbour voxels of a list of voxels

Question

Razakamandimby Tiana 2021-5-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/832898-neighbour-voxels-of-a-list-of-voxels

评论： Razakamandimby Tiana 2021-5-20

I am trying to find all the neighbours voxels of a list of voxels. I already used a function findNeigbours that I downloaded (https://www.mathworks.com/matlabcentral/fileexchange/68549-findneighbours) but it only works for a single voxel and I need to apply the function to a vector.

In my code, I used for loop but it's takiing too long time since the vector list is large. Any solution?

These are the lines I used in my code:

for i = 1:size(vector,1)
    surrounding{i} = findNeighbours(vector(i), [[pxl_x, pxl_y, pxl_z], 26)
end

2 个评论
显示无隐藏无

KSSV 2021-5-18

You can use inbuilt functions to achieve this fast. Show us how your voxels are.

Razakamandimby Tiana 2021-5-19

vector.xlsx

The vector is in the excel file and the size of array in 3dimension is 60 x 60 x 67

请先登录，再进行评论。

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

Answer 1

Matt J 2021-5-18

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/832898-neighbour-voxels-of-a-list-of-voxels#answer_702448

编辑：Matt J 2021-5-18

在 MATLAB Online 中打开

If you have lists of xCoordinates, yCoordinates, and zCoordinates of the voxels, you may simply do,

[dx,dy,dz]=ndgrid([-1,0,1]);
dx=dx(:).'; dy=dy(:).'; dz=dz(:).';
X=xCoordinates(:)+dx; Y=yCoordinates(:)+dy; Z=zCoordinates(:)+dz; %neighbor coordinates

and if the voxel grid has dimensions MxNxP, you might then discard out-of-bounds voxel coordinates by doing,

X(X<1 | X>M)=nan;
Y(Y<1 | Y>N)=nan;
Z(Z<1 | Z>P)=nan;

2 个评论
显示无隐藏无

Razakamandimby Tiana 2021-5-19

Thanks but I have the voxels indexes not their coordinates. But I will try to get the coordinates from the inndexes and see if I get results faster

Matt J 2021-5-19

Just use ind2sub().

请先登录，再进行评论。

Answer 2

Razakamandimby Tiana 2021-5-20

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/832898-neighbour-voxels-of-a-list-of-voxels#answer_704038

在 MATLAB Online 中打开

I could solve this issue. The following code gives all the neighbour voxels (connectivity 26 in 3D) of each voxels in the vector list.

vector % The attached excel file
matrix = [121 121 117];
conn = 26;
neighbor = arrayfun(@(vector) findNeighbours(vector,matrix,conn), vector,'UniformOutput', false)