找回密码
 立即注册

QQ登录

只需一步,快速开始

扫一扫,访问微社区

查看: 17738|回复: 9

[每日一码] 查找离散点边界

[复制链接]

已领礼包: 40个

财富等级: 招财进宝

发表于 2016-4-13 11:54:13 | 显示全部楼层 |阅读模式

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

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

×
利用Auto-cad中的lisp语言找散点边界,适用于凸的,或凹陷不太过度的散点集合

  1. ;;; ========================================================================
  2. ;;; =================
  3. ;;; Title:getarea
  4. ;;;
  5. ;;; Function:  getarea
  6. ;;;
  7. ;;; date:2014/04/15,16,17,18,19,20
  8. ;;;
  9. ;;; By :Chen Meijun
  10. ;;;
  11. ;;; ========================================================================
  12. ;;; =================
  13. (defun c:getarea (/)
  14.   (defun min-of-list (mylist /)
  15.     (setq mymidList nil)
  16.     (setq L (length mylist))
  17.     (setq j 1)
  18.     (setq ax (nth 0 mylist))
  19.     (repeat (- L 1)
  20.       (setq bx (nth j mylist))
  21.       (if (< bx ax)
  22.         (setq ax bx)
  23.       )
  24.       (setq j (1+ j))
  25.     )
  26.     (setq ax ax)
  27.   )                                       ; 自定义函数求表中最小值
  28.   (defun max-of-list (mylist /)
  29.     (setq mymidList nil)
  30.     (setq L (length mylist))
  31.     (setq j 1)
  32.     (setq ax (nth 0 mylist))
  33.     (repeat (- L 1)
  34.       (setq bx (nth j mylist))
  35.       (if (> bx ax)
  36.         (setq ax bx)
  37.       )
  38.       (setq j (1+ j))
  39.     )
  40.     (setq ax ax)
  41.   )                                       ; 自定义函数求表中最大值
  42.   (defun deleteList (myList coord /)
  43.     (setq mymidList nil)
  44.     (setq len (length myList))
  45.     (setq nnn 0)
  46.     (while (< nnn len)
  47.       (if (/= nnn coord)
  48.         (setq mymidList (cons (nth nnn myList) mymidList))
  49.       )
  50.       (setq nnn (1+ nnn))
  51.     )
  52.     (setq re (reverse mymidList))
  53.   )                                       ; 自定义删除表元素

  54.   ;; (setq k (getint  "局部扫描区域半边长,取水平方向最近点距离的1.7倍"))
  55.   (setq ss (ssget))
  56.   (setq n (sslength ss))
  57.   (setq list_pt nil)
  58.   (setq i 0)
  59.   (repeat n
  60.     (setq spt (ssname ss i))
  61.     (setq ept (entget spt))
  62.     (setq l (length ept))
  63.     (setq xyz (cdr (assoc 10 ept)))
  64.     (setq list_pt (cons xyz list_pt))
  65.     (setq i (1+ i))
  66.   )                                       ; 将全部点制表
  67.   (setq n0 (sslength ss))
  68.   (setq listofAx nil)
  69.   (setq i 0)
  70.   (repeat n0
  71.     (setq Ax (car (nth i list_pt)))
  72.     (setq listofAx (cons Ax listofAx))
  73.     (setq i (1+ i))
  74.   )
  75.   (setq Ax1 (min-of-list listofAx))    ; 找到最小x坐标
  76.   (setq Ax2 (max-of-list listofAx))    ; 找到最大x坐标
  77.   (setq j 0)
  78.   (repeat l
  79.     (setq Rx (car (nth j list_pt)))    ; 取表中任意一点
  80.     (if (equal Rx Ax1)
  81.       (progn
  82.         (setq Ay1 (car (cdr (nth j list_pt))))
  83.         (setq Az1 (car (cdr (cdr (nth j list_pt)))))
  84.       )
  85.     )
  86.     (if (equal Rx Ax2)
  87.       (progn
  88.         (setq Ay2 (car (cdr (nth j list_pt))))
  89.         (setq Az2 (car (cdr (cdr (nth j list_pt)))))
  90.       )
  91.     )
  92.     (setq j (1+ j))
  93.   )
  94.   (setq pstart (list Ax1 Ay1 Az1))     ; 找到起始点坐标
  95.   (setq pstart0 pstart)
  96.   (setq Ax0 Ax1)
  97.   (setq Ay0 Ay1)
  98.   (setq Az0 Az1)
  99.   (setq pstart0 (list Ax0 Ay0 Az0))    ; 将起始点赋给pstart0,方便以后循环调
  100.                                        ; 用
  101.   (setq pend (list Ax2 Ay2 Az2))       ; 找到终点坐标
  102.   (setq n1 1)
  103.   (while (> n1 0)                       ; 第一个大循环开始
  104.     (setq rx (+ Ax1 6))
  105.     (setq lx (- Ax1 6))
  106.     (setq uy (+ Ay1 6))
  107.     (setq ly (- Ay1 6))                       ; 确定局部扫描范围
  108.     (setq list_new nil)
  109.     (setq j 0)
  110.     (setq l (length list_pt))
  111.     (repeat l
  112.       (setq pt (nth j list_pt))
  113.       (setq ptx (car pt))
  114.       (setq pty (car (cdr pt)))
  115.       (if (and
  116.             (and
  117.               (< ptx rx)
  118.               (> ptx lx)
  119.             )
  120.             (and
  121.               (< pty uy)
  122.               (> pty ly)
  123.             )
  124.           )
  125.         (setq list_new (cons pt list_new))
  126.       )
  127.       (setq j (1+ j))
  128.     )                                       ; 找到小范围内点
  129.     (setq l (length list_new))
  130.     (setq list_ne nil)
  131.     (setq n 0)
  132.     (repeat l
  133.       (if (equal (nth n list_new) pstart)
  134.         (setq temp n)
  135.       )
  136.       (setq n (1+ n))
  137.     )
  138.     (setq list_new (deletelist list_new temp)) ; 删除倒序排列表中与pstart相?
  139.                                        ; ??
  140.     (setq list_ne list_new)               ; 点倒序排列的表
  141.     (setq list_new (reverse list_new))
  142.     (setq l (length list_new))
  143.     (setq angleAjlist nil)
  144.     (setq angleBjlist nil)
  145.     (setq n 0)
  146.     (repeat l
  147.       (setq A1 (nth n list_new))
  148.       (setq angleBj (angle pstart A1))
  149.       (setq angleBjlist (cons angleBj angleBjlist)) ; B角是倒序的(与对应点?
  150.                                        ; 喾矗?
  151.       (if (> angleBj 1.58)
  152.         (setq angleAj (- 7.85 angleBj))
  153.         (setq angleAj (- 1.57 angleBj))
  154.       )
  155.       (setq angleAjlist (cons angleAj angleAjlist)) ; A角是顺序的(同对应点?
  156.                                        ; ?
  157.       (setq n (1+ n))
  158.     )
  159.     (setq minangleAj (min-of-list angleAjlist))        ; 调用自定义函数求表中最小元
  160.                                        ; 素
  161.     (setq l (length angleAjlist))
  162.     (setq n 0)
  163.     (repeat l
  164.       (setq na (nth n angleAjlist))
  165.       (if (= minangleAj na)
  166.         (setq temp n)

  167.       )
  168.       (setq n (1+ n))
  169.     )
  170.     (setq pa (nth temp list_ne))       ; 找到最小角对应点
  171.     (progn
  172.       (command "line" pstart pa "")
  173.       (setq pstart pa)
  174.       (setq Ax1 (car pstart))
  175.       (setq Ay1 (car (cdr pstart)))    ; (setq pstart(list 100 200 300))
  176.       (setq Az1 (car (cdr (cdr pstart))))

  177.     )
  178.     (if (equal pstart pend)
  179.       (setq n1 -1)
  180.     )
  181.   )                                       ; 第一大循环结束
  182.   (setq pstart pstart0)
  183.   (setq Ax1 Ax0)
  184.   (setq Ay1 Ay0)
  185.   (setq Az1 Az0)
  186.   (setq pend (list Ax2 Ay2 Az2))
  187.   (setq jbmaxangleAj (* 1 pi))
  188.   (setq n2 1)
  189.   (while (> n2 0)                       ; 第二大循环开始
  190.     (setq rx (+ Ax1 6))
  191.     (setq lx (- Ax1 6))
  192.     (setq uy (+ Ay1 6))
  193.     (setq ly (- Ay1 6))                       ; 确定局部扫描范围
  194.     (setq list_new nil)
  195.     (setq j 0)
  196.     (setq l (length list_pt))
  197.     (repeat l
  198.       (setq npt (nth j list_pt))
  199.       (setq nptx (car npt))
  200.       (setq npty (car (cdr npt)))
  201.       (if (and
  202.             (and
  203.               (< nptx rx)
  204.               (> nptx lx)
  205.             )
  206.             (and
  207.               (< npty uy)
  208.               (> npty ly)
  209.             )
  210.           )
  211.         (setq list_new (cons npt list_new))
  212.       )
  213.       (setq j (1+ j))
  214.     )                                       ; 找到小范围内点
  215.     (setq l (length list_new))
  216.     (setq list_ne nil)
  217.     (setq n 0)
  218.     (repeat l
  219.       (if (equal (nth n list_new) pstart)
  220.         (setq temp n)
  221.       )
  222.       (setq n (1+ n))
  223.     )
  224.     (setq list_new (deletelist list_new temp)) ; 删除倒序排列表中与pstart相?
  225.                                        ; ??
  226.     (setq list_ne list_new)               ; 点倒序排列的表
  227.     (setq list_new (reverse list_new))
  228.     (setq l (length list_new))
  229.     (setq angleAjlist nil)
  230.     (setq angleBjlist nil)
  231.     (setq n 0)
  232.     (repeat l
  233.       (setq A1 (nth n list_new))
  234.       (setq angleBj (angle pstart A1))
  235.       (setq angleBjlist (cons angleBj angleBjlist)) ; B角是倒序的(与对应点?
  236.                                        ; 喾矗?


  237.       (if (> angleBj 1.58)
  238.         (setq angleAj (- 7.85 angleBj))
  239.         (setq angleAj (- 1.57 angleBj))
  240.       )
  241.       (setq angleAjlist (cons angleAj angleAjlist)) ; A角是反序的(同对应点?
  242.                                        ; ?
  243.       (setq n (1+ n))
  244.     )


  245.     (setq jbangleAjlist nil)
  246.     (setq l (length list_new))
  247.     (setq n 0)
  248.     (repeat l
  249.       (setq angleAj (nth n angleAjlist))
  250.       (if (and
  251.             (> angleAj (- jbmaxangleAj 1.2))
  252.             (< angleAj (+ jbmaxangleAj 2))
  253.           )
  254.         (setq jbangleAjlist (cons angleAj jbangleAjlist))
  255.       )
  256.       (setq n (1+ n))
  257.     )
  258.     (setq jbmaxangleAj (max-of-list jbangleAjlist)) ; 调用自定义函数求表中最
  259.                                        ; 大元素

  260.     (setq bb jbmaxangleAj)

  261.     ;; (setq maxangleAj(max-of-list
  262.     ;; angleAjlist));;;;;;;;;;;;;;;;;;;;调用自定义函数求表中最大元素
  263.     (setq l (length angleAjlist))
  264.     (setq n 0)
  265.     (repeat l
  266.       (setq na (nth n angleAjlist))
  267.       (if (= jbmaxangleAj na)
  268.         (setq temp n)

  269.       )
  270.       (setq n (1+ n))
  271.     )
  272.     (setq pa (nth temp list_ne))       ; 找到最大角对应点
  273.     (progn
  274.       (command "line" pstart pa "")
  275.       (setq pstart pa)

  276.       (setq Ax1 (car pstart))
  277.       (setq Ay1 (car (cdr pstart)))

  278.       (setq Az1 (car (cdr (cdr pstart))))

  279.     )

  280.     ;; (setq list_pused(cons list_new list_pused))
  281.     (if (equal pstart pend)
  282.       (setq n2 -1)
  283.     )
  284.     (setq jbmaxangleAj bb)

  285.   )                                       ; 第二大循环结束

  286. )


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

已领礼包: 5578个

财富等级: 富甲天下

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

使用道具 举报

已领礼包: 40个

财富等级: 招财进宝

 楼主| 发表于 2016-4-13 13:53:00 | 显示全部楼层
HLCAD 发表于 2016-4-13 13:43
版主,要是有一个将离散点拟合成圆弧的就好了

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

使用道具 举报

已领礼包: 604个

财富等级: 财运亨通

发表于 2016-4-13 15:50:54 | 显示全部楼层
井喷了哇,N版是最近打了鸡血。赞!!!
论坛插件加载方法
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
如果你在论坛求助问题,并且已经从坛友或者管理的回复中解决了问题,请把帖子标题加上【已解决】;
如何回报帮助你解决问题的坛友,一个好办法就是给对方加【D豆】,加分不会扣除自己的积分,做一个热心并受欢迎的人!
回复 支持 反对

使用道具 举报

已领礼包: 1857个

财富等级: 堆金积玉

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

使用道具 举报

已领礼包: 604个

财富等级: 财运亨通

发表于 2016-4-14 21:14:03 来自手机 | 显示全部楼层
newer 发表于 2016-4-13 13:53
那得有规则啊,随机的离散点怎么成圆弧?

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

使用道具 举报

已领礼包: 1857个

财富等级: 堆金积玉

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

使用道具 举报

已领礼包: 1094个

财富等级: 财源广进

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

使用道具 举报

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

本版积分规则

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

GMT+8, 2024-3-29 18:51 , Processed in 0.394957 second(s), 45 queries , Gzip On.

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

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