tout le monde, désolé pour le dérangement, mais je suis tout à fait nouveau avec r face à une difficulté cruciale: je veux créer une carte animée de Russin avec des changements de chômage avec différentes années, comme. Sur l'image, vous pouvez voir les données d'un an
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")
Le résultat que je souhaite obtenir est quelque chose comme une animation ici: http://spatial.ly/2011/02/mapping-londons-population-change-2011-2030/ Cependant, j'ai beaucoup cherché sur Google, lu un certain nombre de thèmes dans http://stackoverflow.com, y compris les éléments suivants: Création d'un film à partir d'une série de tracés dans R , mais ne pouvait toujours pas faire la bonne chose.
Merci d'avance!
Je suis venu avec quelque chose comme ça, quelqu'un peut-il me dire où est l'erreur:
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")
Re The edit: qu'est-ce qui ne va pas avec le code?
—
whuber
Étant donné que votre exemple n'est pas reproductible, il est difficile à résoudre. Quelques choses sautent aux yeux 1) vous appliquez une transformation spatiale dans une boucle, vous le faites donc à plusieurs reprises 2) vous créez un objet appelé "try" qui est également une fonction R 3) vous pouvez parcourir les noms de colonnes réels, c.-à-d. ., pour (i en c ("Var1", "Var2")) la façon dont vous le codez actuellement est très compliquée 4) votre appel à spplot n'est pas correct, vous lui passez un vecteur absurde.
—
Jeffrey Evans
Je suis vraiment désolé d'être si incompréhensible, mais c'est ma première vraie expérience avec R, j'ai ajouté les données dans la question principale, si cela ne vous dérange pas, pouvez-vous suggérer des moyens d'amélioration comme je l'ai vraiment a manqué d'idées
—
Ruvin Rafailov
