处理高色散光谱解读.docx
- 文档编号:11891318
- 上传时间:2023-04-08
- 格式:DOCX
- 页数:59
- 大小:472.44KB
处理高色散光谱解读.docx
《处理高色散光谱解读.docx》由会员分享,可在线阅读,更多相关《处理高色散光谱解读.docx(59页珍藏版)》请在冰豆网上搜索。
处理高色散光谱解读
附件2.IRAF处理高色散光谱
北京师范大学天文系
杨轶0710160009
1.
确定处理数据期间的工作目录,并将文件readlog.cl拷贝到工作目录下;
如:
/home/ngc/w070424echell
2.
启动iraf,打开xgterm和ds9;
打开终端;
bash-3.2$su
密码:
[root@localhostngc]#cdiraf
[root@localhostiraf]#ds9&
[root@localhostiraf]#xgterm-sbr
[root@localhostiraf]#cl
3.
修改login.cl;
在终端中键入:
[root@localhostiraf]#vilogin.cl
键入:
taskreadlog="/home/ngc/w070424echell/readlog.cl"
“”当中应当为之后处理光谱数据的工作目录,并在目录后面加上/readlog.cl
如下:
packageuser
task$adb$bc$cal$cat$comm$cp$csh$date$dbx$df$diff="$foreign"
task$du$find$finger$ftp$grep$lpq$lprm$ls$mail$make="$foreign"
task$man$mon$mv$nm$od$ps$rcp$rlogin$rsh$ruptime="$foreign"
task$rwho$sh$spell$sps$strings$su$telnet$tip$top="$foreign"
task$vi$emacs$w$wc$less$rusers$sync$pwd$gdb="$foreign"
task$xc$mkpkg$generic$rtar$wtar$buglog="$foreign"
#task$fc="$xc-h$*-limfort-lsys-lvops-los"
task$fc=("$"//envget("iraf")//"unix/hlib/fc.csh"//
"-h$*-limfort-lsys-lvops-los")
task$nbugs=("$(setenvEDITOR'buglog-e';"//
"less-Cqm+G"//envget("iraf")//"local/bugs.*)")
task$cls="$clear;ls"
task$clw="$clear;w"
task$pg=("$(less-Cqm$*)")
taskreadlog="/home/ngc/w070424echell/readlog.cl"
加入这一行(不用加$,$代表变量)
if(access("home$loginuser.cl"))
cl<"home$loginuser.cl";
keep
4.
改写log文件;
切换到工作目录下:
[root@localhostiraf]#cd/home/ngc/w070424echell
打开log文件:
[root@localhostw070424echell]#vi20070126log
将log文件改写为如下格式:
24/04/07
25/04/07
w0424001flatflat20:
00:
0018000:
00:
00+00:
00:
002000
w0424002flatflat20:
00:
0018000:
00:
00+00:
00:
002000
w0424003flatflat20:
00:
0018000:
00:
00+00:
00:
002000
w0424004flatflat20:
00:
0018000:
00:
00+00:
00:
002000
w0424005flatflat20:
00:
0018000:
00:
00+00:
00:
002000
w0424006biaszero20:
00:
00000:
00:
00+00:
00:
002000
w0424007biaszero20:
00:
00000:
00:
00+00:
00:
002000
w0424008biaszero20:
00:
00000:
00:
00+00:
00:
002000
w0424009biaszero20:
00:
00000:
00:
00+00:
00:
002000
w0424010biaszero20:
00:
00000:
00:
00+00:
00:
002000
w0424011Tharcomp20:
00:
001000:
00:
00+00:
00:
002000
w0424012Tharcomp20:
00:
002000:
00:
00+00:
00:
002000
w0424013J094543+410417object20:
54:
41360009:
45:
43.66+41:
04:
17.522000
w0424014J094543+410417object22:
01:
26360009:
45:
43.66+41:
04:
17.522000
w0424015J161816+071957object23:
24:
36360016:
18:
16.51+07:
19:
57.562000
w0424016J161816+071957object24:
25:
44360016:
18:
16.51+07:
19:
57.562000
w0424017J161816+071957object25:
32:
10360016:
18:
16.51+07:
19:
57.562000
w0424018J161816+071957object26:
34:
41360016:
18:
16.51+07:
19:
57.562000
w04240195_Serstand27:
44:
3160015:
19:
31.09+01:
45:
01.92004
w0424020a_Boostand27:
57:
411014:
15:
50.83+19:
09:
48.52004
w0424021a_Boostand27:
58:
26514:
15:
50.83+19:
09:
48.52004
w0424022dat_Ophstand28:
01:
4810016:
14:
33.33-03:
42:
15.72004
w0424023dat_Ophstand28:
04:
305016:
14:
33.33-03:
42:
15.72004
w0424024beta_Ophstand28:
07:
485017:
43:
40.21+04:
33:
57.12004
w0424025beta_Ophstand28:
09:
0210017:
43:
40.21+04:
33:
57.12004
w042402631_Aqlstand28:
14:
0190019:
25:
10.7+11:
57:
10.82004
w0424027Tharcomp28:
14:
011000:
00:
00+00:
00:
002000
w0424028Tharcomp28:
14:
01500:
00:
00+00:
00:
002000
w0424029flatflat28:
14:
0118000:
00:
00+00:
00:
002000
w0424030flatflat28:
14:
0118000:
00:
00+00:
00:
002000
w0424031flatflat28:
14:
0118000:
00:
00+00:
00:
002000
w0424032flatflat28:
14:
0118000:
00:
00+00:
00:
002000
w0424033flatflat28:
14:
0118000:
00:
00+00:
00:
002000
w0424034biaszero28:
14:
01000:
00:
00+00:
00:
002000
w0424035biaszero28:
14:
01000:
00:
00+00:
00:
002000
w0424036biaszero28:
14:
01000:
00:
00+00:
00:
002000
w0424037biaszero28:
14:
01000:
00:
00+00:
00:
002000
w0424038biaszero28:
14:
01000:
00:
00+00:
00:
002000
w0424039biaszero28:
14:
01000:
00:
00+00:
00:
002000
w0424040biaszero28:
14:
01000:
00:
00+00:
00:
002000
5.
编辑列表文件;
[root@localhostw070424echell]#ls*.fit>fit1.lst
在终端中键入:
[root@localhostw070424echell]#vifit1.lst
w0424001.fit
w0424002.fit
……
w0424040.fit
再在vi编辑器中键入
:
%s/\.fit//,即去掉文件的后缀名.fit,再键入:
wq,写入并退出,文件变为:
w0424001
w0424002
……
w0424040
键入:
wq,保存并退出;
6.
接下来的操作在xgterm和ds9(saoimage)中进行
将文件格式由*.fits转为*.imh
cl>noao
no>imred
im>spec
sp>tv
tv>ccdred
注意在xgterm中也要切换到工作目录下:
cc>cd/home/ngc/w070424echell
cc>eparfits
IRAF
ImageReductionandAnalysisFacility
PACKAGE=dataio
TASK=rfits
fits_fil=@fit1.lst//.fitdatasource
file_lis=File/extensionslist
iraf_fil=@fit1.lst//.imhIRAFfilename
(make_im=yes)CreateanIRAFimage?
(long_he=no)PrintFITSheadercards?
(short_h=yes)Printshortheader?
(datatyp=)IRAFdatatype
(blank=0.)Blankvalue
(scale=yes)Scalethedata?
(oldiraf=no)UseoldIRAFnameinplaceofiraf_file?
(offset=0)Tapefileoffset
(mode=ql)
---------------------------------------------------------------------------------------
:
g
ls一下可看到每个fits文件对应多出了.imh文件。
这一步操作不必要,但这些fit文件可以留作备份,并移至其他目录下。
剩余的操作对新生成的.imh文件进行;
cc>mkdirfit_backup
cc>cp*.fit./fit_backup
cc>del*.fit
如此工作目录下只剩下.imh文件;
7.
设置天文台信息
cc>epaobserva
IRAF
ImageReductionandAnalysisFacility
PACKAGE=noao
TASK=observatory
command=listCommand(set|list|images)
obsid=BAOObservatorytoset,list,orimagedefault
images=*.imhListofimages
(verbose=no)Verboseoutput?
(observa=BAO)Observatoryidentification
(name=BeijingXingLongObservatory)Observatoryname
(longitu=242.425)Observatorylongitude(degrees)
(latitud=40.385)Observatorylatitude(degrees)
(altitud=950.)Observatoryaltitude(meters)
(timezon=-8.)Observatorytimezone
override=BAOObservatoryidentification
(mode=ql)
---------------------------------------------------------------------------------------
:
g
8.
为.imh加文件头
cc>epareadlog
IRAF
ImageReductionandAnalysisFacility
PACKAGE=user
TASK=readlog
(im=@070424log)Inputimage*
(bf=yes)Ifgetkeywordfromfile?
(keyfile=070424log)filecontainingkeyword*
(format=)Filenameobjnameimtypebtimeexptimeradecepoch
(obj=)Objectnamewhichinputintoimheader
(imtype=)Imagetype
(btime=)TheBeijingtimeofobjectimageobserved
(exptime=0.00)Exposuretimeofobject
(ra=)ObjectRa
(dec=)ObjectDec
(epoch=)Objectepoch
(list1=070424log)*
(mode=ql)
----------------------------------------------------------------------------------------
:
g
带*的项的参数应根据自己的log文件修改
如出现ERRORonline124:
parameter‘type’notfound
readlog(mode=h)
说明图像的类型没有被写入头文件中。
可以用imhead命令查看任一图像的文件头:
cc>imheadw0424013
w0424013[1024,1024][short]:
J094543+410417
Nobadpixels,min=154.,max=2289.
Linestoragemode,physdim[1024,1024],lengthofuserarea1175s.u.
CreatedWed12:
11:
3128-Jul-2010,LastmodifiedWed12:
11:
3128-Jul-2010
Pixelfile"localhost.locald!
/iraf/imdirs/ngc/w0424013.pix"[ok]
BLOCKED=T/CHECKFORPOSSIBLEBLOCKING
OBS-NUM=1/OBSERVATIONSEQUENCENAME
DATE-OBS=04/24/07/DATEOFSTARTOFOBSERVATION
TIME=20:
58:
25/BEIJINGTIME
EXPOSURE=3600.0/EXPOSURETIME(SEC)
SHUTTER='OPEN'/SHUTTERSTATE
RA=18:
58:
45.32/ACTUALRIGHTASCENSION
DEC=+06:
15:
00.83/ACTUALDECLINATION
HA=22:
48:
42.50/ACTUALHOURANGLE
DES-RA=18:
58:
45.27/DESIREDRIGHTASCENSION
DES-DEC=+06:
15:
45.53/DESIREDDECLINATION
EPOCH=2007.33/EPOCHOFDES-RAANDDES-DEC
RECORDS=1457/IMAGESIZEINFITSRECORDS
INSTRUME='PI1k*1kCCD'/THEINSTRUMENT'SNAME
OBSERVER=''/OBSERVER'SNAME
9.
裁剪图像
这一步视图像的质量而定。
可以先打开本底,平场,目标星,比较星,灯谱的图像进行检查,以确定裁剪掉的图像边缘的范围;
cc>epaccdp
IRAF
ImageReductionandAnalysisFacility
PACKAGE=ccdred
TASK=ccdproc
images=w*.imhListofCCDimagestocorrect*裁剪可一次性对所有图像执行
(output=)ListofoutputCCDimages
(ccdtype=)CCDimagetypetocorrect
(max_cac=0)Maximumimagecachingmemory(inMbytes)
(noproc=no)Listprocessingstepsonly?
(fixpix=no)FixbadCCDlinesandcolumns?
(oversca=no)Applyoverscanstripcorrection?
(trim=yes)Trimtheimage?
*裁剪图像
(zerocor=no)Applyzerolevelcorrection?
(darkcor=no)Applydarkcountcorrection?
(flatcor=no)Applyflatfieldcorrection?
(illumco=no)Applyilluminationcorrection?
(fringec=no)Applyfringecorrection?
(readcor=no)Convertzerolevelimagetoreadoutcorrection?
(scancor=no)Convertflatfieldimagetoscancorrection?
(readaxi=line)Readoutaxis(column|line)
(fixfile=)Filedescribingthebadlinesandcolumns
(biassec=[1023:
1024])Overscanstripimagesection
(trimsec=[10:
1018,5:
1021])Trimdatasection
*裁剪尺寸依图像而定
(zero=zero)Zerolevelcalibrationimage
(dark=)Darkcountcalibrationimage
(flat=flatre)Flatfieldimages
(illum=)Illuminationcorrectionimages
(fringe=)Fringecorrectionimages
(minrepl=1.)Minimumflatfieldvalue
(scantyp=shortscan)Scantype(shortscan|longscan)
(nscan=1)Numberofshortscanlines
(interac=no)Fitoverscaninteractively?
(functio=chebyshev)Fittingfunction
(order=1)Numberofpolynomialtermsorsplinepieces
(sample=*)Samplepointstofit
(naverag=1)Numberofsamplepointstocombine
(niterat=1)Numberofrejectioniterations
(low_rej=3.)Lowsigmarejectionfactor
(high_re=3.)Highsigmarejectionfactor
(grow=0.)Rejectiongrowingradius
(mode=ql)
----------------------------------------------------------------------------------------
:
g
10.
合并本底
查看日志文件,确定本底图片的名称;
建立一个文件夹,专门存放本底图片:
cc>mkdirzeroc
cc>mvw0424006.fitszeroc
cc>mvw0424007.fitszeroc
cc>mvw0424008.fitszeroc
cc>mvw0424009.fitszeroc
cc>mvw0424010.fitszeroc
进入存放本底的文件夹进行本底合并;
cc>cdzeroc
cc>epazeroc
IRAF
ImageReductionandAnalysisFacility
PACKAGE=ccdred
TASK=zerocombine
input=*.fitsListofzerolevelimagestocombine
(output=zero)Outputzerolevelname
(combine=average)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 处理 色散 光谱 解读