XDGE几何库应用(30)--多边形最小二乘圆圆心

st788796

相关介绍
http://www.xdcad.net/forum/thread-668040-1-1.html

1. ;;点集二乘圆圆心
   
2. (defun pnts:center2 (pts / lst l)
   
3.   (setq        pts (mapcar '(lambda (x / an r)
   
4.                        (setq an        (angle '(0.0 0.0 0.0) x)
   
5.                              r        (distance '(0.0 0.0) x)
   
6.                        )
   
7.                        (list (* r (cos an)) (* r (sin an)))
   
8.                      )
   
9.                     pts
   
10.             )
    
11.         l   (length pts)
    
12.   )
    
13.   (list        (/ (* 2 (apply '+ (mapcar 'car pts))) l)
    
14.         (/ (* 2 (apply '+ (mapcar 'cadr pts))) l)
    
15.   )
    
16. )
    
17. (defun c:tt (/ e pts pcen pl3d ray3d cc3d an an1 p ptl pt)
    
18.   (if (setq e (car (entsel)))
    
19.     (progn
    
20.       (setq pts          (xdrx_getsamplept e)
    
21.             pcen  (apply 'xdrx_points_centroid pts)
    
22.             pl3d  (xdge::constructor e)
    
23.             ray3d (xdge::constructor "kRay3d");_构造 Ray3d
    
24.             cc3d  (xdge::constructor "kCurveCurveInt3d");_构造相交体
    
25.             an          0.
    
26.             an1          (/ pi 180)
    
27. 
    
28.       )
    
29.       ;;以中心为原点对多边形进行等角度分割采样
    
30.       (repeat 360
    
31.         (xdge::setpropertyvalue ray3d "set" pcen (polar pcen an 1.));_重新设置 Ray3d
    
32.         (xdge::setpropertyvalue
    
33.           cc3d "set" pl3d ray3d        '(0.0 0.0 1.0));_重新设置相交体
    
34.         (if (/= (xdge::getpropertyvalue cc3d "numIntPoints") 0)
    
35.           (setq        ptl (cons (xdge::getpropertyvalue cc3d "intPoint" 0) ptl);_求交点
    
36.                 an  (+ an an1)
    
37.           )
    
38.         )
    
39.       )
    
40.       (xdrx_point_make pcen);_质心
    
41.       (xdge::free)
    
42.       (setq p  (pnts:center2
    
43.                  (mapcar '(lambda (x) (mapcar '- x pcen)) ptl)
    
44.                );_二乘圆圆心
    
45.             pt (xdrx_point_make (mapcar '+ pcen p))
    
46.       )
    
47.       (xdrx_entity_setcolor pt 1)
    
48.     )
    
49.   )
    
50.   (princ)
    
51. )
    

;;点集二乘圆圆心<br /> (defun pnts:center2 (pts / lst l)<br />   (setq        pts (mapcar '(lambda (x / an r)<br />                        (setq an        (angle '(0.0 0.0 0.0) x)<br />                              r        (distance '(0.0 0.0) x)<br />                        )<br />                        (list (* r (cos an)) (* r (sin an)))<br />                      )<br />                     pts<br />             )<br />         l   (length pts)<br />   )<br />   (list        (/ (* 2 (apply '+ (mapcar 'car pts))) l)<br />         (/ (* 2 (apply '+ (mapcar 'cadr pts))) l)<br />   )<br /> )<br /> (defun c:tt (/ e pts pcen pl3d ray3d cc3d an an1 p ptl pt)<br />   (if (setq e (car (entsel)))<br />     (progn<br />       (setq pts          (xdrx_getsamplept e)<br />             pcen  (apply 'xdrx_points_centroid pts)<br />             pl3d  (xdge::constructor e)<br />             ray3d (xdge::constructor "kRay3d");_构造 Ray3d<br />             cc3d  (xdge::constructor "kCurveCurveInt3d");_构造相交体<br />             an          0.<br />             an1          (/ pi 180)<br /> <br />       )<br />       ;;以中心为原点对多边形进行等角度分割采样<br />       (repeat 360<br />         (xdge::setpropertyvalue ray3d "set" pcen (polar pcen an 1.));_重新设置 Ray3d<br />         (xdge::setpropertyvalue<br />           cc3d "set" pl3d ray3d        '(0.0 0.0 1.0));_重新设置相交体<br />         (if (/= (xdge::getpropertyvalue cc3d "numIntPoints") 0)<br />           (setq        ptl (cons (xdge::getpropertyvalue cc3d "intPoint" 0) ptl);_求交点<br />                 an  (+ an an1)<br />           )<br />         )<br />       )<br />       (xdrx_point_make pcen);_质心<br />       (xdge::free)<br />       (setq p  (pnts:center2<br />                  (mapcar '(lambda (x) (mapcar '- x pcen)) ptl)<br />                );_二乘圆圆心<br />             pt (xdrx_point_make (mapcar '+ pcen p))<br />       )<br />       (xdrx_entity_setcolor pt 1)<br />     )<br />   )<br />   (princ)<br /> )<br />

下载地址：晓东文库

zhangq_cai1

理论用于实践~

Lisphk

ST大师，这个半径怎么求？

st788796

Lisphk 发表于 2016-6-6 09:14
ST大师，这个半径怎么求？

1／m ∑√（x-ui）^2+（y-vi）^2
（/ (apply ' + (mapcar '(lambda(x) (distance x pc)) pts)) (length pts)）