MATLAB读取STL文件代码并生成三维模型

function VerMat=ReadSTLFile()
 %读取STL文件代码
filename='model1or.stl'; %指定要打开的文件名
fid=fopen(filename,'r');  %打开文件得到操作柄
VerMat(1,:)=[0 0 0];%初始化保存顶点的数组
k=0;
while(feof(fid)==0)
    line=fgetl(fid);%读取文件的一行
matchev=findstr(line,'vertex');%检查该行中的数据是否包含顶点数据
    if(matchev)
        k=k+1;
        VerMat(k,:)=OpenVertexCode(line);%提取顶点坐标
    end
end
fclose(fid);

function vec=OpenVertexCode(line)
 %输入变量line是一个字符串该字符串中包含一个坐标
 %本程序将这个用字符存储的坐标转换为浮点型数据
 %输出变量是一个数组表示顶点坐标
n=length(line);
k=0;
BlSp(1)=0;%初始化保存空格位置的数组
for i=1:n
    if(line(i)==32)%发现空格
        k=k+1;
        BlSp(k)=i;
    end
end
Pos1=BlSp(1)+1;
Pos2=BlSp(2)-1;
N=Pos2-Pos1+1;
NumCode=line(Pos1:Pos2);%提取表示x坐标的字符
vec(1)=OpenNumCode(NumCode,N);%将字符转换为为浮点型数据
Pos1=BlSp(2)+1;
Pos2=BlSp(3)-1;
N=Pos2-Pos1+1;
NumCode=line(Pos1:Pos2); %提取表示y坐标的字符
vec(2)=OpenNumCode(NumCode,N); %将字符转换为为浮点型数据
Pos1=BlSp(3)+1;
Pos2=n;
N=Pos2-Pos1+1;
NumCode=line(Pos1:Pos2); %提取表示z坐标的字符
vec(3)=OpenNumCode(NumCode,N); %将字符转换为为浮点型数据

function ab=OpenNumCode(NumCode,N)
 %将字符转换为为浮点型数据
is=1;
if(NumCode(1)==45)%发现-号
    is=2;
end
for i=1:N
    if(NumCode(i)==46)%发现小数点
        ip=i;
        break;
    end
end
isf=ip-1;
ise=ip+1;
af=0;
posnum=1/10;
for i=isf:-1:is
    posnum=posnum*10;
    af=af+(NumCode(i)-48)*posnum;
end
bf=0;
posnum=1.0;
for i=ise:N
    posnum=posnum/10;
    bf=bf+(NumCode(i)-48)*posnum;
end
ab=af+bf;
if(is==2)
    ab=-ab;
end

function [pt,pgon]=ConstrucTopology(VerMat)
 %构建拓扑结构
N=length(VerMat(:,1));
Num=N/3;
for i=1:Num
    pgon{i}=[3*(i-1)+1 3*(i-1)+2 3*(i-1)+3];
end
pt=VerMat;

function showpgon(pt,pgon)
 %显示三角网格
view(3);
pat=zeros(1,length(pgon));
hold on
for i=1:length(pgon)
   ptpgoni=pt(pgon{i},:);
   pat(i)=patch(ptpgoni(:,1) ,  ptpgoni(:,2) , ptpgoni(:,3) ,[1.000    0.7000    0.3000], 'facelighting','flat','edgecolor',...
       [0 1.00 1.00], 'edgelighting','flat');
end
hold off
axis equal

function ReadSTLTest()
%读取和显示STL文件
VerMat=ReadSTLFile();
[pt,pgon]=ConstrucTopology(VerMat);
showpgon(pt,pgon);