2 つの .FITS 画像があります。1 つはいくつかの星と銀河の画像です。もう1つは、等高線としてプロットしたい有意性マップです。PyWCSGrid2 はこれを行うための python モジュールですが、しばらくの間、一方を他方に重ねようとしましたが、両方を同時に表示することはできません。これが機能しない理由はありますか?
import matplotlib.pyplot as plt
import sys
import pyfits
import pywcsgrid2
from mpl_toolkits.axes_grid1.axes_divider import make_axes_locatable
imagename=str("example1.fits")
def setup_axes():
ax = pywcsgrid2.subplot(111, header=f_radio[0].header)
return ax
# GET IMAGE1
data = f_radio[0].data #*1000
ax = setup_axes()
# prepare figure & axes
fig = plt.figure(1)
#GET CONTOUR SOURCE:
f_contour_image = pyfits.open("IMAGE2.fits")
data_contour = f_contour_image[0].data
# DRAW CONTOUR
cont = ax.contour(data_contour, [5, 6, 7, 8, 9],
colors=["r","r","r", "r", "r"], alpha=0.5)
# DRAW IMAGE
im = ax.imshow(data, cmap=plt.cm.gray, origin="lower", interpolation="nearest",alpha=1.0)
plt.show()
更新:問題は、2 つの画像の縮尺が異なるため、実際には等高線をプロットしているようですが、左下隅に小さくて waaaay です。どうにかしてそれらを再スケーリングする必要があります。