点在多边形内部

csharp

http://www.cnblogs.com/dwdxdy/p/3230647.html

1. 
   
2. 
   
3.         //  The function will return YES if the point x,y is inside the polygon, or
   
4.         //  NO if it is not.  If the point is exactly on the edge of the polygon,
   
5.         //  then the function may return YES or NO.
   
6.         bool IsPointInPolygon(Point3dCollection pts, Point3d pt)
   
7.         {
   
8.             int i, j;
   
9.             bool c = false;
   
10.             for (i = 0, j = pts.Count - 1; i < pts.Count; j = i++)
    
11.             {
    
12.                 if ((((pts[i].Y <= pt.Y) && (pt.Y < pts[j].Y)) ||
    
13.                     ((pts[j].Y <= pt.Y) && (pt.Y < pts[i].Y)))
    
14.                     && (pt.X < (pts[j].X - pts[i].X) * (pt.Y - pts[i].Y) / (pts[j].Y - pts[i].Y) + pts[i].X))
    
15.                 {
    
16.                     c = !c;
    
17.                 }
    
18.             }
    
19.             return c;
    
20.         }
    

复制代码

csharp

bool IsInsideCurve(Curve curve, Point3d testPt)
        {
            if (!curve.Closed)
            {
                // cannot be inside
                return false;
            }

            Point3d ptOnCurve = curve.GetClosestPointTo(testPt, false);

            if (testPt == ptOnCurve)
            {
                return true;
            }

            if (!curve.IsPlanar)
            {
                return false;
            }

            // check its planar
            Plane plane = curve.GetPlane();

            // make the test ray from the plane
            double epsilon = 2e-6; // ( trust me on this )

            Point3dCollection IntersectionPoints = new Point3dCollection();
            Vector3d normal = plane.Normal;
            Vector3d testVector = normal.GetPerpendicularVector();

            using (Ray ray = new Ray())
            {
                ray.BasePoint = testPt;

                int nGlancingHits = 0, numberOfInters = 0;
                bool bRetryWithOtherRayDirection;

                do
                {
                    bRetryWithOtherRayDirection = false;

                    IntersectionPoints.Clear();

                    ray.UnitDir = testVector;

                    // fire the ray at the curve
                    curve.IntersectWith(ray,
                        Intersect.OnBothOperands,
                        IntersectionPoints,
                        IntPtr.Zero,
                        IntPtr.Zero);

                    numberOfInters = IntersectionPoints.Count;

                    if (numberOfInters == 0)
                    {
                        return false;
                    }

                    nGlancingHits = 0;

                    for (int i = 0; i < IntersectionPoints.Count; ++i)
                    {
                        double hitParam;

                        try
                        {
                            //This try/catch block circumvents an issue with GetParameterAtPoint API
                            hitParam = curve.GetParameterAtPoint(IntersectionPoints[i]);
                        }
                        catch
                        {
                            bRetryWithOtherRayDirection = true;
                            testVector = testVector.RotateBy(5.0 * Math.PI / 180.0, normal);
                            break;
                        }

                        double inParam = hitParam - epsilon;
                        double outParam = hitParam + epsilon;

                        //Loop back inside the curve if param is falling outside of range
                        if (inParam < curve.StartParam)
                            inParam = curve.EndParam - epsilon + (curve.StartParam - inParam);

                        if (outParam > curve.EndParam)
                            outParam = curve.StartParam + epsilon + (curve.EndParam - outParam);

                        IncidenceType inIncidence = CurveIncidence(curve, inParam, testVector, normal);
                        IncidenceType outIncidence = CurveIncidence(curve, outParam, testVector, normal);

                        if (inIncidence == IncidenceType.kIncidenceToFront || outIncidence == IncidenceType.kIncidenceToFront)
                        {
                            bRetryWithOtherRayDirection = true;
                            testVector = testVector.RotateBy(5.0 * Math.PI / 180.0, normal);
                            break;
                        }

                        if ((inIncidence == IncidenceType.kIncidenceToRight && outIncidence == IncidenceType.kIncidenceToLeft) ||
                            (inIncidence == IncidenceType.kIncidenceToLeft && outIncidence == IncidenceType.kIncidenceToRight))
                        {
                            nGlancingHits++;
                        }
                    }
                }
                while (bRetryWithOtherRayDirection);

                return ((numberOfInters + nGlancingHits) % 2

csharp

bool IsInsideCurve(Curve curve, Point3d testPt)
        {
            if (!curve.Closed)
            {
                // cannot be inside
                return false;
            }

            Point3d ptOnCurve = curve.GetClosestPointTo(testPt, false);

            if (testPt == ptOnCurve)
            {
                return true;
            }

            if (!curve.IsPlanar)
            {
                return false;
            }

            // check its planar
            Plane plane = curve.GetPlane();

            // make the test ray from the plane
            double epsilon = 2e-6; // ( trust me on this )

            Point3dCollection IntersectionPoints = new Point3dCollection();
            Vector3d normal = plane.Normal;
            Vector3d testVector = normal.GetPerpendicularVector();

            using (Ray ray = new Ray())
            {
                ray.BasePoint = testPt;

                int nGlancingHits = 0, numberOfInters = 0;
                bool bRetryWithOtherRayDirection;

                do
                {
                    bRetryWithOtherRayDirection = false;

                    IntersectionPoints.Clear();

                    ray.UnitDir = testVector;

                    // fire the ray at the curve
                    curve.IntersectWith(ray,
                        Intersect.OnBothOperands,
                        IntersectionPoints,
                        IntPtr.Zero,
                        IntPtr.Zero);

                    numberOfInters = IntersectionPoints.Count;

                    if (numberOfInters == 0)
                    {
                        return false;
                    }

                    nGlancingHits = 0;

                    for (int i = 0; i < IntersectionPoints.Count; ++i)
                    {
                        double hitParam;

                        try
                        {
                            //This try/catch block circumvents an issue with GetParameterAtPoint API
                            hitParam = curve.GetParameterAtPoint(IntersectionPoints[i]);
                        }
                        catch
                        {
                            bRetryWithOtherRayDirection = true;
                            testVector = testVector.RotateBy(5.0 * Math.PI / 180.0, normal);
                            break;
                        }

                        double inParam = hitParam - epsilon;
                        double outParam = hitParam + epsilon;

                        //Loop back inside the curve if param is falling outside of range
                        if (inParam < curve.StartParam)
                            inParam = curve.EndParam - epsilon + (curve.StartParam - inParam);

                        if (outParam > curve.EndParam)
                            outParam = curve.StartParam + epsilon + (curve.EndParam - outParam);

                        IncidenceType inIncidence = CurveIncidence(curve, inParam, testVector, normal);
                        IncidenceType outIncidence = CurveIncidence(curve, outParam, testVector, normal);

                        if (inIncidence == IncidenceType.kIncidenceToFront || outIncidence == IncidenceType.kIncidenceToFront)
                        {
                            bRetryWithOtherRayDirection = true;
                            testVector = testVector.RotateBy(5.0 * Math.PI / 180.0, normal);
                            break;
                        }

                        if ((inIncidence == IncidenceType.kIncidenceToRight && outIncidence == IncidenceType.kIncidenceToLeft) ||
                            (inIncidence == IncidenceType.kIncidenceToLeft && outIncidence == IncidenceType.kIncidenceToRight))
                        {
                            nGlancingHits++;
                        }
                    }
                }
                while (bRetryWithOtherRayDirection);

                return ((numberOfInters + nGlancingHits) % 2.  

csharp

bool IsInsideCurve(Curve curve, Point3d testPt)
        {
            if (!curve.Closed)
            {
                // cannot be inside
                return false;
            }

            Point3d ptOnCurve = curve.GetClosestPointTo(testPt, false);

            if (testPt == ptOnCurve)
            {
                return true;
            }

            if (!curve.IsPlanar)
            {
                return false;
            }

            // check its planar
            Plane plane = curve.GetPlane();

            // make the test ray from the plane
            double epsilon = 2e-6; // ( trust me on this )

            Point3dCollection IntersectionPoints = new Point3dCollection();
            Vector3d normal = plane.Normal;
            Vector3d testVector = normal.GetPerpendicularVector();

            using (Ray ray = new Ray())
            {
                ray.BasePoint = testPt;

                int nGlancingHits = 0, numberOfInters = 0;
                bool bRetryWithOtherRayDirection;

                do
                {
                    bRetryWithOtherRayDirection = false;

                    IntersectionPoints.Clear();

                    ray.UnitDir = testVector;

                    // fire the ray at the curve
                    curve.IntersectWith(ray,
                        Intersect.OnBothOperands,
                        IntersectionPoints,
                        IntPtr.Zero,
                        IntPtr.Zero);

                    numberOfInters = IntersectionPoints.Count;

                    if (numberOfInters == 0)
                    {
                        return false;
                    }

                    nGlancingHits = 0;

                    for (int i = 0; i < IntersectionPoints.Count; ++i)
                    {
                        double hitParam;

                        try
                        {
                            //This try/catch block circumvents an issue with GetParameterAtPoint API
                            hitParam = curve.GetParameterAtPoint(IntersectionPoints[i]);
                        }
                        catch
                        {
                            bRetryWithOtherRayDirection = true;
                            testVector = testVector.RotateBy(5.0 * Math.PI / 180.0, normal);
                            break;
                        }

                        double inParam = hitParam - epsilon;
                        double outParam = hitParam + epsilon;

                        //Loop back inside the curve if param is falling outside of range
                        if (inParam < curve.StartParam)
                            inParam = curve.EndParam - epsilon + (curve.StartParam - inParam);

                        if (outParam > curve.EndParam)
                            outParam = curve.StartParam + epsilon + (curve.EndParam - outParam);

                        IncidenceType inIncidence = CurveIncidence(curve, inParam, testVector, normal);
                        IncidenceType outIncidence = CurveIncidence(curve, outParam, testVector, normal);

                        if (inIncidence == IncidenceType.kIncidenceToFront || outIncidence == IncidenceType.kIncidenceToFront)
                        {
                            bRetryWithOtherRayDirection = true;
                            testVector = testVector.RotateBy(5.0 * Math.PI / 180.0, normal);
                            break;
                        }

                        if ((inIncidence == IncidenceType.kIncidenceToRight && outIncidence == IncidenceType.kIncidenceToLeft) ||
                            (inIncidence == IncidenceType.kIncidenceToLeft && outIncidence == IncidenceType.kIncidenceToRight))
                        {
                            nGlancingHits++;
                        }
                    }
                }
                while (bRetryWithOtherRayDirection);

                return ((numberOfInters + nGlancingHits) % 2  == 1
}
}