Huidong Tian plot(0,0, type = "n", axes = FALSE, xlab = "", ylab = "") par(usr = c(Lon.lef, Lon.right, Lat.bottom, Lat.top)) map( , add = TRUE)