Voronoi_3D

利用MPT工具箱生成三维Voronoi图并输出拓扑信息到Excel文件。

% 此函数用来输出位于边界面上的点构成的面
function [C] = faceinbound（A）
C=cell（1100）; %存放所有满足要求的面的元胞数组
bound1=0;
bound2=1; %边界为0到1的正方体
[nm]=size（A）;
B0=cell（1m）;%存放位于0边界上的点
for i=1:m
    B0{i}=[];
end
B1=cell（1m）;%存放位于1边界上的点
for i=1:m
    B1{i}=[];
end
for i=1:n
    for j=1:m
        if abs（A（ij）-bound1）<0.00001
            B0{j}=[B0{j};A（i:）];
        end
        if abs（A（ij）-bound2）<0.00001
            B1{j}=[B1{j};A（i:）];
        end
    end
end
k=1;
for i=1:m
    if size（B0{i}1）>=3 %3点以上构成一个面
        C{k}=B0{i};
        k=k+1;
    end
    if size（B1{i}1）>=3
        C{k}=B1{i};
        k=k+1;
    end
end
[nm]=size（C）;
for i=1:m
    if isempty（C（i）） %获取C中非空元素数量
        break
    end
end
C=C（1:i）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        421  2015-01-25 16:48  intersectfordouble.m

     文件     333301  2015-01-25 16:28  test.xls

     文件       2095  2015-01-25 17:00  Voronoin_3D.m

     文件        821  2015-01-25 16:42  faceinbound.m

----------- ---------  ---------- -----  ----

               336638                    4