Miriad memo
Miriad official page at ATNF
MIRIAD for SMA
SMA Lab Instructions in the Second Asian Radio Astronomy School
インストールログ
Miriad official pageのインストールガイドを見ながら。本体の前に、以下をインストール。
- PGPLOT。Linux Wiki PGPLOT に沿って。
- gnu readline。INSTALLファイルの通りに。
コマンドメモ
- fits ファイルを読み込む
- モーメントマップ
- 速度方向スムージング
- cubeの一部抜き出し
- regrid
- Check the start/end velocity of the chunks
- Check the beamsize and fit gaussian to the source
- グリッドをビニングする
- Convolve with a beam (c.f.Image analysis using Miriad)
- ポリゴンで領域を指定し、統計情報を得る
- ポリゴンで領域を指定し、領域情報を書きだす
- 矩形領域のノイズレベルなどを見る
- Primary beam correction
- Plot UV coverage
-
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
-
miriad> imsub in=INFILENAME out=OUTFILENAME region='image(z1,z2)'
-
miriad> regrid in=INFILENAME out=OUTFILENAME axes=1,2 tin=REFIMAGE
-
miriad> uvlist vis=INFILE options=brief,spec
-
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
-
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)'
-
miriad> linmos in=INFILENAME out=OUTFILENAME
- options=sensitivity: Output an image giving the rms noise across the field.
-
miriad> smauvplt vis=INFILENAME axis=uc,vc device=/xw options=nobase
画像表示関係
- コントアマップをカラーpostscriptファイルに描きだす。
- チャネルマップのコントアを引く。
- プロファイルマップ。
- スペクトル表示。
- 画像情報を確認する
- ヘッダ情報を確認する
- ヘッダを書き換える
- 分点変換など
- uv distance vs amplitude plot
- uv distance vs amplitude plot (with S/N)
-
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> itemize in=INFILENAME
-
miriad> puthd in=INFILENAME/HEADER_ITEM value=value
-
miriad> regrid in=INFILENAME out=OUTFILENAME options=equisw
-
miriad> smauvamp vis=INFILENAME axis=uvdis,amp average=MINUTES nbin=No. of bin device=/xw
-
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
- Flag the window edge channels
- Remove the continuum from uv data
- Copy the desired window to new directory
- Check the rest frequency, if necessary
- Put the correct rest frequency, if necessary
- Produce images from visibilities
- CLEAN
- RESTORE
-
miriad> uvflux vis=INFILE line=wide,1,1(for LSB)/2(for USB)
-
miriad> uvflag vis=INFILE select='win(NN)' flagval=f line=chan,Nchan, StartChan, width, step
-
miriad> uvlin vis=INFILE chans=line-free_chan_start,end,start,end out=OUTFILE order=1 mode=line options=nowindow
-
miriad> uvcat vis=INFILE out=OUTFILE select='win(NN)'
-
miriad> prthd in=INFILE
-
miriad> puthd in=INFILE/restfreq value=RESTFREQ
-
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