三维spline如何提取坐标？

tataki

各位好，最近遇到一个问题，在cad中做了两个闭合的三维spline，我是想算一算他们两者最近的距离是多少，当然不用负责的公式就算，就想着用最直观，也是很简单的方法来计算。见附图。
) G; _7 a. r: r# A5 q! Q我的想法是在红色spline上将这条曲线上的每个点计算到蓝色曲线上的每个点距离，也就是一个二次循环，最后找出最小值就行。然而在操作中遇到了两个问题：
1.在spline的图元了保存的控制点坐标怎么提取出来？以前只提过其中的一项，现在面对这几百个数据不知道该怎么提取了。
8 o! u4 d  G* Q$ _5 |2.提取的数据怎么存放呢？我以前是学c++的，在C里面肯定用个数组循环就行了，但是不怎么会lisp，在lisp里面没有数组的概念，因此怎么存大量数据，还有后期的怎么读取都让我很迷惑，请各位懂得朋友帮个忙。在网上也找到了几个提spline坐标的程序，可是里面使用vl那些命令，这个我不知道怎么用，也没法修改成我想要的，哪位高手能帮下我，非常感谢！
; ~1 |8 O3 t9 p9 c9 J! q现附上 这两条曲线的dwg文件。


$ i, r" A/ n9 y& Y0 {8 W! `今早一打开我的帖子，就看到斑竹的答复！很激动，没想到这么快就回复了！
9 w6 |1 U. D) {7 z7 L4 G非常感谢斑竹的帮助！不但给出了源代码，还有详细的注解，斑竹的认真热心和负责让我肃然起敬，以前的我提的几个问题也是斑竹给解决的。刚看了一下发帖时间，半夜02:53....看到这个时间，明白人心理都懂得，再多赞美的词也显得苍白了....
辛苦斑竹了，今天周末请多休息！

woaishuijia

) A! G+ ?- m8 B% ]+ Q  H. d- s9 U

1. 
   
2. (defun c:aaa (/ ss ft fd c1 c2 i j dmin d)
   * I) k% F' [: o0 h" a+ t
3.   (vl-load-com)    ;将 Visual LISP 扩展功能加载到 AutoLISP
   
4.   (setq     ;将一个或多个符号的值设置为相应表达式的值
   
5.     ss (vla-add    ;用 vla- 前缀函数调用 activeX 对象(下面的代码--该方法的第一个参数--指定的选择集集合对象)的 add 方法
   
6.   (vla-get-selectionsets  ;用 vla-get- 前缀函数获取 activeX 对象--当前文档对象--的属性:选择集集合子对象,即使用 add 方法的集合对象
   
7.     (vla-get-activedocument ;用 vla-get- 前缀函数获取 activeX 对象--当前 AutoCAD 应用程序对象--的属性:活动文档子对象
   & \: |/ L+ x: s6 x5 l- ?% @4 T
8.       (vlax-get-acad-object) ;用 vlax-get-acad-object 函数检索当前 AutoCAD 任务中的顶层 AutoCAD 应用程序对象
   
9.     )
   8 q8 U" C/ H. v& r+ A1 I
10.   )
    + x8 P9 m' Y; }: K7 G
11.   "ss"    ;选择集名称
    
12.        )
    ; h, K5 n4 P7 J/ u# Z! `, A0 n0 ~1 B
13.     ft (vlax-make-safearray vlax-vbinteger '(0 . 0))
    2 `# }  A3 _- S& D8 [! {- {
14.      ;定义选择集过滤器的过滤类型数组.用 vlax-make-safearray 函数创建安全数组,整形,一维,一个元素(0 to 0).
    
15.     fd (vlax-make-safearray vlax-vbvariant '(0 . 0))
    
16.      ;定义选择集过滤器的过滤值数组.用 vlax-make-safearray 函数创建安全数组,变体,一维,一个元素(0 to 0).
    ' q, X, W1 N' E" _4 `
17.   )
    
18.   (vlax-safearray-put-element fd 0 "spline" )
    : r) k. t" j6 ]/ q. a! y
19.      ;用 vlax-safearray-put-element 函数为数组 fd 的第一个元素(索引号 0)赋值"spline"
    9 O5 S  @) K1 v* V$ h; Y/ G
20.   (vla-selectonscreen ss ft fd)  ;用 vla- 前缀函数调用 ss 对象的 selectonscreen 方法,由用户在屏幕上选取样条曲线对象
    # O4 Q! t( X$ E2 r; k! W0 c3 w  C
21.   (if (= (vla-get-count ss) 2)  ;如果 ss 集合的 count 属性等于2(即用户选择了两条曲线)则执行下面的代码
    + Q0 k* u( o" f( x3 V8 A
22.     (progn
    . ~6 b# Q5 c$ }- @; ]) |) v/ j
23.       (setq c1 (vla-get-controlpoints ;用 vla-get- 前缀函数获取样条曲线对象的控制点坐标数组并赋值给c1(变量类型为变体,值为双精度数组)
    
24.    (vla-item ss 0) ;用 vla- 前缀函数调用集合对象的 item 方法,获取选择集中的索引号为 0 的样条曲线对象
    
25.         )
    ) y1 h1 z# {% T
26.      c2 (vla-get-controlpoints ;用 vla-get- 前缀函数获取样条曲线对象的控制点坐标数组并赋值给c2(变量类型为变体,值为双精度数组)
    " B, ~' Q+ Z" v0 i
27.    (vla-item ss 1) ;用 vla- 前缀函数调用集合对象的 item 方法,获取选择集中的索引号为 1 的样条曲线对象
    : E6 x$ L* `" B& S0 X9 {
28.         )
    
29.      i  0   ;初始化外层循环变量
    
30.       )
    3 N1 t$ g, Q: ?6 Z: i+ \
31.       (repeat    ;循环,外层
    
32. (/ (1+ (vlax-safearray-get-u-bound
    
33.      ;用 vlax-safearray-get-u-bound 函数获取数组最大下标.第一个参数是数组,第二个参数是维数
    0 ^$ ?& I( k" y/ E6 e
34.    (vlax-variant-value c1);用 vlax-variant-value 函数获取变体变量 c1 的值(双精度数组)
    / t. z3 a# o% J  ~3 n0 H5 [2 {
35.    1
    
36.         )
    8 d+ w" f! ^8 y8 w
37.     )
    * o' }" i! ]8 r: ?! P2 {' W) M$ @
38.     3
    * S' @' \. E! w" E. ^- z4 r. A8 z
39. )    ;用数组最大下标加1再除以3,所得结果做为循环次数--数组中每三个元素为一组代表一个控制点
    : @: H, z& A6 b- [& s& f* v
40.   (setq j 0)   ;初始化内层循环变量
    ; g' m, J7 I- k' g2 ^
41.   (repeat
    ! g7 ^+ C7 t0 B* t
42.     (/ (1+ (vlax-safearray-get-u-bound
    
43.      ;第二条样条曲线控制点坐标数组的最大下标加1再除以3,做为内层循环次数
    5 t( W  c5 d1 m1 b6 B
44.       (vlax-variant-value c2)
    
45.       1
    8 ]- U7 x0 p1 B7 Z- W, n" i% J
46.     )
    0 d2 w" x# i- ^
47.        )
    9 Z- b# Y$ g: D0 t2 \  W
48.        3
    
49.     )
    
50.      (setq d (distance  ;计算第一条曲线索引号为 i 的控制点到第二条曲线索引号为 j 的控制点之间的三维距离
    
51.         (list  ;从 c1 数组中取出第一条曲线索引号为 i 的控制点的点表
    % z7 f, k- o5 t- L, Q; ?5 x1 u
52.    (vlax-safearray-get-element
    : @" `3 g1 S$ C+ K4 s1 @. K
53.      ;用 vlax-safearray-get-element 函数获取数组中的指定的元素
    % `; ]- g# y( p3 U5 q
54.      (vlax-variant-value c1) ;数组
    
55.      (* i 3) ;数组元素下标
    
56.    )
    8 c# \: J" @& C" `, q5 d
57.    (vlax-safearray-get-element
    
58.      (vlax-variant-value c1)
    6 g1 f/ h; c& S) u3 N6 O0 p
59.      (+ (* i 3) 1)
      ]  E9 t+ T* S
60.    )
    
61.    (vlax-safearray-get-element
    
62.      (vlax-variant-value c1)
    8 c* W# w6 ^6 i9 d% ~: H. }
63.      (+ (* i 3) 2)
    
64.    )
    $ P' R7 ~9 z4 U5 ?# M. o: Q5 U% r; f# }
65.         )
    ' X% O* ]+ {/ I$ q: e
66.         (list  ;从 c2 数组中取出第二条曲线索引号为 j 的控制点的点表
    
67.    (vlax-safearray-get-element
    3 O1 P+ Q8 h$ k2 t9 @1 e4 E% h
68.      (vlax-variant-value c2)
    , ?2 R6 ]4 C' t' [
69.      (* j 3)
    
70.    )
    
71.    (vlax-safearray-get-element
    
72.      (vlax-variant-value c2)
    5 i6 z; t; s; f; H. G+ ?: t
73.      (+ (* j 3) 1)
    
74.    )
    
75.    (vlax-safearray-get-element
    4 h' ^1 O8 w. A' f, }
76.      (vlax-variant-value c2)
    
77.      (+ (* j 3) 2)
    
78.    )
    3 }3 {7 p5 Q, N) m$ }4 M: |1 a( P
79.         )
    
80.       )
    - C/ [' ^/ q3 K% R4 J3 ~- s( x: U
81.      )
    7 Y' F# C6 I7 ~! _& w4 Q* H; K  t1 @1 ^
82.      (if dmin   ;检查最小距离是否为 nil (是否首次)
    
83.        (if (> dmin d)  ;不是首次,比较计算所得的距离与已保存的最小距离,如发现更小的距离则保存
    
84.   (setq dmin d)
    
85.        )
    1 C3 T9 l- k! Y' u9 k7 e0 Z
86.        (setq dmin d)  ;首次,把首个距离记录为最小距离
    
87.      )
      D& a7 d/ K0 v* q9 p# R* }. G
88.      (setq j (1+ j))  ;内层循环变量递增
    & M7 u: T) N9 s8 v% P  h
89.   )
    , p1 m3 Y3 A& E# w. o2 C
90.   (setq i (1+ i))  ;外层循环变量递增
    ) [* b; L' P0 Z' n% @- Q
91.       )
    
92.     )
    
93.   )
    
94.   (vla-delete ss)   ;用 vla- 前缀函数调用 ss 集合的 delete 方法,删除用过的选择集
    
95.   (princ dmin)    ;命令行输出计算结果
    . H2 i7 v. G( w7 L) ^
96.   (princ)    ;静默退出
    : e. z' s/ w$ [1 V: y" T; d
97. )
    

复制代码

flyhighone

好资料，收藏了

yoyoho

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

scliukejun

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