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