如何搜索3D体积/表面边界的3D网格坐标索引内部

如何解决如何搜索3D体积/表面边界的3D网格坐标索引内部

该代码的目的是使用给定的数据数组，创建一个等值线0.25。数据是分散的3D数组，其xyz坐标中每个点的值。因此，我们创建了一个3D网格网格阵列，选择了数据表面边界内的网格点，将[data（4，:)]值插值到3D网格阵列中，并获得了这些值的等值面。

我被困在寻找或搜索3-D网格坐标数组索引中所有位于数据点定义的体积内的点的步骤。

是否有任何函数可以通过data（1，:)，data（2，:)，data（3，:)查找表面边界内的xx，yy，zz的索引。那就是使用3D数组找到3-D网格坐标的索引。

主要代码步骤是

◾创建数据点的边界[data（1，：）= x，data（2，：）= y，data（3，：）= z，data（4，：）=值]

◾网格已创建xx yy zz

◾缺少步骤-选择位于数据点边界内部的3-D网格坐标的索引...

◾DATA-创建用于保存插值的3D数组

◾使用Meshgrid DATA创建等值面

% return index for boundary of all points sequence which is inside volume
j = boundary(data(1,:)',data(2,data(3,0.97);
hold on; trisurf(j,data(1,0.97,'FaceColor','red','FaceAlpha',0.1)

% PREPARATION
% ===========
% GENERATE A GRID for x y z co-ordinate
[xx yy zz] = meshgrid(linspace(xx_min,xx_max,20),linspace(yy_min,yy_max,linspace(0,zz_max,20));  % replace,0 1,10 with range of your values
% GENERATE RANDOM DATA for interplot 3D array value v - grid
DATA = zeros(20,20,20);
%Need to interplot the grid points [xx,yy,zz] all 3D grid points- NOT interior points**
DATA= griddata(data(1,data(4,xx,zz); 
% EXTRACT THE  SURFACE
SURF = isosurface(xx,zz,DATA,0.25);  % all x,y,z,v are 3D array 20X20X20
% PLOT THE MASK SURFACE AND DATA
hold on; axis square; axis([0 1 0 1 0 1]); view(3); camlight
p = patch(SURF); isonormals(xx,MASK,p) ;p.FaceColor = 'red'; p.EdgeColor = 'none';

我已将tsearch使用三边界点用于1D数组坐标，但这不适用于3D网格矩阵。例如：

% to located interior points in 3D grid 
xyz3D=[xx(:) yy(:) zz(:)]; % Make 3D grid to a 3D array for using tsearchn
plot3(xyz3D(:,1),xyz3D(:,2),3),'*');view(20,30); % plot 3D grid
tri = delaunayn([data(1,j)' data(2,j)' data(3,j)']); % Generate delaunay triangulization
[tn p] = tsearchn([data(1,j)'],tri,xyz3D) ; % how to use 3D arry for double
tnn=(tn>1);
plot3(xyz3D(tnn,xyz3D(tnn,30); % plot inside grid points

temp= griddata(data(1,data(15,3)); %query interior??
DATA2=reshape(temp,[20,20]);

-此处的tnn索引[delaunayn]并不完全包含所有内部点-有些在边界trisuf之外。

-如何获取边界面内具有某个值（例如0.25）的区域的体积。

enter image description here 图片显示

• 具有网格点的表面边界
• 德拉诺因内部具有[tnn]内点的表面