我尝试使用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()
可以看到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_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
答案 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_ticklabel
和set_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()