# ---------------------------------------
# check if points in polygon-----
points.natbag.csv$pointsIn<-point.in.polygon(points.natbag.csv$longitude,points.natbag.csv$latitude,
polygon.natbag.csv$Longitude,polygon.natbag.csv$Latitude)
# ---------------------------------------
# plot map----------------------------
ggmap(get_map(location = c(lon = mean(polygon.natbag.csv$Longitude), lat = mean(polygon.natbag.csv$Latitude)), zoom = 16, scale = 2, maptype='hybrid'))+
geom_polygon(data=polygon.natbag.csv,aes(x=Longitude,y=Latitude),alpha=0.2,colour="green",fill="green")+
geom_path(data=polygon.natbag.csv,aes(x=Longitude,y=Latitude),colour="green",alpha=0.7,size=0.5)+
geom_point(data = points.natbag.csv, aes(x = longitude, y = latitude,fill=pointsIn), size = 1, shape = 21, alpha = 1)+
guides(fill=FALSE, alpha=FALSE, size=FALSE)
# ---------------------------------------