Again, many many thanks. After following your suggestions (and looking at Mike&#39;s code which is really helpfull) this is what I&#39;m doing but my program crashes:<br><br>// &quot;oPolyDataDel&quot;&nbsp; is a vtkPolyData where I store&nbsp; the output of Delaunay filter .&nbsp; As I set the Aplha value,
<br>// some lines and points appear in the triangulation (I&#39;ve read in the documentation that this&nbsp; is to happen&nbsp; when&nbsp; you&nbsp; set<br>// alpha value) but I only want triangles.&quot;vTriangulos&quot; is an array of ints where I wish to store triangle values consecutively
<br>// So this is what I&#39;ve written:<br><br>&nbsp;&nbsp;&nbsp; int j=0;&nbsp;  <br>&nbsp;&nbsp;&nbsp; int npts,*pts;<br>&nbsp;&nbsp;&nbsp; vtkCellArray* oCellArr= vtkCellArray::New();<br>&nbsp;&nbsp;&nbsp; oCellArr=oPolyDataDel-&gt;GetPolys(); <br>&nbsp;&nbsp;&nbsp; &nbsp;for(int i=0;i&lt;oPolyDataDel-&gt;GetNumberOfPolys();i++)
<br>&nbsp;&nbsp;&nbsp; &nbsp;{<br>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;oCellArr-&gt;GetNextCell(npts,pts);<br>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;vTriangulos[j]= pts[0];<br>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;vTriangulos[j+1]= pts[1];<br>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;vTriangulos[j+2]= pts[2];<br><br>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; j+=3;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;<br>&nbsp;&nbsp;&nbsp; &nbsp;} 
<br><div><div><br>Is that correct? Again, Many thanks<br>&nbsp;</div><br></div><br>