1、CCD狹縫光譜觀測與數據處理北京興隆2.16m光學塑遠鏡，1024X240gam=0865,rdnoise=3.756觀測給出的數據是FIT后綴的，用rfklist-LLnli轉化.然后用readlog.cl寫頭文件信.,dos2unixreadlog.ctaskreadlog=c7user2/hbian/readlog.cl更改observatory.paq/ixafAraf/noao/astutil/人力II入baoextinct.dat（Araf/irafVnoao/lib/onedstsds/）觀測：CCD校準圖象的拍攝：BiasFrames（零秒喑場）10幅DarkFrames（喑場

2、）5幅與當晚最長暴光時間相同（？）Flat-FieldExposure（圓頂平場）510幅保證很高的光電子數（10,000e/pLxel.以上）BvilightFlat（晨昏蒙影犬光平場）35幅用作狹縫照明改正光譜圖象的拍攝：比較光譜（波長定標燈光譜）用作波長定標分光流螢標準星光譜用作流量定標待測星光譜光譜最佳拍攝組合：。一定標燈一標準星一定標燈一待測星一定標燈一。IRAF光譜處理步驟CCD圖彖預處理：主要在noao.Liiired.ccdred包中進行一、設置文件頭信息和處理壞境1確定EIASSEC與TRIMSEC：使用Linexainme或lmplot檢查一幅平場圖彖設置文件頭信息:使用h

3、edit（用readlog批處理）IMAGETYPEzero/dark/flat/object/compEXPTIMEGAINRDNOISEEIASSECTRIMSEC設置處理環境setmstrumentspecphot1）ccdred參數pixeltv=realreal2）ccdproc參數biassec=trimsec=1:1024,3:238（*180:10243:238）4確認文件頭信息被IRAF正確識別：使用ccdhst二、Trim,Overscan和Bias的改正1合并Bias（生成Zero）:使用zerocombmecombine=averagereject=nummxccdty

4、pe=zeronlow=0LlllLgh=1rdnoise=gam=2Tnm、Overscan減Bias：使用ccdprocccdtype=overscan=notnm=vesJzerocor=vesJdarkcor=noflator=nozero=Zero（合并后的Bias）overscan交互式,mteractiveXyes/no）yes:order1或2fq三、Dcuk改正（可忽略丿1、檢查Dark是否顯著，如果不顯著則不作Dark改正:使用lmexamine/implot/imstatistics2、合并Dad合并圓頂（生成Flat）:使用flatconibinecombine=ave

5、ragereject=crrejectccdtype=flatsubsets=yesscale=modenlow=1nhigh=1rdnoise=gain=2、合并晨昏蒙影天光平場（生成Sky）:使用combinecombine=averagereject=crrejectsubsets=yesscale=mode3、圓頂平場的歸一（生成iiFlat）：使用response（noao.twodspec.longslit）calibrate=Flat（Flat為步驟1中合并的平場）normalize=Flat（同上）response=nFlatinteractive=ves交互如下：:funct

6、ionspline3:order25fkjq4、對合并后的天光平場做平場改正：使用ccdprocflatcor=yesflat=nFlat（歸一后的平場）5、檢查狹縫照明的均勻性使用imexamine/iniplot檢查步驟4,生成的天光平場圖彖沿狹縫方向是否均勻，若不均勻性顯著則進行以卜步驟。6、構造狹縫照明函數（生成nsky）：使用Lllumiiiationnbins=5interactive=vesJ交互如下functionorder1sfq7、生成最終的平場校準圖彖（生成pFlat）：使用umrithimanthnFlat*nSkypFlat8、除平場：使用ccdprocflatcor

7、=yesflat=pFlat五、剔除野點（針對光譜圖象）1、修補壞彖元：使用ccdprocfixpix=yesfixfile=badpixel（壞彖元坐標文件）2、剔除宇宙線點：使用cosmicrays/imedit/epixH、光譜的抽取和定標noao.Lmred.kpnoslitnoao.imred.ctioslitnoao丄mred.kpnocoudenao.nvospec一、檢查光譜圖彖1、確定恒星光譜的空間輪廓（狹縫方向）：使用implot確定輪廓人小fXvlimbase-to-base-width2、使用實際光譜分辨率：使用implot確定定標燈光譜譜線寬度fXvlimbase-t

8、o-base-width3、檢查狹縫與CCD行或列的平行性：使用implot對定標燈作多行重疊劃線二、抽取光譜：使用apall1、設置任務參數1）孔徑（apenue）參數line=INDEFnsum=rvLdrhLdth=（空間輪廓寬度）loxer=-O.5*xvidthupper=0.5*widthresize=no2）背景（Background）參數b_sample=b_naver=-100b_function=chebvshevb_order=或23）追跡（Trace）參數t_sumt_setuptjfunct=legendre/spline3t_order=2t_mter=1/2/34

9、）譜求和與減天光參數background=fitweiglit=varianceclean=vesJsaturation=32400readnoise=gain=Lsigma=4hsigirn=45）輸出方式參數format=multispecextras=vesJ2、抽取光譜；運行apaL分強譜與弱譜兩種A、強譜情形1）修改孔徑位置與人小并定孔徑中心findapertures?Yeseditaperture?Yes:upperlowerline:widthcmdIub（進入卜一步）2）修改背景窗II位置與人小并擬合背景sampleordertsf丫q（進入卜一步）3）修改追跡參數并追蹤fun

10、ctionspline3order3dufq（進入卜一步）4）抽取光譜?yesq5）檢查所抽取譜：使用splot三、波長定標1、抽取比較光譜apalcomplout=cobjlref=objlrecen-trace-backmteract-2、確定色散關系：使用identifycoordli=liiielistsJctLohear.datfwidth=（譜線寬度）niterat=1D用光標證認出3、4根線m輸入波長值2）初步擬合色散關系:functioncheb:order2fdjq3）自動尋找譜線14）擬合新的色散關系:functionspline3cheb:order345）檢査結果mis

11、?是否再次找線、擬合？3、確定其他比較光譜的色散關系：使用reidentifyreference=cobjl.msimages=cobj2.msinteract=vesJ4、待測譜與定標譜的匹配：可采取以I、三種方式之一1）直接寫頭文件：使用heditheditobjl.ms.iinhREFSPECl=”cobjl.ms”add+shov+（個比較譜情形）或heditobjl.ms.LinhREFSPECl=cobjla.ms0.5add+ver-show+heditobjl.ms.LinhREFSPECl=cobjlb.ms0.5”add+ver-show+（二個比較譜情形）2）按照觀測時間

12、匹配：使用serjd和refspecsetjd*.imhrefspecobj*.ms.iinhreference=cobj*.ms丄mhselect=uiterpsort=jdconfirm-verb-3）按照列表方式匹配editreftableobjl.mscobjl.msobj2.mscobj2a.mscobj2b.msrefspecobj*.ms.nnhref=cobj*.Lns.imhrefe=reftableselectivevei:bsort=none5、待測譜波長定標：使用dispcor,有四個參數可改變startingwavelengthenduigwavelengthwavelengthperpixelnumberofoutputpoxels四、流量定標1、計算大氣質量并寫入文件頭：使用setairmass2、設置kpnoslit包參數（epakpnoslit）extinct=baoextuict.dat（人/（，消光表，onedstdsjkpnoextmct.dat）cal