Question

我有兴趣一起加入以下省份：伊拉克（Dohuk，Arbil，Kirkuk，As-Sulaymaniyah，Ardebil）以及以下来自土耳其的省份（Erzurum，Elazig，Diyarbakir，Malatya，Sanliurfa，Kars，Gaziantep）和伊朗（科德斯坦）。数据来自GDAM。

library("knitr")
library("pagenum")
library(zoom)
######################################################################################################

##Plot of the middle east
oceanCol = rgb(7,0,30,maxColorValue=255) 
cc<-c("Jordan","Syria","Lebanon","Israel","Palestine","Iraq","Iran","Georgia","Azerbaijan",
      "Turkey","Saudi Arabia","Kuwait","Bahrain","Qatar","United arab Emirates","yemen")

country<-countrycode(cc,"country.name","iso3c")

me<-vector("list",length(country))

for (i in 1:length(country)){ me[[i]]<-getData("GADM", country=country[i],level=1)}
map<-rbind(me[[1]],me[[2]],me[[3]],me[[4]],me[[5]],me[[6]],me[[7]],me[[8]],me[[9]],me[[10]],
            me[[11]],me[[12]],me[[13]],me[[14]],me[[15]],me[[16]])


plot(map,xlim=c(35,60),ylim=c(20,42), col="#ffffff",bg=oceanCol,border="black", 
     main ='State Creation and The Future Boundries Of The Middle East')

text(map, map$NAME_1, cex=0.3)

mtext("Source: TTTTTT July 2017",side=1,line=1)
legend("bottomright", legend=c("©TTTT"),col=white, cex=0.5, bty="n")
addnortharrow(pos="topright",scale=0.5,lwd=1, border="white",cols=c("white","white"), text.col="white")

pagenum(text="TTTT", date=TRUE,
  x=.95, y=.95, just =c('right','top'))

box()
#zm()

作者：王珂，黑龙江民族丛刊HEILONGJIANG

0 个答案: