私はいたるところを見てきましたが、この問題の解決策が見つかりません。かなり簡単なはずですが、見てみましょう。
.FITS 形式のデータ キューブがあり、それを 2D FITS 画像に折りたたむ必要があります。データ キューブには、2 つの空間次元と 1 つのスペクトル/速度次元があります。
立方体にロードしてこれらすべてのレイヤーを平坦化する (つまり、スペクトル/速度軸に沿ってそれらを統合する) 単純な Python ルーチンを探しているだけです。助けてくれてありがとう。
これは少し古い質問ですが、spectral-cubeはこれに対するより良い解決策を提供します。
Teachey の回答に基づく例:
from spectral_cube import SpectralCube
cube = SpectralCube.read(filename)
summed_image = cube.sum(axis=0)
summed_image.hdu.writeto(new_filename)