
#maptoolsf(Windows)
#Qns65Έȏ㍂Ҕ䗦i2006Nj}iRvX}j`
#pdf("C:/GISdata/maptoolsdemo.pdf",family="Japan1")

library(maptools)

#ESRIS{s撬ShapeFileǂݍ
jpn <- read.shape("C:/GISdata/japan_ver61.shp") #ShapeFile=>MapObject

#DBFf[^̎o
jpndata <- jpn$att.data

## Machintosh
#PREF2<-iconv(jpn$att.data$PREF,"shift-jis","utf-8")
#CITY3<-iconv(jpn$att.data$CITY1,"shift-jis","utf-8")
#jpndata<-cbind(jpndata,PREF2)
#jpndata<-cbind(jpndata,CITY3)

#`p|Sf[^̎o
#polylistIuWFNgɕϊ
#polylistɂĂȂƃ|S̓h킯łȂ
jpnpoly <- Map2poly(jpn,region.id=attr(jpn,"region.id"))

#sS_擾
jpnxy <- get.Pcent(jpn)

#Qn̒o
gunmadata <- subset(jpndata,jpndata$PREF=="Qn")
gunmapoly <- subset(jpnpoly,jpndata$PREF=="Qn")

## Macintosh
#gunmadata <- subset(jpndata,jpndata$PREF2=="Qn")
#gunmapoly <- subset(jpnpoly,jpndata$PREF2=="Qn")

#QnsS_擾
#svbgp
x <- subset(jpnxy[,1],jpndata$PREF=="Qn")
y <- subset(jpnxy[,2],jpndata$PREF=="Qn")
## Macintosh
#x <- subset(jpnxy[,1],jpndata$PREF2=="Qn")
#y <- subset(jpnxy[,2],jpndata$PREF2=="Qn")


#Qns65Έȏ㍂Ҕ䗦f[^ǂݍ
aged <- read.delim("C:/GISdata/agedprop.txt")

#Machintosh
#aged <- read.delim("C:/GISdata/agedprop2.txt")

#DBFf[^ɍҔ䗦f[^}[WiJCODEL[ɂāj
gunmadata <- merge(aged,gunmadata,sort=F,by="JCODE")

#}pf[^쐬
DD <- gunmadata$AP2006
cutnum<-16
classes <- cut(DD,seq(min(DD),max(DD),length=cutnum+1),include.lowest=T)
cols <- rev(heat.colors(cutnum))
cols2 <- heat.colors(cutnum)

#n}`
par(cex=0.9)
plot(gunmapoly,col=cols[classes],xlab="",ylab="",axes=F,)
box()

#svbg
text(x+0.0002,y-0.0002,as.character(gunmadata$CITY1),cex=0.7,pos=1,offset=0,col="white")
text(x,y,as.character(gunmadata$CITY1),cex=0.7,pos=1,offset=0,col="black")

##Macintosh
#text(x+0.0002,y-0.0002,as.character(gunmadata$CITY3),cex=0.7,pos=1,offset=0,col="white")
#text(x,y,as.character(gunmadata$CITY3),cex=0.7,pos=1,offset=0,col="black")

#}
legend(max(x)-1.1*(max(x)-min(x)),min(y)+1.37*(max(y)-min(y)),
legend=names(rev(table(classes))),cex=0.7,fill=cols2,text.col='black')

#^Cg
title("Qns65Έȏ㍂Ҕ䗦i2006Nj","oT/QlFV`uRɂیÃf[^͉Kv, http://phi.ypu.jp/msb/")

#dev.off()
