library(ggmap)
library(mapproj)
map <- get_map(location = 'middle east', zoom = 4)
read.table("latlon.txt", header = TRUE) -> tbl
plot<-ggmap(map) + geom_point(aes(x = lon, y = lat), color=tbl$color, size=3, data = tbl) + theme(legend.position = "none")+
  #geom_text(aes(label=pop),data=tbl,hjust=0, vjust=0)
  ggsave(plot=plot,height=7.5,width=11, filename="f4 map.pdf", useDingbats=FALSE)
plot(plot)

map <- get_googlemap(center = 'middle east', zoom = 4,
                 style = 'feature:administrative.country|element:labels|visibility:off')
    plot <- ggmap(map)
    print(plot)