アニメーション パッケージを使用して、マップ上のポイントの「進化する」プロットを生成しようとしています。マップはシェープファイル (readShapeSpatial/readShapeLines 関数) から生成されます。
問題は、for ループでプロットされた場合、結果が加算的であるのに対し、理想的な結果は進化することです。
私が見逃している par() の使用方法はありますか?
私の質問は: ポイント関数からプロットされたポイントだけをクリアし、図全体をクリアせずにシェープファイルを再グラフ化する必要がない方法はありますか?
誰かがコードを見たい場合:
# plotting underlying map
newyork <- readShapeSpatial('nycpolygon.shp')
routes <- readShapeLines('nyc.shp')
par(bg="grey25")
plot(newyork, lwd=2, col ="lightgray")
plot(routes,add=TRUE,lwd=0.1,col="lightslategrey")
# plotting points and save to GIF
ani.options(interval=.05)
saveGIF({
par(bg="grey25")
# Begin loop
for (i in 13:44){
infile <-paste("Week",i,".csv",sep='')
mydata <-read.csv(file = infile, header = TRUE, sep=",")
plotvar <- Var$Para
nclr <- 4
plotclr <-brewer.pal(nclr,"RdPu")
class<- classIntervals(plotvar,nclr,style = "pretty")
colcode <- findColours(class,plotclr)
points(Var$Lon,Var$Lat,col=colcode)
}
})