根据顶点和凸度求弧段圆心及半径

ll_j

1. ;;;根据顶点和凸度求弧段圆心及半径;;;pt1 pt2 顶点坐标，二维点或z坐标为0的三维点
   
2. ;;;bul凸度，按多段线规则，逆时针为正
   
3. ;;;返回由圆心二维点及半径构成的点表
   
4. ;;;例：(ptc '(0.0 0.0) '(1.0 1.0) -1.5) —>((0.291667 0.708333) 0.766032)
   
5. (defun ptc (pt1 pt2 bul / a b h r hh x y)
   
6.   (cond
   
7.     ((= bul 0.0)
   
8.      (princ "\n凸度为0，两点间为直线，无圆心。")
   
9.      (princ)
   
10.     )
    
11.     ((not (and (numberp (car pt1))
    
12.                (numberp (cadr pt1))
    
13.                (numberp (car pt2))
    
14.                (numberp (cadr pt2))
    
15.                (or (not (caddr pt1)) (= (caddr pt1) 0.0))
    
16.                (or (not (caddr pt2)) (= (caddr pt2) 0.0))
    
17.                (< (length pt1) 4)
    
18.                (< (length pt2) 4)
    
19.           )
    
20.      )
    
21.      (princ "\n非法顶点。")
    
22.      (princ)
    
23.     )
    
24.     (t
    
25.      (setq a  (angle pt1 pt2)
    
26.            b  (distance pt1 pt2)
    
27.            h  (- (* b bul 0.5))
    
28.            r  (/ (+ (* b b) (* 4 h h)) 8 h)
    
29.            hh (- r h)
    
30.            x  (+ (/ (+ (car pt1) (car pt2)) 2) (* hh (sin a)))
    
31.            y  (- (/ (+ (cadr pt1) (cadr pt2)) 2) (* hh (cos a)))
    
32.      )
    
33.      (list (list x y) r)
    
34.     )
    
35.   )
    
36. )

;;;根据顶点和凸度求弧段圆心及半径;;;pt1 pt2 顶点坐标，二维点或z坐标为0的三维点<br /> ;;;bul凸度，按多段线规则，逆时针为正<br /> ;;;返回由圆心二维点及半径构成的点表<br /> ;;;例：(ptc '(0.0 0.0) '(1.0 1.0) -1.5) —>((0.291667 0.708333) 0.766032)<br /> (defun ptc (pt1 pt2 bul / a b h r hh x y)<br />   (cond<br />     ((= bul 0.0)<br />      (princ "\n凸度为0，两点间为直线，无圆心。")<br />      (princ)<br />     )<br />     ((not (and (numberp (car pt1))<br />                (numberp (cadr pt1))<br />                (numberp (car pt2))<br />                (numberp (cadr pt2))<br />                (or (not (caddr pt1)) (= (caddr pt1) 0.0))<br />                (or (not (caddr pt2)) (= (caddr pt2) 0.0))<br />                (< (length pt1) 4)<br />                (< (length pt2) 4)<br />           )<br />      )<br />      (princ "\n非法顶点。")<br />      (princ)<br />     )<br />     (t<br />      (setq a  (angle pt1 pt2)<br />            b  (distance pt1 pt2)<br />            h  (- (* b bul 0.5))<br />            r  (/ (+ (* b b) (* 4 h h)) 8 h)<br />            hh (- r h)<br />            x  (+ (/ (+ (car pt1) (car pt2)) 2) (* hh (sin a)))<br />            y  (- (/ (+ (cadr pt1) (cadr pt2)) 2) (* hh (cos a)))<br />      )<br />      (list (list x y) r)<br />     )<br />   )<br /> )

livewin

很给力!谢谢

sicky111

谢谢分享。

yzhyxh

qxlonmsn

感谢大神的代码   非常有用   犇  犇  犇

hh_lj007

正是想要找的方法，谢谢楼主！