It works for my python script:<br><br>from vtk import *<br>s = vtkSphereSource()<br>s.SetPhiResolution(40)<br>s.Update()<br>pd = s.GetOutput()<br>pd.GetNumberOfPoints()<br>loc = vtkPointLocator()<br>loc.SetDataSet(pd)<br>id = vtkIdList()<br>
pt = [0,0,0]<br>loc.FindPointsWithinRadius(5, pt, id)<br>id.GetNumberOfIds()<br>pd.GetNumberOfPoints()<br><br><br>What version of vtk are you using? If you send me a test code that fails with a data set I'll check it out.<br>
<br>Andy<br><br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
Date: Wed, 9 Jun 2010 23:50:08 +0300<br>
From: "A. Burak Yoldemir" <<a href="mailto:yoldemir@gmail.com">yoldemir@gmail.com</a>><br>
Subject: Re: [vtkusers] vtkPointLocator Problem<br>
To: Bill Lorensen <<a href="mailto:bill.lorensen@gmail.com">bill.lorensen@gmail.com</a>>, <a href="mailto:vtkusers@vtk.org">vtkusers@vtk.org</a><br>
Message-ID:<br>
<<a href="mailto:AANLkTikyCw5WIypHmpkjkZgnroGkfmtcij5R7_AaGFa4@mail.gmail.com">AANLkTikyCw5WIypHmpkjkZgnroGkfmtcij5R7_AaGFa4@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Yes, I am aware of that example, but I have no solution to my problem. What<br>
I want is to use a vtkPolyData object as the dataset of a vtkPointLocator.<br>
Here is my problem in a nutshell:<br>
<br>
- I have a number of points in a vtkPoints object.<br>
- Then, I create a vtkPolyData object and use the SetPoints function of it<br>
to define the point coordinates.<br>
- Next, I create a vtkPointLocator object and use the SetDataSet function of<br>
it to build the locator from these points.<br>
- Finally, I use FindPointsWithinRadius function for a point which has many<br>
neighbors (just manually checked), but no neighbors are returned.<br>
- I also tried FindClosestPoint function for each of my points. At each<br>
time, the closest point ID turns out to be 0, however it actually is not.<br>
<br>
What might I be missing?<br>
<br>
<br>
- Ahmet Burak Yoldemir<br>
<br>
<br>
On Wed, Jun 9, 2010 at 1:04 AM, Bill Lorensen <<a href="mailto:bill.lorensen@gmail.com">bill.lorensen@gmail.com</a>>wrote:<br>
<br>
> Have you looked at this example:<br>
> <a href="http://www.vtk.org/Wiki/VTK/Examples/PointLocator_Radius" target="_blank">http://www.vtk.org/Wiki/VTK/Examples/PointLocator_Radius</a><br>
><br>
><br>
> On Tue, Jun 8, 2010 at 4:49 PM, A. Burak Yoldemir <<a href="mailto:yoldemir@gmail.com">yoldemir@gmail.com</a>><br>
> wrote:<br>
> > Hello all,<br>
> ><br>
> > I am having a problem with FindPointsWithinRadius function of<br>
> > vtkPointLocator class. I have a bunch of points in a certain volume, with<br>
> no<br>
> > specific topology. Using the following code snippet, I create a point<br>
> > locator object and specify its data set:<br>
> ><br>
> > vtkPoints* edge_pts = vtkPoints::New();<br>
> ><br>
> > for(int i = 0; i < someLimit; i++)<br>
> ><br>
> > {<br>
> ><br>
> > edge_pts->InsertNextPoint(somePoint);<br>
> ><br>
> > }<br>
> ><br>
> > vtkPolyData* edge_pts_set = vtkPolyData::New();<br>
> ><br>
> > edge_pts_set->SetPoints(edge_pts);<br>
> ><br>
> > edge_pts_set->Update();<br>
> ><br>
> > vtkPointLocator* NeighborLocator = vtkPointLocator::New();<br>
> ><br>
> > NeighborLocator->SetDataSet(edge_pts_set);<br>
> ><br>
> > NeighborLocator->BuildLocator();<br>
> ><br>
> > Then, I use the FindPointsWithinRadius function of NeighborLocator, where<br>
> I<br>
> > (manually) made sure that there are several neighbors of the query point.<br>
> > However, the aforementioned function returns no neighbors. Any ideas on<br>
> > where the problem could be?<br>
> ><br>
> > Regards,<br>
> ><br>
> ><br>
> > Burak<br></blockquote></div><br>