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

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

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

copy


;;;;;;;;;;;;;;;;;;;;;;;;
;;;     角度转换        
;;;;;;;;;;;;;;;;;;;;;;;;
;;;参数说明：
(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 (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坐标吗？为什么要用点表表示？
最后你还是没有说你到底是要做什么

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

发表回复

精选推荐