明经CAD社区

 找回密码
 注册

QQ登录

只需一步,快速开始

搜索
查看: 2888|回复: 9

根据三角网 等高线切横断面

[复制链接]
发表于 2015-7-31 21:50:17 | 显示全部楼层 |阅读模式
论坛上的切断面 程序要么是只能处理直线段 要么处理的数据不能最详细反映地形特征(按距离等分了),有没有根据等高线或者三角网(处理等高线与横断面线 交点 或者三角网跟横断面线交叉点提取值转换为横断面数据的)同时纵断面线支持多段线(多段线含圆弧的)LSP。下面这个是根据三角网或者等高线与横断面线交点生成点位的
  1. ;|;;===========================================================================
  2.       通用函数                                                                ;
  3. 功能:求两个线条对象的交点                                                    ;
  4.       适用对象: Line、Circle、Arc、Ellipse、Polyline、                        ;
  5.       LWPolyline、3dPolyline、Spline                                          ;
  6. 参数:OBJ1   ----对象1                                                        ;
  7.       OBJ2   ----对象2                                                        ;
  8.       Extend ----延伸选项                                                     ;
  9.                  0  acExtendNone                                              ;
  10.                  1  acExtendThisEntity                                        ;
  11.                  2  acExtendOtherEntity                                       ;
  12.                  3  acExtendBoth                                              ;
  13.       ZZZ    ----输出选项                                                     ;
  14.                  "=0"  Z值取0                                                 ;
  15.                  "F1"  取第一个对象上的点                                     ;
  16.                  "F2"  取第二个对象上的点                                     ;
  17.                  "MAX" 取Z值大者                                              ;
  18.                  "MIN" 取Z值小者                                              ;
  19.       Fuzz   ----允许偏差值                                                   ;
  20. 返回:若成功,返回点位表;否则返回nil
  21. 日期:zml84 于2007-11-05
  22. ;;|;
  23. (vl-load-com)
  24. (defun ZL-GETINTERS1 (OBJ1   OBJ2   EXTEND ZZZ  FUZZ /      ENT1
  25.        ENT2   PT10   PT11   PT20  PT21 OBJ11  OBJ22
  26.        ARRAY  LST    LST_PT I  PT PT1    PT2
  27.        Z1     Z2
  28.       )
  29.     ;;0、对参数的格式化处理
  30.     (if (and (= (type EXTEND) 'INT)
  31.       (<= 0 EXTEND 3)
  32. )
  33. ()
  34. (setq EXTEND 0)
  35.     )
  36.     (setq ZZZ (strcase ZZZ))
  37.     ;;======================
  38.     ;;1、获取交点集合>>>>>>>
  39.     (if (and (= (vla-get-objectname OBJ1) "AcDbLine")
  40.       (= (vla-get-objectname OBJ2) "AcDbLine")
  41. )
  42. ;;对直线对象(line) 特别处理
  43. (progn
  44.      (setq ENT1 (entget (vlax-vla-object->ename OBJ1))
  45.     ENT2 (entget (vlax-vla-object->ename OBJ2))
  46.      )
  47.      (setq PT10 (assoc 10 ENT1)
  48.     PT11 (assoc 11 ENT1)
  49.     PT20 (assoc 10 ENT2)
  50.     PT21 (assoc 11 ENT2)
  51.      )
  52.      ;;去除Z坐标
  53.      (setq PT10 (list (cadr PT10) (caddr PT10))
  54.     PT11 (list (cadr PT11) (caddr PT11))
  55.     PT20 (list (cadr PT20) (caddr PT20))
  56.     PT21 (list (cadr PT21) (caddr PT21))
  57.      )
  58.      (setq LST (inters PT10 PT11 PT20 PT21 t))
  59.      (if LST
  60.   (setq LST (append LST '(0)))
  61.      )
  62. )
  63. (progn
  64.      ;;=====================
  65.      ;;复制实体
  66.      (setq OBJ11 (vla-copy OBJ1)
  67.     OBJ22 (vla-copy OBJ2)
  68.      )
  69.      ;;向xy平面投影,将Z坐标改为0
  70.      (TOXY OBJ11)
  71.      (TOXY OBJ22)
  72.      ;;获取交点集合
  73.      (setq ARRAY (vla-intersectwith OBJ11 OBJ22 EXTEND))
  74.      ;;删除复制后的对象
  75.      (vla-delete OBJ11)
  76.      (vla-delete OBJ22)
  77.      ;;由数组转换为表
  78.      (if (and ARRAY
  79.        (> (vlax-safearray-get-u-bound
  80.        (vlax-variant-value ARRAY)
  81.        1
  82.    )
  83.    1
  84.        )
  85.   )
  86.   (progn
  87.       (setq LST (vlax-safearray->list
  88.       (vlax-variant-value ARRAY)
  89.          )
  90.       )
  91.   )
  92.      )
  93. )
  94.     )
  95.     ;;======================
  96.     ;;2、分析整理>>>>>>>
  97.     (setq LST_PT '())
  98.     (if LST
  99. (progn
  100.      (setq I 0)
  101.      (repeat (/ (length LST) 3)
  102.   ;;2.1 获取当前点位
  103.   (setq PT (list (nth I LST)
  104.           (nth (+ 1 I) LST)
  105.           (nth (+ 2 I) LST)
  106.     )
  107.   )
  108.   ;;2.2 获取对象上对应点位
  109.   (setq PT1 (vlax-curve-getclosestpointtoprojection
  110.          OBJ1
  111.          PT
  112.          '(0 0 1)
  113.      )
  114.         PT2 (vlax-curve-getclosestpointtoprojection
  115.          OBJ2
  116.          PT
  117.          '(0 0 1)
  118.      )
  119.   )
  120.   (setq Z1 (caddr PT1)
  121.         Z2 (caddr PT2)
  122.   )

  123.   ;;2.3 效验偏差值
  124.   ;;就是说:过滤:参数中有偏差值选项,却不满足要求的点位
  125.   (if (and FUZZ
  126.     (or (= (type FUZZ) 'REAL)
  127.         (= (type FUZZ) 'INT)
  128.     )
  129.     (not (equal Z1 Z2 FUZZ))
  130.       )
  131.       ;; 空处理
  132.       ()
  133.       ;;2.4 对输出选项的处理
  134.       (progn
  135.    (cond
  136.        ((= ZZZ "F1")
  137.         (setq PT PT1)
  138.        )
  139.        ((= ZZZ "F2")
  140.         (setq PT PT2)
  141.        )
  142.        ((= ZZZ "MAX")
  143.         (if (> Z1 Z2)
  144.      (setq PT PT1)
  145.      (setq PT PT2)
  146.         )
  147.        )
  148.        ((= ZZZ "MIN")
  149.         (if (< Z1 Z2)
  150.      (setq PT PT1)
  151.      (setq PT PT2)
  152.         )
  153.        )
  154.        (t
  155.         (setq PT PT)
  156.        )
  157.    ) ;_结束cond
  158.    (if (member PT LST_PT)
  159.        ()
  160.        (setq LST_PT (cons PT LST_PT))
  161.    )
  162.       ) ;_结束progn
  163.   ) ;_结束if
  164.   (setq I (+ I 3))
  165.      ) ;_结束repeat
  166. ) ;_结束progn
  167.     ) ;_结束if
  168.     ;;3、返回结果>>>>>
  169.     LST_PT
  170. ) ;_结束defun
  171. ;;;============================================================
  172. ;;;功能:曲线实体上每个控制点的z坐标值置为0.0                  
  173. (defun TOXY (OBJ / NAME PT1 TP2)
  174.     ;;取得实体的类型名称
  175.     (setq NAME (vla-get-objectname OBJ))
  176.     (cond
  177. ;;类型1
  178. ;;直线(line)
  179. ((= NAME "AcDbLine")
  180.   ;;取得直线的起终点坐标
  181.   (setq PT1 (vlax-variant-value (vla-get-startpoint OBJ))
  182.         PT2 (vlax-variant-value (vla-get-endpoint OBJ))
  183.   )
  184.   ;;改变z值为0.0
  185.   (vlax-safearray-put-element PT1 2 0.0)
  186.   (vlax-safearray-put-element PT2 2 0.0)
  187.   (vla-put-startpoint OBJ PT1)
  188.   (vla-put-endpoint OBJ PT2)
  189. )
  190. ;;类型2
  191. ;;圆(circle)
  192. ;;圆弧(arc)
  193. ;;椭圆及椭圆弧(ellipse)
  194. ((or (= NAME "AcDbCircle")
  195.       (= NAME "AcDbArc")
  196.       (= NAME "AcDbEllipse")
  197.   )
  198.   ;;取得中心点座标
  199.   (setq PT1 (vlax-variant-value (vla-get-center OBJ)))
  200.   ;;改变中心点座标z值为0.0
  201.   (vlax-safearray-put-element PT1 2 0.0)
  202.   (vla-put-center OBJ PT1)
  203. )
  204. ;;类型3
  205. ;;多段线(polyline、lwpolyline)
  206. ;;拟合的2维多段线(polyline、lwpolyline)
  207. ((or (= NAME "AcDbPolyline")
  208.       (= NAME "AcDb2dPolyline")
  209.   )
  210.   ;;改变标高值为0.0
  211.   (vla-put-elevation OBJ 0.0)
  212. )
  213. ;;类型4
  214. ;;三维多段线(3dpolyline)  
  215. ((= NAME "AcDb3dPolyline")
  216.   ;;取得3维多段线的控制点
  217.   (setq PT1 (vlax-variant-value (vla-get-coordinates OBJ))
  218.         I   0
  219.   )
  220.   (repeat (/ (length (vlax-safearray->list PT1)) 3)
  221.       (vlax-safearray-put-element PT1 (+ I 2) 0.0)
  222.       (setq I (+ I 3))
  223.   )
  224.   (vla-put-coordinates OBJ PT1)
  225. )
  226. ;;类型5
  227. ;;样条曲线(Spline)
  228. ((= NAME "AcDbSpline")
  229.   ;;取得样条曲线的拟合点
  230.   ;;改变每个拟合点的z值为0.0
  231.   (setq PT1 (vlax-variant-value (vla-get-fitpoints OBJ))
  232.         I   0
  233.   )
  234.   (repeat (vla-get-numberoffitpoints OBJ)
  235.       (vlax-safearray-put-element PT1 (+ I 2) 0.0)
  236.       (setq I (+ I 3))
  237.   )
  238.   (vla-put-fitpoints OBJ PT1)
  239.   ;;取得样条曲线的控制点
  240.   ;;改变每个控制点的z值为0.0  
  241.   (setq
  242.       PT2 (vlax-variant-value (vla-get-controlpoints OBJ))
  243.       I  0
  244.   )
  245.   (repeat (vla-get-numberofcontrolpoints OBJ)
  246.       (vlax-safearray-put-element PT2 (+ I 2) 0.0)
  247.       (setq I (+ I 3))
  248.   )
  249.   (vla-put-controlpoints OBJ PT2)
  250. )
  251. (t NIL)
  252.     )
  253. ) ;_结束defun
  254. ;(setq ppzzxx (ZL-GETINTERS (vlax-ename->vla-object(car (entsel) ) ) (vlax-ename->vla-object(car (entsel) ) ) 0 "f2" nil)   )




  255. (DEFUN C:TT()
  256. (vl-load-com)
  257. (setq m_ent1 (car (entsel "\n请选择剖切线: ")))
  258.   (setq m_ent2 (car (entsel "\n请选择一条等高线: ")))
  259. (setq data(entget    m_ent2 ))
  260. (setq layer(assoc 8 data))
  261.   (setq 层名(CDR layer))

  262.   (SETQ SS (SSGET "X" ( list (cons 8  层名)  (cons 0  "LINE,ARC,SPLINE,LWPOLYLINE,POLYLINE"))))
  263. (SSSETFIRST SS)
  264. (setq ii 0)
  265. (repeat (sslength ss);循环选择集长度
  266.    (setq 单个图元   (ssname ss ii))
  267.   
  268. (setq ZL-GETINTERS  (ZL-GETINTERS1 (vlax-ename->vla-object m_ent1) (vlax-ename->vla-object 单个图元) 0 "f2" nil))
  269. ;(setq 高程 (CADDR (vlax-curve-getStartPoint 单个图元))
  270. ;;;      ps2(vlax-curve-getEndPoint m_ent2)
  271.       ;)
  272.    (foreach n ZL-GETINTERS
  273. ;;;     (print n)
  274. (setq xzb (car n))
  275. (setq yzb (cadr n))
  276. (setq 高程 (CADDR n))

  277.   (entmake (list (cons 0  "TEXT") (cons 1 (vl-princ-to-string 高程)) (cons 10 (list  XZB YZB 高程))
  278.                (cons 40 0.5)
  279.                (cons 8 "0000用于图元属性快速输出利用(简称“属性”图层)")
  280.                ))
  281.                    ; 在下面画圆,这样就可以直接用“”读出里程及偏移,而孔号就是该点高程
  282.                    (entmake (list '(0 . "POINT") (cons 10 (list  XZB YZB 高程)) (cons 8 "0000用于图元属性快速输出利用(简称“属性”图层)")))
  283. ;(command "break" 单个图元  (list  XZB YZB 高程) "@")
  284.      ;(command "break" 单个图元 "f" (list  XZB YZB 高程) (list  XZB YZB 高程))
  285.                     ; (command "LINE"  (LIST XZB YZB) '(0 0))

  286.                   

  287.     )

  288.    
  289.   (setq ii (1+ ii))
  290. )

  291.   
  292.   
  293. ;;;  (sssetfirst nil ss)
  294.   )


  295. (defun m_IntersectWith(m_ent1 单个图元 / m_obj1 m_obj2 m_objcopy1 m_objcopy2 m_jdtab m_jdtab1 i)
  296. ; 来源:3楼 [求助]任意两条线的交点坐标-AutoLISP/Visual LISP 编程技术-CAD论坛-明经CAD社区 - Powered by Discuz!
  297. ; http://bbs.mjtd.com/forum.php?mod=viewthread&tid=79868
  298. ;;适用对象: Line、Polyline、LWPolyline、Circle、Arc、Ellipse、3dPolyline、Spline
  299.   ;;支持求空间虚交点,但Z坐标始终为0.0,要求Z坐标,请用(vlax-curve-getClosestPointToProjection)函数
  300.   (setq m_obj1 (vlax-ename->vla-object m_ent1))
  301.   (setq m_obj2 (vlax-ename->vla-object 单个图元))
  302.   
  303.   (setq m_objcopy1 (vla-copy m_obj1));;复制第一条曲线实体
  304.   (setq m_objcopy2 (vla-copy m_obj2));;复制第二条曲线实体
  305.   
  306.   (setq m_objcopy1 (m_ShadowToXY m_objcopy1));;得到投影实体
  307.   (setq m_objcopy2 (m_ShadowToXY m_objcopy2));;得到投影实体
  308.   
  309.   (setq m_jdtab1 (vla-intersectwith m_objcopy1 m_objcopy2 acExtendnone));;得到交点集
  310.   
  311.   (if (> (vlax-safearray-get-u-bound (vlax-variant-value m_jdtab1) 1) 1);;判断有无交点
  312.       (progn
  313. (setq m_jdtab1 (vlax-safearray->list (vlax-variant-value m_jdtab1)));;safearray数组转换为list表
  314. (setq i 0)
  315. (repeat (/ (length m_jdtab1) 3)
  316.    (setq m_jd (list (nth i m_jdtab1) (nth (+ 1 i) m_jdtab1) (nth (+ 2 i) m_jdtab1)));;取得一个交点
  317.    (setq m_jdtab (cons m_jd m_jdtab));;构造交点表((第一个交点) (第二个交点)。。。)
  318.    (setq i (+ 3 i))
  319. )
  320.       )
  321.       (princ"\n两曲线无交点!")
  322.   )
  323.   
  324.   (vla-delete m_objcopy1);;删除复制的第一条曲线实体
  325.   (vla-delete m_objcopy2);;删除复制的第二条曲线实体
  326.   
  327.   (setq m_jdtab m_jdtab);;返回交点表,无交点返回nil
  328. )
  329. (defun m_ShadowToXY(m_obj / m_objname m_pts m_pts1 i)
  330.   ;;对曲线实体m_obj创建一个投影至xy平面的曲线实体,即对曲线实体上每个控制点的z坐标值置为0.0
  331.   ;;输入曲线实体(vla对象)
  332.   ;;返回投影实体(vla对象)
  333.   (setq m_objname (vla-get-objectname m_obj));;取得实体的类型名称
  334. ;;;  (m_princ "\nObjectName:" m_objname)
  335.   (cond
  336.     ((= "AcDbSpline" m_objname);;样条曲线(Spline)
  337.      (setq i 0)
  338.      (setq m_pts (vlax-variant-value (vla-get-fitpoints m_obj)));;取得样条曲线的拟合点
  339.      (setq m_pts1 (vlax-variant-value (vla-get-controlpoints m_obj)));;取得样条曲线的控制点
  340.      (repeat (vla-get-numberoffitpoints m_obj);;循环
  341.        (vlax-safearray-put-element m_pts (+ i 2) 0.0);;改变每个拟合点的z值为0.0
  342.        (setq i (+ i 3))
  343.      )
  344.      (vla-put-fitpoints m_obj m_pts);;更改曲线拟合点属性
  345.      (setq i 0)
  346.      
  347.      (repeat (vla-get-numberofcontrolpoints m_obj);;循环
  348.        (vlax-safearray-put-element m_pts1 (+ i 2) 0.0);;改变每个控制点的z值为0.0
  349.        (setq i (+ i 3))
  350.      )
  351.      (vla-put-controlpoints m_obj m_pts1);;更改曲线控制点属性
  352.     )
  353.    
  354.     ((= "AcDb3dPolyline" m_objname);;三维多段线(3dpolyline)
  355.      (setq i 0)
  356.      (setq m_pts (vlax-variant-value (vla-get-coordinates m_obj)));;取得3维多段线的控制点
  357.      (repeat (/ (length (vlax-safearray->list m_pts)) 3)
  358.        (vlax-safearray-put-element m_pts (+ i 2) 0.0)
  359.        (setq i (+ i 3))
  360.      )
  361.      (vla-put-coordinates m_obj m_pts)
  362.     )
  363.    
  364.     ((= "AcDbLine" m_objname);;直线(line)
  365.      (setq i 0)
  366.      (setq m_pts (vlax-variant-value (vla-get-startpoint m_obj)));;取得直线的起点座标
  367.      (setq m_pts1 (vlax-variant-value (vla-get-endpoint m_obj)));;取得直线的端点座标
  368.      (vlax-safearray-put-element m_pts 2 0.0);;改变起点座标z值为0.0
  369.      (vlax-safearray-put-element m_pts1 2 0.0)
  370.      (vla-put-startpoint m_obj m_pts)
  371.      (vla-put-endpoint m_obj m_pts1)
  372.     )
  373.     ((or (= "AcDbCircle" m_objname);;园(circle)
  374.   (= "AcDbArc" m_objname);;圆弧(arc)
  375.   (= "AcDbEllipse" m_objname);;椭圆及椭圆弧(ellipse)
  376.      )
  377.      (setq m_pts (vlax-variant-value (vla-get-center m_obj)));;取得中心点座标
  378.      (vlax-safearray-put-element m_pts 2 0.0);;改变中心点座标z值为0.0
  379.      (vla-put-center m_obj m_pts)
  380.     )
  381.    
  382.     ((or (= "AcDbPolyline" m_objname);;多段线(polyline、lwpolyline)
  383.   (= "AcDb2dPolyline" m_objname);;拟合的2维多段线(polyline、lwpolyline)
  384.      )
  385.      (vla-put-elevation m_obj 0.0);;改变标高值为0.0
  386.     )
  387.   )
  388.   (setq m_obj m_obj)
  389. )
压缩文件那个EXE能处理纵断面圆弧 和直线 但要分开输入

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

x
发表于 2015-7-31 22:17:12 | 显示全部楼层
支持测绘界人才·~~
发表于 2015-7-31 23:23:42 | 显示全部楼层
依据等高线,交点,
或者依据(平面图中的高程点→三角网),交点,
或者是,等高线+三角网

这里面会有矛盾的地方,(同一个点,针对等高线和三角网的数据应该不一样,有时候或者偏差很大。)
这种情况怎么解决呢?
发表于 2015-8-2 12:13:16 | 显示全部楼层
给你个思路,选择等高线ss,选择剖面线ent,求出ss与ent的交点集pts,将pts按照ent方向排序,随后输出~~~
发表于 2015-8-8 20:23:49 | 显示全部楼层
我根据 ZmL84的思路做了一个用于公路上剖断面的,http://bbs.mjtd.com/thread-110548-1-1.html
发表于 2015-8-9 11:59:30 | 显示全部楼层
,测绘届的人才啊
发表于 2016-2-18 16:30:08 | 显示全部楼层
感谢楼主,看下好用不!!
发表于 2016-2-19 08:46:44 | 显示全部楼层
谢谢分享………………
发表于 2016-2-20 19:54:33 | 显示全部楼层
ding!!!!!!!!!!!!!!!!!!1
发表于 2022-9-28 13:19:04 | 显示全部楼层
支持一下顶一个
您需要登录后才可以回帖 登录 | 注册

本版积分规则

小黑屋|手机版|CAD论坛|CAD教程|CAD下载|联系我们|关于明经|明经通道 ( 粤ICP备05003914号 )  
©2000-2023 明经通道 版权所有 本站代码,在未取得本站及作者授权的情况下,不得用于商业用途

GMT+8, 2024-12-22 23:41 , Processed in 0.506472 second(s), 27 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表