#immoments

cubeデータを様々な方法で軸を減らすことができるcasaのtaskの一つ
- 例えば、cubeを積分強度や速度場、速度分散などにできる
公式ドキュメントはこちら

##指定できるパラメーター

imagename           =         ''        #  Name of the input image
moments             =        [0]        #  List of moments you would like to compute
axis                = 'spectral'        #  The momement axis: ra, dec, lat, long, spectral, or stokes
region              =         ''        #  Region selection. Default is to use the full image.
box                 =         ''        #  Rectangular region(s) to select in direction plane. Default is to use the entire direction plane.
chans               =         ''        #  Channels to use. Default is to use all channels.
stokes              =         ''        #  Stokes planes to use. Default is to use all Stokes planes.
mask                =         ''        #  Mask to use. Default is none.
includepix          =         -1        #  Range of pixel values to include
excludepix          =         -1        #  Range of pixel values to exclude
outfile             =         ''        #  Output image file name (or root for multiple moments)

とりあえずは、

imagename
moments
outfile

を指定すればよい。

momentsは以下を参考に値を指定する。

moments=0   - 積分強度  
moments=1   - 速度場
moments=2   - 速度分散  
moments=3   - 強度の中央値  
moments=4   - coordinateの中央値
moments=5   - スペクトルの平均の標準偏差  
moments=6   - root mean square of the spectrum  
moments=7   - absolute mean deviation of the spectrum  
moments=8   - スペクトルの最大値  
moments=9   - coordinate of the maximum value of the spectrum  
moments=10  - スペクトルの最小値  
moments=11  - coordinate of the minimum value of the spectrum

##参考スクリプト

積分強度を作成し、imageを出力

immoments.py

import os
import pyfits

fname = "oooo.image" or "xxxxx.fits"

fn = os.path.splitext(fname)[0]
outfile = fn+"_mom０"
immoments(imagename = fname,
          moments = [0],
          outfile = outfile+".image")

exportfits(imagename = outfile+".image",
           fitsimage = outfile+".fits",
           dropdeg = True,
           overwrite = True,
           history = False)

その他の関数はこちら