[vtkusers] Intersection of line with 2D PolyData set

Sergei Gasilov sergei.gasilov at gmail.com
Mon Jun 21 08:01:20 EDT 2010

Dear VTK Users,

Is it possible to find all intersetion points of line (ray) with
PolyData set consisting of points and line segments?
My PolyData  set named "ells" is defined by points and lines and it is
couple of embedded ellipses (http://tinypic.com/r/2wdyjqq/6).
Everything is in the z=0 plane.
I have tried to use both vtkOBBtree and vtkModifiedBSPtree, according
to the examples from vtk wiki:

vtkSmartPointer<vtkOBBTree> tree = vtkSmartPointer<vtkOBBTree>::New();
		tree->SetDataSet(ells); //ells is of vtkSmartPointer<vtkPolyData>

The last line leads to the error:  vtkMath::Jacobi: Error extracting
eigenfunctions, and "IntersectWithLine" methos finds no intersection

I noticed also on vtk wiki, that vtkOBBtree is designed to find
"intersection of line with triangles", indeed it works fine with 3D
data set.  Does it mean it is not working with line/polyline cell
types or I'm simply wrong somewhere? Is there solution in vtk for such

With kind regards,

Sergei V. Gasilov

More information about the vtkusers mailing list