使用Matplotlib的半对数图中的纵横比

时间:2014-07-21 13:33:59

标签: python matplotlib aspect-ratio logarithm

当我在matplotlib中绘制函数时,绘图由矩形框起。我希望这个矩形的长度和高度的比率由黄金均值给出,即,dx / dy = 1.618033 ...

如果x和y比例是线性的,我使用谷歌找到了这个解决方案

import numpy as np
import matplotlib.pyplot as pl
golden_mean = (np.sqrt(5)-1.0)/2.0
dy=pl.gca().get_ylim()[1]-pl.gca().get_ylim()[0]
dx=pl.gca().get_xlim()[1]-pl.gca().get_xlim()[0]
pl.gca().set_aspect((dx/dy)*golden_mean,adjustable='box')

如果是对数日志图,我想出了这个解决方案

dy=np.abs(np.log10(pl.gca().get_ylim()[1])-np.log10(pl.gca().get_ylim()[0]))
dx=np.abs(np.log10(pl.gca().get_xlim()[1])-np.log10(pl.gca().get_xlim()[0]))
pl.gca().set_aspect((dx/dy)*golden_mean,adjustable='box')

然而,对于半对数图,当我调用set_aspect时,我得到

UserWarning: aspect is not supported for Axes with xscale=log, yscale=linear

有人能想到解决这个问题的方法吗?

1 个答案:

答案 0 :(得分:7)

最简单的解决方案是记录您的数据,然后使用lin-lin的方法。

然后,您可以标记轴,使其看起来像正常的对数图。

ticks = np.arange(min_logx, max_logx, 1)
ticklabels = [r"$10^{}$".format(tick) for tick in ticks]

pl.yticks(ticks, ticklabels)

如果您的值高于10e9,则需要三对括号,两对用于LaTeX括号,一对用于.format()

ticklabels = [r"$10^{{{}}}$".format(tick) for tick in ticks]

编辑: 如果你还需要0.1ex ... 0.9ex的刻度,你也想要使用次刻度: 它们需要位于log10(1),log10(2),log10(3)......,log10(10),log10(20)......

您可以使用以下命令创建和设置:

minor_ticks = []
for i in range(min_exponent, max_exponent):
    for j in range(2,10):
         minor_ticks.append(i+np.log10(j))


plt.gca().set_yticks(minor_labels, minor=True)