2015年9月25日 星期五

在R語言內劃中國地圖 (Drawing China map in R statistical package)

首先設定好了文件的路徑及要使用的統計包呵~

setwd("路徑, 記得文件夾之間是 / ")

#install.packages("maptools") 要使用這個統計包.
library(maptools)

#install.packages("mapproj")
library(mapproj)

mydat<-readShapePoly("文件名")
plot(mydat) #其實只要輸入這個命令, 就會出到地圖了!

library(ggplot2) #但想圖靚些... 就要用ggplot2統計包

mymap<-ggplot(data=fortify(mydat))+
  geom_polygon(aes(x=long,y=lat,group=id),col="black",fill=I("green"))+
  theme_grey()
print(mymap+coord_map())
print(mymap+coord_map(projection="azequidistant"))

#fortify() 是要將shp地圖文件, 轉回普通的文件啊.
#geom_polygon()是專用來劃地圖文件的, 其中的x是經度, y是緯度, col是地圖邊界的顏色, fill是地圖的底色
#projection="azequidistant" 是使地圖呈球面顯現...


temp<-coordinates(mydat) #加入地圖的經緯度
temp=as.data.frame(temp) #將地圖轉為普通文件格式
temp$NAME=iconv(mydat$name_cn) #在普通文件內, 對應經緯度加入地區名稱
mymap+geom_text(aes(x=V1,y=V2,label=NAME),data=temp)

#上述的都是加上地區名字的命令

    效果也是不錯的,



沒有留言:

張貼留言