精确快速求spline线的最小点及最大点

lgx9612

      之前确因有一贴妙作,发贴不过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. )

(;计算最小点,最大点主程式,spline线<br /> ;ss_minbox_lgx是选集 <br /> ;.......................略<br /> (setq nx_lgx (sslength  ss_minbox_lgx))<br /> (repeat nx_lgx<br /> (setq partname_lgx  (ssname ss_minbox_lgx (setq nx_lgx (1- nx_lgx))))<br /> (setq typepart_lgx (cdr (assoc 0 (entget partname_lgx))))<br /> (if (/= typepart_lgx "SPLINE")<br /> (progn<br /> ;不是spline线的计算<br /> )<br /> (calspline_lgx partname_lgx)<br /> )  <br /> )<br /> )<br /> (defun calspline_lgx(partname_lgx)<br /> (setq pl-temp_min_point nil)<br /> (setq pl-temp_max_point nil)<br /> (setq p_lgx nil)<br /> (setq lowest nil)<br /> (setq bbb (entget partname_lgx))<br /> (while (/= (assoc 10  bbb) nil)<br /> (setq p1_lgx (assoc 10  bbb))<br /> (setq p_lgx (cons (cdr p1_lgx) p_lgx))<br /> (setq bbb (vl-remove p1_lgx bbb))<br /> )<br /> (setq pl-temp_min_point (apply 'mapcar (cons 'min p_lgx)))<br /> (setq pl-temp_max_point (apply 'mapcar (cons 'max p_lgx)))<br /> (setq pl-temp_min_point (list (- (car pl-temp_min_point) 3)(- (cadr pl-temp_min_point) 3)))<br /> (setq pl-temp_max_point (list (+ (car pl-temp_max_point) 3)(+ (cadr pl-temp_max_point) 3)))<br /> (setq vla-object(vlax-ename->vla-object partname_lgx))<br /> (progn<br /> (setq pl_min_point pl-temp_min_point)<br /> (setq pl_max_point pl-temp_max_point)<br /> (setq zzl (/ (- (car pl_max_point)(car pl_min_point)) 10))<br /> (setq  givenPnt (list (car pl_min_point)(cadr pl_max_point)))<br /> (setq type_lgx "y")<br /> (calculate_lgx )<br /> (setq y2 topnt)<br /> )<br /> (progn<br /> (setq pl_min_point pl-temp_min_point)<br /> (setq pl_max_point (list (car pl-temp_max_point)(cadr pl-temp_min_point)))<br /> (setq zzl (/ (- (car pl_max_point)(car pl_min_point)) 10))<br /> (setq  givenPnt (list (car pl_min_point)(cadr pl_max_point)))<br /> (setq type_lgx "y" )<br /> (calculate_lgx)<br /> (setq y1 topnt)<br /> )<br /> (progn<br /> (setq pl_min_point pl-temp_min_point)<br /> (setq pl_max_point (list (car pl-temp_min_point)(cadr pl-temp_max_point)))<br /> (setq zzl (/ (- (cadr pl_max_point)(cadr pl_min_point)) 10))<br /> (setq  givenPnt pl_min_point)<br /> (setq type_lgx "x" )<br /> (calculate_lgx )<br /> (setq x1 topnt)<br /> )<br /> (progn<br /> (setq pl_min_point (list (car pl-temp_max_point)(cadr pl-temp_min_point)))<br /> (setq pl_max_point pl-temp_max_point)<br /> (setq zzl (/ (- (cadr pl_max_point)(cadr pl_min_point)) 10))<br /> (setq  givenPnt pl_min_point)<br /> (setq type_lgx "x" )<br /> (calculate_lgx )<br /> (setq x2 topnt)<br /> )<br /> (setq llpoint (list (car x1) (cadr y1)))<br /> (setq urpoint (list (car x2) (cadr y2)))<br /> )<br /> (defun calculate_lgx()<br /> (setq lowest nil)  <br /> (while (>= zzl 0.02)<br /> (repeat 11<br /> (setq topnt (vlax-curve-getClosestPointTo vla-object givenPnt ))<br /> ;;;(command "line" givenPnt topnt "")<br /> (if (or (>= lowest (distance  givenPnt topnt)) (= lowest nil))<br /> (progn<br /> (setq lowest (distance  givenPnt topnt))<br /> (setq lowestpoint givenPnt)<br /> )<br /> )<br /> (if (= type_lgx "y")<br /> (setq  givenPnt (list (+ (car givenPnt) zzl)(cadr pl_max_point)))<br /> (setq  givenPnt (list (car pl_min_point)(+ (cadr givenPnt) zzl)))<br /> )<br /> )<br /> (if (= type_lgx "y")  <br /> (setq  givenPnt (list (- (car lowestpoint) zzl) (cadr pl_max_point)))<br /> (setq  givenPnt (list (car pl_min_point)  (- (cadr lowestpoint) zzl)))  <br /> )  <br /> (setq zzl (/ zzl 5))<br /> )<br /> (setq topnt (vlax-curve-getClosestPointTo vla-object lowestpoint ));;;givenPnt改为 lowestpoint<br /> (setq lowest nil)  <br /> )

lucas_3333

Gu_xl 发表于 2014-7-17 14:18

风树

什么叫做最小点和最大点啊，有哪方面的应用呢？科普一下吧

自贡黄明儒

什么是最小点？x最小还丫最小。没看明白

lgx9612

不好意思!是最小包容盒子的最小点和最大点.

风树

楼主，你是哪里人啊，写"程式"的人好少哦

自贡黄明儒

lgx9612 发表于 2014-7-16 20:30
不好意思!是最小包容盒子的最小点和最大点.

你用了whi|e，是不是逼近法。G版写的那个更简洁

lgx9612

(;计算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)  
)

风树

能介绍下思路吗？

lgx9612

风树 发表于 2014-7-16 20:47
能介绍下思路吗？

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

风树

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

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