Hi David,<div><br></div><div>Yes indeed, it is a triangle mesh... However as I mentioned vtkDelaunay3D seems to work on all of the other vertices of my mesh, except on the one included the previous e-mail... This seems so odd!</div>
<div><br></div><div>The problem is that I need the 3-D convex hull, since I later want to know if current vertex is enclosed by its neighbors&#39; convex hull... If I use delaunay2D (by projecting my mesh onto a flat surface) would still be possible to check if the point is contained in this convex hull.</div>
<div><br></div><div>Thanks,</div><div>Miguel<br><br><div class="gmail_quote">2011/7/26 David Gobbi <span dir="ltr">&lt;<a href="mailto:david.gobbi@gmail.com">david.gobbi@gmail.com</a>&gt;</span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hi Miguel,<br>
<br>
The vtkDelaunay3D filter is for tetrahedral meshes, but your mesh<br>
looks like a triangle mesh.  So my guess is that you should be using<br>
vtkDelaunay2D instead.  But since Delanay2D is a 2D filter, you might<br>
have to project your mesh onto a flat surface first.<br>
<br>
Also, Delaunay often performs worse on noiseless data.  If you are<br>
working with synthetic data, you will get better results if you<br>
perturb each point with a small, random displacement.  With noiseless<br>
data, it is often the case that three points will lie along a straight<br>
line, and the Delaunay algorithm will consider three such points as<br>
forming a degenerate triangle.  Another problem with noiseless data,<br>
specifically for rectangular grids, is that they can result in<br>
situations where there is more than one valid Delaunay triangulation.<br>
These will cause the algorithm to go into deep recursion.<br>
<br>
 - David<br>
<br>
<br>
2011/7/25 Miguel Sotaquirá &lt;<a href="mailto:msotaquira@gmail.com">msotaquira@gmail.com</a>&gt;:<br>
<div><div></div><div class="h5">&gt; Hi,<br>
&gt; I&#39;m using vtkDelaunay3D to create the local (1-ring neighborhood) convex<br>
&gt; hull for each point in a mesh. For now I&#39;m using a synthetic, perfectly<br>
&gt; regular, noiseless mesh, where each vertex has exactly 6 neighbors.<br>
&gt; I&#39;m having troubles creating this convex hull only on one of the vertices<br>
&gt; (the others are computed correctly). In this particular case the convex hull<br>
&gt; computed does NOT include one of the six neighbors, as seen on the image<br>
&gt; below:<br>
&gt; <a href="http://tinypic.com/r/2944hag/7" target="_blank">http://tinypic.com/r/2944hag/7</a><br>
&gt; the purple point being the center of the neighborhood, the gray spheres its<br>
&gt; neighbors and the gray lines the convex hull obtained from vtkDelaunay3D.<br>
&gt; Here&#39;s the code I&#39;m using (neighbors are introduced as a polydata named<br>
&gt; &quot;pointcloud&quot;):<br>
&gt; vtkSmartPointer&lt;vtkDelaunay3D&gt; delaunay3D =<br>
&gt; vtkSmartPointer&lt;vtkDelaunay3D&gt;::New();<br>
&gt; delaunay3D-&gt;SetInput(pointCloud);<br>
&gt; delaunay3D-&gt;Update();<br>
&gt; during execution I get this warning: &quot;vtkDelaunay3D (0x10345a490): 1<br>
&gt; degenerate triangles encountered, mesh quality suspect&quot;<br>
&gt; How come, since I&#39;m using a synthetic and perfectly regular mesh? As I<br>
&gt; mentioned before, for other points (which have the exact same<br>
&gt; characteristics) the convex hull is computed correctly.... Also, I&#39;ve tried<br>
&gt; using the methods &quot;SetTolerance&quot; and &quot;BoundingTriangulationOn&quot; of<br>
&gt; vtkDelaunay3D but with no success<br>
&gt; Thanks!<br>
&gt; Miguel<br>
</div></div>&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; Please keep messages on-topic and check the VTK FAQ at:<br>
&gt; <a href="http://www.vtk.org/Wiki/VTK_FAQ" target="_blank">http://www.vtk.org/Wiki/VTK_FAQ</a><br>
&gt;<br>
&gt; Follow this link to subscribe/unsubscribe:<br>
&gt; <a href="http://www.vtk.org/mailman/listinfo/vtkusers" target="_blank">http://www.vtk.org/mailman/listinfo/vtkusers</a><br>
&gt;<br>
&gt;<br>
</blockquote></div><br></div>