1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
|
function savebinstl(node, elem, fname, solidname)
%
% savebinstl(node,elem,fname,solidname)
%
% save a tetrahedral mesh to a binary STL (Standard Tessellation Language) file
%
% author: Qianqian Fang, <q.fang at neu.edu>
% date: 2012/02/14
%
% input:
% node: input, surface node list, dimension Nx3
% elem: input, tetrahedral element list; if size(elem,2)==3, it is a surface
% fname: output file name
% solidname: an optional string for the name of the object
%
% -- this function is part of iso2mesh toolbox (http://iso2mesh.sf.net)
%
fid = fopen(fname, 'wb');
if (fid == -1)
error('You do not have permission to save mesh files.');
end
if (nargin < 3)
error('incomplete input');
end
if (isempty(node) || size(node, 2) < 3)
error('invalid node input');
end
hddata = struct('Ver', 1, 'Creator', 'iso2mesh', 'Date', datestr(now));
if (nargin < 4)
solidname = '';
else
hddata.name = solidname;
end
headerstr = savejson('', hddata);
headerstr = regexprep(headerstr, '[\t\n\r]', '');
headerstr = regexprep(headerstr, '": ', '":');
if (length(headerstr) > 80)
headerstr = headerstr(1:80); % this makes the header an invalid json string
else
headerstr(81) = 0;
end
fwrite(fid, headerstr(1:80), 'char');
if (~isempty(elem))
if (size(elem, 2) >= 5)
elem(:, 5:end) = [];
end
if (size(elem, 2) == 4)
elem = meshreorient(node, elem);
elem = volface(elem);
end
ev = surfplane(node, elem);
ev = ev(:, 1:3) ./ repmat(sqrt(sum(ev(:, 1:3) .* ev(:, 1:3), 2)), 1, 3);
len = size(elem, 1);
fwrite(fid, len, 'uint32');
for i = 1:len
fwrite(fid, [ev(i, :)', node(elem(i, :), 1:3)'], 'float32');
fwrite(fid, 0, 'uint16');
end
end
fclose(fid);
|
