hi,<div><br></div><div>Below is the sample code what i was trying.</div><div><br></div><div>int num = polyPd->GetNumberOfPoints(); </div><div>// polyPd is a polydata containing points that forms a polygon</div><div><br>
</div><div>vtkPoints* points = this->del->GetOutput()->GetPoints();</div><div>// del is an instance of vtkDelaunay2d class. Out of this object is used to display</div><div>// scattered points</div><div><br></div>
<div>vtkPoints *newPoints = vtkPoints::New();</div><div>// To store selected points</div><div><br></div><div>double bnds[6], n[3];</div><div><br></div><div>// Get the bounds in 'bnds'</div><div>polyPd->GetPoints()->GetBounds(bnds);</div>
<div><br></div><div>// To get normal direction in 'n'</div><div>vtkPolygon::ComputeNormal(num, </div><div> static_cast<double *>(polyPd->GetPoints()->GetData()->GetVoidPointer(0)), n);</div><div><br>
</div><div>for( int id = 0; id < points->GetNumberOfPoints(); id++ )</div><div>{</div><div> if( vtkPolygon::PointInPolygon( points->GetPoint(id), num, </div><div> static_cast<double *>(polyPd->GetPoints()->GetData()->GetVoidPointer(0)), n) == 1 )</div>
<div> newPoints->InsertNextPoint(points->GetPoint(id));</div><div>}</div><div><br></div><div>// code to create a mapper and actor with </div><div>// the help of newPoints</div><div><br></div><div>I hope there is no difficulty in understanding the other part of code. I got this logic from the mailing list when i searched for PoinInPolygon here</div>
<div><br></div><div><a href="http://public.kitware.com/pipermail/vtkusers/2010-April/107974.html">http://public.kitware.com/pipermail/vtkusers/2010-April/107974.html</a></div><div><br></div><div>Thanks</div><div><br></div>
<div>Warm Regards</div><div>Rakesh Patil</div><div><br></div><div><br><div class="gmail_quote">On Wed, Sep 29, 2010 at 3:49 PM, Jothy <span dir="ltr"><<a href="mailto:jothybasu@gmail.com">jothybasu@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div dir="ltr">Could you please post your code?<br><br>Jothy<br><br><div class="gmail_quote"><div><div></div><div class="h5">
On Wed, Sep 29, 2010 at 11:09 AM, rakesh patil <span dir="ltr"><<a href="mailto:prakeshofficial@gmail.com" target="_blank">prakeshofficial@gmail.com</a>></span> wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204, 204, 204);padding-left:1ex"><div><div></div><div class="h5">Hi,<div><br></div><div>I have been trying to find whether a given point is within a polygon or not. I have been doing this using,</div>
<div><br></div><div>vtkPolygon::PointInPolygon() </div><div><br></div><div>function. But out of 10 runs, in 1 run, i get only few points inside, the polygon. Others just fails and shows 0 points selected.</div>
<div><br></div><div>Can anyone come up with a small example, for this?</div><div><br></div><div>Thanks</div><div><br></div><div>Warm Regards</div><div>Rakesh Patil</div>
<br></div></div>_______________________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the VTK FAQ at: <a href="http://www.vtk.org/Wiki/VTK_FAQ" target="_blank">http://www.vtk.org/Wiki/VTK_FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.vtk.org/mailman/listinfo/vtkusers" target="_blank">http://www.vtk.org/mailman/listinfo/vtkusers</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div dir="ltr">PhD Student<br>Dept. of Medical Physics<br>Clatterbridge Centre for Oncology<br>UK<br></div><br>
</div>
</blockquote></div><br></div>