判断点是否在多边形内的函数

Lispboy

借鉴下老大上面发的论文，论文链接：关于点与简单多边形位置关系的一种算法
给大家写个测试点是否在多边形内的代码，
程序中用到了高飞鸟斑竹的矩阵函数。
用的是射线法。

[pcode=lisp,true]
;;;-----------------------------------------------------------;;
;;; 向量乘标量(系数)                                              ;;
;;; Vector x Scalar - Lee Mac                                      ;;
;;; Args: v - vector in R^n, s - real scalar                      ;;
;;;-----------------------------------------------------------;;
(defun MAT:vxs ( v s )
  (mapcar (function (lambda ( n ) (* n s))) v)
)
(defun GEO:Centroid (Pts / )
  (MAT:vxs (apply 'mapcar (cons '+ pts)) (/ 1.0 (length pts)))
)
;| 参数
   pts:多边形顶点集
   testpoint:测试点
   返回：T（内）NIL（外）
|;
(defun Points:PntIsIn (pts testpoint / centroid i_cnt item p_1 p_2 tmp vp1 x)
  (setq i_cnt nil
        centroid (geo:centroid pts)
        vp1 3e10
        tmp (angle testpoint centroid)
        vp1 (polar testpoint tmp vp1)
        i_cnt nil
  )

  (while (cdr pts)
    (setq p_1 (car pts)
          p_2 (cadr pts)
    )
    (if (setq tmp (inters
                    testpoint
                    vp1
                    p_1
                    p_2
                    't
                  )
        )
      (setq i_cnt (cons tmp i_cnt))
    )
    (setq pts (cdr pts))
  )
  (if i_cnt
    (progn
      (setq tmp (list (car i_cnt)))
      (foreach item (cdr i_cnt)
        (if (apply
              'and
              (mapcar
                '(lambda (x)
                   (not (equal x item 0.001))
                 )
                tmp
              )
            )
          (setq tmp (cons item tmp))
        )
      )
      (setq i_cnt tmp)
    )
  )
  (if (zerop (rem (length i_cnt) 2))
    nil
    t
  )
)
[/pcode]

ScmTools

建议LZ将此贴也发到“LISP/VLISP函数库”板块中, 这里还有另中方法: http://www.xdcad.net/forum/thread-667705-1-1.html

lianghui2018

谢谢楼主分享~~~

skg123

可以做出测试的代码，让大家直接测试

another2121

怎么用？:)