map.germany <- function(daten,
cutpoints=c(0.5,0.6,0.8,1.0,1.2,1.5,2.0), gray=F){
cutpoints <- c(-1e-5, cutpoints, 1e5)
if(gray==T)
farben <- gray((9-as.numeric(cut(daten,cutpoints)))/9)
else{
farben <- (9-as.numeric(cut(daten,cutpoints)))
palette(heat.colors(8))
}
xmin <- 1:length(germany)
xmax <- 1:length(germany)
ymin <- 1:length(germany)
ymax <- 1:length(germany)
for(i in 1:length(germany)){
xmin[i] <- min(germany[[i]][,2],na.rm=T)
xmax[i] <- max(germany[[i]][,2],na.rm=T)
ymin[i] <- min(germany[[i]][,3],na.rm=T)
ymax[i] <- max(germany[[i]][,3],na.rm=T)}
breite <- c(min(xmin),max(xmax))
hoehe <- c(min(ymin),max(ymax))
plot(breite,hoehe,type="n",axes=F, xlab=" ", ylab=" ")
for(k in length(germany):1){
polygon(germany[[k]][,2],germany[[k]][,3],col=farben[k])
}
if(gray==F){
for(i in 1:8){
polygon(c(5800,5800,6100,6100),c(500+350*(i-1),500+350*i,500+350*i,500+350*(i-1)),col=9-i)
}
}
else{
for(i in 1:8){
polygon(c(5800,5800,6100,6100),c(500+350*(i-1),500+350*i,500+350*i,500+350*(i-1)),col=gray((9-i)/9))
}
}
for(i in 2:8){
text(6375,500+350*(i-1),as.character(round(cutpoints[i], 2)),cex=.7,col=rgb(0,0,0))
}
}