-
Notifications
You must be signed in to change notification settings - Fork 3
/
Copy pathvtk_rw.py
144 lines (126 loc) · 5.8 KB
/
vtk_rw.py
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
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
# -*- coding: utf-8 -*-
'''
Functions to read and write vtk files
------------------------------------
* read takes vtk file and returns vertex and face array
* write takes vertex and faces array and optional comment and returns vtk file
* reading and writing of texture not supported yet
TO DO: add reading comments
'''
def read_vtk(file):
'''
Reads ASCII coded vtk files using pandas,
returning vertices, faces and data as three numpy arrays.
Usage:
---------------------
vertices, faces, data = read_vtk('/path/to/vtk/file.vtk')
'''
import pandas as pd
import numpy as np
# read full file while dropping empty lines
vtk_df=pd.read_csv(file, header=None, engine='python')
vtk_df=vtk_df.dropna()
#print vtk_df
# extract number of vertices and faces
number_vertices=int(vtk_df[vtk_df[0].str.contains('POINTS')][0].iloc[0].split()[1])
print("number of vertices:", number_vertices)
number_faces=int(vtk_df[vtk_df[0].str.contains('POLYGONS')][0].iloc[0].split()[1])
print("number of faces:", number_faces)
# read vertices into df and array
start_vertices= (vtk_df[vtk_df[0].str.contains('POINTS')].index.tolist()[0])+1
#print start_vertices
vertex_df=pd.read_csv(file, skiprows=range(start_vertices), nrows=number_vertices, sep='\s*', header=None, engine='python')
#print vertex_df
if np.array(vertex_df).shape[1]==3:
vertex_array=np.array(vertex_df)
# sometimes the vtk format is weird with 9 indices per line, then it has to be reshaped
elif np.array(vertex_df).shape[1]==9:
vertex_df=pd.read_csv(file, skiprows=range(start_vertices), nrows=int(number_vertices), sep='\s*', header=None, engine='python')
vertex_array=np.array(vertex_df.iloc[0:1,0:3])
vertex_array=np.append(vertex_array, vertex_df.iloc[0:1,3:6], axis=0)
vertex_array=np.append(vertex_array, vertex_df.iloc[0:1,6:9], axis=0)
for row in range(1,(int(number_vertices/3)+1)):
for col in [0,3,6]:
vertex_array=np.append(vertex_array, np.array(vertex_df.iloc[row:(row+1),col:(col+3)]),axis=0)
# strip rows containing nans
#vertex_array=vertex_array[ ~np.isnan(vertex_array) ].reshape(number_vertices,3)
#print vertex_array
else:
print "vertex indices out of shape"
# read faces into df and array
start_faces= (vtk_df[vtk_df[0].str.contains('POLYGONS')].index.tolist()[0])+1
#print start_faces
face_df=pd.read_csv(file, skiprows=range(start_faces), nrows=number_faces, sep='\s*', header=None, engine='python')
#print face_df
face_array=np.array(face_df.iloc[:,1:4])
#print face_df.iloc[1,0:4]
# read data into df and array if exists
if vtk_df[vtk_df[0].str.contains('POINT_DATA')].index.tolist()!=[]:
start_data=(vtk_df[vtk_df[0].str.contains('POINT_DATA')].index.tolist()[0])+3
number_data = number_vertices
data_df=pd.read_csv(file, skiprows=range(start_data), nrows=number_data, sep='\s*', header=None, engine='python')
data_array=np.array(data_df)
else:
data_array = np.empty(0)
return vertex_array, face_array, data_array
def write_vtk(filename, vertices, faces, data=None, comment=None):
'''
Creates ASCII coded vtk file from numpy arrays using pandas.
Inputs:
-------
(mandatory)
* filename: str, path to location where vtk file should be stored
* vertices: numpy array with vertex coordinates, shape (n_vertices, 3)
* faces: numpy array with face specifications, shape (n_faces, 3)
(optional)
* data: numpy array with data points, shape (n_vertices, n_datapoints)
NOTE: n_datapoints can be =1 but cannot be skipped (n_vertices,)
* comment: str, is written into the comment section of the vtk file
Usage:
---------------------
write_vtk('/path/to/vtk/file.vtk', v_array, f_array)
'''
import pandas as pd
import numpy as np
# infer number of vertices and faces
number_vertices=vertices.shape[0]
number_faces=faces.shape[0]
if data is not None:
number_data=data.shape[0]
# make header and subheader dataframe
header=['# vtk DataFile Version 3.0',
'%s'%comment,
'ASCII',
'DATASET POLYDATA',
'POINTS %i float'%number_vertices
]
header_df=pd.DataFrame(header)
sub_header=['POLYGONS %i %i'%(number_faces, 4*number_faces)]
sub_header_df=pd.DataFrame(sub_header)
# make dataframe from vertices
vertex_df=pd.DataFrame(vertices)
# make dataframe from faces, appending first row of 3's (indicating the polygons are triangles)
triangles=np.reshape(3*(np.ones(number_faces)), (number_faces,1))
triangles=triangles.astype(int)
faces=faces.astype(int)
faces_df=pd.DataFrame(np.concatenate((triangles,faces),axis=1))
# write dfs to csv
header_df.to_csv(filename, header=None, index=False)
with open(filename, 'a') as f:
vertex_df.to_csv(f, header=False, index=False, float_format='%.3f', sep=' ')
with open(filename, 'a') as f:
sub_header_df.to_csv(f, header=False, index=False)
with open(filename, 'a') as f:
faces_df.to_csv(f, header=False, index=False, float_format='%.0f', sep=' ')
# if there is data append second subheader and data
if data != None:
datapoints=data.shape[1]
sub_header2=['POINT_DATA %i'%(number_data),
'SCALARS EmbedVertex float %i'%(datapoints),
'LOOKUP_TABLE default']
sub_header_df2=pd.DataFrame(sub_header2)
data_df=pd.DataFrame(data)
with open(filename, 'a') as f:
sub_header_df2.to_csv(f, header=False, index=False)
with open(filename, 'a') as f:
data_df.to_csv(f, header=False, index=False, float_format='%.16f', sep=' ')