Miriad memo

Miriad official page at ATNF
SMA Lab Instructions in the Second Asian Radio Astronomy School


Miriad official pageインストールガイドを見ながら。本体の前に、以下をインストール。


  • fits ファイルを読み込む
    • miriad> fits in=INFILENAME op=xyin out=OUTFILENAME
  • モーメントマップ
    • miriad> moment in=INFILENAME out=OUTFILENAME mom=0 region='images(channel1,channel2)'
      mom=2 で出てくるのはσであってFWHMではないことに注意。
  • 速度方向スムージング
    • miriad> hanning in=INFILENAME out=OUTFILENAME width=3
      miriad> imsub in=INFILENAME out=OUTFILENAME incr=1,1,3
  • cubeの一部抜き出し
    • miriad> imsub in=INFILENAME out=OUTFILENAME region='image(z1,z2)'
  • regrid
    • miriad> regrid in=INFILENAME out=OUTFILENAME axes=1,2 tin=REFIMAGE
  • Check the start/end velocity of the chunks
    • miriad> uvlist vis=INFILE options=brief,spec
  • Check the beamsize and fit gaussian to the source
    • miriad> imfit vis=INFILE object=gauss region='arcsec,box(-5,-5,5,5)'
  • グリッドをビニングする
    • miriad> imbin in=INFILENAME bin=1,1,1,1,5,5 out=OUTFILENAME
  • Convolve with a beam (c.f.Image analysis using Miriad)
    • miriad> convol map=INFILENAME out=OUTFILENAME fwhm=HPBW (arcsec) scale=(1.13309*(HPBW/pixsize)^2)^-1
  • ポリゴンで領域を指定し、統計情報を得る
    • miriad> cgcurs in=INFILENAME type=c region='arcsec,box(x1,y1,x2,y2)' slev=a,level levs=3,4,5,6 device=/xw options=stats
  • ポリゴンで領域を指定し、領域情報を書きだす
    • miriad> cgcurs in=INFILENAME type=c region='arcsec,box(x1,y1,x2,y2)' slev=a,level levs=3,4,5,6 device=/xw options=region
      領域ファイルは cgcurs.region という名前で書きだされる。cleanの領域指定などに使える
  • 矩形領域のノイズレベルなどを見る
    • miriad> histo in=INFILENAME region='arcsec,box(x1,y1,x2,y2)'
  • Primary beam correction
    • miriad> linmos in=INFILENAME out=OUTFILENAME
      - options=sensitivity: Output an image giving the rms noise across the field.
  • Plot UV coverage
    • miriad> smauvplt vis=INFILENAME axis=uc,vc device=/xw options=nobase


  • コントアマップをカラーpostscriptファイルに描きだす。
    • miriad> cgdisp in=INFILENAME1,INFILENAME2 type=c,c slev=a,level1,a,level2 levs1=3,6,9,,, levs2=3,6,9,,, device=FILENAME.ps/cps labtyp=hms,dms options=full,blacklab beamtyp=b,l,1
  • チャネルマップのコントアを引く。
    • miriad> cgdisp in=INFILENAME type=[pixel,contour] chan=1,3 slev=a,level device=/xw labtyp=hms,dms options=3value 3format='f5.1' olay=OVERLAY FILENAME
  • プロファイルマップ
    • miriad> cgspec in=INFILENAME1,INFILENAME2 type=pixel,spec device=/xw labtyp=hms,dms options=full olay=grid.txt stick=5,2.0
  • スペクトル表示
    • miriad> imspect in=INFILENAME region='arcsec,box(x1,y1,x2,y2)(z1,z2)' device=/xw yaxis=average hann=1
  • 画像情報を確認する
    • miriad> prthd in=INFILENAME
  • ヘッダ情報を確認する
  • ヘッダを書き換える
    • miriad> puthd in=INFILENAME/HEADER_ITEM value=value
  • 分点変換など
    • miriad> regrid in=INFILENAME out=OUTFILENAME options=equisw
  • uv distance vs amplitude plot
    • miriad> smauvamp vis=INFILENAME axis=uvdis,amp average=MINUTES nbin=No. of bin device=/xw
  • uv distance vs amplitude plot (with S/N)
    • miriad> uvamp vis=INFILENAME bin=No. of bins,width of bins,unit of bins device=/xw

SMA Imaging with Miriad

Standard(?) procedure

  • Check the flux of the primary calibrator
    • miriad> uvflux vis=INFILE line=wide,1,1(for LSB)/2(for USB)
  • Flag the window edge channels
    • miriad> uvflag vis=INFILE select='win(NN)' flagval=f line=chan,Nchan, StartChan, width, step
  • Remove the continuum from uv data
    • miriad> uvlin vis=INFILE chans=line-free_chan_start,end,start,end out=OUTFILE order=1 mode=line options=nowindow
  • Copy the desired window to new directory
    • miriad> uvcat vis=INFILE out=OUTFILE select='win(NN)'
  • Check the rest frequency, if necessary
    • miriad> prthd in=INFILE
  • Put the correct rest frequency, if necessary
    • miriad> puthd in=INFILE/restfreq value=RESTFREQ
  • Produce images from visibilities
    • miriad> invert vis=INFILE map=DIRTYMAP beam=BEAMMAP cell=cellsize imsize=Xsize, Ysize line=vel/chan, Nchan, STARTChan, width, step
      - weighting: default: Unifrom, sup=0: Natural, robust=0.5: Robust
      - options=double,systemp: Double the beam pattern, weighting based on Tsys.
    • miriad> clean map=DIRTYMAP beam=BEAMMAP out=MODELMAP cutoff=3 sigma gain=gain niters=Niters
    • miriad> restor map=DIRTYMAP beam=BEAMMAP model=CELANMAP out=FINALMAP