please include the mailing list in your reply.<br><br><div class="gmail_quote">On Tue, Aug 14, 2012 at 11:26 AM, Agata Krasoń <span dir="ltr">&lt;<a href="mailto:agatakrason@gmail.com" target="_blank">agatakrason@gmail.com</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>Dear Bill,</div><div><br></div><div>Thank You for reply.</div>I have already resolve this problem. <div><br></div>
<div>It was here  :</div><div><br></div><div><div class="im"> for (vtkIdType i = 0; i &lt; number_of_triangles; i++)<br>
    {<br>      vtkIdType a, b, c;<br>      infile &gt;&gt; a &gt;&gt; b &gt;&gt; c;<br>      polys-&gt;InsertNextCell(3);<br></div>    <font color="#ff0000">  polys-&gt;InsertCellPoint(a -1);<br>      polys-&gt;InsertCellPoint(b - 1);<br>

      polys-&gt;InsertCellPoint(c - 1);</font><br><br>    } </div><div><br></div><div><br></div><div>But I have another problem with meshes. Sorry for very basic questions. I am beginig my adventure working  </div><div>with meshes.</div>

<div>I don&#39;t have lots of knowledge  in this area.</div><div>Now I need to display a mesh with attributes in VTK.</div><div>I have already display mesh in vtk, not with attributes.</div><div>Could You help me please ?</div>

<div><br></div><div>My code :</div><div><div>#include &lt;vtkTable.h&gt;</div><div>#include &quot;vtkPolyData.h&quot;</div><div>#include &lt;vtkPointData.h&gt;</div><div>#include &quot;vtkPoints.h&quot;</div><div>#include &quot;vtkCellArray.h&quot;</div>

<div>#include &lt;vtkLine.h&gt;</div><div>#include &lt;vtkFloatArray.h&gt;</div><div>#include &lt;vtkLookupTable.h&gt;</div><div>#include &lt;vtkUnstructuredGrid.h&gt;</div><div><br></div><div>#include &lt;vtkGeometryFilter.h&gt;</div>

<div><br></div><div><br></div><div><br></div><div>#include &lt;vtkParticleReader.h&gt;</div><div>#include &lt;vtkPolyDataReader.h&gt;</div><div>#include &lt;vtkPolyDataWriter.h&gt;</div><div>#include &lt;vtkDelimitedTextReader.h&gt;</div>

<div>#include &lt;vtkUnstructuredGridWriter.h&gt;</div><div>#include &lt;vtkPolyDataWriter.h&gt;</div><div>#include &lt;vtkProperty.h&gt;</div><div><br></div><div>#include &lt;vtkPolyDataMapper.h&gt;</div><div>#include &lt;vtkActor.h&gt;</div>

<div>#include &lt;vtkRenderer.h&gt;</div><div>#include &lt;vtkRenderWindow.h&gt;</div><div>#include &lt;vtkRenderWindowInteractor.h&gt;</div><div class="im"><div><br></div><div>using namespace std;</div><div><br></div><div>
<br></div><div>
int main(int argc, char* argv[])</div><div>{</div><div><br></div></div><div><span style="white-space:pre-wrap">        </span>// file with mesh</div><div class="im"><div>    const char* filename = &quot;heartmesh.txt&quot;;</div>
<div>
<span style="white-space:pre-wrap">        </span>std::ifstream infile(filename);</div><div>    vtkIdType number_of_points, number_of_triangles;</div><div>    infile &gt;&gt; number_of_points &gt;&gt; number_of_triangles;</div>

<div>    vtkPoints* points = vtkPoints::New();</div><div>    points-&gt;SetNumberOfPoints(number_of_points);</div><div><br></div></div><div><span style="white-space:pre-wrap">        </span>//std::cout&lt;&lt;&quot;number_of_points: &quot;&lt;&lt;number_of_points&lt;&lt;std::endl;</div>

<div><span style="white-space:pre-wrap">        </span>//std::cout&lt;&lt;&quot;number_of_triangles: &quot; &lt;&lt;number_of_triangles&lt;&lt;std::endl;</div><div><br></div><div>    for (vtkIdType i = 0; i &lt; number_of_points; i++)</div>
<div class="im">
<div>     {</div><div>        double x, y, z;</div><div>        infile &gt;&gt; x &gt;&gt; y &gt;&gt; z;</div></div><div>        points-&gt;SetPoint(i,x, y, z);</div><div class="im"><div><span style="white-space:pre-wrap">        </span>  </div>

<div>     }</div><div> </div><div>    vtkCellArray* polys = vtkCellArray::New();</div><div><br></div><div>   for (vtkIdType i = 0; i &lt; number_of_triangles; i++)</div><div>    {</div><div>      vtkIdType a, b, c;</div>
<div>
      infile &gt;&gt; a &gt;&gt; b &gt;&gt; c;</div><div>      polys-&gt;InsertNextCell(3);</div></div><div>      polys-&gt;InsertCellPoint(a-1);</div><div>      polys-&gt;InsertCellPoint(b-1);</div><div>      polys-&gt;InsertCellPoint(c-1);</div>

<div><span style="white-space:pre-wrap">        </span> </div><div>    }</div><div><br></div><div><br></div><div>   vtkUnstructuredGrid* grid = vtkUnstructuredGrid::New();</div><div>   grid-&gt;SetPoints(points);</div>
<div>   grid-&gt;SetCells(5, polys);</div><div>   grid-&gt;Update();</div><div> </div><div>   // Here I read a file with attributes scalars</div><div>   vtkDelimitedTextReader *dread = vtkDelimitedTextReader::New();</div>

<div>   dread-&gt;SetFileName(&quot;CSVFile.csv&quot;);</div><div>   dread-&gt;Update();</div><div> </div><div>   vtkTable *tdata = dread-&gt;GetOutput();</div><div><br></div><div>   vtkFloatArray *fd = vtkFloatArray::New();</div>

<div><br></div><div>   vtkIdType nr = tdata-&gt;GetNumberOfRows();</div><div><br></div><div>   fd-&gt;SetNumberOfTuples(nr);</div><div><br></div><div>   for (vtkIdType i=0; i&lt;nr; i++) {</div><div>    std::cout &lt;&lt;  tdata-&gt;GetValue(i, 0).ToFloat() &lt;&lt; &quot;\n&quot;; </div>

<div>    fd-&gt;SetTuple1(i, tdata-&gt;GetValue(i, 0).ToFloat());</div><div>   }</div><div><br></div><div>  vtkPointData *pd = grid-&gt;GetPointData();</div><div>  pd-&gt;SetScalars(fd);</div><div>  pd-&gt;Update();</div>

<div><br></div><div>  // Save in file mesh&amp; attributes</div><div>  vtkUnstructuredGridWriter *writer = vtkUnstructuredGridWriter::New();</div><div>  writer-&gt;SetInput(grid);</div><div>  writer-&gt;SetFileName(&quot;outputMeshWithAttributes.vtk&quot;);</div>

<div>  writer-&gt;Write();</div><div><br></div><div>  vtkGeometryFilter* geometryFilter = vtkGeometryFilter::New();</div><div>  geometryFilter-&gt;SetInput(grid);</div><div>  geometryFilter-&gt;Update();</div><div><br></div>

<div>  vtkPolyData* polydata = geometryFilter-&gt;GetOutput();</div><div><br></div><div>  // Here I need to display mesh with attributes !</div><div>  // Visualization &amp; Display </div><div class="im"><div>  vtkSmartPointer&lt;vtkPolyDataMapper&gt; mapper = vtkSmartPointer&lt;vtkPolyDataMapper&gt;::New();</div>

</div><div>  mapper-&gt;SetInput(polydata);</div><div><br></div><div>  //mapper-&gt;SetLookupTable(lut);</div><div>  mapper-&gt;SetColorModeToMapScalars();</div><div class="im"><div>  vtkSmartPointer&lt;vtkActor&gt; actor = vtkSmartPointer&lt;vtkActor&gt;::New();</div>

<div>  actor-&gt;SetMapper(mapper);</div><div> </div><div><br></div><div>  //Create a renderer, render window, and interactor</div><div>  vtkSmartPointer&lt;vtkRenderer&gt; renderer = vtkSmartPointer&lt;vtkRenderer&gt;::New();</div>

<div>  vtkSmartPointer&lt;vtkRenderWindow&gt; renderWindow = vtkSmartPointer&lt;vtkRenderWindow&gt;::New();</div><div>  renderWindow-&gt;AddRenderer(renderer);</div><div>  vtkSmartPointer&lt;vtkRenderWindowInteractor&gt; renderWindowInteractor = vtkSmartPointer&lt;vtkRenderWindowInteractor&gt;::New();</div>

<div>  renderWindowInteractor-&gt;SetRenderWindow(renderWindow);</div><div><br></div><div>  renderer-&gt;AddActor(actor);</div><div>  renderer-&gt;SetBackground(0, 0, 0); </div><div>  renderWindow-&gt;Render();</div><div>

  renderWindowInteractor-&gt;Start();</div><div><br></div><div><br></div></div><div>  /*polydata-&gt;Delete();</div><div>  geometryFilter-&gt;Delete();</div><div>  writer-&gt;Delete();</div><div>  pd-&gt;Delete();</div><div>
  tdata-&gt;Delete();</div>
<div>  fd-&gt;Delete();</div><div>  dread-&gt;Delete();</div><div>  polys-&gt;Delete();</div><div>  points-&gt;Delete();*/</div><div> </div><div> </div><div> </div><div>  return EXIT_SUCCESS;</div><div>}</div><div><br></div>

<div><br></div></div><div><br></div><div>I received display mesh&amp;attributes only in paraview. I attach a screen.</div><div><br></div><div>agattte</div><div class="HOEnZb"><div class="h5"><div><br></div><div><br><br><div class="gmail_quote">
2012/8/14 Bill Lorensen <span dir="ltr">&lt;<a href="mailto:bill.lorensen@gmail.com" target="_blank">bill.lorensen@gmail.com</a>&gt;</span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Sorry, I cannot read the links to your files.<br><br><div class="gmail_quote">On Tue, Aug 14, 2012 at 8:08 AM, agatte <span dir="ltr">&lt;<a href="mailto:agatakrason@gmail.com" target="_blank">agatakrason@gmail.com</a>&gt;</span> wrote:<br>


<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi ;)<br>
I have a question about a mesh.<br>
I read mesh data from txt file ( I attach this file.)<br>
But I can&#39;t receive good structure of mesh (polydata).<br>
I receive  something like this : (I attach a photo screen )<br>
When I use vtkVertexGlyphFilter. I receive a point cloud with good shape of<br>
heart.<br>
But I need display mesh with triangles/triangleStrips.<br>
Could You look at this code and this file ?<br>
I would appreciate for any help please.<br>
I don&#39;t see an error. What is wrong with this code ?<br>
<br>
My code is here :<br>
#include &lt;vtkVersion.h&gt;<br>
#include &lt;vtkSmartPointer.h&gt;<br>
#include &lt;vtkPolyDataMapper.h&gt;<br>
#include &lt;vtkActor.h&gt;<br>
#include &lt;vtkParticleReader.h&gt;<br>
#include &lt;vtkRenderWindow.h&gt;<br>
#include &lt;vtkRenderWindowInteractor.h&gt;<br>
#include &lt;vtkRenderer.h&gt;<br>
#include &lt;vtkVertexGlyphFilter.h&gt;<br>
#include &lt;sstream&gt;<br>
#include &quot;vtkPolyData.h&quot;<br>
#include &quot;vtkPoints.h&quot;<br>
#include &quot;vtkCellArray.h&quot;<br>
#include &lt;iostream&gt;<br>
#include &lt;fstream&gt;<br>
#include &lt;vtkTriangle.h&gt;<br>
#include &lt;vtkTriangleFilter.h&gt;<br>
#include &lt;vtkContourFilter.h&gt;<br>
#include &lt;vtkDelaunay3D.h&gt;<br>
#include &lt;vtkXMLPolyDataWriter.h&gt;<br>
#include &lt;vtkXMLUnstructuredGridWriter.h&gt;<br>
#include &lt;vtkXMLPolyDataWriter.h&gt;<br>
#include &lt;vtkDataSetSurfaceFilter.h&gt;<br>
#include &lt;vtkPolyDataWriter.h&gt;<br>
#include &lt;vtkExtractEdges.h&gt;<br>
#include &lt;vtkTriangleStrip.h&gt;<br>
#include &lt;vtkLine.h&gt;<br>
#include &lt;vtkXMLPolyDataReader.h&gt;<br>
#include &lt;vtkGeometryFilter.h&gt;<br>
#include &lt;vtkSurfaceReconstructionFilter.h&gt;<br>
#include &lt;vtkDelaunay2D.h&gt;<br>
#include &lt;vtkExtractEdges.h&gt;<br>
#include &lt;vtkFloatArray.h&gt;<br>
#include &lt;vtkLookupTable.h&gt;<br>
#include &lt;vtkStripper.h&gt;<br>
#include &lt;iostream&gt;<br>
#include &lt;vector&gt;<br>
#include &lt;string&gt;<br>
#include &lt;fstream&gt;<br>
#include &quot;vtkProperty.h&quot;<br>
<br>
using namespace std;<br>
<br>
<br>
int main(int argc, char* argv[])<br>
{<br>
<br>
    const char* filename = &quot;heartmesh.txt&quot;;<br>
        std::ifstream infile(filename);<br>
    vtkIdType number_of_points, number_of_triangles;<br>
    infile &gt;&gt; number_of_points &gt;&gt; number_of_triangles;<br>
    vtkPoints* points = vtkPoints::New();<br>
    points-&gt;SetNumberOfPoints(number_of_points);<br>
<br>
        std::cout&lt;&lt;&quot;number_of_points: &quot;&lt;&lt;number_of_points&amp;lt;&amp;lt;std::endl;<br>
        std::cout&amp;lt;&amp;lt;&amp;quot;number_of_triangles: &amp;quot;<br>
&amp;lt;&amp;lt;number_of_triangles&amp;lt;&amp;lt;std::endl;<br>
<br>
    for (vtkIdType i = 0; i &amp;lt; number_of_points; i++)<br>
     {<br>
        double x, y, z;<br>
        infile &gt;&gt; x &gt;&gt; y &gt;&gt; z;<br>
        points-&gt;SetPoint(i, x, y, z);<br>
<br>
     }<br>
<br>
    vtkCellArray* polys = vtkCellArray::New();<br>
<br>
   for (vtkIdType i = 0; i &lt; number_of_triangles; i++)<br>
    {<br>
      vtkIdType a, b, c;<br>
      infile &gt;&gt; a &gt;&gt; b &gt;&gt; c;<br>
      polys-&gt;InsertNextCell(3);<br>
      polys-&gt;InsertCellPoint(a);<br>
      polys-&gt;InsertCellPoint(b);<br>
      polys-&gt;InsertCellPoint(c);<br>
<br>
    }<br>
<br>
  vtkPolyData* polydata = vtkPolyData::New();<br>
  polydata-&gt;SetPoints(points);<br>
  polydata-&gt;SetPolys(polys);<br>
  polydata-&gt;Update();<br>
<br>
  vtkVertexGlyphFilter* glyphFilter = vtkVertexGlyphFilter::New();<br>
  glyphFilter-&gt;SetInputConnection(polydata-&gt;GetProducerPort());<br>
  glyphFilter-&gt;Update();<br>
<br>
  vtkDelaunay2D* delaunay = vtkDelaunay2D::New();<br>
  delaunay-&gt;SetInputConnection(glyphFilter-&gt;GetOutputPort());<br>
  delaunay-&gt;Update();<br>
<br>
  vtkDataSetSurfaceFilter* surfaceFilter = vtkDataSetSurfaceFilter::New();<br>
  surfaceFilter-&gt;SetInputConnection(delaunay-&gt;GetOutputPort());<br>
  surfaceFilter-&gt;Update();<br>
<br>
  vtkTriangleFilter* triangleFilter = vtkTriangleFilter::New();<br>
  triangleFilter-&gt;SetInputConnection(surfaceFilter-&gt;GetOutputPort());<br>
  triangleFilter-&gt;Update();<br>
<br>
  // Visualization<br>
  vtkSmartPointer&lt;vtkPolyDataMapper&gt; mapper =<br>
vtkSmartPointer&lt;vtkPolyDataMapper&gt;::New();<br>
 // mapper-&gt;SetInput(polydata);<br>
  mapper-&gt;SetInputConnection(triangleFilter-&gt;GetOutputPort());<br>
  vtkSmartPointer&lt;vtkActor&gt; actor = vtkSmartPointer&lt;vtkActor&gt;::New();<br>
  actor-&gt;SetMapper(mapper);<br>
<br>
<br>
<br>
  //Create a renderer, render window, and interactor<br>
  vtkSmartPointer&lt;vtkRenderer&gt; renderer =<br>
vtkSmartPointer&lt;vtkRenderer&gt;::New();<br>
  vtkSmartPointer&lt;vtkRenderWindow&gt; renderWindow =<br>
vtkSmartPointer&lt;vtkRenderWindow&gt;::New();<br>
  renderWindow-&gt;AddRenderer(renderer);<br>
  vtkSmartPointer&lt;vtkRenderWindowInteractor&gt; renderWindowInteractor =<br>
vtkSmartPointer&lt;vtkRenderWindowInteractor&gt;::New();<br>
  renderWindowInteractor-&gt;SetRenderWindow(renderWindow);<br>
<br>
  renderer-&gt;AddActor(actor);<br>
  renderer-&gt;SetBackground(0, 0, 0);<br>
  renderWindow-&gt;Render();<br>
  renderWindowInteractor-&gt;Start();<br>
<br>
  polydata-&gt;Delete();<br>
  polys-&gt;Delete();<br>
  points-&gt;Delete();<br>
  glyphFilter-&gt;Delete();<br>
  delaunay-&gt;Delete();<br>
  triangleFilter-&gt;Delete();<br>
  surfaceFilter-&gt;Delete();<br>
<br>
<br>
<br>
<br>
  return EXIT_SUCCESS;<br>
}<br>
<br>
<br>
<br>
agatte<br>
<br>
<a href="http://vtk.1045678.n5.nabble.com/file/n5715236/heartmesh.txt" target="_blank">http://vtk.1045678.n5.nabble.com/file/n5715236/heartmesh.txt</a> heartmesh.txt<br>
<br>
<br>
<br>
<a href="http://vtk.1045678.n5.nabble.com/file/n5715236/Heart_as_Polydata.png" target="_blank">http://vtk.1045678.n5.nabble.com/file/n5715236/Heart_as_Polydata.png</a><br>
<a href="http://vtk.1045678.n5.nabble.com/file/n5715236/Heart_as_PointCloud.png" target="_blank">http://vtk.1045678.n5.nabble.com/file/n5715236/Heart_as_PointCloud.png</a><br>
<br>
<br>
<br><span><font color="#888888">
<br>
--<br>
View this message in context: <a href="http://vtk.1045678.n5.nabble.com/question-about-mesh-tp5715236.html" target="_blank">http://vtk.1045678.n5.nabble.com/question-about-mesh-tp5715236.html</a><br>
Sent from the VTK - Users mailing list archive at Nabble.com.<br>
_______________________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the VTK FAQ at: <a href="http://www.vtk.org/Wiki/VTK_FAQ" target="_blank">http://www.vtk.org/Wiki/VTK_FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.vtk.org/mailman/listinfo/vtkusers" target="_blank">http://www.vtk.org/mailman/listinfo/vtkusers</a><br>
</font></span></blockquote></div><span><font color="#888888"><br><br clear="all"><br>-- <br>Unpaid intern in BillsBasement at noware dot com<br><br>
</font></span></blockquote></div><br></div>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Unpaid intern in BillsBasement at noware dot com<br><br>