判断逆顺时针的方法

Lisphk

已知三点，有什么好的方法判断时针方向吗？

newer

剖砖引玉下，用向量的叉积判断

p2-p1 的向量 叉积 p3-p2的向量，结果>0，逆时针，否则顺时针。

q3_2006

newer 发表于 2014-11-25 12:14
剖砖引玉下，用向量的叉积判断

p2-p1 的向量 叉积 p3-p2的向量，结果>0，逆时针，否则顺时针。

虽然不懂...但用着挺爽....{:soso_e128:}

st788796

用trans 转换也是一个道理

Underway

http://4d-technologies.com/techcenter/
--->Geometry--->
GE_WhatPoly   Find the direction of the polygon (if Cw or CCw)

Lispboy

贴一个纯数学的方法，求N边形面积的（不只是三角形），正值逆时针，负值顺指针

1. (defun XD::Pnts:Area ( pts / pts1 iSum p0 p1 p2 x0 x1 x2 y0 y1 y2) )
   
2.   (setq pts1 pts)
   
3.   (setq iSum 0)
   
4.   (setq p0 (car pts)
   
5.         x0 (car p0)
   
6.         y0 (cadr p0)
   
7.   )
   
8.   (while (cdr pts1)
   
9.     (setq p1 (car pts1)
   
10.           x1 (car p1)
    
11.           y1 (cadr p1)
    
12.           p2 (cadr pts1)
    
13.           y2 (cadr p2)
    
14.           x2 (car p2)
    
15.           iSum (+ iSum (- (* x1 y2) (* x2 y1)))
    
16.           pts1 (cdr pts1)
    
17.     )
    
18.   )
    
19.   (/ (+ iSum (- (* (caar pts1) y0) (* x0 (cadar pts1)))) 2.0)
    
20. )
    

(defun XD::Pnts:Area ( pts / pts1 iSum p0 p1 p2 x0 x1 x2 y0 y1 y2) )<br />   (setq pts1 pts)<br />   (setq iSum 0)<br />   (setq p0 (car pts)<br />         x0 (car p0)<br />         y0 (cadr p0)<br />   )<br />   (while (cdr pts1)<br />     (setq p1 (car pts1)<br />           x1 (car p1)<br />           y1 (cadr p1)<br />           p2 (cadr pts1)<br />           y2 (cadr p2)<br />           x2 (car p2)<br />           iSum (+ iSum (- (* x1 y2) (* x2 y1)))<br />           pts1 (cdr pts1)<br />     )<br />   )<br />   (/ (+ iSum (- (* (caar pts1) y0) (* x0 (cadar pts1)))) 2.0)<br /> )<br />

命令: (XD::Pnts:Area (list p1 p2 p3))
196154.0

命令:
命令: (XD::Pnts:Area (list p1 p3 p2))
-196154.0

st788796

Lispboy 发表于 2014-11-25 19:51
贴一个纯数学的方法，求N边形面积的（不只是三角形），正值逆时针，负值顺指针

如果一个闭合pline圆会是什么情况呢:)

Lispboy

st788796 发表于 2014-11-25 19:56
如果一个闭合pline圆会是什么情况呢

正常，它主要看你给的这些点的顺序。

Lispboy

再提供个方法：

XDRX_API，现成的

命令: (xdrx_points_area p1 p2 p3)
196154.0
命令: (xdrx_points_area p1 p3 p2)
-196154.0

Lispboy

方法还有，使用ARX GE库，直接构造AcGeCircArc2d，用它下面的isClockWise()方法。代码如下：

1. (defun XD::Pnts:IsClockWise (p1 p2 p3 / ge tf)
   
2.   (if (setq ge (xdge::constructor "kCircArc2d" p1 p2 p3))
   
3.     (progn
   
4.       (setq tf (xdge::getpropertyvalue ge "IsClockWise"))
   
5.       (xdge::free ge)
   
6.     )
   
7.   )
   
8.   tf
   
9. )
   

(defun XD::Pnts:IsClockWise (p1 p2 p3 / ge tf)<br />   (if (setq ge (xdge::constructor "kCircArc2d" p1 p2 p3))<br />     (progn<br />       (setq tf (xdge::getpropertyvalue ge "IsClockWise"))<br />       (xdge::free ge)<br />     )<br />   )<br />   tf<br /> )<br />

命令: (XD::PNTS:ISCLOCKWISE p1 p2 p3)
nil

命令: (XD::PNTS:ISCLOCKWISE p1 p3 p2)
T

Lispboy

我发现好多方法啊，如泉涌啊，晓东通用LISP函数库已经有现成的了：

命令: (XD::Pnts:IsClockWise (list p1 p3 p2))
T

Lispboy

还没完，用3点求凸度的方法，判断，+逆时针

1. ;|
   
2. 3点求凸度方法判断逆顺时针，+ 逆时针
   
3. |;
   
4. (defun XD::Pnts:IsClockWise (p1 p2 p3)
   
5.   (setq an (/ (- (angle p2 p3) (angle p1 p2)) 2.0))
   
6.   (< (/ (sin an)(cos an)) 0)
   
7. )
   

;|<br /> 3点求凸度方法判断逆顺时针，+ 逆时针<br /> |;<br /> (defun XD::Pnts:IsClockWise (p1 p2 p3)<br />   (setq an (/ (- (angle p2 p3) (angle p1 p2)) 2.0))<br />   (< (/ (sin an)(cos an)) 0)<br /> )<br />

命令: (XD::PNTS:ISCLOCKWISE p1 p3 p2)
T
命令: (XD::PNTS:ISCLOCKWISE p1 p2 p3)
nil

Lispboy

还有啊，3点求平面法线，Z值+ 逆时针

1. ;|
   
2. 3点求平面法线方法判断，Z值+ 逆时针
   
3. |;
   
4. (defun XD::Pnts:IsClockWise (p1 p2 p3)
   
5.   (setq nrm (xdrx_points_normal p1 p2 p3))
   
6.   (< (last nrm) 0)
   
7. )
   

;|<br /> 3点求平面法线方法判断，Z值+ 逆时针<br /> |;<br /> (defun XD::Pnts:IsClockWise (p1 p2 p3)<br />   (setq nrm (xdrx_points_normal p1 p2 p3))<br />   (< (last nrm) 0)<br /> )<br />

命令: (XD::PNTS:ISCLOCKWISE p1 p2 p3)
nil
命令: (XD::PNTS:ISCLOCKWISE p1 p3 p2)
T

Lispboy

上面用法线法的用ARX 几何库的实现方法

1. (defun XD::Pnts:IsClockWise (p1 p2 p3)
   
2.   (if (setq ge (xdge::constructor "kPlane" p2 p1 p3))   ;;构造AcGePlane几何实体
   
3.      (progn
   
4.        (setq nrm (xdge::getpropertyvalue ge "normal")) ;;Normal方法
   
5.        (xdge::free ge)
   
6.        (< (last nrm) 0.0)
   
7.      )
   
8.   )
   
9. )
   
10. 
    

(defun XD::Pnts:IsClockWise (p1 p2 p3)<br />   (if (setq ge (xdge::constructor "kPlane" p2 p1 p3))   ;;构造AcGePlane几何实体<br />      (progn<br />        (setq nrm (xdge::getpropertyvalue ge "normal")) ;;Normal方法<br />        (xdge::free ge)<br />        (< (last nrm) 0.0)<br />      )<br />   )<br /> )<br /> <br />

csharp

Lispboy 发表于 2014-11-25 20:06
正常，它主要看你给的这些点的顺序。

有点问题，绘制单段ARC的 Pline，然后用 pedit -> close 闭合，这是查看顶点数为 2 ！这个 pline 是顺时针还是逆时针呢？