<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Hi Tim, the input dataset to the probe filter must be constructed from
a space filling cell type (like a VTK_VOXEL). The output from the
contour filter is most likely composed of either polygons or lines
depending on the dimensionality of the input data. These cell types are
not space filling. The space filling cell requirement has to do with
the interpolation algorithm. if the probe doesn't lie inside the data
set then technically speaking interpolation is not possible, in that
case an extrapolation algorithm would be required, but the probe filter
isn't equipped with that, and that my best guess as to why your code
doesn't work.<br>
<br>
see this post
<a class="moz-txt-link-freetext" href="http://public.kitware.com/pipermail/vtkusers/2007-September/092610.html">http://public.kitware.com/pipermail/vtkusers/2007-September/092610.html</a><br>
<br>
<pre class="moz-signature" cols="72">--
Burlen Loring
Information Technologist III
Space Science Center
Institute for the Study of Earth, Oceans, and Space
University of New Hampshire
39 College Road, Durham, NH 03824
Phone: 603-862-1140</pre>
<br>
<br>
Tim Seher wrote:
<blockquote cite="mid200710121338.55671.seher@ipgp.jussieu.fr"
type="cite">
<pre wrap="">Hi!
Sorry, for this second post. Just want to give some additional information.
I have been playing around some more and I still cannot see, why probing a
polydata with a line does not work. All I want to do is to sample a surface
at a given point xy. Any help highly appreciated! Here is, how far I got:
</pre>
<blockquote type="cite">
<pre wrap=""> ### data
reader = vtk.vtkStructuredPointsReader()
reader.SetFileName( infile )
reader.Update()
dims = reader.GetOutput().GetBounds()
### calculate the surfaces
surfaces = vtk.vtkContourFilter()
surfaces.SetInput(reader.GetOutput())
surfaces.GenerateValues(4,2,5)
line = vtk.vtkLineSource()
transform = vtk.vtkTransform()
transform.Translate([0,0,-2])
transform.Scale([100,100,100])
transform.RotateY(90)
tf = vtk.vtkTransformPolyDataFilter()
tf.SetInputConnection(line.GetOutputPort())
tf.SetTransform(transform)
probe = vtk.vtkProbeFilter()
probe.SetInputConnection(tf.GetOutputPort())
probe.SetSource(surfaces.GetOutput())
</pre>
</blockquote>
<pre wrap=""><!---->
</pre>
<blockquote type="cite">
<pre wrap=""> xyplot = vtk.vtkXYPlotActor()
xyplot.AddInput(probe.GetOutput())
ren = vtk.vtkRenderer()
ren.AddActor(xyplot)
ren.SetBackground(0.7, 0.7, 0.7)
win = vtk.vtkRenderWindow()
win.AddRenderer(ren)
win.SetSize(800, 600)
iren = vtk.vtkRenderWindowInteractor()
iren.SetRenderWindow(win)
win.Render()
iren.Initialize()
iren.Start()
</pre>
</blockquote>
<pre wrap=""><!---->
Sorry and thanks a lot! Tim
_______________________________________________
This is the private VTK discussion list.
Please keep messages on-topic. Check the FAQ at: <a class="moz-txt-link-freetext" href="http://www.vtk.org/Wiki/VTK_FAQ">http://www.vtk.org/Wiki/VTK_FAQ</a>
Follow this link to subscribe/unsubscribe:
<a class="moz-txt-link-freetext" href="http://www.vtk.org/mailman/listinfo/vtkusers">http://www.vtk.org/mailman/listinfo/vtkusers</a>
</pre>
</blockquote>
<br>
<br>
<pre class="moz-signature" cols="72">
</pre>
</body>
</html>