[求助]：如何将地形线转化成坐标数据

zhang8755 · 发表于 2006-4-18 23:07:12

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？立即注册

×

断面图中，地形断面线由pline或line构成，如何获得每隔一米的高程，并制作成表。我现在的方法是每隔一米划一条竖线，然后求竖线于地形线的交点，这个方法较笨，速度也很慢。大家帮我想一想还有什么好的方法
如下图所示

xyp1964 · 发表于 2006-4-18 23:34:22

最初由 zhang8755 发布
[B]断面图中，地形断面线由pline或line构成，如何获得每隔一米的高程，并制作成表。我现在的方法是每隔一米划一条竖线，然后求竖线于地形线的交点，这个方法较笨，速度也很慢。大家帮我想一想还有什么好的方法
如下图?.. [/B]

竖线保留：

[FONT=courier new](load "xyp_lib.vlx") ;版本 V.20060413
;|下载和加载通用函数(可在签名栏直接下载后放到搜索路径下)
利用以下任何一种方式（首选第一种）即可加载和运行通用函数内的所有子程序：
★1·在acad.lsp中增加(load"xyp_lib.vlx")
■2·在每个程序内增加(load"xyp_lib.vlx")
■3·在command下，输入(load"xyp_lib.vlx")
■4·在菜单.mnl中增加(load"xyp_lib.vlx")
■5·将xyp_lib.vlx文件直接拽到cad屏幕
[COLOR=red] ★通用函数下载地址：[/COLOR]
[url]http://www.xdcad.net/forum/attachment.php?s=&postid=1606661[/url]
|;
(defun c:test ()
(CMDLA0)
(setvar "osmode" 0)
(setq s1 (car (entsel "\n选择: "))
ptlst (xyp-get-CurveIntersLeng s1 1)
)
(foreach pt ptlst
(xyp-zb pt 500 4);交点坐标标注
(princ "\n")
(princ (cadr pt));交点Y坐标
)
(CMDLA1)
)[/FONT]

zhang8755 · 发表于 2006-4-19 22:00:37

能不能不要竖线，构造一个点沿着线走，每个一段距离记录一组坐标，主要是我需要x坐标为等间距，现在想不出来实现的方法。

xyp1964 · 发表于 2006-4-19 23:53:09

作辅助线的方式：

[FONT=courier new]
(defun c:test ()
(CMDLA0)
(setq s1 (car (entsel "\n选择地形断面: "))
real (UREAL 7 "" "\nx方向间距" real)
pts (xyp-get-CurveStartPoint s1)
pt (getpoint pts "\n方向: ")
dx (- (car pt) (car pts))
ptnlst (LIST)
i -1
)
(if (>= dx 0)
(setq pt (polar pts 0 (* real (setq i (1+ i)))))
(setq pt (polar pts pi (* real (setq i (1+ i)))))
)
(command "line" pt (polar pt (* pi 0.5) real) "")
(setq s2 (entlast))
(while (setq ptn (xyp-get-Inters s1 s2 2))
(setq ptnlst (append ptnlst ptn))
(entdel s2)
(if (>= dx 0)
(setq pt (polar pts 0 (* real (setq i (1+ i)))))
(setq pt (polar pts pi (* real (setq i (1+ i)))))
)
(command "line" pt (polar pt (* pi 0.5) real) "")
(setq s2 (entlast))
)
(entdel s2)
(princ "\n高程点坐标: ")
(foreach pt ptnlst
(princ "\n")
(princ pt)
)
(CMDLA1)
)[/FONT]

[FONT=courier new] 
(defun c:test () 
  (CMDLA0) 
  (setq        s1     (car (entsel "\n选择地形断面: ")) 
        real   (UREAL 7 "" "\nx方向间距" real) 
        pts    (xyp-get-CurveStartPoint s1) 
        pt     (getpoint pts "\n方向: ") 
        dx     (- (car pt) (car pts)) 
        ptnlst (LIST) 
        i      -1 
  ) 
  (if (>= dx 0) 
    (setq pt (polar pts 0 (* real (setq i (1+ i))))) 
    (setq pt (polar pts pi (* real (setq i (1+ i))))) 
  ) 
  (command "line" pt (polar pt (* pi 0.5) real) "") 
  (setq s2 (entlast)) 
  (while (setq ptn (xyp-get-Inters s1 s2 2)) 
    (setq ptnlst (append ptnlst ptn)) 
    (entdel s2) 
    (if        (>= dx 0) 
      (setq pt (polar pts 0 (* real (setq i (1+ i))))) 
      (setq pt (polar pts pi (* real (setq i (1+ i))))) 
    ) 
    (command "line" pt (polar pt (* pi 0.5) real) "") 
    (setq s2 (entlast)) 
  ) 
  (entdel s2) 
  (princ "\n高程点坐标: ") 
  (foreach pt ptnlst 
    (princ "\n") 
    (princ pt) 
  ) 
  (CMDLA1) 
)[/FONT]

hbhu · 发表于 2006-4-20 01:24:12

附件不可以下载的..

xyp1964 · 发表于 2006-4-20 23:24:19

最初由 hbhu 发布
[B]附件不可以下载的.. [/B]

服务器正常可以下载了！

WHGF · 发表于 2006-4-21 07:46:02

可以采用线上加点的方式，每隔固定间隔加一个点，原来线上的点不要，然后输出线坐标。

zhang8755 · 发表于 2006-4-21 18:04:29

xyp版主的程序中
将(command "line" pt (polar pt (* pi 0.5) real) "")
改为
(entmake (list '(0 . "LINE")
(list 10 pt）
(list 11 (polar pt (* pi 0.5) real))
)
)
将
(setq ptnlst (append ptnlst ptn))
改为
(setq ptnlst (cons ptnlst ptn))
效率将大大提高，时间大约是原来的1/10。

xyp1964 · 发表于 2006-4-21 23:53:22

关于command和entmake的运算速度测试（画100000条线）：

[FONT=courier new](defun c:test ()
(CMDLA0)
;;command
(setq i -1)
(xyp-get-StartTime)
(repeat 100000
(setq pt (list (* 10 (setq i (1+ i))) 0))
(command "line" pt (polar pt (* pi 0.5) 200) "")
)
(xyp-get-UsedTime)
;;entmake
(setq i -1)
(xyp-get-StartTime)
(repeat 100000
(setq pt (list (* 10 (setq i (1+ i))) 300))
(entmake (list (cons 0 "LINE")
(cons 10 pt)
(cons 11 (polar pt (* pi 0.5) 200))
)
)
)
(xyp-get-UsedTime)
(CMDLA1)
)[/FONT]

Command: test

用时 : 53.06秒
用时 : 6.17秒

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

[求助] [求助]：如何将地形线转化成坐标数据

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

相关帖子

Re: [求助]：如何将地形线转化成坐标数据