manolas
/
vcglib
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Added PointInsidePolygon function

This commit is contained in:
Nico Pietroni 2012-01-12 00:53:37 +00:00
parent bb7fbdcc88
commit 75aca76eb7
1 changed files with 40 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

40
vcg/space/intersection2.h
View File

@ -253,6 +253,46 @@ namespace vcg {
} }
return false; return false;
} }
template <class ScalarType>
bool PointInsidePolygon(vcg::Point2<ScalarType> p,
const std::vector<vcg::Segment2<ScalarType> > &polygon)
{
int n=polygon.size();
vcg::Box2<ScalarType> BB;
for (int i=0;i<n;i++)
{
BB.Add(polygon[i].P0());
BB.Add(polygon[i].P1());
}
if (!BB.IsIn(p))return false;
ScalarType size=BB.Diag();
///take 4 directions
int inside_test=0;
for (int dir=0;dir<4;dir++)
{
int intersection=0;
vcg::Ray2<ScalarType> r;
vcg::Point2<ScalarType> direct=vcg::Point2<ScalarType>(0,0);
switch (dir)
{
case 0 : direct.X()=1;break;
case 1 : direct.Y()=1;break;
case 2 : direct.X()=-1; break;
default :direct.Y()=-1;
}
r.SetOrigin(p);
r.SetDirection(direct);
for (int i=0;i<n;i++)
{
Point2<ScalarType> p_inters;
if (vcg::RaySegmentIntersection(r,polygon[i],p_inters))intersection++;
}
if ((intersection%2)==1)
inside_test++;
}
return(inside_test>2);
}
//intersection between a circle and a line //intersection between a circle and a line
template<class ScalarType> template<class ScalarType>