[Insight-users] MeshSpatialObject Issues

[Julien Jomier]

Hi Chris,

I guess you are trying to use SOViewer to display a 2D mesh (from the 
file you sent). This is currently not possible, therefore only 3D meshes 
are supported (see example at the end of the email). You can try to 
convert your 2D mesh into a 3D one by setting the Z=0 of your points.

Regarding the SpatialObjectToImageFilter, this filter is looking if the 
point is on the mesh (not inside the mesh). That's why (25,25) is 
returning false and currently this function is not implemented for 
LineCells... You can submit a bug report (www.itk.org/Bug) and we'll try 
to implement it.

Let us know if you still have any issues/concerns.
Thanks,

Julien

----
ObjectType = Scene
NDims = 3
NObjects = 1
ObjectType = Mesh
NDims = 3
ID = 0
TransformMatrix = 1 0 0 0 1 0 0 0 1
Offset = 0 0 0
CenterOfRotation = 0 0 0
ElementSpacing = 1 1 1
PointType = MET_FLOAT
PointDataType = MET_FLOAT
CellDataType = MET_FLOAT
NCellTypes = 2
PointDim = ID x y ...
NPoints = 4
Points =
0 1 1 0
1 100 1 0
2 100 225 0
3 1 225 0
CellType = VRT
NCells = 4
Cells =
4 0
5 1
6 2
7 3
CellType = LNE
NCells = 4
Cells =
0 0 1
1 1 2
2 2 3
3 3 0

Chris McIntosh wrote:
> Hello,
> 
>    I'm trying to use an itk meshSpatialObject as the geometrical
> component to a deformable model. I can deform the mesh without problem, but I'm
> unable to correctly view the mesh using a spatialObject viewer. Also, trying to
> use the spatialObjectToImageFilter results in incorrect output. In fact,
> checking to see if point 25,25 is inline yields the wrong answer as well.
> 
> I've included the meta file below that gets output from the
> spatialObjectWriter.
> 
> Basically, I create a rectangle using 4 points, 4 vertices, and 4 lines. The
> reason I'm using lines is that I don't want to view a surface, I want to
> view a wire-frame mesh. Basically, I'm deforming a spring-mass model and these 
> lines correspond to springs connecting the masses.
> 
> Any help on this issue would be greatly appreciated. Of course, if it will
> help I can provide code chunks.
> 
> Regards,
> 
> Chris McIntosh
> MSc. candidate
> Simon Fraser University
> http://mial.fas.sfu.ca
> 
> 
> -----------------------------------------
> ObjectType = Scene
> NDims = 2
> NObjects = 1
> ObjectType = Mesh
> NDims = 2
> ID = 0
> TransformMatrix = 1 0 0 1
> Offset = 0 0
> CenterOfRotation = 0 0
> ElementSpacing = 1 1
> PointType = MET_FLOAT
> PointDataType = MET_FLOAT
> CellDataType = MET_FLOAT
> NCellTypes = 2
> PointDim = ID x y ...
> NPoints = 4
> Points =
> 0 1 1
> 1 100 1
> 2 100 225
> 3 1 225
> CellType = VRT
> NCells = 4
> Cells =
> 4 0
> 5 1
> 6 2
> 7 3
> CellType = LNE
> NCells = 4
> Cells =
> 0 0 1
> 1 1 2
> 2 2 3
> 3 3 0
> 
> 
> _______________________________________________
> Insight-users mailing list
> Insight-users at itk.org
> http://www.itk.org/mailman/listinfo/insight-users
> 
>