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...
r.
On 04/27/2012 10:27 AM, Xiaogai Li wrote:
> Hi again,
>
> 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 Li<xiaog...(a)gmail.com>
> Date: Fri, Apr 27, 2012 at 9:54 AM
> Subject:
> To: cimr...(a)ntc.zcu.cz
>
>
> Hi Robert,
>
> 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".
>
> atlas_csf-label_legacy_part_works.vtk
>
> 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 Cimrman<cimr...(a)ntc.zcu.cz> wrote:
>
>> Hi,
>>
>> could you send me (off-list) the file? It's a strange error.
>>
>> r.
>>
>>
>> ----- Reply message -----
>> From: "Xiaogai Li"<xiaog...(a)gmail.com>
>> To:<sfepy...(a)googlegroups.com>
>> Subject: About the mesh_tools.smooth
>> Date: Fri, Apr 27, 2012 00:37
>>
>>
>> Hello everyone,
>>
>> 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[0], self.uid[1:] = 0, ii[:-1]
>> IndexError: index out of bounds