Как я могу автоматически установить второстепенные тики на оси журнала? - PullRequest
Новая
       15

Как я могу автоматически установить второстепенные тики на оси журнала?

0 голосов
/ 28 июня 2018

Я ищу что-то вроде AutoMinorLocator для логарифмической шкалы.

Я знаю, что второстепенные тики по умолчанию включены, однако вызов метода .hist() удаляет их: 

>>> ax.get_xticks(minor=True)
array([2.e-02, 3.e-02, 4.e-02, 5.e-02, 6.e-02, 7.e-02, 8.e-02, 9.e-02,
       2.e-01, 3.e-01, 4.e-01, 5.e-01, 6.e-01, 7.e-01, 8.e-01, 9.e-01,
       2.e+00, 3.e+00, 4.e+00, 5.e+00, 6.e+00, 7.e+00, 8.e+00, 9.e+00,
       2.e+01, 3.e+01, 4.e+01, 5.e+01, 6.e+01, 7.e+01, 8.e+01, 9.e+01,
       2.e+02, 3.e+02, 4.e+02, 5.e+02, 6.e+02, 7.e+02, 8.e+02, 9.e+02,
       2.e+03, 3.e+03, 4.e+03, 5.e+03, 6.e+03, 7.e+03, 8.e+03, 9.e+03]) 
>>> ax.hist([], bins=[1e-7, 1e-3])
(array([0.]), array([1.e-07, 1.e-03]), <a list of 1 Patch objects>)
>>> ax.get_xticks(minor=True)
array([], dtype=float64)

Следующий код: 

major = ax.get_xticks(minor=False)
ax.set_xticks(np.concatenate([np.linspace(0.2, 0.9, 8) * major.min()] +\
                             [np.linspace(1, 9, 9) * m
                              for m in major]),
              minor=True)

хорошо работает для мелких тиков, но разрушает автоматические ограничения ...

1 Ответ

0 голосов
/ 28 июня 2018

Лучшее, что у меня есть, это: 

start, end = ax01.get_xlim()
startLg = int(np.ceil(np.log10(start)))
endLg = int(np.floor(np.log10(end)))
major = np.logspace(startLg, endLg, endLg - startLg + 1)
ax.set_xticks([x for x in chain(np.linspace(0.2, 0.9, 8) * major.min(),
                                *[np.linspace(1, 9, 9) * m for m in major])
               if start <= x <= end],
              minor=True)

но результат не выглядит потрясающе.

Кажется, что Matplotlib автоматически решает скрыть мелкие тики, так как основные тики пропускаются каждые второе десятилетие: 

>>> ax.hist([], bins=[1e-7, 1e-3])
(array([0.]), array([1.e-07, 1.e-03]), <a list of 1 Patch objects>)
>>> ax.get_xticks(minor=False)
array([1.e-10, 1.e-08, 1.e-06, 1.e-04, 1.e-02, 1.e+00])
...