ok, I think I have a suspect: the working mesh has element type 12 (VTK_HEXAHEDRON), while the two that do not work have 11 (VTK_VOXEL). I have updated SfePy to support pixels and voxels - check out the current git master. The sinus.vtk seems to work, about the other file I am not sure, as it takes a long time, but it can be read...
On 04/27/2012 10:27 AM, Xiaogai Li wrote:
Here is one more model with the same problem. This geometry is much simpler and only a few elements. Thanks. BR, XG
---------- Forwarded message ---------- From: Xiaogai Lixiaog...@gmail.com Date: Fri, Apr 27, 2012 at 9:54 AM Subject: To: cimr...@ntc.zcu.cz
Thanks and please find the attached file.
smooth_voxel_mesh.py is the python script I used to run the smoothing.
atlas_csf-label_legacy_notwork.vtk is the file that with the problem "IndexError: index out of bounds".
is the file that is part of the whole image, and this image works with the python script.
So I wonder what is the problem in file atlas_csf-label_legacy_notwork.vtk.
Thanks a lot for your time. BR Xiaogai Li
On Fri, Apr 27, 2012 at 9:08 AM, Robert Cimrmancimr...@ntc.zcu.cz wrote:
could you send me (off-list) the file? It's a strange error.
I met a strange error when using smooth_mesh. It seems that I could not create the object Domain, however, I have verified that the voxel-mesh was quite nice.
coordssmooth = smooth_mesh(m, 4, 0.6307, -0.6347, None, True, False) File "/usr/local/lib/python2.7/dist-packages/sfepy/mesh/mesh_tools.py", line 505, in smooth_mesh domain = Domain('mesh', mesh) File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/domain.py", line 223, in __init__ self.setup_facets(verbose=verbose) File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/domain.py", line 375, in setup_facets obj.setup_unique() File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/facets.py", line 281, in setup_unique self.uid, self.uid[1:] = 0, ii[:-1] IndexError: index out of bounds