[Insight-users] Force field & mesh3D?

[Luis Ibanez]

Hi Marie,


Thanks for your clarification.

This is good news for you because deforming the Mesh
under a final displacement field is much easier than
running an iterative deformable model.


Your request was logged as Bug # 1333
http://www.itk.org/Bug/bug.php?op=show&bugid=1333&pos=0

and in response we added a new class

    itkWarpMeshFilter

in the directory

    Insight/Code/BasicFilters

This class takes two inputs:

   1) an itk::Mesh to be deformed
   2) an itk::Image of vectors

The image of vectors represents the deformation field.

the coordinates of the mesh are modified by mapping them
using the vectors of the deformation field.

For an example on how to use this class you may want to
look at the test


   Insight/Testing/Code/BasicFilters/
                         itkWarpMeshFilterTest.cxx


Please let us know if you encounter any problems,


    Thanks


       Luis



-------------------------
Danoush Marie Jo wrote:

> Hi luis,
> 
> Thanks for your answer!!
> 
>  
> 
> In fact I have the next problem that you were supposed!!
> 
> First, I’ve an initial segmentation of an image (my object 3D). In the 
> other hand, I’ve a final deformation field resulting from register the 
> image against another image. This field, I get it from 
> [registration->GetLastTransformParameters()]
> 
> Note that, I’ve many files of deformation for each frame time!!    
> 
>  
>E
> Actually, my problem is to deform this mesh by used the deformation 
> field (from the image registration)?? Could you help me by an example 
> code for doing that?
> 
>  
> 
> Many many thanks in advance
> 
> Marie
> 
>  
> 
> 
> 
> */Luis Ibanez <luis.ibanez at kitware.com>/* wrote:
> 
> 
>     Hi Marie,
> 
> 
>     You are right,
> 
>     You can use the itkDeformableMesh3DFilter for the purpose of
>     deforming your Mesh under the influence of a deformation field.
> 
>     Note that you also have the option of using the
> 
>     itk::DeformableSimplexMesh3DBalloonForceFilter
> 
>     http://www.itk.org/Insight/Doxygen/html/classitk_1_1DeformableSimplexMesh3DBalloonForceFilter.html
> 
>     contributed by Thomas Boettger, and further refined
>     by Leila Baghdadi.
> 
> 
> 
>     In fact, you may find useful to use directly the GUI
>     application in
> 
> 
>     InsightApplications/
>     DeformableModelSimplexMesh
> 
>     Where this filter is combined with the visualization and
>     GUI elements that make easier to use it.
> 
> 
>     On a second thought...
>     Note that it is different to deform the Mesh using a *final*
>     deformation field resulting from the image registration than!
>     to use a *force field* like the one resulting form a gradient
>     computation.
> 
>     If you actually have a *final* deformation field, your problem
>     then seems to be similar to an atlas based segmentation where
>     you have an initial segmentation (your mesh) of an image, then
>     you register the image against another image, and finally you
>     use the deformation field (from the image registration) in order
>     to map the initial mesh.
> 
>     If this is what you are doing, then you don't need the Deformable
>     filter at all, you can just use the deformation field as a transform
>     for Warping the mesh. For that purpose you can use the filter:
> 
>     http://www.itk.org/Insight/Doxygen/html/classitk_1_1TransformMeshFilter.html
> 
> 
> 
>     Please let us know which one of the two scenarios above are you
>     dealing with.
> 
> 
>     Thanks
> 
> 
>     Luis
> 
> 
> 
> 
>     ------------
>      > Hi,
>      >
>      > I have an initial mesh (from vtkPolyDATA).
>      >
>      > I characterized the displacements field from deformed medical
>     images (itkRegistration4).
>      >
>      > Therefore, I'd like to use forces derived from the deformed image
>     (last resultants) to deform the mesh, so the
>     itkDeformableMesh3DFilter seems to be the answer to this problem. Is
>     it?
>      >
>      > Could any body use this itk tool to do something like this?
>      >
>      > Is somebody could help me by an example code for doing this
>      > purpose?
>      >
>      > Many thanks in advance
>      >
>      > Marie.
>      >
>      >
>      >
> 
> 
> 
> ------------------------------------------------------------------------
> Créez gratuitement votre Yahoo! Mail avec *100 Mo de stockage !*
> Créez votre Yahoo! Mail 
> <http://fr.rd.yahoo.com/mail/taglines/*http://fr.rd.yahoo.com/evt=25917/*http://fr.rd.yahoo.com/mail/mail_taglines_100/default/*http://fr.benefits.yahoo.com/>
> 
> *Le nouveau Yahoo! Messenger est arrivé !* Découvrez toutes les 
> nouveautés pour dialoguer instantanément avec vos amis. Téléchargez 
> GRATUITEMENT ici ! 
> <http://fr.rd.yahoo.com/mail/taglines/*http://fr.rd.yahoo.com/evt=26111/*http://fr.rd.yahoo.com/messenger/mail_taglines/default/*http://fr.messenger.yahoo.com> 
>