Alle, entschuldigen Sie die Störung, aber ich bin ziemlich neu mit einer entscheidenden Schwierigkeit: Ich möchte eine animierte Karte von Russin mit Veränderungen der Arbeitslosigkeit mit unterschiedlichen Jahren erstellen, wie z. Auf dem Bild sehen Sie Daten für ein Jahr
require(sp)
require(maptools)
require(RColorBrewer)
require(rgdal)
rus<-url("http://www.filefactory.com/file/4h1hb5c1cw7r/n/RUS_adm1_RData")
print(load(rus))
unempl <- read.delim2(file="C:\\unempl1.txt", header = TRUE,
sep = ";",quote = "", dec=",", stringsAsFactors=F)
gadm_names <-gadm$NAME_1
total <- length(gadm_names)
pb <- txtProgressBar(min = 0, max = total, style = 3)
order <- vector()
for (i in 1:total){
order[i] <- agrep(gadm_names[i], unempl$region,
max.distance = 0.2)[1]
setTxtProgressBar(pb, i) # update progress bar
}
col_no <- as.factor(as.numeric(cut(unempl$data[order],
c(0,2.5,5,7.5,10,15,100))))
levels(col_no) <- c("<2,5%", "2,5-5%", "5-7,5%",
"7,5-10%", "10-15%", ">15%")
gadm$col_no <- col_no
myPalette<-brewer.pal(6,"Purples")
proj4.str <- CRS("+init=epsg:3413 +lon_0=105")
gadm.prj <- spTransform(gadm, proj4.str)
spplot(gadm.prj, "col_no", col=grey(.9), col.regions=myPalette,
main="Unemployment in Russia by region")
Das Ergebnis, das ich bereit bin zu bekommen, ist so etwas wie eine Animation hier: http://spatial.ly/2011/02/mapping-londons-population-change-2011-2030/ Ich habe jedoch viel gegoogelt und eine Reihe von Themen gelesen in http://stackoverflow.com, einschließlich der folgenden Informationen : Erstellen eines Films aus einer Reihe von Handlungen in R , aber immer noch nicht das Richtige.
Vielen Dank im Voraus!
Ich habe mir so etwas ausgedacht. Kann mir bitte jemand sagen, wo der Fehler liegt:
require(animation)
require(sp)
require(RColorBrewer)
require(classInt)
require(rgdal)
rus<-url("http://www.filefactory.com/file/4h1hb5c1cw7r/n/RUS_adm1_RData")
print(load(rus))
unempl1 <- read.delim2(file="C:\\unempl11.txt", header = TRUE,
sep = ";",quote = "", dec=",", stringsAsFactors=F)
unempl2<- read.delim2(file="C:\\unempl12.txt", header = TRUE,
sep = ";",quote = "", dec=",", stringsAsFactors=F)
gadm_names <-gadm$NAME_1
total <- length(gadm_names)
pb <- txtProgressBar(min = 0, max = total, style = 3)
order <- vector()
for (i in 1:total){
order[i] <- agrep(gadm_names[i], unempl1$region,
max.distance = 0.2)[1]
setTxtProgressBar(pb, i) # update progress bar
}
for (l in 1:total){
order[l] <- agrep(gadm_names[l], unempl2$region,
max.distance = 0.2)[1]
setTxtProgressBar(pb, i) # update progress bar
}
col_no_1 <- as.factor(as.numeric(cut(unempl1$data[order],
c(0,2.5,5,7.5,10,15,100))))
col_no_2<- as.factor(as.numeric(cut(unempl2$data[order],
c(0,2.5,5,7.5,10,15,100))))
saveHTML(
for(k in 1:2) {
try<-get(paste("col_no_", k, sep = ""))
levels(try) <- c("<2,5%", "2,5-5%", "5-7,5%",
"7,5-10%", "10-15%", ">15%")
gadm$col_no <- try
myPalette<-brewer.pal(6,"Purples")
proj4.str <- CRS("+init=epsg:3413 +lon_0=105")
gadm.prj <- spTransform(gadm, proj4.str)
spplot(gadm.prj, "col_no", col=grey(.9), col.regions=myPalette,
main="Unemployment in Russia by region")
},img.name = "map", htmlfile = "unrus2.html")
Betreff Die Bearbeitung: Was läuft mit dem Code falsch?
—
whuber
Da Ihr Beispiel nicht reproduzierbar ist, ist die Fehlerbehebung schwierig. Ein paar Dinge springen heraus 1) Sie wenden eine räumliche Transformation in einer Schleife an, also tun Sie dies wiederholt 2) Sie erstellen ein Objekt namens "try", das auch eine R-Funktion ist 3) Sie könnten durch tatsächliche Spaltennamen iterieren, d. H. ., für (i in c ("Var1", "Var2")) ist die Art und Weise, wie Sie es derzeit codiert haben, sehr kompliziert. 4) Ihr Aufruf von spplot ist nicht korrekt, Sie übergeben ihm einen Unsinnvektor.
—
Jeffrey Evans
Es tut mir wirklich leid, dass ich so ein Unverständnis habe, aber dies ist meine erste echte Erfahrung mit R überhaupt. Ich habe die Daten in der Hauptfrage hinzugefügt. Wenn es Sie nicht stört, können Sie bitte Verbesserungsmöglichkeiten vorschlagen, wie ich es wirklich bin lief von Ideen
—
Ruvin Rafailov
