はFITSはGalSim

Question

セットアップ

• のMacOS 10.12.6
• のPython 2.7
• Galsim 1.4.4
• numpyの1.13.1
• Astropy 1.3.3でファイルを読み取ることができません。
• Fitsio

ゴール

私の目標は、シミュレートされた幅広い分野HSTの画像（12288個の* 12288ピクセル、.fits形式でおよそ600メガバイト）を取ると、それにWFIRST検出機能を追加することです。フィットファイルを読み込む

問題

残念ながら、私が直面している困難。
主な質問は、ガルシムで画像を操作できるように画像を調整する方法です。私はgalsimを試してみました

を試してみた

は、メソッド 'galsim.fits.read（ 『fine_name.fitsを』）' にフィットし、私が持っている：私は試してみました

galsim.Image(bounds=galsim.BoundsI(xmin=1, xmax=12288, ymin=1, ymax=12288), array= 
    array([[ 0., 0., 0., ..., 0., 0., 0.], 
[ 0., 0., 0., ..., 0., 0., 0.], 
[ 0., 0., 0., ..., 0., 0., 0.], 
..., 
[ 0., 0., 0., ..., 0., 0., 0.], 
[ 0., 0., 0., ..., 0., 0., 0.], 
[ 0., 0., 0., ..., 0., 0., 0.]], dtype=float32), wcs=galsim.PixelScale(1.0)) 

fitsio.read(“file_name.fits”)でも、私は 'galsim.Convolve'を通してPSFを畳み込むことができず、検出器の機能を追加することもできません。

-Thanks

Answer 1

あなたがgalsim.fits.readのdocstringを読めば、あなたはあなたの実験で発見したとして、それは、galsim.Imageを返すことがわかります。 GSObjectを作成するには、galsim.InterpolatedImageをインスタンス化する必要があります。これにより、畳み込みや他の操作を行うことができます。 galsim.InterpolatedImageのドキュメントストリングを読んだり、galsim demo11.pyを見て、これが実際に行われているかの例を見ることができます（イメージファイルをInterpolatedImageに読み込んで、何か他のものと畳み込む）。

検出機能doはGSObjectsではなく画像に追加されます。この例については、GalSimのdemo13.pyを参照してください。

GalSimの基本的な使用方法は、クイックリファレンスガイドhttps://github.com/GalSim-developers/GalSim/blob/releases/1.4/doc/GalSim_Quick_Reference.pdfとGalSimリポジトリのデモで説明しています。