如何找到 3D 网格所有内部边缘的坐标？

Question

我有一张指定尺寸的3D图片，代码如下：

% input image
imageSizeY = 30;
imageSizeX = 20;
imageSizeZ = 5;

% obtain size of internal vertices
ny = imageSizeY - 1;
nx = imageSizeX - 1;
nz = imageSizeZ - 1;

% generate 3D mesh (image coordinate system)
[columnsInImage, rowsInImage, pagesInImage] = meshgrid(1:nx, 1:ny, 1:nz);

% obtain coordinates of vertices
vertexCoords = [rowsInImage(:), columnsInImage(:), pagesInImage(:)];

% obtain count of vertices (from top to bottom)
countVertex= ny*(vertexCoords(:,2) - 1) + vertexCoords(:,1) + ny*nx*(vertexCoords(:,3)-1);

在附图中，内部顶点标记为点，边缘标记为x

我已经能够找到图像所有内部顶点（从上到下）的坐标，如代码所示。然而，我需要找到所有边缘的坐标，但我不确定我是否确切知道如何实现这一目标。

Recent figure

Answer 1

因为它是一个规则的网格，所以我对坐标使用了逻辑运算来找到它们的最小值和最大值：

% deal vertexes into xyz coordinates
[y,x,z] = deal(vertexCoords(:,1),vertexCoords(:,2),vertexCoords(:,3));
nv = numel(x); % number of vertexes
allV = (1:nv).'; % indexes of all vertexes
% compute cube edges indexes => edge = [point1,point2]
edgeIdxs1 = [allV(y < ny),allV(y < ny) + 1];
edgeIdxs2 = [allV(x < nx),allV(x < nx) + ny];
edgeIdxs3 = [allV(z < nz),allV(z < nz) + ny*nx];
edgeIdxs = [edgeIdxs1;edgeIdxs2;edgeIdxs3];
% xyz coordinates of the edges
xEdges = [x(edgeIdxs(:,1)),x(edgeIdxs(:,2))];
yEdges = [y(edgeIdxs(:,1)),y(edgeIdxs(:,2))];
zEdges = [z(edgeIdxs(:,1)),z(edgeIdxs(:,2))];
% "outer" edges (as defined by you) 
outerEdges = all(xEdges == nx,2) | all(xEdges == 1,2) ...
    | all(yEdges == ny,2) | all(yEdges == 1,2) ...
    | all(zEdges == nz,2) | all(zEdges == 1,2);
% "inner" edges indexes
innerEdges = ~outerEdges;
innerEdgeIdxs = edgeIdxs(innerEdges,:);
% xyz coordinates of the inner edges (to plot x on their middle)
xInnerEdges = [x(innerEdgeIdxs(:,1)),x(innerEdgeIdxs(:,2))];
yInnerEdges = [y(innerEdgeIdxs(:,1)),y(innerEdgeIdxs(:,2))];
zInnerEdges = [z(innerEdgeIdxs(:,1)),z(innerEdgeIdxs(:,2))];
% "outer" vertexes indexes
outerVertexes = (x == nx | x == 1) | (y == ny | y == 1) | (z == nz | z == 1);
% "inner" vertexes indexes (to plot o on them)
innerVertexes = ~outerVertexes;
% plot all edges
plot3([x(edgeIdxs(:,1)),x(edgeIdxs(:,2))]',...
    [y(edgeIdxs(:,1)),y(edgeIdxs(:,2))]',...
    [z(edgeIdxs(:,1)),z(edgeIdxs(:,2))]','k');
hold on
% plot x's
plot3(mean(xInnerEdges,2),mean(yInnerEdges,2),mean(zInnerEdges,2),'xr');
% plot o's
plot3(x(innerVertexes),y(innerVertexes),z(innerVertexes),'oc');
axis equal
xlabel('x');ylabel('y');zlabel('z');

这是一个示例：

% input image
imageSizeY = 5;
imageSizeX = 5;
imageSizeZ = 5;

如何找到 3D 网格所有内部边缘的坐标？

How to find the coordinates of all internal edges of a 3D mesh?

matlab

image-processing

3d

coordinates

edges