Hi Vincent,<div><br></div><div>This is great news, just downloaded the software. It would be great sending me the code snippet for converting the mesh.</div><div>What code do you use for surface extraction?</div><div><br></div>
<div>Thx,</div><div><br></div><div>S.<br><br><div class="gmail_quote">On Tue, Feb 23, 2010 at 10:39 PM, Vincent Magnotta <span dir="ltr">&lt;<a href="mailto:vincent-magnotta@uiowa.edu">vincent-magnotta@uiowa.edu</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Just to follow-up to Luis&#39;s e-mail.<br>
<br>
There is another solution to the mesh visualization question. What we have<br>
done in some of our code is to have a parallel copy of the mesh in VTK that<br>
we update the point locations based on the results of the FE analysis. This<br>
can then be written to disk and visualized and animated in Paraview.<br>
<br>
The software that was mentioned for mesh generation (IA-FEMesh) is both a 3D<br>
Slicer module and has a stand-alone version<br>
(<a href="http://www.ccad.uiowa.edu/mimx" target="_blank">http://www.ccad.uiowa.edu/mimx</a>). Given a surface, you can generate a<br>
hexahedral mesh and assign material properties. The mesh can be saved in VTK<br>
format where the material properties are stored as Cell data. You will need<br>
to write some code to create an ITK FE mesh from the resulting VTK mesh. I<br>
have an example if that will be helpful. The ITK FEM framework is an area<br>
where some refactoring is needed. Hopefully you will some work in this area<br>
with the push towards ITK4.<br>
<br>
Vince<br>
<div><div></div><div class="h5"><br>
<br>
On 2/23/10 4:23 PM, &quot;Luis Ibanez&quot; &lt;<a href="mailto:luis.ibanez@kitware.com">luis.ibanez@kitware.com</a>&gt; wrote:<br>
<br>
&gt; Hi Steve,<br>
&gt;<br>
&gt; Option (A)<br>
&gt;<br>
&gt; The original FEM code in ITK has a compile flag:<br>
&gt;<br>
&gt;           FEM_BUILD_VISUALIZATION<br>
&gt;<br>
&gt; but, given that Visualization was not in the scope of<br>
&gt; ITK, this code has probably never been exercised.    :-(<br>
&gt;<br>
&gt;<br>
&gt; Option (B)<br>
&gt;<br>
&gt; You could use the filters that overlap a grid on an<br>
&gt; image, and resample that grid image through the<br>
&gt; deformation field computed by the FEM registration<br>
&gt; method.<br>
&gt;<br>
&gt; See<br>
&gt;              Insight/Code/Review/<br>
&gt;                            itkGridImageSource.h<br>
&gt;<br>
&gt; and its explanation in the Insight Journal paper<br>
&gt;<br>
&gt;             &quot;Gridding Graphic Graticules&quot;<br>
&gt;             Tustison N., Avants B., Gee J.<br>
&gt;             University of Pennsylvania<br>
&gt;             <a href="http://hdl.handle.net/1926/475" target="_blank">http://hdl.handle.net/1926/475</a><br>
&gt; <a href="http://www.insight-journal.org/browse/publication/140" target="_blank">http://www.insight-journal.org/browse/publication/140</a><br>
&gt;<br>
&gt;<br>
&gt; -----<br>
&gt;<br>
&gt; Regarding the process of creating a customized<br>
&gt; Mesh:<br>
&gt;<br>
&gt;         We don&#39;t have such functionality in ITK.<br>
&gt;<br>
&gt; Mesh Generation is outside of the scope of ITK.<br>
&gt;<br>
&gt; You will find very useful to look at:<br>
&gt;<br>
&gt; <a href="http://www.na-mic.org/Wiki/index.php/NA-MIC_NCBC_Collaboration:Automated_FE_Me" target="_blank">http://www.na-mic.org/Wiki/index.php/NA-MIC_NCBC_Collaboration:Automated_FE_Me</a><br>
&gt; sh_Development<br>
&gt;<br>
&gt; as a possible solution.<br>
&gt;<br>
&gt; This is an Open Source application based<br>
&gt; on ITK and VTK<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;      Regards,<br>
&gt;<br>
&gt;<br>
&gt;                Luis<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; -------------------------------------------------------<br>
&gt; On Mon, Feb 22, 2010 at 9:52 PM, Steve Lancey &lt;<a href="mailto:steve.lancey@gmail.com">steve.lancey@gmail.com</a>&gt; wrote:<br>
&gt;&gt; Hi,<br>
&gt;&gt;<br>
&gt;&gt; I have two questions:<br>
&gt;&gt;<br>
&gt;&gt; 1) I would like to store the generated mesh after every iteration of the<br>
&gt;&gt; energy minimization procedure and visualize them later (as a binary images<br>
&gt;&gt; or paths).<br>
&gt;&gt; What is the best way to achieve this?<br>
&gt;&gt;<br>
&gt;&gt; 2) What is the best tool to create a mesh manually, including the assignment<br>
&gt;&gt; of mesh element properties?<br>
&gt;&gt;<br>
&gt;&gt; Thanks for taking your time,<br>
&gt;&gt;<br>
&gt;&gt; Steve Lancey<br>
&gt;&gt;<br>
&gt;&gt; _____________________________________<br>
&gt;&gt; Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
&gt;&gt;<br>
&gt;&gt; Visit other Kitware open-source projects at<br>
&gt;&gt; <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
&gt;&gt;<br>
&gt;&gt; Kitware offers ITK Training Courses, for more information visit:<br>
&gt;&gt; <a href="http://www.kitware.com/products/protraining.html" target="_blank">http://www.kitware.com/products/protraining.html</a><br>
&gt;&gt;<br>
&gt;&gt; Please keep messages on-topic and check the ITK FAQ at:<br>
&gt;&gt; <a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
&gt;&gt;<br>
&gt;&gt; Follow this link to subscribe/unsubscribe:<br>
&gt;&gt; <a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt; _____________________________________<br>
&gt; Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
&gt;<br>
&gt; Visit other Kitware open-source projects at<br>
&gt; <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
&gt;<br>
&gt; Kitware offers ITK Training Courses, for more information visit:<br>
&gt; <a href="http://www.kitware.com/products/protraining.html" target="_blank">http://www.kitware.com/products/protraining.html</a><br>
&gt;<br>
&gt; Please keep messages on-topic and check the ITK FAQ at:<br>
&gt; <a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
&gt;<br>
&gt; Follow this link to subscribe/unsubscribe:<br>
&gt; <a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
<br>
</div></div>----------------------<br>
Associate Professor<br>
Department of Radiology<br>
0453-D JCP<br>
200 Hawkins Drive<br>
Iowa City, IA 52242<br>
E-mail: <a href="mailto:vincent-magnotta@uiowa.edu">vincent-magnotta@uiowa.edu</a><br>
Phone: 319-356-8255 Fax: 319-353-6275<br>
Website: <a href="http://www.radiology.uiowa.edu" target="_blank">http://www.radiology.uiowa.edu</a><br>
<br>
<br>
<br>
</blockquote></div><br></div>