R/inpoly.R

Defines functions `inpoly`

`inpoly` <-
  function(x,y, POK )
  {
    ##  
    kin = fields::in.poly(cbind(x, y) ,cbind(POK$x, y =POK$y) )

    G = rep(0,length(x))
    G[kin] = 1

    return(G)
  }

Try the GEOmap package in your browser

Any scripts or data that you put into this service are public.

GEOmap documentation built on Sept. 1, 2023, 5:09 p.m.