function?[x?y?z?varargout]?=?stlread（filename）
%?This?function?reads?an?STL?file?in?binary?format?into?matrixes?X?Y?and
%?Z?and?C.??C?is?optional?and?contains?color?rgb?data?in?5?bits.??
%
%?USAGE:?[x?y?z?c]?=?stlread（filename）;
%
%?To?plot?use?patch（xyzc）?or?patch（xyz）
%
%?Written?by?Doron?Harlev

if?nargout>4
????error（‘Too?many?output?arguments‘）
end
use_color=（nargout==4）;

fid=fopen（filename?‘r‘）;?%Open?the?file?assumes?STL?Binary?format.
if?fid?==?-1?
????error（‘File?could?not?be?opened?check?name?or?path.‘）
end

ftitle=fread（fid80‘uchar=>schar‘）;?%?Read?file?title
num_facet=fread（fid1‘int32‘）;?%?Read?number?of?Facets

fprintf（‘\ntitle:?%s\n‘?char（ftitle‘））;
fprintf（‘Num?Facets:?%d\n‘?num_facet）;

%?Preallocate?memory?to?save?running?time
x=zeros（3num_facet）;?y=zeros（3num_facet）;?z=zeros（3num_facet）;
if?use_color
????c=uint8（zeros（3num_facet））;
end

h?=?waitbar（0‘Please?wait...‘）;
for?i=1:num_facet
????norm=fread（fid3‘float32‘）;?%?normal?coordinates?ignored?for?now
????ver1=fread（fid3‘float32‘）;?%?vertex?1
????ver2=fread（fid3‘float32‘）;?%?vertex?2
????ver3=fread（fid3‘float32‘）;?%?vertex?3
????col=fread（fid1‘uint16‘）;?%?color?bytes
????if?（bitget（col16）==1?&?use_color）
????????r=bitshift（bitand（2^16-1?col）-10）;
????????g=bitshift（bitand（2^11-1?col）-5）;
????????b=bitand（2^6-1?col）;
????????c（:i）=[r;?g;?b];
????end
????x（:i）=[ver1（1）;?ver2（1）;?ver3（1）];?%?convert?to?matlab?“patch“?compatible?format
????y（:i）=[ver1（2）;?ver2（2）;?ver3（2）];
????z（:i）=[ver1（3）;?ver2（3）;?ver3（3）];
????if?mod（ifloor（num_facet/10））==0
????????waitbar（i/num_faceth）;
????end
end
if?use_color
????varargout（1）={c};
end
fclose（fid）;
close（h）;

%?For?more?information?http://rpdrc.ic.polyu.edu.hk/old_files/stl_binary_format.htm

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1840??2018-07-11?20:18??stlread.m