首先設定好了文件的路徑及要使用的統計包呵~
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)
#上述的都是加上地區名字的命令
效果也是不錯的,
沒有留言:
張貼留言