当在lattice::xyplot
和trellis.device
中使用layout
时,你可以在一个PDF中获得几个页面:
trellis.device(pdf, file="myfile.pdf")
data(mtcars)
xyplot(hp~mpg|gear, data=mtcars, layout=c(1, 1))
dev.off()
我想使用相同的方法来获得一个multi-page SVG。我尝试过cairo::svg
,包gridSVG
、SVGAnnotation
和RSVGTipsDevice
都没有成功:只保存了trellis
对象的最后一页。
有没有使用R代码的解决方案?
谢谢!
发布于 2011-11-12 18:20:04
我最终决定直接请教grid和gridSVG的创建者Paul Murrell。他亲切地提供了一些很好的建议和代码,用类似于@mbq的建议的方法完全解决了这个问题。我修改了他的代码,编写了这个函数:
library(gridSVG)
library(XML)
animateTrellis <- function(object, file='animatedSVG.svg',
duration=.1, step=2, show=TRUE){
nLayers <- dim(object)
stopifnot(nLayers>1)
for (i in seq_len(nLayers)){
p <- object[i]
label <- p$condlevels[[1]][i]
##Create intermediate SVG files
g <- grid.grabExpr(print(p, prefix=label),
name=label)
if (i==1){ ## First frame
ga <- animateGrob(g, group=TRUE,
visibility="hidden",
duration=duration, begin=step)
} else if (i==nLayers){ ##Last Frame
gg <- garnishGrob(g, visibility='hidden')
ga <- animateGrob(gg, group=TRUE,
visibility="visible",
duration=duration, begin=step*(i-1))
} else { ##any frame
gg <- garnishGrob(g, visibility='hidden')
gaV <- animateGrob(gg, group=TRUE,
visibility="visible",
duration=duration, begin=step*(i-1))
ga <- animateGrob(gaV, group=TRUE,
visibility="hidden",
duration=duration, begin=step*i)
}
grid.newpage()
grid.draw(ga)
fich <- tempfile(fileext='.svg')
gridToSVG(fich)
## Combine all
if (i==1) {
svgTop <- xmlParse(fich)
nodeTop <- getNodeSet(svgTop,
"//svg:g[@id='gridSVG']",
c(svg="http://www.w3.org/2000/svg"))[[1]]
} else {
svgChildren <- xmlParse(fich)
node <- getNodeSet(svgChildren,
"//svg:g[@id='gridSVG']/*",
c(svg="http://www.w3.org/2000/svg"))
addChildren(nodeTop, node)
}
unlink(fich)
}
saveXML(svgTop, file=file)
dev.off()
if (show) browseURL(file)
invisible(svgTop)
}
然后我可以生成带有动画的SVG文件:
p <- xyplot(Sepal.Length~Petal.Length|Species, data=iris, layout=c(1, 1))
animateTrellis(p, file='iris.svg')
https://stackoverflow.com/questions/7554162
复制