2

我正在尝试根据随时间逐步建立的散点图创建动画。用例是我有一个包含大约 200 万个点的数据库，每个点都有时间戳，并且想要生成显示特定日期或之前的所有点的帧。

在不保存图像的情况下，我可以通过首先调用 来做到这一点plot()，然后有一个 for 循环，使用该points()函数递增地绘制每一天的数据。

当我尝试使用下面的代码保存图像时，我收到“尚未调用 plot.new”的错误。据我了解，dev.off()需要保存图像，但这也会关闭正在绘制的设备。有没有办法解决这个问题？由于数据的大小，不得不为每一帧重新绘制数据并不是一个很好的选择。

plot(info$lon, info$lat, xlim=c(0,30), ylim=c(30,60))
for (i in c(1:length(allDates))){
    filename=paste(sprintf('%05d', i), ".png", sep="")  
    png(filename=fileName)
    # (code that gets the data for a particular date via a database query)
    points(info$lon, info$lat, cex=0.1)
    dev.off()
}

更新： @roman-lustrik 的评论ggsave()是我正在寻找的内容，结果如下：

plotObj = ggplot(...) + geom_point() + xlim(...) + ylim(...)
for (i in c(1:length(allDates))){
    filename=paste(sprintf('%05d', i), ".png", sep="")  
    # (code that gets the data for a particular date via a database query)
    plotObj = plotObj + geom_point(data=info, aes(x=lon, y=lat), size=0.5)
    print(plotObj)
    ggsave(filename=filename, width=6, height=6)
}

但是，这仍然有点慢，所以我目前快速渲染图像的解决方案是使用与原始代码类似的代码，但我只是plot()用来渲染具有单个日期数据的帧（使用透明背景）。为了逐步堆叠图像，我使用了一个 bash 脚本，该脚本使用 imagemagickconvert -composite命令将两个图像混合在一起。然后将此混合图像与下一个日期的图像混合，依此类推，直到最终图像显示所有数据：

#!/bin/bash
for i in $files
do
    convert $prevFile $i -composite ./stackedImages/$i
    prevFile=./stackedImages/$i
done
4

1 回答 1

1

如果我理解了，您想要获得几个具有不同点数的 png 文件，并且第一个点由plot(info$lon, info$lat, xlim=c(0,30), ylim=c(30,60)).

你可以这样做：

temp1 <- info$lon
temp2 <- info$lat
for (i in c(1:length(allDates))){
        filename=paste(sprintf('%05d', i), ".png", sep="")  
        png(filename=fileName)
        plot(temp1, temp2, xlim=c(0,30), ylim=c(30,60))
        # (code that gets the data for a particular date via a database query)
        points(info$lon, info$lat, cex=0.1)
        dev.off()
        temp1 <- c(temp1,info$lon)
        temp2 <- c(temp2,info$lat)
}
于 2012-07-23T12:03:02.133 回答