不确定拟合点数的样条曲线的绘制

island1000 · 发表于 2003-12-3 18:17:00

本人在解决不确定拟合点数的样条曲线的绘制时，用了一个土办法。具体是，先根据终点和起点画样条曲线。然后将样条曲线的序列数据a1取出，加入剩下不确定数量的点列。当然点列应该是(11 .. .. ..)的样式。然后(entmod a1)。这样就能画了。但是，当我移动坐标系后，出现了问题。请高手们帮我看看：（程序有点乱，请多包涵）

;;;;;;;;;;;;;;;;;;;;;;;;
;;; 角度转换
;;;;;;;;;;;;;;;;;;;;;;;;
;;;参数说明：
(defun Degrees->Radians(numberOfDegrees)
(* pi (/ numberOfDegrees 180.0))
)
;;;;;;;;;;;;;;;;;;;;;;;;
;;;二维点转换为三维点
;;;;;;;;;;;;;;;;;;;;;;;;
;;;参数说明：
(defun VD:2dPoint->3dPoint(2dPt)
(list (float (car 2dPt)) (float (cdr 2dPt)) 0.0)
)
;;;;;;;;;;;;;;;;;;;;;;;;
;;; 三维点转换为二维点
;;;;;;;;;;;;;;;;;;;;;;;;
;;;参数说明：
(defun VD:3dPoint->2dPoint(3dPt)
(list (float (car 3dPt))(float (cadr 3dPt)))
)
;;;;;;;;;;;;;;;;;;;;;;;;
;;; 绘制箭头
;;;;;;;;;;;;;;;;;;;;;;;;
;;;参数说明：
(defun vd:arrow
(pt ang ratio / ang1
angleBackward angleUpward angleDownward
arrowLength arrowWidth halfWidth pt1
p1 p2 s1
)
(setq ang1 (Degrees->Radians ang)
angleBackward (+ ang1 pi)
angleUpward (+ ang1 (/ pi 2.0))
angleDownward (+ angleUpward pi)
arrowLength (* 2.5 ratio)
arrowWidth (* 0.83 ratio)
halfWidth (/ arrowWidth 2.0)
pt1 (polar pt angleBackward arrowLength)
p1 (polar pt1 angleUpward halfWidth)
p2 (polar pt1 angleDownward halfWidth)
)
(command ".line" pt p1 p2 "c" "")
(setq s1 (ssget "C" pt p1))
(command ".hatch" "solid" s1 "")
(princ)
)
;;;;;;;;;;;;;;;;;;;;;;;;
;;;绘制坐轴及横向分格线
;;;;;;;;;;;;;;;;;;;;;;;;
;;;参数说明：
(defun vd:coordinate_drawing (p1 p2 space ratio /
nStart i nTotal a
xEnd yEnd original
tp tp1 textHeight
s1
)
(setq
nStart (fix (/ (cadr p1) space))
i nStart
nTotal (- (fix (/ (cadr p2) space)) nStart)
a (* ratio 20)
Xend (list (+ (car p1) a) (* nStart space) 0.0)
yend (list 0.0 (+ (* a 0.5) (cadr p2)) 0.0)
original (list 0.0 (cadr xend) 0.0)
a (* ratio 3.0)
tp (polar xend (* pi 1.3) (* 1.5 a))
textHeight (* ratio 3.0)
)
(vd:arrow xend 0 ratio)
(vd:arrow yend 90 ratio)
(command "line" original yend "")
(command "line" original xend "")
(setq s1 (ssget "L"))
(command "text" tp textHeight 0.0 "X" "")
(setq tp (polar yend (* pi 1.25) (* 1.414 a)))
(command "text" tp textHeight 0.0 "Y" "")
(setq tp (polar original pi a))
(if (= i 0)
(progn
(command "text" tp textHeight 0.0 "A" "")
(setq tp1
(polar
original
0.0
(+ (distance Xend original) a)
)
)
(command "text" tp1 textHeight 0.0 "A" "")
)
)
(repeat nTotal
(setq i (1+ i))
(command "offset" space s1 yend "")
(setq s1 (ssget "L"))
(cond
((= i 1) (setq str "B"))
((= i 2) (setq str "C"))
((= i 3) (setq str "D"))
((= i 4) (setq str "E"))
((= i 5) (setq str "F"))
((= i 6) (setq str "G"))
((= i 7) (setq str "H"))
((= i 8) (setq str "I"))
((= i 9) (setq str "J"))
((= i 10) (setq str "K"))
((= i 11) (setq str "L"))
((= i 12) (setq str "M"))
((= i 13) (setq str "N"))
((= i 14) (setq str "O"))
((= i 15) (setq str "P"))
((= i 16) (setq str "Q"))
((= i 17) (setq str "R"))
((= i 18) (setq str "S"))
((= i 19) (setq str "T"))
((= i 20) (setq str "U"))
((= i 21) (setq str "V"))
((= i 22) (setq str "W"))
)
(setq tp (polar tp (/ pi 2.0) space)
tp1 (polar tp1 (/ pi 2.0) space))
(command "text" tp textHeight 0.0 str "")
(command "text" tp1 textHeight 0.0 str "")
)
(princ)
)
;;;;;;;;;;;;;;;;;;;;;;;;
;;; 断面数据输入
;;;;;;;;;;;;;;;;;;;;;;;;
(defun vd:dataInput(/ fname f i blst a lst n str)
(setq fname (getfiled "Select a data file" "" "txt" 8)
f (open fname "r")
i 0
Blst nil
a 0.0
)
(while (setq str (read-line f))
(setq lst nil
n (read str)
)
(if (< n 10)
(setq str (substr str 3))
(setq str (substr str 4))
)
(setq flag (read str)
str (substr str 3)
)
(while (read str)
(setq a (read str)
i (+ 2 (strlen (rtos a 2)))
str (substr str i)
lst (cons a lst)
)
)
(setq lst (reverse lst)
lst (cons n (cons flag lst))
Blst (append Blst (list lst))
n (1+ n)
)
)
(close f)
Blst
)
;;;;;;;;;;;;;;;;;;;;;;;;
;;; 断面数据处理
;;;;;;;;;;;;;;;;;;;;;;;;
(defun VD:dataProcessing(lst space / n i a b c d l1 l2 templist)
(setq n (length lst)
i (/ n 2)
)
(repeat i
(setq a (car lst)
b (cadr lst)
l1 (cdr a)
l2 (cdr b)
c (VD:dp1 l1 space)
d (VD:dp1 l2 space)
templist (cons (car a)
(list
c
d
)
)
)
(setq lst (append (cddr lst) (list templist))
)
)
lst
)
(defun VD:dp1 (lst space / n j flag startPt endPt y temp)
(setq n (length lst)
flag (car lst)
startPt (cons (nth 2 lst) (cadr lst))
endPt (cons (cadr (reverse lst))
(last lst)
)
lst (reverse (cddr (reverse (cdddr lst))))
y (*
(fix (/ (cdr startPt) space))
space
)
j (- n 5)
)
(repeat j
(setq y (+ y space)
temp (cons (car lst) y)
lst (cdr lst)
lst (append lst (list temp))
)
)
(setq lst (cons startPt lst)
lst (reverse(cons endPt (reverse lst)))
lst (cons flag lst)
)
lst
)
;;;;;;;;;;;;;;;;;;;;;;;;
;;;二维点转换为三维点
;;;;;;;;;;;;;;;;;;;;;;;;
;;;参数说明：
(defun VD:2dPoints->3dPoints (plst / n 2dPt)
(setq n (length plst))
(repeat n
(setq plst (append plst
(list (list 11
(float (car (car plst)))
(float (cdr (car plst)))
0.0
)
)
)
plst (cdr plst)
)
)
plst
)
;;;;;;;;;;;;;;;;;;;;;;;;
;;; 断面轮廓线绘制
;;;;;;;;;;;;;;;;;;;;;;;;
(defun VD:draw (ptLst / dLst p1 p2 a1)
(setq dLst nil
dLst (VD:2dPoints->3dPoints (cdr (cadr ptLst)))
p1 (cdr (car dLst))
p2 (cdr (last dLst))
dLst (cdr (vl-remove (last dLst) dLst))
)
(princ (reverse dLst))
(command ".spline" p2 p1 "" "" "")
(setq a1 (reverse (entget (ssname (ssget "L") 0))))
(princ a1)
(setq a1 (reverse (append (list (car a1))
dLst
(cdr a1)
)
)
)
(setq a1 (entmod a1))
(princ a1)
(setq dLst (VD:2dPoints->3dPoints (cdr (caddr ptLst)))
p1 (cdr(car dLst))
p2 (cdr(last dLst))
dLst (cdr (vl-remove (last dLst) dLst))
)
(princ (reverse dLst))
(command ".spline" p2 p1 "" "" "")
(setq a1 (reverse (entget (ssname (ssget "L") 0))))
(princ a1)
(setq a1 (reverse (append (list (car a1))
dLst
(cdr a1)
)
)
)
(setq a1 (entmod a1))
(princ a1)
(setq p1 (VD:2dPoint->3dPoint (car(cdr (cadr ptLst))))
p2 (VD:2dPoint->3dPoint (car(cdr (caddr ptLst))))
)
(command ".line" p1 p2 "")
(setq p1 (VD:2dPoint->3dPoint (last (cdr(cadr ptLst))))
p2 (VD:2dPoint->3dPoint (last (cdr(caddr ptLst))))
)
(command ".line" p1 p2 "")
(setq p1 (VD:2dPoint->3dPoint (car(cdr (cadr ptLst))))
p2 (VD:2dPoint->3dPoint (last (cdr(caddr ptLst))))
)
(VD:coordinate_drawing p1 p2 20.0 2)
)
;;;;;;;;;;;;;;;;;;;;;;;;
;;; 断面尺寸标注
;;;;;;;;;;;;;;;;;;;;;;;;
(defun VD:ddim (ptLst ratio / pt tp ptlst1 i)
(setq ptlst1 (cdr (cadr ptLst))
i (length ptlst1)
textHeight (* ratio 3.0)
)
(repeat i
(setq
pt (vd:2dPoint->3dPoint (car ptlst1))
tp (polar (list 10.0 (cadr pt) 0.0) (/ pi 2) (* ratio 1.0))
ptlst1 (cdr ptlst1)
str (rtos (car pt) 2)
)
(if (= (rem (cadr pt) 20.0) 0.0)
(progn
(command "text" tp textHeight 0.0 str "")
(vd:arrow pt 0 ratio)
)
(progn
(setq p1 (list 0.0 (cadr pt) 0.0)
p2 (polar p1 0 (car pt))
)
(command "dim1" "Hor" p1 pt p2 "")
(setq p1 (list (car pt) 0.0 0.0)
p2 (polar p1 (/ pi 2) (cadr pt))
)
(command "dim1" "Ver" p1 pt p2 "")
)
)
)
(setq ptlst1 (cdr (caddr ptLst))
i (length ptlst1)
)
(repeat i
(setq
pt (vd:2dPoint->3dPoint (car ptlst1))
tp (polar (list 280 (cadr pt) 0.0) (/ pi 2) (* ratio 1.0))
ptlst1 (cdr ptlst1)
str (rtos (car pt) 2)
)
(if (= (rem (cadr pt) 20.0) 0.0)
(progn
(command "text" tp textHeight 0.0 str "")
(vd:arrow pt 0 ratio)
)
(progn
(setq p1 (list (car pt) 0.0 0.0)
p2 (polar p2 180 30.0)
)
(command "dim1" "Ver" p1 pt p2 "")
(setq p1 (list 0.0 (cadr pt) 0.0)
p2 (polar p1 0 (car pt))
)
(command "dim1" "Hor" p1 pt p2 "")
)
)
)
)
;;;;;;;;;;;;;;;;;;;;;;;;
;;; 主程序
;;;;;;;;;;;;;;;;;;;;;;;;
(defun c:VD()
(setq scmdecho (getvar "cmdecho"))
(setvar "cmdecho" 0)
(setq ptLst nil
ptLst (VD:dataInput)
ptLst (VD:dataProcessing ptLst 20.0)
n (length ptLst)
i 1
ratio 2
a (* ratio 50.0)
)
(repeat n
(setq pptLst (assoc i ptLst)
i (1+ i)
)
;;; 确定作图区域（左上角点与右下角点）
(setq leftTop (vd:2dPoint->3dPoint (last (caddr pptLst)))
rightBottom (vd:2dPoint->3dPoint (cadr (caddr pptLst)))
rightBottom (list (+ (car rightBottom) a) 0.0 0.0)
leftTop (list 0.0 (+ (* a 0.5) (cadr leftTop)) 0.0)
newOriginal (list (car rightBottom) 0.0 0.0)
)
(VD:draw pptLst)
(VD:ddim pptLst 2.0)
(command "ucs" "m" newOriginal)
)
(command "ucs" "")
;;; (VD:dataOut )
(setvar "cmdecho" scmdecho)
)

;;;;;;;;;;;;;;;;;;;;;;;; 
;;;     角度转换         
;;;;;;;;;;;;;;;;;;;;;;;; 
;;;参数说明： 
(defun Degrees->Radians(numberOfDegrees) 
    (* pi (/ numberOfDegrees 180.0)) 
) 
;;;;;;;;;;;;;;;;;;;;;;;; 
;;;二维点转换为三维点    
;;;;;;;;;;;;;;;;;;;;;;;; 
;;;参数说明： 
(defun VD:2dPoint->3dPoint(2dPt) 
    (list (float (car 2dPt)) (float (cdr 2dPt)) 0.0) 
) 
;;;;;;;;;;;;;;;;;;;;;;;; 
;;;  三维点转换为二维点 
;;;;;;;;;;;;;;;;;;;;;;;; 
;;;参数说明： 
(defun VD:3dPoint->2dPoint(3dPt) 
    (list (float (car 3dPt))(float (cadr 3dPt))) 
) 
;;;;;;;;;;;;;;;;;;;;;;;; 
;;;      绘制箭头       
;;;;;;;;;;;;;;;;;;;;;;;; 
;;;参数说明： 
(defun vd:arrow 
               (pt ang               ratio           /               ang1 
                   angleBackward           angleUpward angleDownward 
                   arrowLength arrowWidth  halfWidth   pt1 
                   p1               p2           s1 
                  ) 
  (setq        ang1              (Degrees->Radians ang) 
        angleBackward (+ ang1 pi) 
        angleUpward   (+ ang1 (/ pi 2.0)) 
        angleDownward (+ angleUpward pi) 
        arrowLength   (* 2.5 ratio) 
        arrowWidth    (* 0.83 ratio) 
        halfWidth     (/ arrowWidth 2.0) 
        pt1              (polar pt angleBackward arrowLength) 
        p1              (polar pt1 angleUpward halfWidth) 
        p2              (polar pt1 angleDownward halfWidth) 
  ) 
  (command ".line" pt p1 p2 "c" "") 
  (setq s1 (ssget "C" pt p1)) 
  (command ".hatch" "solid" s1 "") 
  (princ) 
) 
;;;;;;;;;;;;;;;;;;;;;;;; 
;;;绘制坐轴及横向分格线 
;;;;;;;;;;;;;;;;;;;;;;;; 
;;;参数说明： 
(defun vd:coordinate_drawing (p1 p2         space         ratio         / 
                                 nStart         i         nTotal         a 
                                 xEnd         yEnd         original 
                                 tp         tp1         textHeight 
                                 s1 
                                ) 
  (setq 
    nStart     (fix (/ (cadr p1) space)) 
    i               nStart 
    nTotal     (- (fix (/ (cadr p2) space)) nStart) 
    a               (* ratio 20) 
    Xend       (list (+ (car p1) a) (* nStart space) 0.0) 
    yend       (list 0.0 (+ (* a 0.5) (cadr p2)) 0.0) 
    original   (list 0.0 (cadr xend) 0.0) 
    a               (* ratio 3.0) 
    tp               (polar xend (* pi 1.3) (* 1.5 a)) 
    textHeight (* ratio 3.0) 
  ) 
  (vd:arrow xend 0 ratio) 
  (vd:arrow yend 90 ratio) 
  (command "line" original yend "") 
  (command "line" original xend "") 
  (setq s1 (ssget "L")) 
  (command "text" tp textHeight 0.0 "X" "") 
  (setq tp (polar yend (* pi 1.25) (* 1.414 a))) 
  (command "text" tp textHeight 0.0 "Y" "") 
  (setq tp (polar original pi a)) 
  (if (= i 0) 
   (progn 
     (command "text" tp textHeight 0.0 "A" "") 
     (setq tp1 
            (polar 
              original 
              0.0 
              (+ (distance Xend original) a) 
            ) 
     ) 
     (command "text" tp1 textHeight 0.0 "A" "") 
   ) 
 ) 
  (repeat nTotal 
    (setq i (1+ i)) 
    (command "offset" space s1 yend "") 
    (setq s1 (ssget "L")) 
    (cond 
      ((= i 1) (setq str "B")) 
      ((= i 2) (setq str "C")) 
      ((= i 3) (setq str "D")) 
      ((= i 4) (setq str "E")) 
      ((= i 5) (setq str "F")) 
      ((= i 6) (setq str "G")) 
      ((= i 7) (setq str "H")) 
      ((= i 8) (setq str "I")) 
      ((= i 9) (setq str "J")) 
      ((= i 10) (setq str "K")) 
      ((= i 11) (setq str "L")) 
      ((= i 12) (setq str "M")) 
      ((= i 13) (setq str "N")) 
      ((= i 14) (setq str "O")) 
      ((= i 15) (setq str "P")) 
      ((= i 16) (setq str "Q")) 
      ((= i 17) (setq str "R")) 
      ((= i 18) (setq str "S")) 
      ((= i 19) (setq str "T")) 
      ((= i 20) (setq str "U")) 
      ((= i 21) (setq str "V")) 
      ((= i 22) (setq str "W")) 
    ) 
    (setq tp (polar tp (/ pi 2.0) space) 
          tp1 (polar tp1 (/ pi 2.0) space)) 
    (command "text" tp textHeight 0.0 str "") 
    (command "text" tp1 textHeight 0.0 str "") 
 
  ) 
  (princ) 
) 
;;;;;;;;;;;;;;;;;;;;;;;; 
;;;    断面数据输入      
;;;;;;;;;;;;;;;;;;;;;;;; 
(defun vd:dataInput(/ fname f i blst a lst n str) 
  (setq        fname (getfiled "Select a data file" "" "txt" 8) 
        f     (open fname "r") 
        i     0 
        Blst  nil 
        a     0.0 
  ) 
  (while (setq str (read-line f)) 
   (setq lst nil 
         n   (read str) 
   ) 
   (if (< n 10) 
     (setq str (substr str 3)) 
     (setq str (substr str 4)) 
   ) 
    (setq flag (read str) 
          str  (substr str 3) 
    ) 
    (while (read str) 
      (setq a        (read str) 
            i        (+ 2 (strlen (rtos a 2))) 
            str        (substr str i) 
            lst        (cons a lst) 
      ) 
    ) 
    (setq lst  (reverse lst) 
          lst  (cons n (cons flag lst)) 
          Blst (append Blst (list lst)) 
          n    (1+ n) 
    ) 
  ) 
  (close f) 
  Blst 
) 
;;;;;;;;;;;;;;;;;;;;;;;; 
;;;    断面数据处理      
;;;;;;;;;;;;;;;;;;;;;;;; 
(defun VD:dataProcessing(lst space / n i a b c d l1 l2 templist) 
  (setq        n (length lst) 
        i (/ n 2) 
  ) 
(repeat        i 
  (setq        a         (car lst) 
        b         (cadr lst) 
        l1         (cdr a) 
        l2         (cdr b) 
        c         (VD:dp1 l1 space) 
        d         (VD:dp1 l2 space) 
        templist (cons (car a) 
                       (list 
                         c 
                         d 
                       ) 
                 ) 
  ) 
  (setq        lst (append (cddr lst) (list templist)) 
  ) 
) 
  lst 
) 
(defun VD:dp1 (lst space / n j flag startPt endPt y temp) 
  (setq        n        (length lst) 
        flag        (car lst) 
        startPt        (cons (nth 2 lst) (cadr lst)) 
        endPt        (cons (cadr (reverse lst)) 
                      (last lst)                      
                ) 
        lst        (reverse (cddr (reverse (cdddr lst)))) 
        y        (* 
                   (fix (/ (cdr startPt) space)) 
                   space 
                ) 
        j        (- n 5) 
  ) 
  (repeat j 
    (setq y    (+ y space) 
          temp (cons (car lst) y) 
          lst  (cdr lst) 
          lst  (append lst (list temp)) 
    ) 
  ) 
  (setq        lst   (cons startPt lst) 
        lst   (reverse(cons endPt (reverse lst))) 
        lst   (cons flag lst) 
  ) 
  lst 
) 
;;;;;;;;;;;;;;;;;;;;;;;; 
;;;二维点转换为三维点    
;;;;;;;;;;;;;;;;;;;;;;;; 
;;;参数说明： 
(defun VD:2dPoints->3dPoints (plst / n 2dPt) 
 (setq n     (length plst)) 
  (repeat n 
    (setq plst (append plst 
                       (list (list 11 
                                   (float (car (car plst))) 
                                   (float (cdr (car plst))) 
                                   0.0 
                             ) 
                       ) 
               ) 
          plst (cdr plst) 
    ) 
  ) 
  plst 
) 
;;;;;;;;;;;;;;;;;;;;;;;; 
;;;   断面轮廓线绘制    
;;;;;;;;;;;;;;;;;;;;;;;; 
(defun VD:draw (ptLst / dLst p1 p2 a1) 
  (setq        dLst nil 
        dLst (VD:2dPoints->3dPoints (cdr (cadr ptLst))) 
        p1   (cdr (car dLst)) 
        p2   (cdr (last dLst)) 
        dLst (cdr (vl-remove (last dLst) dLst)) 
  ) 
  (princ (reverse dLst)) 
  (command ".spline" p2 p1 "" "" "") 
  (setq a1 (reverse (entget (ssname (ssget "L") 0)))) 
  (princ a1) 
  (setq a1 (reverse (append (list (car a1)) 
                            dLst 
                            (cdr a1) 
                    ) 
           ) 
  ) 
  (setq a1 (entmod a1)) 
  (princ a1) 
  (setq        dLst    (VD:2dPoints->3dPoints (cdr (caddr ptLst))) 
        p1        (cdr(car dLst)) 
        p2        (cdr(last dLst)) 
        dLst    (cdr (vl-remove (last dLst) dLst)) 
  ) 
  (princ (reverse dLst)) 
  (command ".spline" p2 p1 "" "" "") 
  (setq        a1 (reverse (entget (ssname (ssget "L") 0)))) 
  (princ a1) 
  (setq        a1 (reverse (append (list (car a1)) 
                            dLst 
                            (cdr a1) 
                    ) 
           ) 
  ) 
  (setq a1 (entmod a1)) 
  (princ a1) 
  (setq        p1 (VD:2dPoint->3dPoint (car(cdr (cadr ptLst)))) 
        p2 (VD:2dPoint->3dPoint (car(cdr (caddr ptLst)))) 
  ) 
  (command ".line" p1 p2 "") 
  (setq        p1 (VD:2dPoint->3dPoint (last (cdr(cadr ptLst)))) 
        p2 (VD:2dPoint->3dPoint (last (cdr(caddr ptLst)))) 
  ) 
  (command ".line" p1 p2 "") 
  (setq        p1 (VD:2dPoint->3dPoint (car(cdr (cadr ptLst)))) 
        p2 (VD:2dPoint->3dPoint (last (cdr(caddr ptLst)))) 
  ) 
 
  (VD:coordinate_drawing p1  p2 20.0 2) 
 
) 
;;;;;;;;;;;;;;;;;;;;;;;; 
;;;   断面尺寸标注       
;;;;;;;;;;;;;;;;;;;;;;;; 
(defun VD:ddim (ptLst ratio / pt tp ptlst1 i) 
  (setq        ptlst1           (cdr (cadr ptLst)) 
        i           (length ptlst1) 
        textHeight (* ratio 3.0) 
  ) 
  (repeat i 
    (setq 
      pt     (vd:2dPoint->3dPoint (car ptlst1)) 
      tp     (polar (list 10.0 (cadr pt) 0.0) (/ pi 2) (* ratio 1.0)) 
      ptlst1 (cdr ptlst1) 
      str    (rtos (car pt) 2) 
    ) 
    (if        (= (rem (cadr pt) 20.0) 0.0) 
      (progn 
        (command "text" tp textHeight 0.0 str "") 
        (vd:arrow pt 0 ratio) 
      ) 
      (progn 
        (setq p1 (list 0.0 (cadr pt) 0.0) 
              p2 (polar p1 0 (car pt)) 
        ) 
        (command "dim1" "Hor" p1 pt p2 "") 
        (setq p1 (list (car pt) 0.0 0.0) 
              p2 (polar p1 (/ pi 2) (cadr pt)) 
        ) 
        (command "dim1" "Ver" p1 pt p2 "") 
      ) 
    ) 
  ) 
  (setq        ptlst1 (cdr (caddr ptLst)) 
        i      (length ptlst1) 
  ) 
  (repeat i 
    (setq 
      pt     (vd:2dPoint->3dPoint (car ptlst1)) 
      tp     (polar (list 280 (cadr pt) 0.0) (/ pi 2) (* ratio 1.0)) 
      ptlst1 (cdr ptlst1) 
      str    (rtos (car pt) 2) 
    ) 
    (if        (= (rem (cadr pt) 20.0) 0.0) 
      (progn 
        (command "text" tp textHeight 0.0 str "") 
        (vd:arrow pt 0 ratio) 
      ) 
      (progn 
        (setq p1 (list (car pt) 0.0 0.0) 
              p2 (polar p2 180 30.0) 
        ) 
        (command "dim1" "Ver" p1 pt p2 "") 
        (setq p1 (list 0.0 (cadr pt) 0.0) 
              p2 (polar p1 0 (car pt)) 
        ) 
        (command "dim1" "Hor" p1 pt p2 "") 
      ) 
    ) 
  ) 
) 
;;;;;;;;;;;;;;;;;;;;;;;; 
;;;        主程序       
;;;;;;;;;;;;;;;;;;;;;;;; 
(defun c:VD() 
    (setq scmdecho (getvar "cmdecho")) 
    (setvar "cmdecho" 0) 
    (setq ptLst        nil 
          ptLst        (VD:dataInput) 
          ptLst        (VD:dataProcessing ptLst 20.0) 
          n        (length ptLst) 
          i        1 
          ratio 2 
          a        (* ratio 50.0) 
 
    ) 
  (repeat n 
    (setq pptLst (assoc i ptLst) 
          i         (1+ i) 
    ) 
;;;    确定作图区域（左上角点与右下角点） 
    (setq leftTop     (vd:2dPoint->3dPoint (last (caddr pptLst))) 
          rightBottom (vd:2dPoint->3dPoint (cadr (caddr pptLst))) 
          rightBottom (list (+ (car rightBottom) a) 0.0 0.0) 
          leftTop     (list 0.0 (+ (* a 0.5) (cadr leftTop)) 0.0) 
          newOriginal (list (car rightBottom) 0.0 0.0) 
    ) 
    (VD:draw pptLst) 
    (VD:ddim pptLst 2.0) 
    (command "ucs" "m" newOriginal) 
  ) 
  (command "ucs" "") 
;;;    (VD:dataOut ) 
    (setvar "cmdecho" scmdecho) 
)

数据文件如下：

meflying · 发表于 2003-12-3 18:36:00

你的数据是什么格式，，，
看完你的数据处理以后，看那些结果，不知道你到底是要做什么？
能否详细说说你具体是要做什么？

island1000 · 发表于 2003-12-4 10:19:00

是纯文本格式的,一个断面有两条样条曲线,每行中的数据是一条样条曲线的数据,第一个数据表示第几个断面.第二个数据是样条曲线所在侧的识别号.这里只有正和反两侧.0表示正,1表示反.剩下的全是曲线上的点的数据.这里有个特殊情况，就是中间各个点的y间距是相等的，在程序中用了space变量就是这个间距。比如第一组，起点是（269.51 0.0）终点是（122.14 251.66），中间的数据全部是各个点的坐标的数据。读入进行数据处理以后得到的数据序列如下

当然，这里我用了3组同样的断面数据，因为这些数据还是我自己整理出来的，总共有20多个，我只是想测试一下，应该完全可以的。
不好意思，太麻烦你了，我先谢谢了！

meflying · 发表于 2003-12-4 10:52:00

还是不明白。。。
你说的断面是什么的断面？
为什么起点是取一个数据和0，而终点又是在数据中直接取两个，到底哪个数据是X坐标，哪个数据是Y坐标？
处理后的数据中0、20、40。。。是什么东西？Y坐标吗？为什么要用点表表示？
最后你还是没有说你到底是要做什么

账号		自动登录	找回密码
密码			注册

不确定拟合点数的样条曲线的绘制

本帖子中包含更多资源

本帖子中包含更多资源