Hi David,<br><br>So I don't use Paraview that often, usually to quick prototype something like this project. I've created a sphere, and run that glyph through PolyDataNormals so that I get point normals. I must be doing something incorrect, because when I try to run PointSetSurfaceReconstructionExample on my file, I get: <br>
<br>vtkPointSetSurfaceReconstruction (0x8a5f268): The input to vtkPointSetSurfaceReconstruction must have point normals. <br><br>I've also run the Input.vtp file from the source folder as input, which does run successfully, but the output file (attached) is unreadable by VMTK's vmtkimagereader script, which outputs the following error:<br>
<br>Reading VTK XML image file.<br>ERROR: In /usr/lib/vtk/VTK/IO/vtkXMLReader.cxx, line 601<br>vtkXMLImageDataReader (0x9d733e0): Cannot find ImageData element in file.<br><br>Since VMTK is my destination, this amounts to 'game over' for me unless I can determine which ImageData element isn't in the Input.vti image. Any suggestions would be appreciated!<br>
<br><div class="gmail_quote">On Mon, Apr 19, 2010 at 12:17 PM, David Welch <span dir="ltr"><<a href="mailto:dmwelch@engineering.uiowa.edu" target="_blank">dmwelch@engineering.uiowa.edu</a>></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;">I'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'll upgrade to the nightly build and post the results on the midas journal page.<br>
<br>Thanks!<div><div></div><div><br><br><div class="gmail_quote">On Mon, Apr 19, 2010 at 11:54 AM, David Doria <span dir="ltr"><<a href="mailto:daviddoria%2Bvtk@gmail.com" target="_blank">daviddoria+vtk@gmail.com</a>></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><div class="gmail_quote">On Mon, Apr 19, 2010 at 12:39 PM, David Welch <span dir="ltr"><<a href="mailto:dmwelch@engineering.uiowa.edu" target="_blank">dmwelch@engineering.uiowa.edu</a>></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'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'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'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></div></div><div><div></div><div>-- <br>David Welch<br>Graduate Student<br>Dept. of Biomedical Engineering<br>University of Iowa<br>Lab: (319) 335-5279<br>
</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>