I&#39;ve tried to compile the code, and it seems that one needs to have a version of VTK greater than the latest release to run this without compilation errors.  Specifically, you use vtkMath::Subtract, which is not included in 5.4.2.  I&#39;ll upgrade to the nightly build and post the results on the midas journal page.<br>
<br>Thanks!<br><br><div class="gmail_quote">On Mon, Apr 19, 2010 at 11:54 AM, David Doria <span dir="ltr">&lt;<a href="mailto:daviddoria%2Bvtk@gmail.com">daviddoria+vtk@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div class="im"><div class="gmail_quote">On Mon, Apr 19, 2010 at 12:39 PM, David Welch <span dir="ltr">&lt;<a href="mailto:dmwelch@engineering.uiowa.edu" target="_blank">dmwelch@engineering.uiowa.edu</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">

Thank you!   I&#39;ll let you know how it works out.  How can I add this to my Paraview setup?<div><div></div><div></div></div></blockquote></div><div><br></div></div>Let&#39;s keep the discussion on the mailing list so everyone can input and benefit.<div>

<br><div>As it was part of a bigger submission (surface reconstruction), I didn&#39;t wrap vtkVoxelizePolyData directly. I just tried to quickly make it into a Paraview plugin, but it seems I have separated the code a bit differently than I remember.</div>

<div><br></div><div>vtkVoxelizePolyData simply creates the vtkImageData with reasonable bounds over the polydata. Computing the distance to the surface was actually done in vtkPointSetSurfaceReconstrion, and then I used a vtkContourFilter to extract the 0 level set (the surface). You can definitely reuse most of this code - you just have to check the sign on the SignedDistanceToPlane function and mark the point as inside or outside.</div>

<div><br></div><div>You could also ust vtkVoxelizePolyData to make the grid, then use something like this:</div><div><a href="http://www.vtk.org/Wiki/VTK/Examples/PointInsideObject" target="_blank">http://www.vtk.org/Wiki/VTK/Examples/PointInsideObject</a></div>

<div><br></div><div>on every point to do the inside/outside marking.</div>
<div><br></div><div>Good luck,<br><font color="#888888"><br>David</font></div>
</div>
</blockquote></div><br><br clear="all"><br>-- <br>David Welch<br>Graduate Student<br>Dept. of Biomedical Engineering<br>University of Iowa<br>Lab: (319) 335-5279<br>