[求助]请问怎么计算弧线对应的圆心角？？？

361553894

怎么获取objType为"arc" 的弧线的圆心角？？？
目前做法是
Angle = axlDistance(car(seg->startEnd), cadr(seg->startEnd))/axlDistance(car(seg->startEnd), seg->xy)*180.0/3.1415926
但是怎么分辨小于半圆还是大于半圆？

361553894

deargds 发表于 2016-4-27 11:03
试下这个函数，可以区分优弧

1. procedure( wlGetSegAngle(@optional (seg nil))
   
2. prog( (angle angle1 angle2)
   
3. unless(seg
   
4. axlSetFindFilter(?enabled '("noall" "linesegs" "clinesegs") ?onButtons '("noall" "linesegs" "clinesegs"))
   
5. seg = car(axlGetSelSet(axlSingleSelectPoint()))
   
6. )
   
7. unless(seg return(nil))
   
8. when(seg->objType == "line"
   
9. angle = wlGetVectorAngle(car(seg->startEnd) cadr(seg->startEnd))
   
10. )
    
11. when(seg->objType == "arc"
    
12. angle1 = wlGetVectorAngle(seg->xy car(seg->startEnd))
    
13. angle2 = wlGetVectorAngle(seg->xy cadr(seg->startEnd))
    
14. angle = angle2 - angle1
    
15. if(seg->isClockwise then
    
16. when(angle > 0 angle = angle - 360)
    
17. else
    
18. when(angle < 0 angle = angle + 360)
    
19. )
    
20. )
    
21. return(angle)
    
22. ))
    
23. 
    
24. procedure( wlGetVectorAngle(start end)
    
25. let( (slope angle pi)
    
26. pi = 3.1415926
    
27. slope = axlLineSlope(list(start end))
    
28. if(slope then
    
29. angle = atan(slope) * 180.0 / pi
    
30. when(car(end) < car(start)
    
31. angle = angle + 180.0
    
32. )
    
33. when(angle < 0 angle = angle + 360.0)
    
34. else
    
35. if(cadr(end) > cadr(start) then
    
36. angle = 90.0
    
37. else
    
38. angle = 270.0
    
39. )
    
40. )
    
41. angle
    
42. ))

复制代码

这是我写的获取角度的函数，直线直接获取起点到终点的角度，顺时针弧线获取负的旋转角度，逆时针弧线获取正的旋转角度

deargds

试下这个函数，可以区分优弧

1. defun( GetARC_Angle (db)
   
2. prog((r xy start startX startY endX endY end len pi)
   
3.         len = nil
   
4.         pi = 3.1415926
   
5.         when(db->objType == "arc"
   
6.                 r = db->radius
   
7.                 xy = db->xy
   
8.                 start = nth(0 db->startEnd)
   
9.                 end = nth(1 db->startEnd)
   
10.        
    
11.                 startX = xCoord(start)-xCoord(xy)
    
12.                 startY = yCoord(start)-yCoord(xy)
    
13.                 endX = xCoord(end)-xCoord(xy)
    
14.                 endY = yCoord(end)-yCoord(xy)
    
15.                        
    
16.                 startR = GetPointRadian(startX startY)
    
17.                 endR = GetPointRadian(endX endY)
    
18. 
    
19.                 if(db->isClockwise then
    
20.                         if(startR>endR then
    
21.                                 len = (startR-endR)*r
    
22.                         else
    
23.                                 len = (2*3.1415926+startR-endR)*r
    
24.                         )
    
25.                 else
    
26.                         if(startR>endR then
    
27.                                 len = (2*3.1415926+endR-startR)*r
    
28.                         else
    
29.                                 len = (endR-startR)*r
    
30.                         )
    
31.                 )
    
32.         )
    
33.         return((len*180)/(3.1415926*r))
    
34. ))
    
35. 
    
36. defun( GetPointRadian (x y)
    
37. prog((a pi)
    
38.         pi = 3.1415926
    
39.         cond(
    
40.                 (x==0
    
41.                         cond(
    
42.                                 (y==0 a=0)
    
43.                                 (y>0 a=pi/2.0)
    
44.                                 (t a=3.0*pi/2)
    
45.                         )
    
46.                 )
    
47. 
    
48.                 (x>0
    
49.                         cond(
    
50.                                 (y==0 a=0)
    
51.                                 (y>0 a=atan(y/x))
    
52.                                 (t a=2.0*pi+atan(y/x))
    
53.                         )
    
54.                 )
    
55. 
    
56.                 (t
    
57.                         cond(
    
58.                                 (y==0 a=pi)
    
59.                                 (y>0 a=pi+atan(y/x))
    
60.                                 (t a=pi+atan(y/x))
    
61.                         )
    
62.                 )
    
63.         )
    
64.         return(a)
    
65. )
    
66. )
    
67. ;deargds

复制代码