0

我想知道是否有人知道如何使用 python 包 matplotlib.imshow 绘制适合图像，并使用相应的世界坐标系值，甚至可能是右升或赤纬作为 x 和 y 值而不是物理像素值，类似于底部本页图： http: //astroplotlib.stsci.edu/page_images.htm

不幸的是，提供的脚本在 IDL 中......我还不精通......

如果我概述了我的 gridspec 布局可能会有所帮助：

fig = pyplot.figure(figsize=(11,11))

gridspec_layout = gridspec.GridSpec(3,3)
gridspec_layout.update(hspace=0.0, wspace=0.0)

hdulist_org_M33_UVM2 = fits.open('myfits.fits')
wcs = WCS(hdulist_org_M33_UVM2[0].header)

pyplot_2 = fig.add_subplot(gridspec_layout[2])

ax = WCSAxes(fig, [0.1, 0.1, 0.8, 0.8], wcs=wcs)

pyplot_2.add_axes(ax)

但没有运气。

非常感谢。

4

1 回答 1

0

一种解决方案可能是使用subplot参数 toFITSFigure，并从您的 gridspec 中获取边界。

大致如下：

from matplotlib import pyplot
from matplotlib import gridspec
import aplpy

fig = pyplot.figure(figsize=(11, 11))
gridspec_layout = gridspec.GridSpec(3, 3)
gridspec_layout.update(hspace=0.0, wspace=0.0)
axes = fig.add_subplot(gridspec_layout[2])
m33 = aplpy.FITSFigure('wfpcii.fits', figure=fig,
                       subplot=list(gridspec_layout[2].get_position(fig).bounds),
                       # dimensions and slices are not normally necessary;
                       # my test-figure has 3 axes
                       dimensions=[0, 1], slices=[0])
print(dir(gridspec_layout[2]))
print(gridspec_layout[2].get_position(fig).bounds)
m33.show_colorscale()
pyplot.show()

不是很漂亮，但它会工作。如果我遇到一种将 aFITSFigure直接附加到轴的更简单方法，我将修改此答案或添加一个新答案。

于 2016-06-13T07:55:33.607 回答