(XD::Pnts:Area)纯数学方法求N边形面积（符号判断顺逆时针）

Lispboy

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. )
    
21. 
    

(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 /> <br />

st788796

这个早发过的

Lispboy

st788796 发表于 2014-11-25 19:59
这个早发过的

发过了吗，我也感觉面熟，老了{:soso_e127:}