1options("rgdal_show_exportToProj4_warnings"="none") 2library(sp) 3 4# open polygon: 5print(point.in.polygon(1:10,1:10,c(3,5,5,3),c(3,3,5,5))) 6# closed polygon: 7print(point.in.polygon(1:10,rep(4,10),c(3,5,5,3,3),c(3,3,5,5,3))) 8