在R的网格图形中保留纵横比

Question

要使用graphics我所说的低级包装设备绘制一个高度比其宽度大2倍的"交叉"矩形:

xlim <- c(0, 500)
ylim <- c(0, 1000)
plot.new()
plot.window(xlim, ylim, asp=1)
rect(xlim[1], ylim[1], xlim[2], ylim[2])
lines(c(xlim[1], xlim[2]), c(ylim[1], ylim[2]))
lines(c(xlim[1], xlim[2]), c(ylim[2], ylim[1]))

该图有一个很好的功能:保留纵横比,这样如果我改变绘图窗口的大小,我会得到相同的高宽比.

如何获得与grid图形相同的结果？

Answer 1

您应该创建一个使用方形标准化父坐标的视口，请参阅?unit：

"snpc"：（...）这对于制作占视口一定比例但必须是正方形（或具有固定长宽比）的东西很有用。

这是代码：

library('grid')
xlim <- c(0, 500)
ylim <- c(0, 1000)
grid.newpage() # like plot.new()
pushViewport(viewport( # like plot.window()
        x=0.5, y=0.5, # a centered viewport
        width=unit(min(1,diff(xlim)/diff(ylim)), "snpc"), # aspect ratio preserved
        height=unit(min(1,diff(ylim)/diff(xlim)), "snpc"),
        xscale=xlim, # cf. xlim
        yscale=ylim  # cf. ylim
))
# some drawings:
grid.rect(xlim[1], ylim[1], xlim[2], ylim[2], just=c(0, 0), default.units="native")
grid.lines(xlim, ylim, default.units="native")
grid.lines(xlim, rev(ylim), default.units="native")

default.units例如中的参数强制grid.rect绘图函数使用本机 ( xscale/ yscale) 视口坐标。 just=c(0, 0)表示xlim[1], ylim[1]表示矩形的左下节点。