设为首页收藏本站CAD论坛CAD教程CAD下载明经币充值
切换到宽版

明经CAD社区

 找回密码
 注册

QQ登录

只需一步，快速开始

搜索
明经CAD社区»社区 › 开发篇 › AutoLISP/Visual LISP 编程技术 › 精确快速求spline线的最小点及最大点
12345下一页
返回列表 发新帖
查看: 6017|回复: 45

[讨论] 精确快速求spline线的最小点及最大点

  [复制链接]
lgx9612
发表于 2014-7-16 20:21 | 显示全部楼层 |阅读模式
本帖最后由 lgx9612 于 2014-7-27 13:56 编辑

      之前确因有一贴妙作,发贴不过10贴的不给看,故在不了解版规的情况下发了不少图片充数,确有炫耀和广告之嫌.在此实在是不好意思,虽写出几个小程序但本人笨拙,写的代码量很多,就不一一展示了.对本人来说计算最小点,最大点虽教程到处都可找到,但spline线对我来说是小有困难,费了我很大的心思.下面就spline线的最小点及最大点给出如下代码,请大家批评指正.
  2. (;计算最小点,最大点主程式,spline线
  3. ;ss_minbox_lgx是选集
  4. ;.......................略
  5. (setq nx_lgx (sslength  ss_minbox_lgx))
  6. (repeat nx_lgx
  7. (setq partname_lgx  (ssname ss_minbox_lgx (setq nx_lgx (1- nx_lgx))))
  8. (setq typepart_lgx (cdr (assoc 0 (entget partname_lgx))))
  9. (if (/= typepart_lgx "SPLINE")
  10. (progn
  11. ;不是spline线的计算
  12. )
  13. (calspline_lgx partname_lgx)
  14. )  
  15. )
  16. )
  17. (defun calspline_lgx(partname_lgx)
  18. (setq pl-temp_min_point nil)
  19. (setq pl-temp_max_point nil)
  20. (setq p_lgx nil)
  21. (setq lowest nil)
  22. (setq bbb (entget partname_lgx))
  23. (while (/= (assoc 10  bbb) nil)
  24. (setq p1_lgx (assoc 10  bbb))
  25. (setq p_lgx (cons (cdr p1_lgx) p_lgx))
  26. (setq bbb (vl-remove p1_lgx bbb))
  27. )
  28. (setq pl-temp_min_point (apply 'mapcar (cons 'min p_lgx)))
  29. (setq pl-temp_max_point (apply 'mapcar (cons 'max p_lgx)))
  30. (setq pl-temp_min_point (list (- (car pl-temp_min_point) 3)(- (cadr pl-temp_min_point) 3)))
  31. (setq pl-temp_max_point (list (+ (car pl-temp_max_point) 3)(+ (cadr pl-temp_max_point) 3)))
  32. (setq vla-object(vlax-ename->vla-object partname_lgx))
  33. (progn
  34. (setq pl_min_point pl-temp_min_point)
  35. (setq pl_max_point pl-temp_max_point)
  36. (setq zzl (/ (- (car pl_max_point)(car pl_min_point)) 10))
  37. (setq  givenPnt (list (car pl_min_point)(cadr pl_max_point)))
  38. (setq type_lgx "y")
  39. (calculate_lgx )
  40. (setq y2 topnt)
  41. )
  42. (progn
  43. (setq pl_min_point pl-temp_min_point)
  44. (setq pl_max_point (list (car pl-temp_max_point)(cadr pl-temp_min_point)))
  45. (setq zzl (/ (- (car pl_max_point)(car pl_min_point)) 10))
  46. (setq  givenPnt (list (car pl_min_point)(cadr pl_max_point)))
  47. (setq type_lgx "y" )
  48. (calculate_lgx)
  49. (setq y1 topnt)
  50. )
  51. (progn
  52. (setq pl_min_point pl-temp_min_point)
  53. (setq pl_max_point (list (car pl-temp_min_point)(cadr pl-temp_max_point)))
  54. (setq zzl (/ (- (cadr pl_max_point)(cadr pl_min_point)) 10))
  55. (setq  givenPnt pl_min_point)
  56. (setq type_lgx "x" )
  57. (calculate_lgx )
  58. (setq x1 topnt)
  59. )
  60. (progn
  61. (setq pl_min_point (list (car pl-temp_max_point)(cadr pl-temp_min_point)))
  62. (setq pl_max_point pl-temp_max_point)
  63. (setq zzl (/ (- (cadr pl_max_point)(cadr pl_min_point)) 10))
  64. (setq  givenPnt pl_min_point)
  65. (setq type_lgx "x" )
  66. (calculate_lgx )
  67. (setq x2 topnt)
  68. )
  69. (setq llpoint (list (car x1) (cadr y1)))
  70. (setq urpoint (list (car x2) (cadr y2)))
  71. )
  72. (defun calculate_lgx()
  73. (setq lowest nil)  
  74. (while (>= zzl 0.02)
  75. (repeat 11
  76. (setq topnt (vlax-curve-getClosestPointTo vla-object givenPnt ))
  77. ;;;(command "line" givenPnt topnt "")
  78. (if (or (>= lowest (distance  givenPnt topnt)) (= lowest nil))
  79. (progn
  80. (setq lowest (distance  givenPnt topnt))
  81. (setq lowestpoint givenPnt)
  82. )
  83. )
  84. (if (= type_lgx "y")
  85. (setq  givenPnt (list (+ (car givenPnt) zzl)(cadr pl_max_point)))
  86. (setq  givenPnt (list (car pl_min_point)(+ (cadr givenPnt) zzl)))
  87. )
  88. )
  89. (if (= type_lgx "y")  
  90. (setq  givenPnt (list (- (car lowestpoint) zzl) (cadr pl_max_point)))
  91. (setq  givenPnt (list (car pl_min_point)  (- (cadr lowestpoint) zzl)))  
  92. )  
  93. (setq zzl (/ zzl 5))
  94. )
  95. (setq topnt (vlax-curve-getClosestPointTo vla-object lowestpoint ));;;givenPnt改为 lowestpoint
  96. (setq lowest nil)  
  97. )
回复

使用道具 举报

lucas_3333
发表于 2014-9-9 10:25 | 显示全部楼层
本帖最后由 lucas_3333 于 2016-12-3 14:40 编辑
Gu_xl 发表于 2014-7-17 14:18

回复 支持 反对

使用道具 举报

风树
发表于 2014-7-16 20:25 | 显示全部楼层
什么叫做最小点和最大点啊，有哪方面的应用呢？科普一下吧
回复 支持 反对

使用道具 举报

自贡黄明儒
发表于 2014-7-16 20:28 来自手机 | 显示全部楼层
什么是最小点？x最小还丫最小。没看明白
回复 支持 反对

使用道具 举报

lgx9612
 楼主| 发表于 2014-7-16 20:30 | 显示全部楼层
不好意思!是最小包容盒子的最小点和最大点.
回复 支持 反对

使用道具 举报

风树
发表于 2014-7-16 20:40 | 显示全部楼层
楼主，你是哪里人啊，写"程式"的人好少哦
回复 支持 反对

使用道具 举报

自贡黄明儒
发表于 2014-7-16 20:41 来自手机 | 显示全部楼层
lgx9612 发表于 2014-7-16 20:30
不好意思!是最小包容盒子的最小点和最大点.

你用了whi|e，是不是逼近法。G版写的那个更简洁
回复 支持 反对

使用道具 举报

lgx9612
 楼主| 发表于 2014-7-16 20:45 | 显示全部楼层
本帖最后由 lgx9612 于 2014-7-27 13:59 编辑

(;计算spline线最小包容盒子的左下角点和右上角点
;ss_minbox_lgx是选集
;.....................省略
(setq nx_lgx (sslength  ss_minbox_lgx))
(repeat nx_lgx
(setq partname_lgx  (ssname ss_minbox_lgx (setq nx_lgx (1- nx_lgx))))
(setq typepart_lgx (cdr (assoc 0 (entget partname_lgx))));;;加上这行,这行少写了
(if (/= typepart_lgx "SPLINE")
(progn
;不是spline线的计算
)
(calspline_lgx partname_lgx)
)  
)
)
(defun calspline_lgx(partname_lgx)
(setq pl-temp_min_point nil)
(setq pl-temp_max_point nil)
(setq p_lgx nil)
(setq lowest nil)
(setq bbb (entget partname_lgx))
(while (/= (assoc 10  bbb) nil)
(setq p1_lgx (assoc 10  bbb))
(setq p_lgx (cons (cdr p1_lgx) p_lgx))
(setq bbb (vl-remove p1_lgx bbb))
)
(setq pl-temp_min_point (apply 'mapcar (cons 'min p_lgx)))
(setq pl-temp_max_point (apply 'mapcar (cons 'max p_lgx)))
(setq pl-temp_min_point (list (- (car pl-temp_min_point) 3)(- (cadr pl-temp_min_point) 3)))
(setq pl-temp_max_point (list (+ (car pl-temp_max_point) 3)(+ (cadr pl-temp_max_point) 3)))
(setq vla-object(vlax-ename->vla-object partname_lgx))
(progn
(setq pl_min_point pl-temp_min_point)
(setq pl_max_point pl-temp_max_point)
(setq zzl (/ (- (car pl_max_point)(car pl_min_point)) 10))
(setq  givenPnt (list (car pl_min_point)(cadr pl_max_point)))
(setq type_lgx "y")
(calculate_lgx )
(setq y2 topnt)
)
(progn
(setq pl_min_point pl-temp_min_point)
(setq pl_max_point (list (car pl-temp_max_point)(cadr pl-temp_min_point)))
(setq zzl (/ (- (car pl_max_point)(car pl_min_point)) 10))
(setq  givenPnt (list (car pl_min_point)(cadr pl_max_point)))
(setq type_lgx "y" )
(calculate_lgx)
(setq y1 topnt)
)
(progn
(setq pl_min_point pl-temp_min_point)
(setq pl_max_point (list (car pl-temp_min_point)(cadr pl-temp_max_point)))
(setq zzl (/ (- (cadr pl_max_point)(cadr pl_min_point)) 10))
(setq  givenPnt pl_min_point)
(setq type_lgx "x" )
(calculate_lgx )
(setq x1 topnt)
)
(progn
(setq pl_min_point (list (car pl-temp_max_point)(cadr pl-temp_min_point)))
(setq pl_max_point pl-temp_max_point)
(setq zzl (/ (- (cadr pl_max_point)(cadr pl_min_point)) 10))
(setq  givenPnt pl_min_point)
(setq type_lgx "x" )
(calculate_lgx )
(setq x2 topnt)
)
(setq llpoint (list (car x1) (cadr y1)))
(setq urpoint (list (car x2) (cadr y2)))
)
(defun calculate_lgx()
(setq lowest nil)  
(while (>= zzl 0.02)
(repeat 11
(setq topnt (vlax-curve-getClosestPointTo vla-object givenPnt ))
;;;(command "line" givenPnt topnt "")
(if (or (>= lowest (distance  givenPnt topnt)) (= lowest nil))
(progn
(setq lowest (distance  givenPnt topnt))
(setq lowestpoint givenPnt)
)
)
(if (= type_lgx "y")
(setq  givenPnt (list (+ (car givenPnt) zzl)(cadr pl_max_point)))
(setq  givenPnt (list (car pl_min_point)(+ (cadr givenPnt) zzl)))
)
)
(if (= type_lgx "y")  
(setq  givenPnt (list (- (car lowestpoint) zzl) (cadr pl_max_point)))
(setq  givenPnt (list (car pl_min_point)  (- (cadr lowestpoint) zzl)))  
)  
(setq zzl (/ zzl 5))
)
(setq topnt (vlax-curve-getClosestPointTo vla-object lowestpoint ))
(setq lowest nil)  
)
回复 支持 反对

使用道具 举报

风树
发表于 2014-7-16 20:47 | 显示全部楼层
能介绍下思路吗？
回复 支持 反对

使用道具 举报

lgx9612
 楼主| 发表于 2014-7-16 21:17 | 显示全部楼层
风树 发表于 2014-7-16 20:47
能介绍下思路吗？

      通过曲线自身的参数求出能够框住曲线的外框来,再用外框边上取点求出距离曲线最近点,四边用同样的方法分别求出最大x,y最小x,y即可.
回复 支持 反对

使用道具 举报

风树
发表于 2014-7-16 21:33 | 显示全部楼层
lgx9612 发表于 2014-7-16 21:17
通过曲线自身的参数求出能够框住曲线的外框来,再用外框边上取点求出距离曲线最近点,四边用同样的方 ...

很好，
另外的这个方法也许可行：
原曲线端点+极值点排序后确定y方向上的最大最小值，
原曲线旋转90度后用同样的方式确定原来x方向上的最大最小值

回复 支持 反对

使用道具 举报

12345下一页
返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 注册

本版积分规则

小黑屋|手机版|CAD论坛|CAD教程|CAD下载|联系我们|关于明经|明经通道 ( 粤ICP备05003914号 )  
©2000-2023 明经通道 版权所有 本站代码，在未取得本站及作者授权的情况下，不得用于商业用途

GMT+8, 2024-4-20 06:18 , Processed in 5.769517 second(s), 28 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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