无聊，出来求助下，各位高手请解答下在下的问题吧

邹锋

在下在写一个LISP程序时遇到个问题，就是在多线段中，已知圆弧两点坐标与凸度值，求出这个圆弧的半径与圆心坐标，要用LISP代码写，求代码，，



以下是我在网络上搜索的用C语言写 的

1. （1）AutoCAD中约定：凸度为0是直线顶点，它与下一个顶点连接为一直线；凸度不为0是圆弧顶点，它与下一个顶点连接为一圆弧；凸度值为负表示顺时针圆弧，凸度值为正表示逆时针圆弧；凸度绝对值小于1表示圆弧包角小于180°，凸度绝对值大于1表示圆弧包角大于180°。凸度与圆弧包角的关系是：圆弧包角= 4×arctan|凸度值|。 void lwpolylineToArc(CPoint3d BeginPoint,CPoint3d EndPoint,double u,CPoint3d &CenterPoint)
   
2.   {
   
3.         double centerAngle;//包角
   
4.         centerAngle=4*atan(abs(u));
   
5.         centerAngel=centerAngel/pi;
   
6. 
   
7.         double x1,x2,y1,y2;//圆弧起始点和终止点
   
8.        x1=BeginPoint.x;
   
9.       x2=EndPoint.x;
   
10.       y1=BeginPoint.y;
    
11.       y2=EndPoint.y;
    
12. 
    
13.      double L; //弦长
    
14.      L=sqrt(pow((x1-x2),2)+pow((y1-y2),2));
    
15.   
    
16.         double R;//圆弧半径
    
17.         R=0.5*L/sin(0.5*centerAngle);
    
18. 
    
19.         //已知圆上两点和半径，求圆心坐标
    
20.         double h;//圆心到弦的距离
    
21.         h=sqrt(R*R-L*L/4);
    
22. 
    
23.         double k;//起始点和终止点连线的中垂线斜率
    
24.         double xc,yc;//圆心坐标
    
25.         double xa,ya; //起始点和终止点连线的中点横纵坐标
    
26.         xa=0.5*(x1+x2);
    
27.         ya=0.5*(y1+y2);
    
28. 
    
29.         //弦的方向角（0-2PI之）
    
30. 
    
31.          double angle;//起点到终点的弦向量与x正方向之间的倾斜角
    
32.          angle=acos((x2-x1)/sqrt(pow(x2-x1,2)+pow(y2-y1,2)));
    
33. 
    
34.          double amass; //弦向量与X轴正向单位向量的叉积
    
35.          amass = y1-y2;//由（由(x2-x1)*0-1*(y2-y1)）得到
    
36.          
    
37.          if (amass<0)
    
38.          {  angle=-angle;
    
39.           angle=2*PI+angle;
    
40.         }
    
41.       
    
42.    
    
43.         
    
44.       double DirectionAngel;//弦中点到圆心的直线向量的方向角（0-2PI之间）
    
45.       if ((u>0 && centerAngle<PI)||(u<0 && centerAngle>PI))
    
46.       DirectionAngel=angle+PI/2;
    
47.       if((u<0 && centerAngle<PI)||(u>0 && centerAngle>PI))
    
48.       DirectionAngel=angle-PI/2;
    
49.       if (DirectionAngel>2*PI)
    
50.        DirectionAngel= DirectionAngel-2*PI;
    
51. 
    
52.        double d;//圆心到弦的距离
    
53.       d=sqrt(R*R-L*L/4);
    
54.       if (DirectionAngle=0)
    
55.      {
    
56.       xc=xa+d;
    
57.       yc=ya;
    
58.     }
    
59.      else if(DirectionAngle=PI/2)
    
60.    {
    
61.      xc=xa;
    
62.      yc=ya+d;
    
63.    }
    
64.    else if (DirectionAngle=PI)
    
65.    {
    
66.     xc=xa-d;
    
67.     yc=xa;
    
68.    }
    
69.    else if (DirectionAngle=PI+PI/2)
    
70.    {
    
71.     xc=xa;
    
72.     yc=xa-d;
    
73.    }
    
74.    else
    
75.    {
    
76.     double nslope,k;//nslope 为弦的斜率，K为弦中垂线的斜率
    
77.     double nAngle;//中垂线的倾斜角；
    
78.     double X,Y; //圆心相对于弦中心点的坐标偏移量
    
79.    
    
80.     nslope = (y2 - y1) / (x2-x1);
    
81.     k = -1 / nslope;
    
82.        nAngle = atan(k) ;
    
83.        X = cos(nAngle) * d;
    
84.        Y = sin(nAngle) * d;
    
85. 
    
86.        if (DirectionAngle > PI / 2 && DirectionAngle < PI )
    
87.     {X = -X;
    
88.         Y = -Y;
    
89.     }
    
90.        if (DirectionAngle > PI && DirectionAngle < (PI + PI / 2) )
    
91.     {
    
92.           X = -X;
    
93.           Y = -Y;
    
94.     }
    
95. 
    
96.        xc=xa + X;
    
97.        yc=ya+ Y;
    
98. 
    
99.       CenterPoint.x=xc;
    
100.       CenterPoint.y=yc;
     
101.       CenterPoint.z=0.0;

复制代码

Andyhon

Ref:
http://bbs.mjtd.com/thread-89561-1-1.html
==>bulge

nzl1116

(setq Ang (* (atan tudu) 2.0)
        Dist  (* (distance Pnt1 Pnt2) 0.5)
        R     (/ Dist (sin Ang))
        D1   (vlax-curve-getDistAtPoint curve Pnt1)
        D2   (vlax-curve-getDistAtPoint curve Pnt2)
        D3   (* (+ D1 D2) 0.5)
        Pnt3 (vlax-curve-getPointAtDist curve D3)
        Mid  (polar Pnt1 (angle Pnt1 Pnt2) Dist)
        Cen (polar Pnt3 (angle Pnt3 Mid) R)
)

邹锋

nzl1116 发表于 2012-4-8 09:21
(setq Ang (* (atan tudu) 2.0)
        Dist  (* (distance Pnt1 Pnt2) 0.5)
        R     (/ Dist (si ...

(defun c:test()
  (setq pnt1(getpoint))
  (setq pnt2(getpoint))
  (setq tudu -0.414214);;画个4分之1圆
(setq Ang (* (atan tudu) 2.0)
      Dist  (* (distance Pnt1 Pnt2) 0.5)
      R     (/ Dist (sin Ang))
      D1   (vlax-curve-getDistAtPoint curve Pnt1)
      D2   (vlax-curve-getDistAtPoint curve Pnt2)
      D3   (* (+ D1 D2) 0.5)
      Pnt3 (vlax-curve-getPointAtDist curve D3)
      Mid  (polar Pnt1 (angle Pnt1 Pnt2) Dist)
      Cen (polar Pnt3 (angle Pnt3 Mid) R)
)
  (prin1 cen))

无效！！！

Gu_xl

1. ;;参数 pt = 圆弧起点 np = 圆弧终点 gx = 弓弦比
   
2. ;;返回表 '(圆心 半径)
   
3. (defun get_Center_radius (pt np gx / pa bj xc gg rr cp mdp)
   
4.     (setq pa (angle pt np)
   
5.           bj (* (atan (abs gx)) 4)
   
6.           xc (* 0.5 (distance pt np))
   
7.            gg (* gx xc)
   
8.           )
   
9.     (if (/= 0.0 gg)
   
10.       (progn
    
11.         (setq  rr (/ (+ (* xc xc) (* gg gg)) (* 2 gg)))
    
12.           
    
13.     (setq cp (polar pt pa xc)
    
14.           mdp cp
    
15.           cp (polar cp (+ pa pi2) (- rr gg))
    
16.           )
    
17.     (list cp rr)
    
18.         )
    
19.       (list pt 0.0)
    
20.       )
    
21.     )
    

;;参数 pt = 圆弧起点 np = 圆弧终点 gx = 弓弦比 <br /> ;;返回表 '(圆心 半径)<br /> (defun get_Center_radius (pt np gx / pa bj xc gg rr cp mdp)<br />     (setq pa (angle pt np)<br />           bj (* (atan (abs gx)) 4)<br />           xc (* 0.5 (distance pt np))<br />            gg (* gx xc)<br />           )<br />     (if (/= 0.0 gg)<br />       (progn<br />         (setq  rr (/ (+ (* xc xc) (* gg gg)) (* 2 gg)))<br />           <br />     (setq cp (polar pt pa xc)<br />           mdp cp<br />           cp (polar cp (+ pa pi2) (- rr gg))<br />           )<br />     (list cp rr)<br />         )<br />       (list pt 0.0)<br />       )<br />     )<br />

nzl1116

邹锋 发表于 2012-4-8 20:42
(defun c:test()
  (setq pnt1(getpoint))
  (setq pnt2(getpoint))

其中的变量curve指的是多段线

邹锋

nzl1116 发表于 2012-4-9 08:28
其中的变量curve指的是多段线

谢谢了，我昨晚想到了，还要选择多线段再转对象类型就是curve了

Andyhon

Ref:
http://www.google.com/search?as_q=&as_epq=revpoly.Lsp