如何用MATLAB读stl并显示点云文件

function [VertexData,FVCD,isBinary]=stl2matlab(stlfile)
% STL2MATLAB reads STL-file, ASCII or binary format, into Matlab
%
% Usage:
%
% [VertexData, FVCD]=stl2matlab(stlfile)
%
% Input:
%
% stlfile - STL-file
%
% Output:
%
% VertexData - Cellarray containing matrices with vertices
% FVCD - FaceVertexColorData
% isBinary - Flag with STL-file information
%
% To plot the STL-object use plotSTL(VertexData) or plotSTL(VertexData,FVCD)
%
% m-file can be downloaded at
% http://www.mathworks.com/matlabcentral/fileexchange/13253-iges-toolbox
%
% written by Per Bergstrm --
%

if nargin ==
    warning('No STL-file is specified');
end

VertexData=cell(,);

try

    % Try to read an STL ASCII file
    [VertexData{}, VertexData{}, VertexData{}, FVCD] = stlAread(stlfile);
    isBinary=false;

catch

    try

        % Try to read an STL binary file
        [VertexData{}, VertexData{}, VertexData{}, FVCD] = stlBread(stlfile);
        isBinary=true;

    catch

        error('File could not be read!')

    end

end

function [X, Y, Z, FVCD] = stlAread(stlfile)
% Reads an STL ASCII file

fid=fopen(stlfile,'r');

fileTitle=sscanf(fgetl(fid),'%*s %s');

vnum=;
fclr=;
testASCII=true;
lineCount=;

while feof(fid) ==
    stlLine=fgetl(fid);
    keyWord=sscanf(stlLine,'%s');
    if strncmpi(keyWord,'c',) == ;
        fclr=sscanf(stlLine,'%*s %f %f %f');
    elseif strncmpi(keyWord,'v',) == ;
        vnum=vnum+;
        vertex(:,vnum)=sscanf(stlLine,'%*s %f %f %f');
        clr(:,vnum)=fclr;
    elseif testASCII
        lineCount=lineCount+;
        if lineCount>
            if vnum>
                testASCII=false;
            else
                error('File is not an STL ASCII file!')
            end
        end
    end
end

X=[vertex(,::end);vertex(,::end);vertex(,::end)];
Y=[vertex(,::end);vertex(,::end);vertex(,::end)];
Z=[vertex(,::end);vertex(,::end);vertex(,::end)];

FVCD=reshape(clr,,vnum/);

fclose(fid);

function [X, Y, Z, FVCD] = stlBread(stlfile)
% Reads an STL binary file

fid=fopen(stlfile,'r');

fileTitle=fread(fid,,'uchar=>schar');
fnum=fread(fid,,'int32');

X=zeros(,fnum);
Y=zeros(,fnum);
Z=zeros(,fnum);

FVCD=uint8(zeros(,fnum));

for i=:fnum,
    normal=fread(fid,,'float32');
    vertex1=fread(fid,,'float32');
    vertex2=fread(fid,,'float32');
    vertex3=fread(fid,,'float32');
    clr=fread(fid,,'uint16');
    if bitget(clr,)==
        rd=bitshift(bitand(,clr),-);
        grn=bitshift(bitand(,clr),-);
        bl=bitand(,clr);
        FVCD(:,i)=[rd;grn;bl];
    end
    X(:,i)=[vertex1(); vertex2(); vertex3()];
    Y(:,i)=[vertex1(); vertex2(); vertex3()];
    Z(:,i)=[vertex1(); vertex2(); vertex3()];
end

fclose(fid);

把这个函数放到Matlab中去，然后

b = figure();
model = stl2matlab('rawmodel-983.stl');
patch(model{},model{},model{},'b');
axis equal;
view(,);
mesh(x, y,z);
title('topographie map')

具体什么原理不太懂，反正这样真的可以

就是图丑了一点2333333333333333333