Scripts: triangles.py

#
# Access vertices and triangles of volume surface mesh.
#
from VolumeViewer import volume_list
for v in volume_list():
    for p in v.surface_piece_list:
        # vertices is N by 3 numpy array of floats (x,y,z) vertex positions
        # triangles is T by 3 numpy array of indices into vertices giving 3 triangle vertices.
        vertices, triangles = p.geometry

        print('%s map, %d vertices, %d triangles' % (v.name, len(vertices), len(triangles)))
        print ('vertex coordinates:')
        for x,y,z in vertices[:5,:]:
            print ('  %8.5g %8.5g %8.5g' % (x,y,z))
        print('...')
        print('triangle vertex indices:')
        for i0,i1,i2 in triangles[:5,:]:
            print ('  %6d %6d %6d' % (i0,i1,i2))
        print('...')