<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="margin:0 0 0 .8ex;border-left:1px #ccc solid;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>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">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><br>David</div>
</div>