子图投影会覆盖tick_params

时间:2018-02-05 11:00:48

标签: python matplotlib astropy

我尝试使用fit标头提供的坐标绘制图像,并通过astropy的WCS获得。

from astropy import wcs
import numpy as np
import matplotlib.pyplot as plt
hdr = r[1].header #r ist the repsective fits files, I copied content "w" at the end of the page
w = wcs.WCS(hdr)
ax = plt.subplot(projection=w)
ax.imshow(np.ones((100,100)),origin='lower')

ax.tick_params(axis='x', which='major', labelsize='large',width=25)
ax.tick_params(axis='y', which='major', labelsize='large')
plt.show()

plot with preoction on

可以看到tick_params被忽略了。

wenn我也这样做,但关闭了预测,即:

ax = plt.subplot()
ax.imshow(np.ones((100,100)),origin='lower')

ax.tick_params(axis='x', which='major', labelsize='large',width=25)
ax.tick_params(axis='y', which='major', labelsize='large')
plt.show()

tick_parmas works again

Tick_params正在重新开始工作。

知道这里可能出现什么问题吗?

WCS是:

 print(w)
 WCS Keywords

 Number of WCS axes: 2
 CTYPE : 'RA---TAN'  'DEC--TAN'  
 CRVAL : 266.41798186205955  -29.006968367892327  
 CRPIX : 248.5  340.0  
 NAXIS : 497  680

1 个答案:

答案 0 :(得分:2)

WCS投影完全取代matplotlib轴,请参阅Ticks, tick labels, and grid lines。因此,您不能再使用matplotlib方法,或者至少您不能指望它们对实际绘图有任何影响。

相反,它需要使用WCS的方法。 所以,如果

ax = plt.subplot(projection=wcs)

WCSAxesSubplot,您可以将x轴设为ax.coords[0],将y轴设为ax.coords[1]。然后你可以设置ticklabel大小

ax.coords[0].set_ticklabel(size="large")

,刻度线宽度为

ax.coords[0].set_ticks(width=25)

set_ticklabelset_ticks两种方法是astropy.visualization.wcsaxes.coordinate_helpers.CoordinateHelper类的方法。我不确定是否有一些完整的可用方法参考,但您可以随时查看the source code以查看公开的方法。

一些完整的示例(基于文档中的one of the examples):

import matplotlib.pyplot as plt

from astropy.wcs import WCS
from astropy.io import fits
from astropy.utils.data import get_pkg_data_filename

filename = get_pkg_data_filename('galactic_center/gc_msx_e.fits')

hdu = fits.open(filename)[0]
wcs = WCS(hdu.header)

ax = plt.subplot(projection=wcs)

ax.imshow(hdu.data, vmin=-2.e-5, vmax=2.e-4, origin='lower')

ax.coords.grid(True, color='white', ls='solid')
ax.coords[0].set_axislabel('Galactic Longitude')
ax.coords[1].set_axislabel('Galactic Latitude')

ax.coords[0].set_ticks(width=25)
ax.coords[0].set_ticklabel(size="large")

plt.show()

enter image description here