powered by
shadeborder(reg, lat, lon, col = 0, col.names = c("lon", "lat"))
lat
lon
lines
geocontour
Proj