三维spline如何提取坐标？

tataki

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

/ ]7 Q" z7 X. t* c- G& w1 k
今早一打开我的帖子，就看到斑竹的答复！很激动，没想到这么快就回复了！
非常感谢斑竹的帮助！不但给出了源代码，还有详细的注解，斑竹的认真热心和负责让我肃然起敬，以前的我提的几个问题也是斑竹给解决的。刚看了一下发帖时间，半夜02:53....看到这个时间，明白人心理都懂得，再多赞美的词也显得苍白了....
辛苦斑竹了，今天周末请多休息！

woaishuijia

1. 
   # I& m! H( y: t5 j# `: E% y
2. (defun c:aaa (/ ss ft fd c1 c2 i j dmin d)
   2 B( W/ u* D. D  T; U
3.   (vl-load-com)    ;将 Visual LISP 扩展功能加载到 AutoLISP
   
4.   (setq     ;将一个或多个符号的值设置为相应表达式的值
   2 H: {6 [& F; r7 S6 V
5.     ss (vla-add    ;用 vla- 前缀函数调用 activeX 对象(下面的代码--该方法的第一个参数--指定的选择集集合对象)的 add 方法
   % D8 P+ K, j, Y3 `3 c% L0 @) @
6.   (vla-get-selectionsets  ;用 vla-get- 前缀函数获取 activeX 对象--当前文档对象--的属性:选择集集合子对象,即使用 add 方法的集合对象
   
7.     (vla-get-activedocument ;用 vla-get- 前缀函数获取 activeX 对象--当前 AutoCAD 应用程序对象--的属性:活动文档子对象
   9 j, M% T# T$ y: Z& e2 Q+ X% A
8.       (vlax-get-acad-object) ;用 vlax-get-acad-object 函数检索当前 AutoCAD 任务中的顶层 AutoCAD 应用程序对象
   
9.     )
   
10.   )
    
11.   "ss"    ;选择集名称
    
12.        )
    
13.     ft (vlax-make-safearray vlax-vbinteger '(0 . 0))
    2 Q" X% ?/ O$ w2 O, F) Z  U
14.      ;定义选择集过滤器的过滤类型数组.用 vlax-make-safearray 函数创建安全数组,整形,一维,一个元素(0 to 0).
    3 e( \) T4 P2 U% t5 a
15.     fd (vlax-make-safearray vlax-vbvariant '(0 . 0))
    
16.      ;定义选择集过滤器的过滤值数组.用 vlax-make-safearray 函数创建安全数组,变体,一维,一个元素(0 to 0).
    
17.   )
    $ {7 z3 G" I$ {9 I, f; [$ @
18.   (vlax-safearray-put-element fd 0 "spline" )
    4 f! D6 p& z$ p& _; q( y2 P, A
19.      ;用 vlax-safearray-put-element 函数为数组 fd 的第一个元素(索引号 0)赋值"spline"
    9 ]. j, G9 C6 g* w0 |* U3 o3 b
20.   (vla-selectonscreen ss ft fd)  ;用 vla- 前缀函数调用 ss 对象的 selectonscreen 方法,由用户在屏幕上选取样条曲线对象
    
21.   (if (= (vla-get-count ss) 2)  ;如果 ss 集合的 count 属性等于2(即用户选择了两条曲线)则执行下面的代码
    
22.     (progn
    
23.       (setq c1 (vla-get-controlpoints ;用 vla-get- 前缀函数获取样条曲线对象的控制点坐标数组并赋值给c1(变量类型为变体,值为双精度数组)
    
24.    (vla-item ss 0) ;用 vla- 前缀函数调用集合对象的 item 方法,获取选择集中的索引号为 0 的样条曲线对象
    ' C: y7 Q$ L* ^2 O" {
25.         )
    ; i/ u- l* [, H/ s( w% r& |
26.      c2 (vla-get-controlpoints ;用 vla-get- 前缀函数获取样条曲线对象的控制点坐标数组并赋值给c2(变量类型为变体,值为双精度数组)
    
27.    (vla-item ss 1) ;用 vla- 前缀函数调用集合对象的 item 方法,获取选择集中的索引号为 1 的样条曲线对象
    / {& W: A& h% n9 w; l2 j
28.         )
    % a5 M  ^, Y. H
29.      i  0   ;初始化外层循环变量
    
30.       )
    
31.       (repeat    ;循环,外层
    
32. (/ (1+ (vlax-safearray-get-u-bound
    9 B+ X, B3 a1 n6 P' q
33.      ;用 vlax-safearray-get-u-bound 函数获取数组最大下标.第一个参数是数组,第二个参数是维数
    & ]/ v3 ~- k/ J
34.    (vlax-variant-value c1);用 vlax-variant-value 函数获取变体变量 c1 的值(双精度数组)
    2 w- K# U0 r, m/ Y. }7 F+ w2 A
35.    1
    2 P3 R3 n/ p# p9 e* ?5 J
36.         )
    
37.     )
    # l3 ^2 h* X: G' O( d" V( ]* b0 N
38.     3
    
39. )    ;用数组最大下标加1再除以3,所得结果做为循环次数--数组中每三个元素为一组代表一个控制点
    
40.   (setq j 0)   ;初始化内层循环变量
    3 U* d6 o3 B7 z2 N" u
41.   (repeat
    
42.     (/ (1+ (vlax-safearray-get-u-bound
    
43.      ;第二条样条曲线控制点坐标数组的最大下标加1再除以3,做为内层循环次数
    5 C' x$ o; @, k0 }) E' a( S
44.       (vlax-variant-value c2)
    9 R/ E' q! Q8 G7 B5 x
45.       1
    
46.     )
    
47.        )
    8 @2 o$ k0 v2 r) |3 S/ ?+ O* x
48.        3
    
49.     )
    
50.      (setq d (distance  ;计算第一条曲线索引号为 i 的控制点到第二条曲线索引号为 j 的控制点之间的三维距离
    4 A( }: T+ F* w( c: e1 s7 R0 D
51.         (list  ;从 c1 数组中取出第一条曲线索引号为 i 的控制点的点表
    
52.    (vlax-safearray-get-element
    7 j' f5 D" G2 c$ q( ~5 p7 V0 R1 ^
53.      ;用 vlax-safearray-get-element 函数获取数组中的指定的元素
    ) c; W5 V7 v+ M8 c* R; M; `
54.      (vlax-variant-value c1) ;数组
    8 p" I0 ^) M, ?. E. k  N
55.      (* i 3) ;数组元素下标
    
56.    )
    
57.    (vlax-safearray-get-element
    
58.      (vlax-variant-value c1)
    
59.      (+ (* i 3) 1)
    
60.    )
    
61.    (vlax-safearray-get-element
    
62.      (vlax-variant-value c1)
    
63.      (+ (* i 3) 2)
    6 T6 b3 P, z! j" G
64.    )
    
65.         )
    1 v6 O$ b! D7 Y  P" Z
66.         (list  ;从 c2 数组中取出第二条曲线索引号为 j 的控制点的点表
    $ I# |: I7 n+ q2 b/ Q: d. y
67.    (vlax-safearray-get-element
    
68.      (vlax-variant-value c2)
    
69.      (* j 3)
    
70.    )
    
71.    (vlax-safearray-get-element
    8 L, b& g0 S( G" ~0 H
72.      (vlax-variant-value c2)
    ! ~. T3 p3 s# q
73.      (+ (* j 3) 1)
    
74.    )
    * H+ O3 c4 w( [7 N8 N
75.    (vlax-safearray-get-element
    * F( {- g8 s, q, ^
76.      (vlax-variant-value c2)
    
77.      (+ (* j 3) 2)
    
78.    )
    ! B) E5 ]7 \% r# @& T) q! T4 z
79.         )
    3 l5 _1 z6 D5 P! n
80.       )
    8 g9 Y+ L! C2 H# P+ n/ G
81.      )
    3 {# m( S5 ?- ~; o- D# U( B" Y
82.      (if dmin   ;检查最小距离是否为 nil (是否首次)
    
83.        (if (> dmin d)  ;不是首次,比较计算所得的距离与已保存的最小距离,如发现更小的距离则保存
    5 n  v: Z2 l" b% H4 }
84.   (setq dmin d)
    ( D$ o3 F* Q/ w  f/ S$ d' m7 R
85.        )
    + ~: P& @$ v. _& n: R% L+ S/ \
86.        (setq dmin d)  ;首次,把首个距离记录为最小距离
    8 q. v/ I( u0 Y6 n
87.      )
    2 a7 H- _: P5 _9 I  }
88.      (setq j (1+ j))  ;内层循环变量递增
    
89.   )
    
90.   (setq i (1+ i))  ;外层循环变量递增
      ?4 D& b. T1 f8 l
91.       )
    5 I: o4 q# K3 N: E8 j
92.     )
    ) f$ l  p: e6 v5 _8 ?# {- A+ g& t% e
93.   )
    / C1 r- u& o' c: k2 l
94.   (vla-delete ss)   ;用 vla- 前缀函数调用 ss 集合的 delete 方法,删除用过的选择集
    3 f& c9 c9 ^- Z- z, s' W, g' Y
95.   (princ dmin)    ;命令行输出计算结果
    
96.   (princ)    ;静默退出
    ! L& Q% D% U2 L3 E8 h6 x% T" k1 `, s1 F
97. )
    

复制代码

flyhighone

好资料，收藏了

yoyoho

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

scliukejun

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