设为首页收藏本站

晓东CAD家园-论坛

 找回密码
 立即注册

QQ登录

只需一步,快速开始

扫一扫,访问微社区

查看: 1471|回复: 3

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

[复制链接]

点击这里给我发消息

已领礼包: 1259个

财富等级: 财源广进

发表于 2014-11-6 08:40:13 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
相关介绍
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.       (repeat 360
  30.         (xdge::setpropertyvalue ray3d "set" pcen (polar pcen an 1.));_重新设置 Ray3d
  31.         (xdge::setpropertyvalue
  32.           cc3d "set" pl3d ray3d        '(0.0 0.0 1.0));_重新设置相交体
  33.         (if (/= (xdge::getpropertyvalue cc3d "numIntPoints") 0)
  34.           (setq        ptl (cons (xdge::getpropertyvalue cc3d "intPoint" 0) ptl);_求交点
  35.                 an  (+ an an1)
  36.           )
  37.         )
  38.       )
  39.       (xdrx_point_make pcen);_质心
  40.       (xdge::free)
  41.       (setq p  (pnts:center2
  42.                  (mapcar '(lambda (x) (mapcar '- x pcen)) ptl)
  43.                );_二乘圆圆心
  44.             pt (xdrx_point_make (mapcar '+ pcen p))
  45.       )
  46.       (xdrx_entity_setcolor pt 1)
  47.     )
  48.   )
  49.   (princ)
  50. )


下载地址:晓东文库


本帖被以下淘专辑推荐:

论坛插件加载方法
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
如果你在论坛求助问题,并且已经从坛友或者管理的回复中解决了问题,请把帖子标题加上【已解决】;
如何回报帮助你解决问题的坛友,一个好办法就是给对方加【D豆】,加分不会扣除自己的积分,做一个热心并受欢迎的人!

已领礼包: 849个

财富等级: 财运亨通

发表于 2014-11-6 10:10:09 | 显示全部楼层
理论用于实践~
论坛插件加载方法
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
如果你在论坛求助问题,并且已经从坛友或者管理的回复中解决了问题,请把帖子标题加上【已解决】;
如何回报帮助你解决问题的坛友,一个好办法就是给对方加【D豆】,加分不会扣除自己的积分,做一个热心并受欢迎的人!
回复 支持 反对

使用道具 举报

已领礼包: 19个

财富等级: 恭喜发财

发表于 2016-6-6 09:14:36 | 显示全部楼层
ST大师,这个半径怎么求?
论坛插件加载方法
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
如果你在论坛求助问题,并且已经从坛友或者管理的回复中解决了问题,请把帖子标题加上【已解决】;
如何回报帮助你解决问题的坛友,一个好办法就是给对方加【D豆】,加分不会扣除自己的积分,做一个热心并受欢迎的人!
回复 支持 反对

使用道具 举报

点击这里给我发消息

已领礼包: 1259个

财富等级: 财源广进

 楼主| 发表于 2016-6-6 23:34:57 来自手机 | 显示全部楼层
本帖最后由 st788796 于 2016-6-7 07:38 编辑
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))
论坛插件加载方法
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
如果你在论坛求助问题,并且已经从坛友或者管理的回复中解决了问题,请把帖子标题加上【已解决】;
如何回报帮助你解决问题的坛友,一个好办法就是给对方加【D豆】,加分不会扣除自己的积分,做一个热心并受欢迎的人!
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|申请友链|小黑屋|手机版|Archiver|辽公网安备|晓东CAD家园 ( 辽ICP备15016793号 )  

GMT+8, 2020-5-29 02:45 , Processed in 0.073833 second(s), 38 queries , Gzip On, WinCache On.

Powered by Discuz! X3.4

Copyright © 2001-2020, Tencent Cloud.

快速回复 返回顶部 返回列表