Как определить направление вращения замкнутой полилинии?
Чтобы определить направление вращения полилинии по часовой или против часовой стрелки можно пройтись по всем сегментам полилинии и проверить лежит ли он справа или слева. Следующий код делает именно это. Он подсчитывает сумму числа поворотов для всех сегментов. Повороты возможны внутри сегментов (если они дуговые) или между последовательными сегментами. Результат работы функции – число поворотов. Если поворот налево – число положительное, если поворот направо – отрицательное. Заметим, что алгоритм не определяет, является ли полилиния самопересекающейся и может даже давать ошибочные результаты если полилиния действительно самопересекающаяся.
- Adesk::Boolean polyWindingNumber( AcDbPolyline*pPoly, int&winding )
- {
- if( !pPoly->isClosed() )
- {
- return Adesk::kFalse;
- }
- double turn = 0;
- int nSegs = pPoly->numVerts();
- AcGeVector3dArray startTans(nSegs);
- AcGeVector3dArray endTans(nSegs);
- double TWOPI = 6.28318530718;
- double PIBYTWO = 1.570796326795;
- for( int i=0;i<nSegs; i++ )
- {
- if( pPoly->segType(i) == AcDbPolyline::kArc )
- {
- AcGeCircArc2d arc;
- pPoly->getArcSegAt(i, arc );
- AcGeVector2d startTan;
- AcGeVector2d endTan;
- AcGeVector2dArray startDerivs;
- arc.evalPoint( arc.startAng(), 1, startDerivs );
- startTan = startDerivs[0];
- AcGeVector2dArray endDerivs;
- arc.evalPoint( arc.endAng(), 1, endDerivs);
- endTan = endDerivs[0];
- startTans.append(AcGeVector3d(startTan.x, startTan.y,0.0));
- endTans.append(AcGeVector3d(endTan.x, endTan.y, 0.0));
- double ang = arc.endAng() - arc.startAng();
- turn += arc.isClockWise() ? -ang : ang;
- }
- else if(pPoly->segType(i) == AcDbPolyline::kLine )
- {
- AcGeLineSeg2d line;
- pPoly->getLineSegAt(i, line );
- AcGeVector2d tan2d = line.endPoint() -line.startPoint();
- AcGeVector3d tan = AcGeVector3d( tan2d.x, tan2d.y, 0.0);
- startTans.append(tan);
- endTans.append(tan);
- }
- }
- nSegs = startTans.length();
- for(int i=0;i<nSegs; i++ )
- {
- double angle = endTans[i].angleTo( startTans[(i+1)%nSegs] );
- AcGeVector3d norm = endTans[i].crossProduct(
- startTans[(i+1)%nSegs] );
- angle = norm.z > 0 ? angle: -angle;
- turn += angle;
- }
- turn = turn / TWOPI;
- double lower = floor( turn );
- double tol = 1e-6;
- if( (turn - lower ) < tol)
- winding = (int)lower;
- else if( (( lower + 1 ) - turn ) < tol )
- winding = (int)(lower + 1);
- else
- return Adesk::kFalse;
- return Adesk::kTrue;
- }
- static void AsdkUtilsPolyWind(void)
- {
- ads_name eName;
- ads_point pt;
- if( RTNORM !=
- acedEntSel(L"\nВыберите полилинию для проверки: ", eName, pt) )
- return;
- AcDbObjectId id;
- acdbGetObjectId( id, eName );
- AcDbPolyline*pPoly;
- acdbOpenObject(pPoly, id, AcDb::kForRead );
- if(pPoly == NULL )
- return;
- int winding =0;
- if( polyWindingNumber(pPoly, winding ) )
- {
- acutPrintf(L"\nЧисло поворотов полигона %d", winding );
- }
- else
- {
- acutPrintf(L"\nНевозможно вычислить число поворотов полилинии ");
- }
- pPoly->close();
- }
Источник: http://adndevblog.typepad.com/autocad/2013/01/how-to-test-the-winding-turn-of-a-closed-polyline.html
Обсуждение: http://adn-cis.org/forum/index.php?topic=387
Опубликовано 09.12.2013