如何计算曲线交点？

晨曦CAD

给定输入一个曲线图元curve（*LINE，ARC均可）和两个点pp1、pp2，

假设过两点的线段与曲线有交点，大佬有没有什么好的办法计算出交点？

附一个我现有的方案:
弊端：需要创建一个LINE图元

2. (defun lineinsect(curve pp1 pp2 / obj1 obj2 objen intp xintp closedpoint)
   
3.     ;closedpoint返回计算在点列表pbll中距离点pa最近的点
   
4.     (defun closedpoint(pa pbll / padistlst cx)
   
5.       (setq padistlst(mapcar'(lambda(cx)(distance pa cx))pbll))
   
6.       (nth (car(HG::LIST::INDEX (HG::LIST::MIN padistlst)padistlst)) pbll)
   
7.       )
   
8.       
   
9.     (setq obj1(vlax-ename->vla-object curve))
   
10.     ；创建pp1、pp2为起点终点的LINE
    
11.     (setq objen(HG::ACT::PUT_ENAME "LINE" (LIST (trans2wcs pp1) (trans2wcs pp2))))
    
12.     (setq obj2(vlax-ename->vla-object objen))
    
13.     (setq intp(vla-intersectwith obj1 obj2 acExtendNone))
    
14.     (ENTDEL objen)
    
15.     (setq xintp(vlax-safearray->list (vlax-variant-value intp)))
    
16.     (if(> (length xintp)3)
    
17.       (trans2ucs(closedpoint pp1 (HG::LIST::DIVIDE3 xintp)));辅助函数HG::LIST::DIVIDE3将列表三三分割
    
18.       (trans2ucs xintp);辅助函数trans2ucs将点转化为UCS坐标系
    
19.     )
    
20.   )

(defun lineinsect(curve pp1 pp2 / obj1 obj2 objen intp xintp closedpoint)<br />     ;closedpoint返回计算在点列表pbll中距离点pa最近的点<br />     (defun closedpoint(pa pbll / padistlst cx)<br />       (setq padistlst(mapcar'(lambda(cx)(distance pa cx))pbll))<br />       (nth (car(HG::LIST::INDEX (HG::LIST::MIN padistlst)padistlst)) pbll)<br />       )<br />       <br />     (setq obj1(vlax-ename->vla-object curve))<br />     ；创建pp1、pp2为起点终点的LINE<br />     (setq objen(HG::ACT::PUT_ENAME "LINE" (LIST (trans2wcs pp1) (trans2wcs pp2))))<br />     (setq obj2(vlax-ename->vla-object objen))<br />     (setq intp(vla-intersectwith obj1 obj2 acExtendNone))<br />     (ENTDEL objen)<br />     (setq xintp(vlax-safearray->list (vlax-variant-value intp)))<br />     (if(> (length xintp)3)<br />       (trans2ucs(closedpoint pp1 (HG::LIST::DIVIDE3 xintp)));辅助函数HG::LIST::DIVIDE3将列表三三分割<br />       (trans2ucs xintp);辅助函数trans2ucs将点转化为UCS坐标系<br />     )<br />   )

liuhe

http://bbs.mjtd.com/forum.php?mo ... hlight=%BD%BB%B5%E3    大哥，有专门的函数求交点Intersectwith

晨曦CAD

liuhe 发表于 2023-1-4 14:20
http://bbs.mjtd.com/forum.php?mo ... hlight=%BD%BB%B5%E3    大哥，有专门的函数求交点Intersectwith

我给的代码例子是用Intersectwith的，有没有不用Intersectwith的方法，因为只有一条曲线图元，不能直接满足Intersectwith的输入参数要求

飞雪神光

获取 曲线图元 点表 用inters 循环求交点

wzg356

飞雪神光 发表于 2023-1-4 17:24
获取 曲线图元 点表 用inters 循环求交点

真曲线如圆弧是不行的

飞雪神光

wzg356 发表于 2023-1-4 17:36
真曲线如圆弧是不行的

根据精度需求等分曲线 但是速度会慢 还有ssnamex 这个也可以取得交点 也是精度上也有问题

20060510412

该用辅助线的时候，就得用辅助线。

纯算法的效率，有可能还不如使用辅助线。

xyp1964

1. ;; (lineinsect (setq s0(car(entsel"\n选择: "))) (setq p0 (getpoint  "\n基点<退出>: ")) (setq p1 (getpoint p0 "\n基点<退出>: ")))
   
2. (defun lineinsect (curve p1 p2 / LST S1)   
   
3.   (defun ListDiv (lst nn / a tmp)
   
4.     ;; (ListDiv '(1 2 3 4 5 6 7 8 9 ) 3)
   
5.     (while (>= (length lst) nn)
   
6.       (setq tmp '())
   
7.       (repeat nn (setq tmp (cons (car lst) tmp) lst (cdr lst)))
   
8.       (setq a (cons (reverse tmp) a))
   
9.     )
   
10.     (if lst (setq a (cons lst a)))
    
11.     (reverse a)
    
12.   )
    
13.   (command "line" "non" p1 "non" p2 "")
    
14.   (setq lst (vlax-invoke(vlax-ename->vla-object curve)'intersectwith(vlax-ename->vla-object (setq s1 (entlast)))0))
    
15.   (ENTDEL s1)
    
16.   (ListDiv lst 3)
    
17. )
    

;; (lineinsect (setq s0(car(entsel"\n选择: "))) (setq p0 (getpoint  "\n基点<退出>: ")) (setq p1 (getpoint p0 "\n基点<退出>: ")))<br /> (defun lineinsect (curve p1 p2 / LST S1)    <br />   (defun ListDiv (lst nn / a tmp)<br />     ;; (ListDiv '(1 2 3 4 5 6 7 8 9 ) 3)<br />     (while (>= (length lst) nn)<br />       (setq tmp '())<br />       (repeat nn (setq tmp (cons (car lst) tmp) lst (cdr lst)))<br />       (setq a (cons (reverse tmp) a))<br />     )<br />     (if lst (setq a (cons lst a)))<br />     (reverse a)<br />   )<br />   (command "line" "non" p1 "non" p2 "")<br />   (setq lst (vlax-invoke(vlax-ename->vla-object curve)'intersectwith(vlax-ename->vla-object (setq s1 (entlast)))0))<br />   (ENTDEL s1)<br />   (ListDiv lst 3)<br /> )<br />

magicheno

xyp1964 发表于 2023-1-7 15:16

感谢院长分享~