Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://lnfm1.sai.msu.ru/~kupry/SOFTWARE/prog_calib_12.3.pdf Äàòà èçìåíåíèÿ: Wed Mar 14 20:52:03 2012 Äàòà èíäåêñèðîâàíèÿ: Mon Oct 1 23:03:10 2012 Êîäèðîâêà:

Step-by-step nÀvod na redukci dat pozorovÀnÌ protuberancÌ na MFS Verze kalibracnÌho programu: 12.3 pro windows Verze nÀvodu: 1.1
(autor - Pavol Schwartz, pelozena do rustiny ­ Juri Kupjakov)

MFS : 12.3 ( Windows XP) :1.1
( - , ­ ) Vsechny adresÀe a programy pro kalibraci a uklÀdÀnÌ dat jsou v adresÀi c:\MFS_zpracovani. V podaresÀi software jsou IDL programy a soubory, kterÈ tyto programy pouzÌvajÌ. PodadresÀ tmp_in slouzÌ jako Çlozist docasnÙch soubor. V podadresÀi raw_data jsou uklÀdÀny surovÀ data pejmenovanÀ do formÀtu, kterÙ je potebnÙ pro sprÀvnou funkci kalibracnÌho software. , c:\MFS_zpracovani. software IDL , . tmp_in . raw_data . 1) ZkopÌrujeme fits soubory slit-jaw obrÀzk vsech pozorovÀnÌ z danÈho dne do podadresÀe tmp_in. 1) fts SJ tmp_in. 2) V IDL zkompilujeme a spustÌme program renameMFSfiles na pejmenovÀnÌ soubor na nÀzvy, jakÈ jsou vyzadovÀny kalibracnÌm software. Nejprve program pozaduje vybr adresÀe, kde se nachÀzejÌ soubory, kterÈ se majÌ pejmenovat vybereme .tmp_in. 2) IDL renameMFSfiles . . tmp_in. 3) Z prvnÌho menu vybereme moznost 1 ,,take all files ..." 3) 1 ,,take all files ..." 4) V dalsÌm menu vybereme moznost 1 ,,slit-jaw images" 4) 1 ,,slit-jaw images" 5) ZadÀme datum pozorovÀnÌ ve formÀtu YYMMDD; tedy napÌklad pro 18. kvtna 2011 zadÀme 110518. 5) YYMMDD ( , 18. 2011 110518).


6) VyskocÌ dialogovÈ okno pro vÙbr adresÀe, kde se ulozÌ pejmenovanÈ soubory . vybereme podadresÀ raw_data. 6) , . raw_data. 7) Pokud pejmenovÀnÌ probhlo bez problÈm, vsechny soubory v tmp_in vymazeme. 7) , tmp_in . 8) Po pekonvertovÀnÌ z BMP, vsechny snÌmky H spekter z danÈho dne ulozÌme do podadresÀe tmp_in. DbÀme na to, aby jejich vertikÀlnÌ orientace byla shodnÀ s orientacÌ slit-jaw obrÀzk. 8) *.bmp H-alpha *.fts ( bmp2fits_v2.0.pro), tmp_in ( , SJ H-alpha). 9) V IDL spustÌme program renameMFSfiles a pouzÌvÀme ho stejn, jako v pÌpad soubor slitjaw obrÀzk (body 2-6) s tÌm rozdÌlem, ze z druhÈho menu (bod 4) vybereme. moznost 2 ,,H_alpha spectrum". Jsou pejmenovÀny vsechny soubory vrÀtan dark-frames a flat-frames. Pro uklÀdÀnÌ pejmenovanÙch soubor vybereme podadresÀ raw_data, stejn jako to bylo u slit-jaw soubor. 9) IDL renameMFSfiles SJ (26). (4) 2 ,,H_alpha spectrum". ( dark-frames flat-frames), raw_data. 10)Po ÇspsnÈm pejmenovÀnÌ soubor obrÀzk H spektra, obsah podadresÀe tmp_in vymazeme. Kdyz se te podÌvÀme do podadresÀe raw_data uvidÌme tam soubory slit-jaw obrÀzk, kterÙch nÀzvy zacÌnajÌ ,,SJ_". Budou tam taky soubory spekter vrÀtan dark- a flat-frames, kterÙch nÀzvy zacÌnajÌ ,,Halp_". 10) H, tmp_in. raw_data slit-jaw , ,,SJ_" ( dark- flat-frames) ,,Halp_". 11)V adresÀi c:\MFS_zpracovani vytvoÌme podadresÀ (dÀle jej budeme nazÙvat ,,vÙstupnÌ adresÀ"), kterÈho nÀzev bude datum pozorovÀnÌ ve formÀtu YYMMDD. Do tohoto adresÀe budeme uklÀdat zkalibrovanÀ data. 11) c:\MFS_zpracovani ( ) YYMMDD. 12) SpustÌme SolarSoft pomocÌ sswidl.bat, kterÙ se nachÀzÌ na plose. Ten nastavÌ potebnÈ parametry a automaticky spustÌ IDL 12) SolarSoft sswidl.bat, . IDL . 13)V IDL zkompilujeme a spustime program mfsff, pomoci kterÈho zkonstruujeme flatfield flatfield. 13) IDL mfsff,


flatfield. 14)Na prvnÌ otÀzku ,,Are you running sswidl" odpovÌm y, ale v podstat stacÌ stisknout enter, to y je tam defaultnÌ. V prvnÌm menu zvolÌme moznost 2 ,,construct flatfield" 14) ,,Are you running sswidl" Y, enter (Y -- , y ­ ). 2 ,,construct flatfield". 15)Pokud na dalsÌ otazku ,,Take all files ending 'D.fts' from a specific directory?" odpovÌme y, tak nabhne okno, ve kterÈm vybereme adresÀ, z kterÈho budou vybrÀny vsechny soubory dark-frames (jejich nÀzvy koncÌ na D.fts). Tedy vybereme podadresÀ raw_data. Pokud odpovÌme n, nabhne okno, kde mzeme vybrat jednotlivÈ dark-frame soubory tak, ze stlacÌme ctrl a mysÌ na n klikÀme. 15) ,,Take all files ending 'D.fts' from a specific directory?" y, dark-frames (D.fts). raw_data. n, ( ctrl , . 16)DalsÌ otÀzka se tÙka vÙbru flat-frame soubor. Podobn, jak to bylo pro dark-frames v bod 15, i tady, pokud odpovÌme y, mzeme vybrat adresÀ (raw_data), ze kterÈho budou vybrÀny vsechny soubory koncÌcÌ na F.fts. Pokud odpovÌme n, flat-frame soubory vybÌrÀme manuÀln. 16) flat-frame . dark-frames 15, .. y, (raw_data), F.fts. n, flat-frame . 17)Oteve se nÀm okno se spektrem zprmrovanÙm z vybrÀnÙch flat-frame- s odectenÙm prmrnÙm dark-frame. V okn ,,Xloadct" nastavÌme ,,Gamma Correction" slider doprava tak, aby bylo dobe vidt horizontÀlnÌ pÀsy ve spektru. V okn ,,Xloadct" likneme na tlacÌtko ,,Done". 17) flat-frame dark-frame. ,,Xloadct" ,,Gamma Correction" . ,,Xloadct" ,,Done". 18)Stiskneme enter, potom klikneme postupn na oba horizontÀlnÌ vlasy. 18) enter, - . 19)Stiskneme enter. Te jsme v rezimu oznacovÀnÌ sklonu horizontÀlnÌch pÀs ve spektru. Stiskneme enter, a mzeme oznacit jeden horizontÀlnÌ pÀs v levÈ cÀsti spektra levÙm tlacÌtkem mysi. Potom oznacÌme pozici tohoto pÀsu v pravÈ cÀsti spektra levÙm tlacÌtkem mysi. DalsÌ pÀsy mzeme takto oznacit vzdy po stisknutÌ enter. Pokud po stisknutÌ enter klikneme do okna se spektrem pravÙm tlacÌtkem mysi, vyjdeme z rezimu oznacovÀnÌ vodorovnÌch pÀs a program pokracuje automaticky vÙpoctem Çhlu sklonu spektra a spektrum o tento sklon koriguje. 19) enter. . . enter . . enter. enter, ,


. 20)Te se bude urcovat zakivenÌ spektra. Stiskneme enter a levÙm tlacÌtkem mysi oznacÌme levÙ dolnÌ a nÀsledn pravÙ hornÌ roh obdÈlnÌkovÈ oblasti, kterÀ by mla obsahovat dobe viditelnou Çzkou spektrÀlnÌ caru, ze kterÈ bude program schopen urcit automaticky jejÌ zakivenÌ. Po vÙbru oblasti, program nakreslÌ graf relativnÌho posunu spektrÀlnÌ cÀry ve smru disperze v pixlech v zÀvislosti na poloze na strbine. Protoze posun je urcen v celÙch pixlech, je poteba to vyhladit (anglicky smooth), abychom to mohli fitovat. ZadÀm tzv. smooth factor, mÀm odzkousenÈ, ze 7 stacÌ, aby se to dalo dobe nafitovat polynomem 3. stupn. 20) . enter , . , . .. (ang. smooth), fitovat ( ). .. smooth factor, 7, 3 . 21)Pokud jsme spokojenÌ s vyhlazenÌm, odpovÌme y a zacÌname s fitovÀnÌm. Te mzeme zadat startovacÌ hodnoty koeficient polynomu. MÀm odzkousenÈ, ze z defaultnÌch startovacÌch hodnot to vzdy dobe konvergovalo, tedy stacÌ stisknout 4 krÀt enter. Te se bude program ptÀt u kazdÈho koeficientu, jestli ho chceme zahrnout do fitovÀnÌ, anebo ho fixovat. Jelikoz chceme zahrnout vsechny koeficienty, znova stiskneme 4 krÀt enter. 21) , y . . , ( ) 4 enter. . , 4 enter. 22)SpustÌ se automatickÈ fitovÀnÌ a po jeho skoncenÌ (zlomek sekundy), program oteve okno z cervenÙmi bodmi a zelenÙm fitem. Pokud je fit dobry (coz bylo zatÌm vzdy), nenÌ poteba pokracovat v iteracÌch . stiskneme n. Program koriguje zakivenÌ ve flat-frame , chvÌli to trvÀ. 22) ( ), . ( ), , n. flat-frame ( ). 23)Az program napÌse ,,Done" a ukÀze flat-frame korigovanÙ o zakivenÌ, stiskneme enter a program nÀm po chvili zobrazÌ tzv. hard and soft flat-field. 23) ,,Done" flat-frame enter .. hard soft flat-field. 24)Po dalsÌm stisknutÌ enter, se ukÀze jeste tzv. slit flat. A po dalsim stisku enter se ukÀzou uz jenom dva flatfieldy: StatickÙ, kterÙ bude pro vsechny spektra ten samÙ a tzv. slit flat, kterÙ mze bÙt rzn vertikÀln posunutÙ pro jednotlivÀ spektra. 24) enter .. slit flat. enter


2 flatfiel(): , slit flat, . 25)Program se zeptÀ, jestli chceme ulozit vÙslednÈ flatfieldovÈ matice, odpovÌme y. Na otÀzku, jestli chceme zmnit defaultnÌ adresÀ c:\Program Files\IDL64 na uklÀdÀnÌ soubor, odpovÌme y a vybereme ,,vÙstupnÌ adresÀ", aby se do nj ten idl-save soubor s flatfieldovÙmi maticemi ulozil. Hned potom zadÀme nÀzev souboru bez koncovky, do kterÈho se flat-fieldovÈ matice ulozÌ. K nÀzvu souboru se automaticky pidÀ koncovka .idl. Chtl bych, aby se ten soubor s flat-fieldama jmenoval FF_YYMMDD.idl, kde YY je rok, MM msÌc a DD je den. Tedy pro pozorovÀnÌ z 18.5.2011 by se soubor s flat-fieldama ml jmenovat FF_110518.idl 25) , flatfield , y. , c:\Program Files\IDL64 , y « » , idl-save . . idl. FF_YYMMDD.idl, YY , MM DD . 18.5.2011 flat-field() FF_110518.idl. 26)Po konstrukci flatfieldu mzeme zacÌt redukovat spektra pozorovÀnÌ protuberance. Znovu spustÌme program mfsff a v prvnÌm menu vybereme moznost 3 ,,use flatfield". 26) flatfield . mfsff a 3 ,,use flatfield". 27)Vybereme idl-save soubor s flatfieldama . ml by bÙt ve ,,vÙstupnÌm" adresÀi. 27) idl-save ( 110607). 28)Vybereme s podadresÀe raw_data fits soubor se spektrem, kterÈ chceme zredukovat. 28) raw_data fits , . 29)V okn ,,Xloadct" posuneme slider Gamma Correction do leva tak, aby bylo dobe vidt oba horizontÀlnÌ vlasy. Klikneme na tlacÌtko ,,Done". 29) ,,Xloadct" Gamma Correction , . ,,Done". 30)Stiskneme enter a klikneme postupn na oba horizontÀlnÌ vlasy. 30) enter . 31)Stiskneme enter a cekÀm az program zkoriguje spektrum na sklon a zakivenÌ . trvÀ to nekolik desÌtek sekund. 31) enter ( ). 32)OznacÌme oblast na disku, kterÀ se vyuzije na urcenÌ vertikÀlnÌho posunu slit-matice. PrvnÌm


kliknutÌm levÈho tlacÌtka mysi oznacÌme jednu hranici a druhÙm kliknutÌm druhou hranici tÈto oblasti. Je nezÀdoucÌ, aby touto oblastÌ prochÀzel kterÙkoliv z horizontÀlnÌch vlas. 32) , slit- . , . - . 33)Program nakreslÌ graf, kde jsou vykresleny pmrnÈ intenzity podÈl strbiny z vybranÈ oblasti a vertikÀlnÌ prbeh slit-flat matice. JejÌ posun jiz program urcil automaticky, mzeme jej jeste manuÀln zmnit, pokud na otÀzku ,,Are you satisfied?" odpvÌme n. Pokud tedy odpovÌme n, mzeme mnit posun . odpovÌme znova n a mzeme zadat novÑu hodnotu posunu. Po kazdÈ odpovedi ,,n", zadÀvÀme novou hodnotu posunu. Kdyz jsme jiz spokojeni, odpovÌme y. Program se zeptÀ ,,Apply the drift?" jestli uplatnit na slit flatfield posun. Pokud se nÀm podailo najÌt takovÙ posun, ze prbh prmrnÈ intenzity podÈl strbiny dobe koreloval s vertikÀlnÌm prbhem slit-flat matice, odpovÌme y. Pokud ne, odpovÌme n; skutecn i tato moznost se mi stala napr. u pozorovÀnÌ z 18.4., kde byly spektra dost zasumnÈ. 33) , slit-flat . , , , ,,Are you satisfied?" n. n ­ . ,,n" ­ . y. ,,Apply the drift?" slit flatfield ? , slitflat , y. , n ( 18.4, ). 34)CekÀme, az program vypocte vÙslednÙ flatfield a uplatnÌ ho na spektrum. Az napÌse ,,Done" je to hotovo a program oteve ti okna: v jednom je spektrum ped flatfieldovanÌm, ve druhÈm spektrum po uplatnenÌ flatfieldu a ve tetÌm samotnÌ flatfield. 34) flatfield . ,,Done" ­ 3 : flatfiel(), flatfiel() flatfield. 35)Stiskneme enter zobrazÌ se prmrnÙ profil z disku a porovnÀvacÌ ,,sample" profil. Pokud. je pozorovanÙ profil obrÀcenÙ ve vlnovÈ dÈlce oproti ,,sample" profilu, odpovÌme y a pozorovanÙ profil a spektrum se pevrÀtÌ ve vlnovÈ dÈlce. Pokud je po pevrÀcenÌ pozorovanÙ profil sprÀvn, na otÀzku ,,Are you satisfied?" odpovÌme y. 35) enter, ,,sample" . ,,sample" , y, . , ,,Are you satisfied?" y. 36)Na otÀzku ,,Estimate dispersion and create wvl vector?" odpovÌme y. OznacÌme nkolik spektrÀlnÌch car levÙm tlacÌtkem mysi v okn ,,avg observed spectrum". OznacovanÌ ukoncÌme kliknutÌm do okna pravÙm tlacÌtkem mysi. Stiskneme enter, a ty samÈ spektrÀlnÌ cÀry oznacÌme v okn ,,sample avg disc spectrum". Nabhne okno s kalibracnÌ pimkou pro vlnovou dÈlku. 36) ,,Estimate dispersion and create wvl vector?" y.


,,avg observed spectrum". . enter ,,sample avg disc spectrum". . 37)Stiskneme enter a na otÀzku ,,Make absolute calibration ..." odpovÌme y. Stiskneme znovu enter a v okn ,,flatfielded raw spectra" oznacÌme levÙm tlacÌtkem mysi ob hranice disku ve spektru. 37) enter ,,Make absolute calibration ..." y. enter ,,flatfielded raw spectra" . 38)Nabhne okno, kde vybereme fits soubor se slit-jaw komplementÀrnÌm ke zpracovÀvanÈmu spektru. 38) slit-jaw fits . 39)V okn ,,slit-jaw image" oznacÌme mysÌ postupne oba prsecnÌky strbiny s horizontÀlnÌmi vlasami. 39) ,,slit-jaw image" . 40)Stiskneme enter, a v okn ,,slit-jaw image" oznacÌme levÙm tlacÌtkem mysi nejmÈn ctyy body na limbu. NenÌ zas dobrÈ oznacit tch bod vÌce nez pt az sest. OznacovÀnÌ ukoncÌme kliknutÌm do okna pravÙm tlacÌtkem mysi a program nafituje limb cÀstÌ kruznice. Pokud program fitne limb dobe, odpovÌme y na otÀzku ,,Are you satisfied?". Kdyz je fitovÀnÌ limbu spatn, odpovÌme n a vrÀtime se k bodu 39 a vyzkousÌme jinÈ rozlozenÌ bod na limbu. 40) enter ,,slit-jaw image" 4 . - 5 6. . , y ,,Are you satisfied?". , n 39 . 41)Klikneme nkde do cÀsti prinÀlezÌcÌ disku ve spektru v okn ,,flatfielded raw spectra ­ real relative intensities". Profil z tÈto pozice se pouzije na kalibraci pomocÌ profilu Davida. NeklikÀme pÌlis blÌzko limbu, protoze pro tyto pozice nemÀme Davidv profil. 41) , , ,,flatfielded raw spectra ­ real relative intensities". Davida. , .. Davida. 42)Stiskneme enter, a oznacÌme dva body nad limbem, kterÈ budou slouzit na odhad rozptÙlenÈho svtla mimo disku (off-limb scattered light). 42) enter 2 , (off-limb scattered light). 43)Stiskneme enter, vyskocÌ okno s profilem z pozice na disku vybranÈ v bod 41, klikneme na


sted tohto profilu. 43) enter 41, . 44)Stiskneme enter, a oznacÌme na profilu Çsek v kontinuu (mimo jakÙchkoliv spektrÀlnÌch car) nalevo od H, stiskneme enter a oznacÌme Çsek v kontinuu napravo od H. Tyto dva Çseky budou slouzit na odhad off-limb rozptÙlenÈho svtla. Stiskneme enter a oznacÌme oblast H cÀry. 44) enter ( ) H, enter H. off-limb . enter H. 45)Stiskneme enter a vyskocÌ okno ,,tabulated profile", kde je zobrazen cerven Davidv profil pro ve vybranÈ pozici na disku a profil z tÈto pozice. KlikÀnÌm levÙm tlacÌtkem mysi vybÌrÀme body v profilu, kterÈ budou pouzity na konstrukci kalibracnÌ pÌmky. VÙbr bod ukoncÌme kliknutÌm do okna pravÙm tlacÌtkem mysi. Nabhne okno kalibrace jest pe odstrannÌm rozptÙlenÈho svtla a zapocÌtanÌm Zadigovho filtru. KalibracnÌ pÌmka neprochÀzÌ bodem [0,0], coz je zpsobeho rozptÙlenÙm svetlem . Z tÈto kalibracnÌ kivky program vypocÌtÀ intenzitu rozptÙlenÈho svtla na disku. 45) enter ,,tabulated profile", Davida . . . Zadigovho . [0,0] . . 46)Na otÀzku ,,Is a filter put at the disk ..." odpovÌme y. Pokud byl filter nasazen na limb pesn, odpovÌme i na dalsÌ otÀzku y. Pokud ne, odpovÌme na ni n a nÀsledn jsme vyzvÀni programem, aby jsme ob hranice vÙskitu filtru oznacili v spektru levÙm tlacÌtkem mysi. 46) Is a filter put at the disk ..." y. , y. , n . 47)Z menu ,,Filter transmitivity" vybereme propustnost pouzitÈho Zadigovho filtru, tedy pro 12.7% zadÀme moznost 2. ZobrazÌ se nÀm spektrum v reÀlnÙch intenzitÀch po zapocÌtÀnÌ propustnosti filtr. Protuberanci jiz v spektru nevidÌme, protoze intenzita v H je omnoho nizsÌ nez skutecnÀ intenzita disku. Pokud chceme vidt disk i protuberanci, odpovÌme y na otÀzku ,,View well both filtered ...?" a disk se zobrazÌ v modrÈ a protuberance v cervenÈ barv. Pokud jsme spokojenÌ s nastavenÌm pozic disku a filtru odpovÌme y na otÀzku ,,Are you satisfied?" Pokud odpovÌme n, vrÀtÌme se do bodu 37. 47) ,,Filter transmitivity" Zadigovho , 12.7% 2. . .. H . , y ,,View well both filtered ...?" , . y ,,Are you satisfied?". , n 37.


48)Stiskneme enter, VypÌse se kalibracnÌ koeficient, kterÙ byl uz vypocten po odcÌtÀnÌ rozptÙlenÈho svtla na disku a zahrnutÌ propustnosti Zadigovho filtru. 48) enter, , Zadigovho . 49)Na otÀzku ,,Take into account variations ...?" . jestli chceme vzÌt do Çvahy zmny kalibracnÌho koeficientu s vlnovou dÈlkou, odpovÌme y. Vyskoci okno, kde program nafitoval odchÙlky kalibrovanÈho vybranÈho profilu od Davidova profilu kvadratickou funkcÌ a vypÌse jejÌ parametry fitu s odchÙlkami v procentech. Pokud je fit dobrÙ odpovÌme y. Ale pokud fit nenÌ dobrÙ, nebo odchÙlka parametru a je velkÀ (nad 50%), odpovÌme n a program udlÀ lineÀrnÌ fit. Pokud je fit dobrÙ, odpovÌme y. Pokud ani lineÀrnÌ fit nenÌ dobrÙ, nebo odchÙlka parametru a je vtsÌ nez 50%, odpovÌme n a tÌm pÀdem se zmny kalibracnÌho koeficientu s vlnovou dÈlkou neberou v Çvahu. Z mÈ zkusenosti, lineÀrnÌ fit je dostacujÌcÌ, v nkterÙch pÌpadech byla dokonce zÀvislost kalibracnÈho koeficientu na vlnovÈ dÈlce zanedbatelnÀ. 49) ,,Take into account variations ...?" y. - Davidova . , y. , 50% n . , y. 50% n . , . 50)Pokud vsechny korekce a kalibrace probhly v poÀdku, mzeme zkalibrovanÈ spektrum ulozit. Na dotaz ,,Save the flatfielded spektrum into file?" odpovÌme y 50) , . ,,Save the flatfielded spektrum into file?" y. 51)I na dalsi otÀzku odpovÌme y a vyskocÌ okno, ve kterÈm vybereme, kde se mÀ zredukovanÈ spektrum a slit-jaw ulozit vybereme .,,vÙstupnÌ adresÀ". Vybereme moznost 2 ,,fits files with standard header" 51) y. slit-jaw. 2 ,,fits files with standard header". 52)Do ,,vÙstupnÌho adresÀe" se ulozili ti fits soubory: slit-jaw se souadnicovÙm systÈmem a pozicemi prsecnÌk horizontÀlnÌch vlas se strbinou, flatfieldovanÈ spektrum s intenzitami v count-ech a flatfieldovanÈ spektrum s intenzitmi v erg/cm2/s/sr/Hz. 52) 3 fits : slit-jaw , flatfield() flatfield() [erg/cm2 /s/sr/Hz]. 53)Pro kalibraci dalsÌho spektra z danÈho pozorovÀnÌ, cteme tento nÀvod od bodu 26


53) , 26. S pÌpadnÙmi dotazmi a/anebo problÈmami piste na pschwartz@astro.sk S pozdravem Pavol Schwartz pschwartz@astro.sk