三维spline如何提取坐标？

tataki

+ S* Y: W, i/ n# o& H: a: T
0 E) m. `+ Y0 W! P+ N8 B, n- V各位好，最近遇到一个问题，在cad中做了两个闭合的三维spline，我是想算一算他们两者最近的距离是多少，当然不用负责的公式就算，就想着用最直观，也是很简单的方法来计算。见附图。
我的想法是在红色spline上将这条曲线上的每个点计算到蓝色曲线上的每个点距离，也就是一个二次循环，最后找出最小值就行。然而在操作中遇到了两个问题：
" N- f) g# J0 Q% K! q  V1.在spline的图元了保存的控制点坐标怎么提取出来？以前只提过其中的一项，现在面对这几百个数据不知道该怎么提取了。
+ b* x, l: @; X( T" w2.提取的数据怎么存放呢？我以前是学c++的，在C里面肯定用个数组循环就行了，但是不怎么会lisp，在lisp里面没有数组的概念，因此怎么存大量数据，还有后期的怎么读取都让我很迷惑，请各位懂得朋友帮个忙。在网上也找到了几个提spline坐标的程序，可是里面使用vl那些命令，这个我不知道怎么用，也没法修改成我想要的，哪位高手能帮下我，非常感谢！
现附上 这两条曲线的dwg文件。


7 W) l4 `" ~+ V, @9 C今早一打开我的帖子，就看到斑竹的答复！很激动，没想到这么快就回复了！
非常感谢斑竹的帮助！不但给出了源代码，还有详细的注解，斑竹的认真热心和负责让我肃然起敬，以前的我提的几个问题也是斑竹给解决的。刚看了一下发帖时间，半夜02:53....看到这个时间，明白人心理都懂得，再多赞美的词也显得苍白了....
5 p% Y1 h9 t) v' f: r0 l# N辛苦斑竹了，今天周末请多休息！

woaishuijia

1. 
   " b! P3 _  ~9 g
2. (defun c:aaa (/ ss ft fd c1 c2 i j dmin d)
   ( T4 W# V4 Z  d$ M, q: N8 R
3.   (vl-load-com)    ;将 Visual LISP 扩展功能加载到 AutoLISP
   
4.   (setq     ;将一个或多个符号的值设置为相应表达式的值
   
5.     ss (vla-add    ;用 vla- 前缀函数调用 activeX 对象(下面的代码--该方法的第一个参数--指定的选择集集合对象)的 add 方法
   ! o- p! E% _, [9 X) A- s
6.   (vla-get-selectionsets  ;用 vla-get- 前缀函数获取 activeX 对象--当前文档对象--的属性:选择集集合子对象,即使用 add 方法的集合对象
   0 [' U  `: k9 j  h5 [+ d
7.     (vla-get-activedocument ;用 vla-get- 前缀函数获取 activeX 对象--当前 AutoCAD 应用程序对象--的属性:活动文档子对象
   ! b' W7 c' _6 l/ d
8.       (vlax-get-acad-object) ;用 vlax-get-acad-object 函数检索当前 AutoCAD 任务中的顶层 AutoCAD 应用程序对象
   
9.     )
   
10.   )
    $ j! G/ P9 p9 ^5 q) ?
11.   "ss"    ;选择集名称
    3 ]6 p& U# s& i& C9 D9 Y! J
12.        )
    
13.     ft (vlax-make-safearray vlax-vbinteger '(0 . 0))
    
14.      ;定义选择集过滤器的过滤类型数组.用 vlax-make-safearray 函数创建安全数组,整形,一维,一个元素(0 to 0).
    . z. L* w' [' j* o7 ]7 E/ w9 i3 A: `0 Q
15.     fd (vlax-make-safearray vlax-vbvariant '(0 . 0))
    
16.      ;定义选择集过滤器的过滤值数组.用 vlax-make-safearray 函数创建安全数组,变体,一维,一个元素(0 to 0).
    % Q( v" [/ l5 ^/ i3 j6 X" ^% u& t/ I
17.   )
    
18.   (vlax-safearray-put-element fd 0 "spline" )
    
19.      ;用 vlax-safearray-put-element 函数为数组 fd 的第一个元素(索引号 0)赋值"spline"
    4 H0 k2 J- c0 W8 l# s3 c% _
20.   (vla-selectonscreen ss ft fd)  ;用 vla- 前缀函数调用 ss 对象的 selectonscreen 方法,由用户在屏幕上选取样条曲线对象
    . Q( G+ m: s, f* d  h) a; u
21.   (if (= (vla-get-count ss) 2)  ;如果 ss 集合的 count 属性等于2(即用户选择了两条曲线)则执行下面的代码
    
22.     (progn
    3 r3 y& x2 S5 T/ o8 l; n
23.       (setq c1 (vla-get-controlpoints ;用 vla-get- 前缀函数获取样条曲线对象的控制点坐标数组并赋值给c1(变量类型为变体,值为双精度数组)
    " d  L/ I: \  P. P
24.    (vla-item ss 0) ;用 vla- 前缀函数调用集合对象的 item 方法,获取选择集中的索引号为 0 的样条曲线对象
    
25.         )
    
26.      c2 (vla-get-controlpoints ;用 vla-get- 前缀函数获取样条曲线对象的控制点坐标数组并赋值给c2(变量类型为变体,值为双精度数组)
    
27.    (vla-item ss 1) ;用 vla- 前缀函数调用集合对象的 item 方法,获取选择集中的索引号为 1 的样条曲线对象
    
28.         )
    
29.      i  0   ;初始化外层循环变量
    1 j$ c* H/ P) ^8 C2 d5 k' b+ i
30.       )
    
31.       (repeat    ;循环,外层
    
32. (/ (1+ (vlax-safearray-get-u-bound
    2 H: @3 f9 h, x4 F' [3 [& w' n$ E
33.      ;用 vlax-safearray-get-u-bound 函数获取数组最大下标.第一个参数是数组,第二个参数是维数
    % x6 M  \8 f2 {0 s# X! s
34.    (vlax-variant-value c1);用 vlax-variant-value 函数获取变体变量 c1 的值(双精度数组)
    
35.    1
    
36.         )
    7 u0 I3 ^5 D# g6 {3 J
37.     )
    : a5 A1 q  _3 }$ [$ Z
38.     3
    1 }2 \6 G5 u$ m" a8 M8 o
39. )    ;用数组最大下标加1再除以3,所得结果做为循环次数--数组中每三个元素为一组代表一个控制点
      S  c4 S& e9 y( V4 P# T7 q1 t- q4 G
40.   (setq j 0)   ;初始化内层循环变量
    " O9 g$ Z! H" s3 F. o% ^, g; ~
41.   (repeat
    
42.     (/ (1+ (vlax-safearray-get-u-bound
    6 K+ \6 {1 J( j9 S. O* ?) l7 @
43.      ;第二条样条曲线控制点坐标数组的最大下标加1再除以3,做为内层循环次数
    8 B% s/ x  I- r2 w6 O, p
44.       (vlax-variant-value c2)
    
45.       1
    & U8 a; i* n. _1 x4 C* s
46.     )
    7 Y% j! z3 Z' w- |
47.        )
    
48.        3
    0 G' m1 @% r9 M3 J
49.     )
    ; ^9 M; J8 e# ^7 R' U, i3 e# o: m
50.      (setq d (distance  ;计算第一条曲线索引号为 i 的控制点到第二条曲线索引号为 j 的控制点之间的三维距离
    
51.         (list  ;从 c1 数组中取出第一条曲线索引号为 i 的控制点的点表
    
52.    (vlax-safearray-get-element
    
53.      ;用 vlax-safearray-get-element 函数获取数组中的指定的元素
    + H/ `  G; l& n* u) n3 m# q! v
54.      (vlax-variant-value c1) ;数组
    
55.      (* i 3) ;数组元素下标
    " [, S& n' i0 _7 k
56.    )
    
57.    (vlax-safearray-get-element
      s! ~; d' b4 o8 j& Q/ Y" L
58.      (vlax-variant-value c1)
    
59.      (+ (* i 3) 1)
    4 c# h, x0 v7 X4 Q6 I, ^
60.    )
    + @9 ^, ^( E; n/ _
61.    (vlax-safearray-get-element
    
62.      (vlax-variant-value c1)
    % _8 q# E& i) u" [7 C  w3 i
63.      (+ (* i 3) 2)
    , _. M5 X7 f, J) l, l
64.    )
    ) b* y7 b; n6 |3 K& `
65.         )
    
66.         (list  ;从 c2 数组中取出第二条曲线索引号为 j 的控制点的点表
    
67.    (vlax-safearray-get-element
    
68.      (vlax-variant-value c2)
    
69.      (* j 3)
    
70.    )
    ( E: C# d7 J0 x# E0 i. m: u" P9 L
71.    (vlax-safearray-get-element
    0 ~+ Z6 a7 w, S+ \
72.      (vlax-variant-value c2)
    0 p) H3 p5 F8 a
73.      (+ (* j 3) 1)
    1 a3 x; a' e7 T; N. I7 ?
74.    )
    
75.    (vlax-safearray-get-element
    
76.      (vlax-variant-value c2)
    
77.      (+ (* j 3) 2)
    / n. T8 H* |) N7 D! B) P1 c( l
78.    )
    
79.         )
    ) C) N" q) Y2 d( @) [% V6 ~
80.       )
    
81.      )
    2 h4 k7 ^" w4 d+ \( j
82.      (if dmin   ;检查最小距离是否为 nil (是否首次)
    & F3 L( i" w# P% z
83.        (if (> dmin d)  ;不是首次,比较计算所得的距离与已保存的最小距离,如发现更小的距离则保存
    ; m2 _( ]* ~; J% v
84.   (setq dmin d)
    * W( Y& E$ M: i+ G
85.        )
    
86.        (setq dmin d)  ;首次,把首个距离记录为最小距离
    
87.      )
    
88.      (setq j (1+ j))  ;内层循环变量递增
    1 `! u: r! f! D/ G1 F
89.   )
    6 D' F9 \* p% u! k& l3 B
90.   (setq i (1+ i))  ;外层循环变量递增
    . b  N9 D$ @2 d3 |8 m2 a* }
91.       )
    
92.     )
    
93.   )
    
94.   (vla-delete ss)   ;用 vla- 前缀函数调用 ss 集合的 delete 方法,删除用过的选择集
    % m0 i& r8 M4 j' c
95.   (princ dmin)    ;命令行输出计算结果
    5 s: |1 d3 W  z# P6 _8 R8 O
96.   (princ)    ;静默退出
    
97. )
    

复制代码

flyhighone

好资料，收藏了

yoyoho

感谢版主分享程序,收藏学习了!

scliukejun

楼上好强！不过，如果自已动手，去用VBA吧，应该比LISP好用