Hello Kevin,<br><br>Thanks a lot for the help. I have one more question. <br><br>If I have a point coordinate is it possible to find its pointid and cell id of the cell in which the point lies? If yes then how?<br><br>Amitesh
<br><br><div><span class="gmail_quote">On 8/18/06, <b class="gmail_sendername">Kevin  H. Hobbs</b> &lt;<a href="mailto:hobbsk@ohiou.edu">hobbsk@ohiou.edu</a>&gt; wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
On Fri, 2006-08-18 at 00:42 -0500, Amitesh Agarwal wrote:<br>&gt; Hi all,<br>&gt;<br>&gt; I have read an unstructured grid and have converted the grid in<br>&gt; polydata. I have been trying to access individual cell information and
<br>&gt; point information serially in the polydata but to no avail. The<br>&gt; problem is that all of these have some kind of vtkid assigned to each<br>&gt; of the cell and point. I have grown immensly frustrated trying to find
<br>&gt; a way to access a list of those vtkids so that I could be able to use<br>&gt; them to access the individual cells and points. Any ideas how to do<br>&gt; it? Any help will be greatly appreciated.<br>&gt;<br>&gt; Thanks
<br>&gt;<br>&gt; Amitesh<br><br><br><br><br>// Make a reader.<br>vtkXMLPolyDataReader * reader = vtkXMLPolyDataReader::New();<br>reader-&gt;SetFileName( in_poly_file_name );<br>reader-&gt;Update();<br><br>in_poly = reader-&gt;GetOutput();
<br>in_poly-&gt;BuildLinks();<br><br>// Get the cells and points.<br>vtkCellArray * in_lines&nbsp;&nbsp;= in_poly-&gt;GetLines();<br>vtkPoints&nbsp;&nbsp;&nbsp;&nbsp;* in_points = in_poly-&gt;GetPoints();<br><br>// How many are there?<br>vtkIdType num_cells&nbsp;&nbsp;= in_lines-&gt;GetNumberOfCells();
<br>vtkIdType num_points = in_points-&gt;GetNumberOfPoints();<br><br>// Visit the points.<br>vtkIdType point;<br>for (point = 0; point &lt; num_points; ++point )<br>{<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;// What is the point's location?<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;double * pcoord;
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;pcoord = in_poly-&gt;GetPoint( point );<br><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;std::cout &lt;&lt; pcord[0] &lt;&lt; &quot; &quot;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&lt;&lt; pcord[1] &lt;&lt; &quot; &quot;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&lt;&lt; pcord[2] &lt;&lt; std::endl;<br>
}<br><br><br></blockquote></div><br>