测井曲线的读取-问题何在

mogm127

(defun C:WTCJ ()
  (command "osnap" "off" )
  (COMMAND "insert" "c:\\Temp\\AB线CD段EF隧道DZ-01-023孔.dwg" "0,0,0" "1" "1" "0")
  (COMMAND "_explode" "all")
  (COMMAND "_zoom" "e")
  (princ "\n读取综合测井文件数据，绘制曲线。")
  (if (setq FILE (getfiled "选择dat文件"
      ""
      "dat"
      4
   )
      )
    (progn
      ;;读模式打开lsp文件
      (setq A1 (open FILE "r"))
      ;;逐行读入
      (while (/= A1 nil)
 (setq STR (repeat 2 (read-line F1)))
 (if(= F1 "自然电位")
 (progn
        (setq min1 (getreal "请输入自然电位最小值"))
        (setq max1 (getreal "请输入自然电位最大值"))
        (rtos (setq mid1 (/ max1 4)));转成字符串
        (rtos (setq mid2 (/ max1 2)))
        (rtos (setq mid3 (* (/ max1 4) 3)))
        (setq mid (strcat mid1 "   " mid2 "   " mid3))
        (command "text" "s" "BG_ST" "j" "bc" "110.012,2.375,0" "2.5" "90d" mid)
        (setq STR (repeat 4 (read-line F1)));钻孔编号
        (setq LST STR)
        (command "text" "s" "BG_ST" "j" "bl" "-140.524,42.492,0" "3.6" "90d" LST);写图签
        (setq len (strlen STR))
             (setq LST1 (strcat(substr STR 1 1) "  "))
             (setq len (- len 1))
             (if (/= len 0)
                (progn
;                      (setq n 1)
                      (setq n (+ n 1))
                      (setq LST1 (strcat LST1 (substr LST1 n 1) "  "))
                      (setq len (- len 1))
                )
             )
        (command "text" "s" "BG_ST" "j" "bl" "90.960,48.634,0" "5.0" "90d" LST1);写标题
        (setq STR (repeat 3 (read-line F1)));套管长度
        (setq LST2 (strcat STR "m"))
        (command "text" "s" "BG_ST" "j" "bl" "217.39,32.308,0" "2.5" "90d" LST2)
        (atoi LST2)
        (setq LST3 (* LST2 2));1:500故乘2
       
        (setq pa1 (polar '(30.0 0.0) (* pi 1.5) LST3))
        (setq pa2 (polar '(80.0 0.0) (* pi 1.5) LST3))
        (setq pa3 (polar '(32.382 -2.608) (* pi 1.5) (- LST3 2)))
        (setq pa4 (polar '(32.154 -2.608) (* pi 1.5) (- LST3 2)))
        (command "pline" '(32.382 -2.608) pa3 "")
        (command "pline" '(32.154 -2.608) pa4 "")
        (command "pline" pa1 pa2 "")
        (command "text" "s" "BG_ST" "j" "bc" "55,-116,0" "2.0" "90d" LST3)
        (command "text" "s" "BG_ST" "j" "bc" "70,-116,0" "2.0" "90d" "套管")
        (setq STR (read-line F1));测试日期
        (command "text" "s" "BG_ST" "j" "mc" "-58.524,11.258,0" "2.4" "90d" STR)
        (setq STR (repeat 2 (read-line F1)));钻孔标高
        (setq LST4 (strcat STR "m"))
        (command "text" "s" "BG_ST" "j" "bl" "103.604,32.308,0" "2.5" "90d" LST4)
        (setq STR (repeat 140 (read-line F1)));测井数据152行
        (atof STR);转成实数
        (list STR);转到串行
        (setq data1 (car STR))
        (setq data2 (cadr STR))
        (setq data11 (atof data1))
        (setq data21 (atof data2))
        (setq Y (* data11 -2));1:500计算Y坐标
        (if(< data21 min1) (setq X 90))
        (if(> data21 max1) (setq X 130))
        (if(and (> data21 min1) (< data21 max1))
             (setq X (+ (* (/ (- data21 min1) (- max1 min1)) 40) 90));min1<data21<max1时，X=[(data21-min1)/(max1-min1)]*40+90
             (setq G1 '(X,Y))
         )   
             (while (/= STR nil)
                    (line-line())
                    (command "layer" "s" "图层8" "")
                    (command "pline" G1 G2 "w" "0.0" "0.0")
                    (setq G1 G2)
             )
             (defun line-line (STR)
             (progn
                   (setq STR (read-line F1))
                   (atof STR)
                   (list STR)
                   (setq data1 (car STR))
                   (setq data2 (cadr STR))
                   (setq data11 (atof data1))
                   (setq data21 (atof data2))
                   (setq Y (* data11 -2))
                   (if(< data21 min1) (setq X 90))
                   (if(> data21 max1) (setq X 130))
                   (if(and (> data21 min1) (< data21 max1))
                       (setq X (+ (* (/ (- data21 min1) (- max1 min1)) 40) 90))
                       (setq G2 '(X,Y))
                   )
                  
              )
              )
     ))   
   )          
      ))
     (close A1)
     (princ "\nOK!")
     (prompt "\n********* << C:WTCJ >> *********")
)
      
 问题：函数的功能是读测井曲线，并画在CAD图上。为什么会是参数类型错误呢？要怎么改？

mogm127

高手何在?进来帮帮忙吧.

xianaihua

源程序乱得一团糟，硬的头皮改了一下

2. (defun C:WTCJ ()
3.   (defun line-line ()
4.     (setq STR (read-line F1))
5.     (setq lst (strParse STR " ")) ;转换成列表
6.     (setq data1 (car lst))
7.     (setq data2 (cadr lst))
8.     (setq data11 (atof data1))
9.     (setq data21 (atof data2))
10.     (setq Y (* data11 -2))
11.     (if (< data21 min1)
12.       (setq X 90)
13.     )
14.     (if (> data21 max1)
15.       (setq X 130)
16.     )
17.     (if (and (> data21 min1) (< data21 max1))
18.       (setq
19.         X
20.   (+ (* (/ (- data21 min1) (- max1 min1)) 40) 90)
21.       )
22.     )
23.     (setq G2 (list X Y))
24.   )
25.   
26.   (command "osnap" "off")
27.   (COMMAND "insert"
28.     "C:\\Downloads\\AB线CD段EF隧道DZ-01-023孔\\AB线CD段EF隧道DZ-01-023孔.dwg"
29.     "0,0,0"       "1"
30.     "1"        "0"
31.    )
32.   (COMMAND "_explode" "all")
33.   (COMMAND "_zoom" "e")
34.   (princ "\n读取综合测井文件数据，绘制曲线。")
35.   (if (setq FILE (getfiled "选择dat文件"
36.       ""
37.       "dat"
38.       4
39.    )
40.       )
41.     (progn
42.       ;;读模式打开lsp文件
43.       (setq F1 (open FILE "r"))
44.       ;;逐行读入
45.       (while (/= A1 nil)
46. (setq STR (repeat 2 (read-line F1)))
47. (if (= STR "自然电位")
48.    (progn
49.      (setq min1 (getreal "请输入自然电位最小值"))
50.      (setq max1 (getreal "请输入自然电位最大值"))
51.      (setq mid (strcat (rtos (setq mid1 (/ max1 4))) ;转成字符串
52.          "  "
53.          (rtos (setq mid2 (/ max1 2)))
54.          "  "
55.          (rtos (setq mid3 (* (/ max1 4) 3)))
56.         )
57.      )
58.      (command "text"  "s"      "BG_ST"  "j"
59.        "bc"  "110.012,2.375,0"  "2.5"
60.        "0"  mid
61.       )
62.      (setq STR (repeat 4 (read-line F1))) ;钻孔编号
63.      (setq LST STR)
64.      (command "text"  "s"      "BG_ST"  "j"
65.        "bl"  "-140.524,42.492,0"  "3.6"
66.        "0"  LST
67.       )   ;写图签
68. ;;;     (setq len (strlen STR))
69. ;;;     (setq LST1 (strcat (substr STR 1 1) "  "))
70. ;;;     (setq len (- len 1))
71. ;;;     (if (/= len 0)
72. ;;;       (progn
73. ;;;     ;                      (setq n 1)
74. ;;;  (setq n (+ n 1))
75. ;;;  (setq LST1 (strcat LST1 (substr LST1 n 1) "  "))
76. ;;;  (setq len (- len 1))
77. ;;;       )
78. ;;;     )
79.      (command "text"  "s"      "BG_ST"  "j"
80.        "bl"  "90.960,48.634,0"  "5.0"
81.        "0"  LST
82.       )   ;写标题
83.      (setq STR (repeat 3 (read-line F1))) ;套管长度
84.      (setq LST2 (strcat STR "m"))
85.      (command "text"  "s"      "BG_ST"  "j"
86.        "bl"  "217.39,32.308,0"  "2.5"
87.        "0"  LST2
88.       )
89.      
90.      (setq LST3 (* (atof  STR) 2)) ;1:500故乘2
91.      (setq pa1 (polar '(30.0 0.0) (* pi 1.5) LST3))
92.      (setq pa2 (polar '(80.0 0.0) (* pi 1.5) LST3))
93.      (setq pa3 (polar '(32.382 -2.608) (* pi 1.5) (- LST3 2)))
94.      (setq pa4 (polar '(32.154 -2.608) (* pi 1.5) (- LST3 2)))
95.      (command "pline" '(32.382 -2.608) pa3 "")
96.      (command "pline" '(32.154 -2.608) pa4 "")
97.      (command "pline" pa1 pa2 "")
98.      (command "text" "s" "BG_ST" "j" "bc" "55,-116,0" "2.0" "0"
99.        LST3)
100.      (command "text" "s" "BG_ST" "j" "bc" "70,-116,0" "2.0" "0"
101.        "套管")
102.      (setq STR (read-line F1)) ;测试日期
103.      (command "text"  "s"      "BG_ST"  "j"
104.        "mc"  "-58.524,11.258,0"  "2.4"
105.        "0"  STR
106.       )
107.      (setq STR (repeat 2 (read-line F1))) ;钻孔标高
108.      (setq LST4 (strcat STR "m"))
109.      (command "text"  "s"      "BG_ST"  "j"
110.        "bl"  "103.604,32.308,0"  "2.5"
111.        "0"  LST4
112.       )
113.      (setq STR (repeat 140 (read-line F1))) ;测井数据152行
114.      (setq lst (strParse STR " "));转换成列表
115.      (setq data1 (car lst))
116.      (setq data2 (cadr lst))
117.      (setq data11 (atof data1))
118.      (setq data21 (atof data2))
119.      (setq Y (* data11 -2)) ;1:500计算Y坐标
120.      (if (< data21 min1)
121.        (setq X 90)
122.      )
123.      (if (> data21 max1)
124.        (setq X 130)
125.      )
126.      (if (and (> data21 min1) (< data21 max1))
127.        (setq X (+ (* (/ (- data21 min1) (- max1 min1)) 40) 90))
128.      ;min1<data21<max1时，X=[(data21-min1)/(max1-min1)]*40+90
129.      )
130.      (setq G1 (list X Y))
131.      (while (/= STR nil)
132.        (line-line)
133.        (command "layer" "s" "图层8" "")
134.        (command "pline" G1 G2 "w" "0.0" "0.0" "")
135.        (setq G1 G2)
136.      )
137.      
138.    )
139. )
140.       )
141.     )
142.   )
143.   (close A1)
144.   (princ "\nOK!")
145. )
146. ;;;将具有分隔符的字符串解析为列表
147. (defun strParse (Str Delimiter / SearchStr StringLen return n char)
148.   (setq SearchStr Str)
149.   (setq StringLen (strlen SearchStr))
150.   (setq return '())
151.   (while (> StringLen 0)
152.     (setq n 1)
153.     (setq char (substr SearchStr 1 1))
154.     (while (and (/= char Delimiter) (/= char ""))
155.       (setq n (1+ n))
156.       (setq char (substr SearchStr n 1))
157.     ) ;_ end of while
158.     (setq return (cons (substr SearchStr 1 (1- n)) return))
159.     (setq SearchStr (substr SearchStr (1+ n) StringLen))
160.     (setq StringLen (strlen SearchStr))
161.   ) ;_ end of while
162.   (reverse return)
163. ) ;_ end of defun

(defun C:WTCJ ()   (defun line-line ()     (setq STR (read-line F1))     (setq lst (strParse STR " ")) ;转换成列表     (setq data1 (car lst))     (setq data2 (cadr lst))     (setq data11 (atof data1))     (setq data21 (atof data2))     (setq Y (* data11 -2))     (if (< data21 min1)       (setq X 90)     )     (if (> data21 max1)       (setq X 130)     )     (if (and (> data21 min1) (< data21 max1))       (setq         X   (+ (* (/ (- data21 min1) (- max1 min1)) 40) 90)       )     )     (setq G2 (list X Y))   )      (command "osnap" "off")   (COMMAND "insert"     "C:\\Downloads\\AB线CD段EF隧道DZ-01-023孔\\AB线CD段EF隧道DZ-01-023孔.dwg"     "0,0,0"       "1"     "1"        "0"    )   (COMMAND "_explode" "all")   (COMMAND "_zoom" "e")   (princ "\n读取综合测井文件数据，绘制曲线。")   (if (setq FILE (getfiled "选择dat文件"       ""       "dat"       4    )       )     (progn       ;;读模式打开lsp文件       (setq F1 (open FILE "r"))       ;;逐行读入       (while (/= A1 nil) (setq STR (repeat 2 (read-line F1))) (if (= STR "自然电位")    (progn      (setq min1 (getreal "请输入自然电位最小值"))      (setq max1 (getreal "请输入自然电位最大值"))      (setq mid (strcat (rtos (setq mid1 (/ max1 4))) ;转成字符串          "  "          (rtos (setq mid2 (/ max1 2)))          "  "          (rtos (setq mid3 (* (/ max1 4) 3)))         )      )      (command "text"  "s"      "BG_ST"  "j"        "bc"  "110.012,2.375,0"  "2.5"        "0"  mid       )      (setq STR (repeat 4 (read-line F1))) ;钻孔编号      (setq LST STR)      (command "text"  "s"      "BG_ST"  "j"        "bl"  "-140.524,42.492,0"  "3.6"        "0"  LST       )   ;写图签 ;;;     (setq len (strlen STR)) ;;;     (setq LST1 (strcat (substr STR 1 1) "  ")) ;;;     (setq len (- len 1)) ;;;     (if (/= len 0) ;;;       (progn ;;;     ;                      (setq n 1) ;;;  (setq n (+ n 1)) ;;;  (setq LST1 (strcat LST1 (substr LST1 n 1) "  ")) ;;;  (setq len (- len 1)) ;;;       ) ;;;     )      (command "text"  "s"      "BG_ST"  "j"        "bl"  "90.960,48.634,0"  "5.0"        "0"  LST       )   ;写标题      (setq STR (repeat 3 (read-line F1))) ;套管长度      (setq LST2 (strcat STR "m"))      (command "text"  "s"      "BG_ST"  "j"        "bl"  "217.39,32.308,0"  "2.5"        "0"  LST2       )            (setq LST3 (* (atof  STR) 2)) ;1:500故乘2      (setq pa1 (polar '(30.0 0.0) (* pi 1.5) LST3))      (setq pa2 (polar '(80.0 0.0) (* pi 1.5) LST3))      (setq pa3 (polar '(32.382 -2.608) (* pi 1.5) (- LST3 2)))      (setq pa4 (polar '(32.154 -2.608) (* pi 1.5) (- LST3 2)))      (command "pline" '(32.382 -2.608) pa3 "")      (command "pline" '(32.154 -2.608) pa4 "")      (command "pline" pa1 pa2 "")      (command "text" "s" "BG_ST" "j" "bc" "55,-116,0" "2.0" "0"        LST3)      (command "text" "s" "BG_ST" "j" "bc" "70,-116,0" "2.0" "0"        "套管")      (setq STR (read-line F1)) ;测试日期      (command "text"  "s"      "BG_ST"  "j"        "mc"  "-58.524,11.258,0"  "2.4"        "0"  STR       )      (setq STR (repeat 2 (read-line F1))) ;钻孔标高      (setq LST4 (strcat STR "m"))      (command "text"  "s"      "BG_ST"  "j"        "bl"  "103.604,32.308,0"  "2.5"        "0"  LST4       )      (setq STR (repeat 140 (read-line F1))) ;测井数据152行      (setq lst (strParse STR " "));转换成列表      (setq data1 (car lst))      (setq data2 (cadr lst))      (setq data11 (atof data1))      (setq data21 (atof data2))      (setq Y (* data11 -2)) ;1:500计算Y坐标      (if (< data21 min1)        (setq X 90)      )      (if (> data21 max1)        (setq X 130)      )      (if (and (> data21 min1) (< data21 max1))        (setq X (+ (* (/ (- data21 min1) (- max1 min1)) 40) 90))      ;min1<data21<max1时，X=[(data21-min1)/(max1-min1)]*40+90      )      (setq G1 (list X Y))      (while (/= STR nil)        (line-line)        (command "layer" "s" "图层8" "")        (command "pline" G1 G2 "w" "0.0" "0.0" "")        (setq G1 G2)      )          ) )       )     )   )   (close A1)   (princ "\nOK!") ) ;;;将具有分隔符的字符串解析为列表 (defun strParse (Str Delimiter / SearchStr StringLen return n char)   (setq SearchStr Str)   (setq StringLen (strlen SearchStr))   (setq return '())   (while (> StringLen 0)     (setq n 1)     (setq char (substr SearchStr 1 1))     (while (and (/= char Delimiter) (/= char ""))       (setq n (1+ n))       (setq char (substr SearchStr n 1))     ) ;_ end of while     (setq return (cons (substr SearchStr 1 (1- n)) return))     (setq SearchStr (substr SearchStr (1+ n) StringLen))     (setq StringLen (strlen SearchStr))   ) ;_ end of while   (reverse return) ) ;_ end of defun

xianaihua

mogm127

谢谢楼上的，有空来长沙，请你喝茶。