.NET

[#1]

Hello. A query, any of you know any method or function that tells me the orientation of a polyline, if it is clockwise or counter-clockwise

Thanks

[#2]

Hi,

Search google for Signed Area, depending on convention, if the signed area is positive then you have a ccw, if not cw, and if is 0 then is a degenerated polyline (just a line).

Gaston Nunez

[#3]

Try Clockwise function written by _gile from this page:

http://www.acadnetwork.com/topic-250.msg447.html#msg447

[#4 · Accepted solution]

Many thanks:

Using the formula for finding areas using determinants i was able to develop a simple function.

Greetings and thanks

Private Function SentidoPoli(ByVal pto3DPoli As Point3dCollection) As Boolean
        Dim r As Integer = pto3DPoli.Count
        Dim Matris(r, 2) As Double
        Dim rrIz As ArrayList = New ArrayList
        Dim rrDe As ArrayList = New ArrayList
        Dim sumaIz As Double = 0
        Dim sumaDe As Double = 0
        Dim area As Double = 0
        For i As Integer = 0 To r - 1
            Matris(i, 0) = pto3DPoli.Item(i).X
            Matris(i, 1) = pto3DPoli.Item(i).Y
            Matris(i, 2) = pto3DPoli.Item(i).Z
        Next
        For i As Integer = 0 To r - 2
            rrIz.Add(Matris(i, 0) * Matris(i + 1, 1))
            rrDe.Add(Matris(i + 1, 0) * Matris(i, 1))
        Next
        For i As Integer = 0 To rrIz.Count - 1
            sumaIz += rrIz.Item(i)
            sumaDe += rrDe.Item(i)
        Next
        area = 0.5 * ((sumaIz) - (sumaDe))
        If area > 0 Then
            Return True
        Else
            Return False
        End If
    End Function

[#5]

How about comparing the start and end points for only the first 2 line segments?

[#6]

Are you sure that's a good idea?

[#7]

Hmmm, you made me doubt my suggestion. Let's see.

I'm writing this whiling analizing the cases, so I won't know my final conclusion until the end of the post.

The succession of polyline segments are created in 1 of 4 cases. I'm attaching the image describing those 4 cases.

Note: case 4 is the most comon case.

In all cases the point for the end of segment 1 is the same as the start of segment 2, so comparing them won't give us any information. That is why I suggested comparing the start point of segment 1 and the end point of segment 2.

Let's analize the cases:

Let's call the coordinates Seg1Start.X, Seg1Start.Y, Seg2Start.X and Seg2Start.Y.

Case 1:

• The Seg1Start.X is higher than Seg2End.X and Seg1Start.Y lower than Seg2End.Y, that means Seg1Start is to the right and below of Seg2End. Seg2 can be as long as you want. The angle can go from 0 to 180 degrees. If the angle is 0 and if Seg2End and Seg1Start are the same, the result of comparing them would be ambiguous, something that would have to be controlled. If the angle is 180. Seg2End would be left of Seg1Start, making the polyline clockwise. Anything out of that angle range would make Case1 turn into Case3. Case 1 is clockwise.

Case 2:

• Anything out of angle range turns it into Case4. Overlapping of Seg1Start and Seg2End points must be controlled for ambiguity. Case2 is counter-clockwise.

Case 3:

• Anything out of angle range turns it into Case1. Overlapping of Seg1Start and Seg2End points must be controlled for ambiguity. Case3 is counter-clockwise.

Case 4:

• Anything out of angle range turns it into Case2. Overlapping of Seg1Start and Seg2End points must be controlled for ambiguity. Case4 is clockwise.

I might be missing something (validation for weird cases). It's just a quick look at the problem.

Can you take a look at my post about hatches?

I always create threads and no one replies with ideas or advise @.@My hatch post

[#8]

First two segment are idetical...

[#9]

Well, we'd have to establish the criteria to clasify a polygon as clockwise or counter-clockwise.

I thought it depended on the order of the vertices of the polygon.

It's a matter that requires a bit of reading time I don't have right now.

You are right, that's another case. I think I was just considering quadrilaterals XDDD

I found this: http://stackoverflow.com/questions/1165647/how-to-determine-if-a-list-of-polygon-points-are-in-clock...

and this: http://debian.fmi.uni-sofia.bg/~sergei/cgsr/docs/clockwise.htm

[#10 · Accepted solution]

Hi,

As said gasty1001, you can use the algebraic area (signed area) with non self crossed polylines.

Here's a little static class defining extension methods (works with arc segment polyline too).

    public static class AlgebraicArea
    {
        public static double GetArea(Point2d pt1, Point2d pt2, Point2d pt3)
        {
            return (((pt2.X - pt1.X) * (pt3.Y - pt1.Y)) -
                        ((pt3.X - pt1.X) * (pt2.Y - pt1.Y))) / 2.0;
        }

        public static double GetArea(this CircularArc2d arc)
        {
            double rad = arc.Radius;
            double ang = arc.IsClockWise ?
                arc.StartAngle - arc.EndAngle :
                arc.EndAngle - arc.StartAngle;
            return rad * rad * (ang - Math.Sin(ang)) / 2.0;
        }

        public static double GetArea(this Polyline pline)
        {
            CircularArc2d arc = new CircularArc2d();
            double area = 0.0;
            int last = pline.NumberOfVertices - 1;
            Point2d p0 = pline.GetPoint2dAt(0);

            if (pline.GetBulgeAt(0) != 0.0)
            {
                area += pline.GetArcSegment2dAt(0).GetArea();
            }
            for (int i = 1; i < last; i++)
            {
                area += GetArea(p0, pline.GetPoint2dAt(i), pline.GetPoint2dAt(i + 1));
                if (pline.GetBulgeAt(i) != 0.0)
                {
                    area += pline.GetArcSegment2dAt(i).GetArea(); ;
                }
            }
            if ((pline.GetBulgeAt(last) != 0.0) && pline.Closed)
            {
                area += pline.GetArcSegment2dAt(last).GetArea();
            }
            return area;
        }
    }

Using example:

        [CommandMethod("Test")]
        public void ClockwisePolyline()
        {
            Document doc = AcAp.DocumentManager.MdiActiveDocument;
            Database db = HostApplicationServices.WorkingDatabase;
            Editor ed = doc.Editor;
            using (Transaction tr = db.TransactionManager.StartTransaction())
            {
                BlockTableRecord btr =
                    (BlockTableRecord)db.CurrentSpaceId.GetObject(OpenMode.ForWrite);
                PromptEntityOptions peo = new PromptEntityOptions("\nSelect a polyline: ");
                peo.SetRejectMessage("Selected object is not a polyline.");
                peo.AddAllowedClass(typeof(Polyline), false);
                while (true)
                {
                    PromptEntityResult per = ed.GetEntity(peo);
                    if (per.Status != PromptStatus.OK) break; ;
                    Polyline pline = (Polyline)per.ObjectId.GetObject(OpenMode.ForRead);
                    double area = pline.GetArea();
                    Application.ShowAlertDialog(
                        string.Format("{0}\nArea = {1}", area < 0 ? "CW" : "CCW", area));
                }
                tr.Commit();
            }
        }