<div dir="ltr"><div class="gmail_extra"><div><div>Mmh thanks. I think I'll stick with the probe filter.</div></div><div><br></div><div>I've been able to get  this example working</div><div><a href="http://www.itk.org/Wiki/VTK/Examples/Cxx/PolyData/InterpolateMeshOnGrid">http://www.itk.org/Wiki/VTK/Examples/Cxx/PolyData/InterpolateMeshOnGrid</a><br></div><div>if I have a real-world pt cloud, I'm a bit unclear how I should make the grid. <br></div><div><br></div><div>If xll and yll are the lower x and y coordinates of an axis aligned bounding box</div><div><br></div><div><div>size_t gridSize = 100;</div><div>    for ( size_t x = 0; x < gridSize; x++ )</div><div>    {</div><div>        for ( size_t y = 0; y < gridSize; y++ )</div><div>        {</div><div>            gridPoints->InsertNextPoint ( x+xll, y+yll, 0);</div><div>        }</div><div>    }</div></div><div><br></div><div>does not seem to produce any output.</div><div><br></div><div>Any suggestions?</div><div><br></div><div class="gmail_quote">On Mon, Sep 15, 2014 at 4:39 PM, 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:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Another filter that might work is vtkShepardMethod, which takes<br>
unstructured points (e.g. an arbitrary set of samples) as input and<br>
produces a vtkImageData as output.<br>
<br>
The difference between vtkProbeFilter and vtkShepardMethod is<br>
that the former performs the interpolation based on the connectivity<br>
of the points (i.e. based on the cells), while the latter ignores the<br>
connectivity and does the interpolation purely based on distance.<br>
<span class=""><font color="#888888"><br>
 - David<br>
</font></span><div class=""><div class="h5"><br>
On Mon, Sep 15, 2014 at 4:23 PM, Chris Marsh <<a href="mailto:chris.marsh@usask.ca">chris.marsh@usask.ca</a>> wrote:<br>
> Thanks David, this looks perfect.<br>
><br>
><br>
> On Mon, Sep 15, 2014 at 4:10 PM, David Gobbi <<a href="mailto:david.gobbi@gmail.com">david.gobbi@gmail.com</a>> wrote:<br>
>><br>
>> Hi Chris,<br>
>><br>
>> This sounds like a job for vtkProbeFilter.  The probe filter takes two<br>
>> inputs,<br>
>> the "input" and the "source".  The idea is that it resamples the scalars<br>
>> from<br>
>> the "source" onto the geometry of the "input".  So, in this case, the<br>
>> "source"<br>
>> is your unstructured grid, and the "input" is any vtkImageData that has<br>
>> the<br>
>> origin, spacing, and extent that you want to use for the regridding.  The<br>
>> output of vtkProbeFilter will be an image that has the same geometry as<br>
>> the input, but that has pixel values that have been interpolated from your<br>
>> unstructured grid.  You should be able to find several examples on the<br>
>> wiki.<br>
>><br>
>>  - David<br>
>><br>
>><br>
>> On Mon, Sep 15, 2014 at 3:59 PM, chrism <<a href="mailto:chris.marsh@usask.ca">chris.marsh@usask.ca</a>> wrote:<br>
>> > Is it possible to (easily) convert an unstructured grid to a structured<br>
>> > grid<br>
>> > (Image data)?<br>
>> ><br>
>> > Is unstructured -> polydata -> image data  the only way? If so, how<br>
>> > example<br>
>> > (<a href="http://www.vtk.org/Wiki/VTK/Examples/Cxx/PolyData/DataSetSurfaceFilter" target="_blank">http://www.vtk.org/Wiki/VTK/Examples/Cxx/PolyData/DataSetSurfaceFilter</a><br>
>> > ??)<br>
>> ><br>
>> > Ultimately, the unstructured grid represents a topography, and I am<br>
>> > looking<br>
>> > to 'rasterize' it for use in a GIS.<br>
><br>
><br>
</div></div></blockquote></div><br></div></div>