paraview vtk文件乱码_Matlab处理数据导出Paraview可读的vtk文件（一）

谁践踏了优雅 2022-12-31 15:22 67阅读 0赞

1 function vtkwrite( filename,dataType,varargin )2 %VTKWRITE Writes 3D Matlab array into VTK file format.3 % vtkwrite(filename,'structured\_grid',x,y,z,'vectors',title,u,v,w) writes4 % a structured 3D vector data into VTK file, with name specified by the string

5 %filename. (u,v,w) are the vector components at the points (x,y,z). x,y,z6 % should be 3-D matrices like those generated by meshgrid, where

7 % point(ijk) isspecified by x(i,j,k), y(i,j,k) and z(i,j,k).8 %The matrices x,y,z,u,v,w must all be the same size and contain9 % corrresponding position and vector component. The stringtitle specifies10 %the name of the vector field to be saved.11 %

12 % vtkwrite(filename,'structured\_grid',x,y,z,'scalars',title,r) writes a 3D13 % scalar data into VTK file whose name is specified by the string

14 % filename. r isthe scalar value at the points (x,y,z). The matrices15 %x,y,z,r must all be the same size and contain the corresponding position16 %and scalar values.17 %

18 % vtkwrite(filename,'structured\_grid',x,y,z,'vectors',title,u,v,w,'scalars',19 %title2,r) writes a 3D structured grid that contains both vector and scalar values.20 %x,y,z,u,v,w,r must all be the same size and contain the corresponding21 %positon, vector and scalar values.22 %

23 % vtkwrite(filename, 'structured\_points', title, m) saves matrix m (could24 % be 1D, 2D or 3D array) into vtk asstructured points.25 %

26 % vtkwrite(filename, 'structured\_points', title, m, 'spacing', sx, sy, sz)27 % allows user to specify spacing. (default: 1, 1, 1). This isthe aspect28 %ratio of a single voxel.29 %

30 % vtkwrite(filename, 'structured\_points', title, m, 'origin', ox, oy, oz)31 % allows user to speicify origin of dataset. (default: 0, 0, 0).32 %

33 % vtkwrite(filename,'unstructured\_grid',x,y,z,'vectors',title,u,v,w,'scalars',34 %title2,r) writes a 3D unstructured grid that contains both vector and scalar values.35 %x,y,z,u,v,w,r must all be the same size and contain the corresponding36 %positon, vector and scalar values.37 %

38 % vtkwrite(filename, 'polydata', 'lines', x, y, z) exports a 3D line where

39 %x,y,z are coordinates of the points that make the line. x, y, z are40 % vectors containing the coordinates of points of the line, wherepoint(n)41 % isspecified by x(n), y(n) and z(n).42 %

43 % vtkwrite(filename,'polydata','lines',x,y,z,'Precision',n) allows you to44 % specify precision of the exported number up to n digits after decimal

45 % point. Default precision is 3digits.46 %

47 % vtkwrite(filename,'polydata','triangle',x,y,z,tri) exports a list of48 % triangles where x,y,z are the coordinates of the points and tri isan49 % m\*3matrix whose rows denote the points of the individual triangles.50 %

51 % vtkwrite(filename,'polydata','tetrahedron',x,y,z,tetra) exports a list52 % of tetrahedrons wherex,y,z are the coordinates of the points53 % and tetra is an m\*4matrix whose rows denote the points of individual54 %tetrahedrons.55 %

56 % vtkwrite('execute','polydata','lines',x,y,z) will save data with default

57 % filename ''matlab\_export.vtk'and automatically loads data into

58 %ParaView.59 %

60 % Version 2.3

61 % Copyright, Chaoyuan Yeh, 2016

62 % Codes are modified from William Thielicke and David Gingras's submission.

64 if strcmpi(filename,'execute'), filename = 'matlab\_export.vtk'; end65 fid = fopen(filename, 'w');66 %VTK files contain five major parts67 % 1. VTK DataFile Version68 fprintf(fid, '\# vtk DataFile Version 2.0\\n');69 % 2. Title70 fprintf(fid, 'VTK from Matlab\\n');71

73 binaryflag = any(strcmpi(varargin, 'BINARY'));74 if any(strcmpi(varargin, 'PRECISION'))75 precision = num2str(varargin\{find(strcmpi(vin, 'PRECISION'))+1\});76 else

77 precision = '2';78 end79

80 switchupper(dataType)81 case 'STRUCTURED\_POINTS'

82 title = varargin\{1\};83 m = varargin\{2\};84 if any(strcmpi(varargin, 'spacing'))85 sx = varargin\{find(strcmpi(varargin, 'spacing'))+1\};86 sy = varargin\{find(strcmpi(varargin, 'spacing'))+2\};87 sz = varargin\{find(strcmpi(varargin, 'spacing'))+3\};88 else

89 sx = 1;90 sy = 1;91 sz = 1;92 end93 if any(strcmpi(varargin, 'origin'))94 ox = varargin\{find(strcmpi(varargin, 'origin'))+1\};95 oy = varargin\{find(strcmpi(varargin, 'origin'))+2\};96 oz = varargin\{find(strcmpi(varargin, 'origin'))+3\};97 else

98 ox = 0;99 oy = 0;100 oz = 0;101 end102 \[nx, ny, nz\] =size(m);103 setdataformat(fid, binaryflag);104

105 fprintf(fid, 'DATASET STRUCTURED\_POINTS\\n');106 fprintf(fid, 'DIMENSIONS %d %d %d\\n', nx, ny, nz);107 fprintf(fid, \['SPACING', num2str(sx), ' ', num2str(sy), ' ',...108 num2str(sz), '\\n'\]);109 fprintf(fid, \['ORIGIN', num2str(ox), ' ', num2str(oy), ' ',...110 num2str(oz), '\\n'\]);111 fprintf(fid, 'POINT\_DATA %d\\n', nx\*ny\*nz);112 fprintf(fid, \['SCALARS', title, 'float 1\\n'\]);113 fprintf(fid,'LOOKUP\_TABLE default\\n');114 if ~binaryflag115 spec = \['%0.', precision, 'f'\];116 fprintf(fid, spec, m(:)');

117 else

118 fwrite(fid, m(:)','float','b');

119 end120

121 case \{'STRUCTURED\_GRID','UNSTRUCTURED\_GRID'\}122 % 3. The format data proper is saved in(ASCII or Binary). Use123 % fprintf to write data in the case of ASCII and fwrite forbinary.124 if numel(varargin)<6, error('Not enough input arguments'); end125 setdataformat(fid, binaryflag);126 % fprintf(fid, 'BINARY\\n');127 x = varargin\{1\};128 y = varargin\{2\};129 z = varargin\{3\};130 if sum(size(x)==size(y) & size(y)==size(z))~=length(size(x))131 error('Input dimesions do not match')132 end133 n\_elements =numel(x);134 % 4. Type of Dataset ( can be STRUCTURED\_POINTS, STRUCTURED\_GRID,135 %UNSTRUCTURED\_GRID, POLYDATA, RECTILINEAR\_GRID or FIELD )136 %This part, dataset structure, begins with a line containing the137 % keyword 'DATASET'followed by a keyword describing the type of dataset.138 %Then the geomettry part describes geometry and topology of the dataset.139 if strcmpi(dataType,'STRUCTURED\_GRID')140 fprintf(fid, 'DATASET STRUCTURED\_GRID\\n');141 fprintf(fid, 'DIMENSIONS %d %d %d\\n', size(x,1), size(x,2), size(x,3));142 else

143 fprintf(fid, 'DATASET UNSTRUCTURED\_GRID\\n');144 end145 fprintf(fid, \['POINTS' num2str(n\_elements) 'float\\n'\]);146 output = \[x(:)'; y(:)'; z(:)'\];

147

148 if ~binaryflag149 spec = \['%0.', precision, 'f'\];150 fprintf(fid, spec, output);151 else

152 fwrite(fid, output, 'float', 'b');153 end154 % 5.This final part describe the dataset attributes and begins with the155 % keywords 'POINT\_DATA' or 'CELL\_DATA', followed by an integer number156 % specifying the number of points of cells. Other keyword/data combination157 %then define the actual dataset attribute values.158 fprintf(fid, \['\\nPOINT\_DATA'num2str(n\_elements)\]);159 %Parse remaining argument.160 vidx = find(strcmpi(varargin,'VECTORS'));161 sidx = find(strcmpi(varargin,'SCALARS'));162 if vidx~=0

163 for ii = 1:length(vidx)164 title = varargin\{vidx(ii)+1\};165 %Data enteries begin with a keyword specifying data type166 %and numeric format.167 fprintf(fid, \['\\nVECTORS', title,'float\\n'\]);168 output = \[varargin\{ vidx(ii) + 2 \}(:)';...

169 varargin\{ vidx(ii) + 3 \}(:)';...

170 varargin\{ vidx(ii) + 4 \}(:)'\];

171

172 if ~binaryflag173 spec = \['%0.', precision, 'f'\];174 fprintf(fid, spec, output);175 else

176 fwrite(fid, output, 'float', 'b');177 end178 % fwrite(fid, \[reshape(varargin\{vidx(ii)+2\},1,n\_elements);...179 % reshape(varargin\{vidx(ii)+3\},1,n\_elements);...180 % reshape(varargin\{vidx(ii)+4\},1,n\_elements)\],'float','b');181 end182 end183 if sidx~=0

184 for ii = 1:length(sidx)185 title = varargin\{sidx(ii)+1\};186 fprintf(fid, \['\\nSCALARS', title,'float\\n'\]);187 fprintf(fid, 'LOOKUP\_TABLE default\\n');188 if ~binaryflag189 spec = \['%0.', precision, 'f'\];190 fprintf(fid, spec, varargin\{ sidx(ii) + 2\});191 else

192 fwrite(fid, varargin\{ sidx(ii) + 2\}, 'float', 'b');193 end194 % fwrite(fid, reshape(varargin\{sidx(ii)+2\},1,n\_elements),'float','b');195 end196 end197

198 case 'POLYDATA'

199

200 fprintf(fid, 'ASCII\\n');201 if numel(varargin)<4, error('Not enough input arguments'); end202 x = varargin\{2\}(:);203 y = varargin\{3\}(:);204 z = varargin\{4\}(:);205 if numel(varargin)<4, error('Not enough input arguments'); end206 if sum(size(x)==size(y) & size(y)==size(z))~=length(size(x))207 error('Input dimesions do not match')208 end209 n\_elements =numel(x);210 fprintf(fid, 'DATASET POLYDATA\\n');211 if mod(n\_elements,3)==1

212 x(n\_elements+1:n\_elements+2,1)=\[0;0\];213 y(n\_elements+1:n\_elements+2,1)=\[0;0\];214 z(n\_elements+1:n\_elements+2,1)=\[0;0\];215 elseif mod(n\_elements,3)==2

216 x(n\_elements+1,1)=0;217 y(n\_elements+1,1)=0;218 z(n\_elements+1,1)=0;219 end220 nbpoint =numel(x);221 fprintf(fid, \['POINTS' num2str(nbpoint) 'float\\n'\]);222

223 spec = \[repmat(\['%0.', precision, 'f'\], 1, 9), '\\n'\];224

225 output = \[x(1:3:end-2), y(1:3:end-2), z(1:3:end-2),...226 x(2:3:end-1), y(2:3:end-1), z(2:3:end-1),...227 x(3:3:end), y(3:3:end), z(3:3:end)\]';

228

229 fprintf(fid, spec, output);230

231 switch upper(varargin\{1\})232 case 'LINES'

233 if mod(n\_elements,2)==0

234 nbLine = 2\*n\_elements-2;235 else

236 nbLine = 2\*(n\_elements-1);237 end238 conn1 = zeros(nbLine,1);239 conn2 = zeros(nbLine,1);240 conn2(1:nbLine/2) = 1:nbLine/2;241 conn1(1:nbLine/2) = conn2(1:nbLine/2)-1;242 conn1(nbLine/2+1:end) = 1:nbLine/2;243 conn2(nbLine/2+1:end) = conn1(nbLine/2+1:end)-1;244 fprintf(fid,'\\nLINES %d %d\\n',nbLine,3\*nbLine);245 fprintf(fid,'2 %d %d\\n',\[conn1';conn2'\]);246 case 'TRIANGLE'

247 ntri = length(varargin\{5\});248 fprintf(fid,'\\nPOLYGONS %d %d\\n',ntri,4\*ntri);249 fprintf(fid,'3 %d %d %d\\n',(varargin\{5\}-1)');

250 case 'TETRAHEDRON'

251 ntetra = length(varargin\{5\});252 fprintf(fid,'\\nPOLYGONS %d %d\\n',ntetra,5\*ntetra);253 fprintf(fid,'4 %d %d %d %d\\n',(varargin\{5\}-1)');

254 end255 end256 fclose(fid);257 if strcmpi(filename,'matlab\_export.vtk')258 switchcomputer259 case \{'PCWIN','PCWIN64'\}260 !paraview.exe --data='matlab\_export.vtk' &

261 % Exclamation point character isa shell escape, the rest of the262 %input line will be sent to operating system. It can not take263 % variables, though. The & at the end of line will returncontrol to264 % Matlab even when the outside process isstill running.265 case \{'GLNXA64','MACI64'\}266 !paraview --data='matlab\_export.vtk' &

267 end268 end269 end270

271 function setdataformat(fid, binaryflag)272

273 if ~binaryflag274 fprintf(fid, 'ASCII\\n');275 else

276 fprintf(fid, 'BINARY\\n');277 end278 end