绘制地图

在R中画地图先从简单的maps包开始。

library("maps")

在这个maps包中有一些数据集，用命令data(package=”maps”)，可以看到如下数据：

canada.cities          Database of Canadian cities 
county.fips            FIPS county codes for US County Map 
countyMapEnv           United States County Map 
franceMapEnv           France Map 
italyMapEnv            Italy Map 
nzMapEnv               New Zealand Basic Map 
ozone                  Sample datasets 
state.carto.center     United States State Population Cartogram Map 
state.cartoMapEnv      United States State Population Cartogram Map 
state.fips             FIPS state codes for US 48 State Map 
state.vbm.center       United States State Visibility Base Map 
state.vbmMapEnv        United States State Visibility Base Map 
stateMapEnv            United States State Boundaries Map 
unemp                  Sample datasets 
us.cities              Database of US cities 
usaMapEnv              United States Coast Map 
votes.repub            Sample datasets 
world.cities           Database of world cities 
world2MapEnv           Pacific Centric Low resolution World Map 
worldMapEnv            Low resolution World Map

 

画出美国各州的边界：

map('state', panel.first=grid())

axis(1,lwd=0)

axis(2,lwd=0)

axis(3,lwd=0)

axis(4,lwd=0)

box()

标出美国的城市：

data(us.cities)

points(us.cities$long, us.cities$lat, pch=19)

标出人口大于一百万的城市

cities <- us.cities[us.cities$pop > 1000000,]

points(cities$long, cities$lat, pch=19)

text(cities$long, cities$lat, cities$name, col="blue", cex=0.8, pos=1)

 

强大的ggmap中的get_map()可以方便获得各种地理或遥感图像数据，可惜现在的google被屏蔽，试不了了。

下面这个是用VPN后调用get_map取得的美国地图图像。

 

另外也可以用OpenStreetMap包。

library("OpenStreetMap")

map <- openmap(c(50,-130), c(20,-60),type="apple-iphoto")

plot(map)