<div dir="ltr">Hi all,<div><br></div><div>Thanks for persevering with me. I got it working correctly with vtkPlaneSource.</div><div><br></div><div>The problem ended up being a PEBKAC with how I was using the Delaunay2D class: When converting from CGAL to vtkUnstructuredgrid, you need to index the face values and vertices separately while maintaining connectivity. Because you are defining the connectivity, this works. However! with vtkDelaunay2D, you can't do this. Rather, for each vertex that is added to 'points', the face data must be added 3-times. Which this isn't quite perfect, it should work fine, unless there is a way of using an unstructured mesh and a probe filter? <a href="http://www.vtk.org/Wiki/VTK/Examples/Cxx/PolyData/GeometryFilter">http://www.vtk.org/Wiki/VTK/Examples/Cxx/PolyData/GeometryFilter</a>   suggests you can? So I will have to play around with that more.</div><div><br></div><div>Thanks again, much appreciated.</div></div><div class="gmail_extra"><br clear="all"><div><div dir="ltr"><span style="color:rgb(102,102,102)">Chris Marsh</span><br style="color:rgb(102,102,102)"><span style="color:rgb(102,102,102)">PhD Student</span><br style="color:rgb(102,102,102)"><span style="color:rgb(102,102,102)"></span><a href="http://chrismarsh.ca/" style="color:rgb(102,102,102)" target="_blank">chrismarsh.ca</a><br style="color:rgb(102,102,102)"><br style="color:rgb(102,102,102)"><span style="color:rgb(102,102,102)">13 Kirk Hall</span><br style="color:rgb(102,102,102)"><span style="color:rgb(102,102,102)">University of Saskatchewan</span><br style="color:rgb(102,102,102)"></div></div>
<br><div class="gmail_quote">On Wed, Sep 17, 2014 at 11:00 AM, David Gobbi <span dir="ltr"><<a href="mailto:david.gobbi@gmail.com" target="_blank">david.gobbi@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">On Wed, Sep 17, 2014 at 10:45 AM, Chris Marsh <<a href="mailto:chris.marsh@usask.ca">chris.marsh@usask.ca</a>> wrote:<br>
> Hi David,<br>
><br>
> Thanks for the reply, this is slowly starting to make sense. That call did<br>
> in fact fix the bounds issues. However, the output still looks like it did<br>
> in the above screenshot.<br>
><br>
>>The plane is already composed of cells, so you should not run it<br>
> through Delaunay.<br>
><br>
> This is what perplexed me about the above linked example.<br>
><br>
> What is the best way to write the probe results to a file then?<br>
<br>
</span>At the beginning, I was hinting that you should probe with a<br>
vtkImageData, which you could have written out as an image<br>
(either as raw data, or as a png, etc).  But now that you have<br>
the probed data as a polydata plane, you'll have to get the<br>
scalars that are associated with the points that make up the<br>
polydata.  GetPointData()->GetScalars() will return an array<br>
that contains all the scalars.  See the VTK text for an in-depth<br>
description of how data is stored in a VTK data set.<br>
<span class="HOEnZb"><font color="#888888"><br>
 - David<br>
</font></span></blockquote></div><br></div>