第 4 章 データ整約実習 27
4.2.2 A-B
A B imarith
“ab”
−−> ! s e d ’ s /\( .∗ \) / ab\1 / ’ c r o b j 1 a . l s t > abobj1a . l s t
−−> i m a r i t h @crobj1a . l s t − @crobj1b . l s t @abobj1a . l s t
4.2.3 ドームフラットの作成とフラットフィールド
imcombine
i n p u t = L i s t o f images t o combine
output = L i s t o f output images
( h e a d e r s = ” ” ) L i s t o f h e a d e r f i l e s ( o p t i o n a l ) ( bpmasks = ” ” ) L i s t o f bad p i x e l masks ( o p t i o n a l ) ( r e j m a s k s = ” ” ) L i s t o f r e j e c t i o n masks ( o p t i o n a l )
( n r e j m a s k s = ” ” ) L i s t o f number r e j e c t e d masks ( o p t i o n a l ) ( expmasks = ” ” ) L i s t o f e x p o s u r e masks ( o p t i o n a l )
( s i gm as = ” ” ) L i s t o f sigma images ( o p t i o n a l ) ( imcmb = ” $ I ” ) Keyword f o r IMCMB keywords
4.2. 33 ( l o g f i l e = ”STDOUT” ) Log f i l e
( combine = ” median ” ) Type o f combine o p e r a t i o n ( r e j e c t = ” s i g c l i p ” ) Type o f r e j e c t i o n
( p r o j e c t = no ) P r o j e c t h i g h e s t d i m e n s i on o f i n p u t images ? ( o u t t y p e = ” r e a l ” ) Output image p i x e l d a t a t y p e
( o u t l i m i t s = ” ” ) Output l i m i t s ( x1 x2 y1 y2 . . . ) ( o f f s e t s = ” none ” ) Input image o f f s e t s
( masktype = ” none ” ) Mask type ( maskvalue = ” 0 ” ) Mask v a l u e
( blank = 0 . 0 ) Value i f t h e r e a r e no p i x e l s ( s c a l e = ” e x p o s u r e ” ) Image s c a l i n g
( z e r o = ” none ” ) Image z e r o p o i n t o f f s e t ( w e i g h t = ” none ” ) Image w e i g h t s
( s t a t s e c = ” ” ) Image s e c t i o n f o r computing s t a t i s t i c s ( expname = ”EXPTIME” ) Image h e a d e r e x p o s u r e time keyword ( l t h r e s h o l d = INDEF) Lower t h r e s h o l d
( h t h r e s h o l d = INDEF) Upper t h r e s h o l d
( nlow = 1 ) minmax : Number o f low p i x e l s t o r e j e c t ( nhigh = 1 ) minmax : Number o f h i g h p i x e l s t o r e j e c t
( nkeep = 1 ) Minimum t o keep ( pos ) o r maximum t o r e j e c t ( neg ) ( m c l i p = y e s ) Use median i n sigma c l i p p i n g a l g o r i t h m s ?
( l s i g m a = 3 . 0 ) Lower sigma c l i p p i n g f a c t o r ( hsigma = 3 . 0 ) Upper sigma c l i p p i n g f a c t o r
( r d n o i s e = ” 0 . ” ) c c d c l i p : CCD r e a d o u t n o i s e ( e l e c t r o n s ) ( g a i n = ” 1 . ” ) c c d c l i p : CCD g a i n ( e l e c t r o n s /DN)
( s n o i s e = ” 0 . ” ) c c d c l i p : S e n s i t i v i t y n o i s e ( f r a c t i o n )
( s i g s c a l e = 0 . 1 ) T o l e r a n c e f o r sigma c l i p p i n g s c a l i n g c o r r e c t i o n s ( p c l i p = −0.5) p c l i p : P e r c e n t i l e c l i p p i n g parameter
( grow = 0 . 0 ) Radius ( p i x e l s ) f o r n e i g h b o r r e j e c t i o n imcombine
−−> imcombine @ f l a t 1 . l s t HK500 CDFN2 Domeflat1 . f i t s
1
−−> i m a r i t h HK500 CDFN2 Domeflat1 . f i t s / 1 0 0 0 0 . HK500 CDFN2 Domeflat1 . f i t s
34 4 fixpix
−−> f i x p i x HK500 CDFN2 Domeflat1 . f i t s mdpdb$bpm/ nlbpm1 FF64r . f i t s A-B
“fl”
−−> ! s e d ’ s /\( .∗ \) / f l\1 / ’ abobj1a . l s t > f l o b j 1 a . l s t
−−> i m a r i t h @abobj1a . l s t / HK500 CDFN2 Domeflat1 . f i t s @ f l o b j 1 a . l s t
4.2.4 ゆがみ補正
geotran
i n p u t = Input data
output = Output data
d a t a b a s e = Name o f GEOMAP d a t a b a s e f i l e t r a n s f o r m s = Names o f c o o r d i n a t e t r a n s f o r m s i n
d a t a b a s e f i l e
( geometry = ” g e o m e t r i c ” ) T r a n s f o r m a t i o n type ( l i n e a r , g e o m e t r i c ) ( x i n = INDEF) X o r i g i n o f i n p u t frame i n p i x e l s ( y i n = INDEF) Y o r i g i n o f i n p u t frame i n p i x e l s ( x s h i f t = INDEF) X o r i g i n s h i f t i n p i x e l s
( y s h i f t = INDEF) Y o r i g i n s h i f t i n p i x e l s ( xout = INDEF) X o r i g i n o f output frame i n
r e f e r e n c e u n i t s
( yout = INDEF) Y o r i g i n o f output frame i n r e f e r e n c e u n i t s
( xmag = INDEF) X s c a l e o f i n p u t p i c t u r e i n p i x e l s p e r r e f e r e n c e u n i t
( ymag = INDEF) Y s c a l e o f i n p u t p i c t u r e i n p i x e l s p e r r e f e r e n c e u n i t
( x r o t a t i o n = INDEF) X a x i s r o t a t i o n i n d e g r e e s ( y r o t a t i o n = INDEF) Y a x i s r o t a t i o n i n d e g r e e s
( xmin = INDEF) Minimum r e f e r e n c e x v a l u e o f output p i c t u r e
( xmax = INDEF) Maximum r e f e r e n c e x v a l u e o f output p i c t u r e
( ymin = INDEF) Minimum r e f e r e n c e y v a l u e o f output p i c t u r e
4.3. 35 ( ymax = INDEF) Maximum r e f e r e n c e y v a l u e o f output
p i c t u r e
( x s c a l e = 1 . 0 ) X s c a l e o f output p i c t u r e i n r e f e r e n c e u n i t s p e r p i x e l
( y s c a l e = 1 . 0 ) Y s c a l e o f output p i c t u r e i n r e f e r e n c e u n i t s p e r p i x e l
( n c o l s = INDEF) Number o f columns i n t h e output p i c t u r e ( n l i n e s = INDEF) Number o f l i n e s i n t h e output p i c t u r e ( xsample = 1 . 0 ) C o o r d i n a t e s u r f a c e s am p l i n g i n t e r v a l
i n x
( ysample = 1 . 0 ) C o o r d i n a t e s u r f a c e s am p l i n g i n t e r v a l i n y
( i n t e r p o l a n t = ” l i n e a r ” ) I n t e r p o l a n t
( boundary = ” c o n s t a n t ” ) Boundary e x t e n s i o n ( n e a r e s t , c o n s t a n t , r e f l e c t , wrap )
( c o n s t a n t = 0 . 0 ) Constant boundary e x t e n s i o n ( f l u x c o n s e r v e = y e s ) P r e s e r v e image f l u x ?
( n x b l o c k = 5 1 2 ) X d i m e n s i on o f working b l o c k s i z e i n p i x e l s
( n y b l o c k = 5 1 2 ) Y d i m e n s i on o f working b l o c k s i z e i n p i x e l s
( v e r b o s e = y e s ) P r i n t mes s ages about t h e p r o g r e s s o f t h e t a s k
geotran “gc”
−−> ! s e d ’ s /\( .∗ \) / gc\1 / ’ f l o b j 1 a . l s t > g c o b j 1 a . l s t
−−> g e o t r a n @ f l o b j 1 a . l s t @gcobj1a . l s t
mdpdb$geomap/ m c s d i s t c r r 1 f e b 0 7 n e w . dbs m c s d i s t c r r 1 f e b 0 7 n e w . gmp
geotran MDPDB
mdpdb$geomap/mcsdistcrr1 feb07new.dbs mcsdistcrr1 feb07new.gmp 1 md-pdb$geomap/mcsdistcrr2 feb07new.dbs mcsdistcrr2 feb07new.gmp 2 database transform 3
4.3 データの切り出し
3 MCSRED
36 4
MCSMDP maskplot MDP
−−> maskplot CDFN MASK02 . mdp image=gcflabcrMCSA00057147 . f i t s raw+
MODS11-0390 imcopy
−−> ! s e d ’ s /\( .∗ \)\. f i t s /\1 MODS11−0390\. f i t s / ’ g c o b j 1 a . l s t >
gcMODS11−0390. l s t
−−> ! s e d ’ s /\( .∗ \) /\1 [∗, 1 7 5 5 : 1 8 4 0 ] / ’ g c o b j 1 a . l s t > c u t . l s t
−−> imcopy @cut . l s t @gcMODS11−0390. l s t
4.4 個別処理
4.4.1 波長較正
X OH A-B
−−> ! s e d ’ s /\( .∗ \) / f l s k y\1 / ’ c r o b j 1 a . l s t > f l s k y 1 a . l s t
−−> i m a r i t h @crobj1a . l s t / HK500 CDFN2 Domeflat1 . f i t s @ f l s k y 1 a . l s t
−−> ! s e d ’ s /\( .∗ \) / gc\1 / ’ f l s k y 1 a . l s t > gc s k y 1a . l s t
−−> g e o t r a n @ f l s k y 1 a . l s t @gcsky1a . l s t
mdpdb$geomap/ m c s d i s t c r r 2 f e b 0 7 n e w . dbs m c s d i s t c r r 2 f e b 0 7 n e w . gmp
−−> ! s e d ’ s /\( .∗ \)\. f i t s /\1 MODS11−0390\. f i t s / ’ gc s k y 1a . l s t
> gcskyMODS11−0390. l s t
−−> ! s e d ’ s /\( .∗ \) /\1 [∗, 1 7 5 5 : 1 8 4 0 ] / ’ gc s k y 1a . l s t > c u t . l s t
−−> imcopy @cut . l s t @gcskyMODS11−0390. l s t
identify identify
X
images = Images c o n t a i n i n g f e a t u r e s t o be i d e n t i f i e d
4.4. 37 c r v a l = Approximate c o o r d i n a t e ( a t r e f e r e n c e
p i x e l )
c d e l t = Approximate d i s p e r s i o n
( s e c t i o n = ” middle l i n e ” ) S e c t i o n t o apply t o two d i m e n s i o n a l images
( d a t a b a s e = ” d a t a b a s e ” ) Database i n which t o r e c o r d f e a t u r e data
( c o o r d l i s t = ” mdpdb$ohlist / l i s t N S H K 500 ” ) User c o o r d i n a t e l i s t ( u n i t s = ” ” ) C o o r d i n a t e u n i t s
( nsum = ” 2 0 ” ) Number o f l i n e s / columns / bands t o sum
i n 2D images
( match = −3.0) C o o r d i n a t e l i s t matching l i m i t ( m a x f e a t u r e s = 5 0 ) Maximum number o f f e a t u r e s f o r
a u t o m a t i c i d e n t i f i c a t i o n
( zwidth = 1 0 0 . 0 ) Zoom graph width i n u s e r u n i t s ( f t y p e = ” e m i s s i o n ” ) F e a t u r e type
( f w i d t h = 8 . 0 ) F e a t u r e width i n p i x e l s ( c r a d i u s = 5 . 0 ) C e n t e r i n g r a d i u s i n p i x e l s ( t h r e s h o l d = 0 . 0 ) F e a t u r e t h r e s h o l d f o r c e n t e r i n g
( minsep = 2 . 0 ) Minimum p i x e l s e p a r a t i o n ( f u n c t i o n = ” chebyshev ” ) C o o r d i n a t e f u n c t i o n
( o r d e r = 4 ) Order o f c o o r d i n a t e f u n c t i o n ( sample = ”∗” ) C o o r d i n a t e sample r e g i o n s ( n i t e r a t e = 1 0 ) R e j e c t i o n i t e r a t i o n s ( l o w r e j e c t = 3 . 0 ) Lower r e j e c t i o n sigma ( h i g h r e j e c t = 3 . 0 ) Upper r e j e c t i o n sigma
( grow = 0 . 0 ) R e j e c t i o n growing r a d i u s
( a u t o w r i t e = no ) A u t o m a t i c a l l y w r i t e t o d a t a b a s e ( g r a p h i c s = ” s t d g r a p h ” ) G r aphics output d e v i c e
( c u r s o r = ” ” ) G r aphics c u r s o r i n p u t
( a i d p a r s = ” ” ) Automatic i d e n t i f i c a t i o n a l g o r i t h m p a r a m e t e r s
section nsum coordlist
OH Rousselot et al. (2000) [3] HK500
MCSMDP 4.1
identify
−−> i d e n t i f y gcflskycrMCSA00057147 MODS11−0390. f i t s
38 4
1400016000180002000022000Angstrom 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 x1e+4
14356.719 14518.961
14698.437
14887.699
15055.55 15187.14 15240.954 15287.789 15332.403 15432.156 15540.328 15656.963
15833.272 15972.596 16030.832 16128.609 16235.377
16502.365
16692.38
16903.679 16955.078 17008.757 17123.659 17247.926 17386.696 17449.967
17653.222
17880.299 17993.962 18118.494
19250.306 19350.119
19771.862
20008.163
20275.839 20339.497 20412.68 20499.364 20563.548
20729.015
20909.569
21115.857 21176.557 21249.592
21507.308
21711.17 21802.312 21873.518 21955.638 22052.366 22124.875
22311.799
22517.969
OH line list for HK500 grism (mdpdb$ohlist/list_NS_HK500)
4.1: HK500 OH
4.4. 39
4.2 X
20
4.1
4.1: identify
m d c w f l
q identify
? b t
x X
y Y
?
“m” 4.1 5
“m” 4
5 “f” 4.3
residual 3 Chebychev
5
identify iraf
4.2
4H 4 1
X
5order=4 iraf order order-1
40 4
4.2: identify
4.2: identify
h,i,j,k,l f d u c s z
:order n n-1
:function
:niterate sigma rejection :low reject sigma rejection :high reject sigma rejection
q
?
4.4. 41
4.3: identify
“q”
“l”
“f” HK500
7.72 [˚A/pixel] RMS ˚A “q”
“q” identify
Write f e a t u r e data t o t h e d a t a b a s e ( y e s ) ?
return database
(database/idgcflskycrMCSA00057147 MODS11-0390)
reidentify identify
reidentify
42 4 r e f e r e n c e = R e f e r e n c e image
images = Images t o be r e i d e n t i f i e d answer = ” y e s ” F i t d i s p e r s i o n f u n c t i o n
i n t e r a c t i v e l y ?
c r v a l = Approximate c o o r d i n a t e ( a t r e f e r e n c e p i x e l )
c d e l t = Approximate d i s p e r s i o n ( i n t e r a c t i v e = ”no ” ) I n t e r a c t i v e f i t t i n g ?
( s e c t i o n = ” middle l i n e ” ) S e c t i o n t o apply t o two d i m e n s i o n a l images
( newaps = y e s ) R e i d e n t i f y a p e r t u r e s i n images not i n r e f e r e n c e ?
( o v e r r i d e = y e s ) O v e r r i d e p r e v i o u s s o l u t i o n s ? ( r e f i t = y e s ) R e f i t c o o r d i n a t e f u n c t i o n ? ( t r a c e = y e s ) Trace r e f e r e n c e image ?
( s t e p = ” 2 0 ” ) Step i n l i n e s / columns / bands f o r t r a c i n g an image
( nsum = ” 2 0 ” ) Number o f l i n e s / columns / bands t o sum ( s h i f t = ” 0 . ” ) S h i f t t o add t o r e f e r e n c e f e a t u r e s
(INDEF t o s e a r c h ) ( s e a r c h = 0 . 0 ) S e a r c h r a d i u s
( n l o s t = 1 0 0 ) Maximum number o f f e a t u r e s which may be l o s t
( c r a d i u s = 5 . 0 ) C e n t e r i n g r a d i u s
( t h r e s h o l d = 0 . 0 ) F e a t u r e t h r e s h o l d f o r c e n t e r i n g ( a d d f e a t u r e s = no ) Add f e a t u r e s from a l i n e l i s t ?
( c o o r d l i s t = ” mdpdb$ohlist / l i s t N S H K 500 ” ) User c o o r d i n a t e l i s t ( match = −3.0) C o o r d i n a t e l i s t matching l i m i t ( m a x f e a t u r e s = 5 0 ) Maximum number o f f e a t u r e s f o r
a u t o m a t i c i d e n t i f i c a t i o n ( minsep = 2 . 0 ) Minimum p i x e l s e p a r a t i o n ( d a t a b a s e = ” d a t a b a s e ” ) Database
( l o g f i l e s = ” l o g f i l e ” ) L i s t o f l o g f i l e s
( p l o t f i l e = ” ” ) P l o t f i l e f o r r e s i d u a l s ( v e r b o s e = y e s ) Verbose output ?
( g r a p h i c s = ” s t d g r a p h ” ) G r aphics output d e v i c e ( c u r s o r = ” ” ) G r aphics c u r s o r i n p u t
4.4. 43
( a i d p a r s = ” ” ) Automatic i d e n t i f i c a t i o n a l g o r i t h m p a r a m e t e r s
−−> r e i d e n t i f y gcflskycrMCSA00057147 MODS11−0390. f i t s gcflskycrMCSA00057147 MODS11−0390. f i t s
20
6 X,Y
2
X Y
fitcoords fitcoords
images = Images whose c o o r d i n a t e s a r e t o be f i t ( f i t n a m e = ” ” ) Name f o r c o o r d i n a t e f i t i n t h e
d a t a b a s e
( i n t e r a c t i v e = y e s ) F i t c o o r d i n a t e s i n t e r a c t i v e l y ? ( combine = no ) Combine i n p u t c o o r d i n a t e s f o r a
s i n g l e f i t ? ( d a t a b a s e = ” d a t a b a s e ” ) Database
( d e l e t i o n s = ” d e l e t i o n s . db ” ) D e l e t i o n l i s t f i l e ( not used i f n u l l ) ( f u n c t i o n = ” chebyshev ” ) Type o f f i t t i n g f u n c t i o n
( x o r d e r = 4 ) X o r d e r o f f i t t i n g f u n c t i o n ( y o r d e r = 3 ) Y o r d e r o f f i t t i n g f u n c t i o n ( l o g f i l e s = ”STDOUT, l o g f i l e ” ) Log f i l e s
( p l o t f i l e = ” p l o t f i l e ” ) P l o t l o g f i l e
( g r a p h i c s = ” s t d g r a p h ” ) G r aphics output d e v i c e ( c u r s o r = ” ” ) G r aphics c u r s o r i n p u t function xorder identify, reidentify
−−> f i t c o o r d s gcflskycrMCSA00057147 MODS11−0390
identify
fit-coords 4.3 fitcoords X Y
6interactive=yes
44 4
4.3: fitcoords
x[key],y[key] X Y key
r f
:order n n-1
:function
q fitcoords
?
[key]
x X
y Y
z s
r (s-z)
“xxyy” “r” X Y
identify reidentify
X X x Y (“xxyr”)
xorder RMS
2 2
X
y Y “xyyr” Y
q fitcoords
1 (MCSA00057147)
MOIRCS 3
7
1
fitcoords transform transform
i n p u t = Input images
output = Output images
f i t n a m e s = Names o f c o o r d i n a t e f i t s i n t h e d a t a b a s e
( minput = ” ” ) Input masks
7MCSMDP
4.4. 45
4.4: 2 :
( moutput = ” ” ) Output masks
( d a t a b a s e = ” d a t a b a s e ” ) I d e n t i f y d a t a b a s e ( i n t e r p t y p e = ” l i n e a r ” ) I n t e r p o l a t i o n type
( x1 = INDEF) Output s t a r t i n g x c o o r d i n a t e ( x2 = INDEF) Output e n d i n g x c o o r d i n a t e ( dx = INDEF) Output X p i x e l i n t e r v a l ( nx = INDEF) Number o f output x p i x e l s ( x l o g = no ) L o g a r i t h m i c x c o o r d i n a t e ?
( y1 = INDEF) Output s t a r t i n g y c o o r d i n a t e ( y2 = INDEF) Output e n d i n g y c o o r d i n a t e ( dy = INDEF) Output Y p i x e l i n t e r v a l ( ny = INDEF) Number o f output y p i x e l s ( y l o g = no ) L o g a r i t h m i c y c o o r d i n a t e ? ( f l u x = y e s ) Conserve f l u x p e r p i x e l ? ( blank = INDEF) Value f o r out o f r a n g e p i x e l s ( l o g f i l e s = ”STDOUT, l o g f i l e ” ) L i s t o f l o g f i l e s
transform “tr”
−−> ! s e d ’ s /\( .∗ \) / t r\1 / ’ gcMODS11−0390. l s t > trMODS11−0390. l s t
−−> t r a n s f o r m @gcMODS11−0390. l s t @trMODS11−0390. l s t gcflskycrMCSA00057147 MODS11−0390
transform ds9 4.4
Y
transform FITS
ds9 information panel
−−> mdpdisplay gcflabcrMCSA00057147 MODS11−0390. f i t s frame=1
−−> mdpdisplay trgcflabcrMCSA00057147 MODS11−0390. f i t s frame=2
46 4
4.4.2 残りスカイ引き
A-B A-B
transform Y
Y
background background
“bg”
i n p u t = Input images t o be background
s u b t r a c t e d
output = Output background s u b t r a c t e d images
( a x i s = 2 ) Axis a l o n g which background i s f i t and s u b t r a c t e d
( i n t e r a c t i v e = y e s ) S e t f i t t i n g p a r a m e t e r s i n t e r a c t i v e l y ? ( sample = ”∗” ) Sample o f p o i n t s t o u s e i n f i t
( n a v e r a g e = 1 ) Number o f p o i n t s i n sample a v e r a g i n g ( f u n c t i o n = ” chebyshev ” ) F i t t i n g f u n c t i o n
( o r d e r = 3 ) Order o f f i t t i n g f u n c t i o n ( l o w r e j e c t = 3 . 0 ) Low r e j e c t i o n i n sigma o f f i t ( h i g h r e j e c t = 3 . 0 ) High r e j e c t i o n i n sigma o f f i t
( n i t e r a t e = 3 ) Number o f r e j e c t i o n i t e r a t i o n s ( grow = 0 . 0 ) R e j e c t i o n growing r a d i u s
( g r a p h i c s = ” s t d g r a p h ” ) G r aphics output d e v i c e ( c u r s o r = ” ” ) G r aphics c u r s o r i n p u t
−−> ! s e d ’ s /\( .∗ \) / bg\1 / ’ trMODS11−0390. l s t > bgMODS11−0390. l s t
−−> background @trMODS11−0390. l s t @bgMODS11−0390. l s t Fit column
X ds9
1
background identify
2 (order=3)
“s”
4.5. 47 1
“q” Fit column column
ds9
4.5 足し合わせ
4 A-B
A B
A B 8
B A
−−> h s e l e c t @bgMODS11−0390. l s t ” $I , K DITWID” y e s bgtrgcflabcrMCSA00057147 MODS11−0390. f i t s 3 . 0 0 0 bgtrgcflabcrMCSA00057151 MODS11−0390. f i t s 3 . 0 0 0 bgtrgcflabcrMCSA00057159 MODS11−0390. f i t s 3 . 0 0 0 bgtrgcflabcrMCSA00057163 MODS11−0390. f i t s 3 . 0 0 0
3 MOIRCS 0.117
arcsec/pixel 3.0/0.117∼26pixel imshift
−−> ! s e d ’ s /\( .∗ \) / sh\1 / ’ bgMODS11−0390. l s t > shMODS11−0390. l s t
−−> i m s h i f t @bgMODS11−0390. l s t @shMODS11−0390. l s t 0 26 imarith
−−> ! s e d ’ s /\( .∗ \) / ng\1 / ’ shMODS11−0390. l s t > ngMODS11−0390. l s t
−−> i m a r i t h @shMODS11−0390. l s t ∗ −1 @ngMODS11−0390. l s t imcombine
median imcombine
i n p u t = ”” L i s t o f images t o combine
output = ”” L i s t o f output images
( h e a d e r s = ” ” ) L i s t o f h e a d e r f i l e s ( o p t i o n a l ) ( bpmasks = ” ” ) L i s t o f bad p i x e l masks ( o p t i o n a l ) ( r e j m a s k s = ” ” ) L i s t o f r e j e c t i o n masks ( o p t i o n a l ) ( n r e j m a s k s = ” ” ) L i s t o f number r e j e c t e d masks
48 4 ( o p t i o n a l )
( expmasks = ” ” ) L i s t o f e x p o s u r e masks ( o p t i o n a l ) ( s i gm as = ” ” ) L i s t o f sigma images ( o p t i o n a l )
( imcmb = ” $ I ” ) Keyword f o r IMCMB keywords ( l o g f i l e = ”STDOUT” ) Log f i l e
( combine = ” median ” ) Type o f combine o p e r a t i o n ( r e j e c t = ” s i g c l i p ” ) Type o f r e j e c t i o n
( p r o j e c t = no ) P r o j e c t h i g h e s t d i m e n s i on o f i n p u t images ?
( o u t t y p e = ” r e a l ” ) Output image p i x e l d a t a t y p e ( o u t l i m i t s = ” ” ) Output l i m i t s ( x1 x2 y1 y2 . . . )
( o f f s e t s = ” none ” ) Input image o f f s e t s ( masktype = ” none ” ) Mask type
( maskvalue = ” 0 ” ) Mask v a l u e
( blank = 0 . 0 ) Value i f t h e r e a r e no p i x e l s ( s c a l e = ” e x p o s u r e ” ) Image s c a l i n g
( z e r o = ” none ” ) Image z e r o p o i n t o f f s e t ( w e i g h t = ” e x p o s u r e ” ) Image w e i g h t s
( s t a t s e c = ” ” ) Image s e c t i o n f o r computing s t a t i s t i c s ( expname = ”EXPTIME” ) Image h e a d e r e x p o s u r e time keyword ( l t h r e s h o l d = INDEF) Lower t h r e s h o l d
( h t h r e s h o l d = INDEF) Upper t h r e s h o l d
( nlow = 1 ) minmax : Number o f low p i x e l s t o r e j e c t ( nhigh = 1 ) minmax : Number o f h i g h p i x e l s t o r e j e c t ( nkeep = 1 ) Minimum t o keep ( pos ) o r maximum t o
r e j e c t ( neg )
( m c l i p = y e s ) Use median i n sigma c l i p p i n g a l g o r i t h m s ? ( l s i g m a = 3 . 0 ) Lower sigma c l i p p i n g f a c t o r
( hsigma = 3 . 0 ) Upper sigma c l i p p i n g f a c t o r
( r d n o i s e = ” 0 . ” ) c c d c l i p : CCD r e a d o u t n o i s e ( e l e c t r o n s ) ( g a i n = ” 1 . ” ) c c d c l i p : CCD g a i n ( e l e c t r o n s /DN)
( s n o i s e = ” 0 . ” ) c c d c l i p : S e n s i t i v i t y n o i s e ( f r a c t i o n ) ( s i g s c a l e = 0 . 1 ) T o l e r a n c e f o r sigma c l i p p i n g s c a l i n g
c o r r e c t i o n s
( p c l i p = −0.5) p c l i p : P e r c e n t i l e c l i p p i n g parameter ( grow = 0 . 0 ) Radius ( p i x e l s ) f o r n e i g h b o r r e j e c t i o n ( mode = ” a l ” )
4.6. 49
−−> imcombine @bgMODS11−0390. l s t , @ngMODS11−0390. l s t HK500 MODS11−0390
4.6 フラックス較正
1 1 ADU
(Fλ[erg sec−1cm−2˚A−1])
Nobs(λ) Fλ,int R(λ)
Nobs(λ) =R(λ)×Fλ,int (4.1)
A Castelli (2004) [1]
8
2
2 HK500 M53735.fits
2 apall apall
2 1 2
‘yes’
−−> a p a l l HK500 M53735 . f i t s
R e c e n t e r a p e r t u r e s f o r HK500 M53735? ( ’ yes ’ ) : R e s i z e a p e r t u r e s f o r HK500 M53735? ( ’ yes ’ ) : Edit a p e r t u r e s f o r HK500 M53735? ( ’ yes ’ ) :
8
50 4
4.5: apall ( :
: :
: )
4.5
“b”
4.5 “s”
“q” “q”
Trace a p e r t u r e s f o r HK500 M53735?
F i t t r a c e d p o s i t i o n s f o r HK500 M53735 i n t e r a c t i v e l y ? F i t c u r v e t o a p e r t u r e 1 o f HK500 M53735 i n t e r a c t i v e l y
4.5
“s” “:order”
Write a p e r t u r e s f o r HK500 M53735 t o d a t a b a s e ? E x t r a c t a p e r t u r e s p e c t r a f o r HK500 M53735?
Review e x t r a c t e d s p e c t r a from HK500 M53735?
4.5 HK500 M53735.ms.fits
N(λ)
4.6. 51
R(λ) MCSMDP rescurve
s t d d a t a = ”” one−D s t a n d a r d s t a r data
stdmag = Vega magnitude o f t h e s t a n d a r d a t m a g f i l t e r output = ”” output r e s p o n s e c u r v e
( m a g f i l t e r = ” m d p d b $ f i l t e r /2 mass K . f i t s ” ) t r a n s m i s s i o n c u r v e o f f i l t e r f o r stdmag ( moddata = ” mdpdb$standard /A0V g50 . s p c ” ) s t a n d a r d s t a r model spectrum
( vegamodel = ” m d p d b $ f i l t e r / a l p h a l y r s t i s 0 0 3 . f i t s ” ) Vega model spectrum
(A0V g50.spc) K Vega
Vega Vega
MDPDB
HIPPARCOS A0 M53735 K 2MASS
8.856
−−> r e s c u r v e HK500 M53735 . ms . f i t s 8 . 8 5 6 resc CDFN2 HK500 . f i t s
resc CDFN2 HK500.fits R(λ) 9
2 R(λ) 2
MCSMDP mdpfcalib
1 (=
−−> m d p f c a l i b HK500 MODS11−0390. f i t s resc CDFN2 HK500 . f i t s HK500 MODS11−0390 f l . f i t s
HK500 MODS11-0390fl.fits [erg sec−1cm−2˚A−1]
9 FITS ds9
53
第 5 章 データ解析実習
54 5
5.1 赤方偏移を求める
z∼2
(Yoshikawa et al. 2010 [5]) 2
2 ds9 A-B
2 Y
−−> mdpdisplay HK500 MODS11−0390 f l . f i t s frame=1
splot 2 1
splot
−−> s p l o t HK500 MODS11−0390 f l . f i t s
splot 5.1 1
“:nsum”
ds9 2
5 7
(“w”) “a”
“k”
“k”
Gaussian FWHM
HK500 1.3µm 2.5µm z∼2
6000˚A
5.2 AGN 1 2
1
5.2. 55 5.1: splot
a
“a”
w
k Gaussian
“k”
s
q fitcoords
?
:nsum n n
5.2:
[˚A]
[SII] 6716, 6731 1:1
Hα 6563
[NII] 6548, 6583 Hα 1:2
[OIII] 4959, 5007 1:3
Hβ 4861
5.2 赤方偏移とフラックスから光度を求める
splot Hα
(F) (L)
L=F×4πd2L(z) (5.1)
dL(z)
MCSMDP
2 1
−−> cosmology 2 . 0
Luminosity D i s t a n c e (m) : 4 . 8 0 0 e+26
Angular Diameter D i s t a n c e (m) : 5 . 3 3 3 e+25 Look−Back Time ( yr ) : 1 . 0 2 4 e+10
56 5 参考・輝線の光度について
Hα
アパーチャ
1
銀河系での星間減光
NED Extinction Calculator2
E(B−V) = 0.011
目標の銀河での星間減光
2http://nedwww.ipac.caltech.edu/forms/calculator.html
57
参考文献
[1] Castelli, F., Kurucz, R. L. 2004, arXiv:astro-ph/0405087
[2] McLean, 1997, Electronic Imaging in Astronomy; Detectors and Instrumentation (John Wiley & Sons Ltd., 1997)
[3] Rousselot et al. 2000, A&A, 354, 1134 [4] Suzuki, R., et al. 2008, PASJ, 60, 1347 [5] Yoshikawa, T., et al. 2010, ApJ, 718, 112
59
付 録 A MCSMDP による自動リダク
ション
60 A MCSMDP
A.1 全体の流れと MCSMDP の特徴
MCSMDP
3.1 mdpproc,
mdpcombine
MCSMDP S/N
宇宙線・バッドピクセル処理 A-B
スペクトルの切り出し MDP
波長較正:各フレームでのOH夜光同定 1 OH
-波長較正:自動OH夜光同定 1 OH
スカイフレームの作成
multi-extention FITS (MEF) 2 extention
エラーの推定 Poisson
3 extention
SPECFITによる輝線フィット STSDAS specfit [NII]-Hα
HK500 Gaussian
Poisson MOIRCS
Poisson
σi =
√g×xi
g ×1
r = 1 r
rxi
g, (A.1)
xi g e−ADU−1 r response curve
ADU erg−1sec cm2˚A 1
σobj=
s2Σobjσ2i texp
, (A.2)
texp
√2 A-B
A.2. 61
A.2 実際の処理
A.2.1 フラットフレーム
mkdflat infile outfile
bpm (MDPDB )
offfile
normilize imstat
−−> m k d f l a t @ f l a t 1 . l s t HK500 CDFN2 Domeflat1 . f i t s mdpdb$bpm/ nlbpm1 FF64r . f i t s
A.2.2 テンプレートファイル作り
mktemplate
pyraf
mktemplate MCSMDP
mktemplate.conf
mktemplate
−−> copy MCSMDP$doc/ mktemplate . c o n f .
−−> ! v i mktemplate . c o n f
−−> mktemplate mktemplate . c o n f > l o g . py
mktemplate.conf bash
“=” “$”
## p a r a m e t e r s
c o o r d l i s t =”mdpdb\$ o h l i s t / l i s t N S H K 500 ” m d p f i l e=”CDFN MASK02 . mdp”
p r e f i x =”HK500”
# c a l i b r a t i o n data f o r c h i p 1
bpm1=”mdpdb\$bpm/ nlbpm1 FF64r . f i t s ”
62 A MCSMDP f l a t 1 =”HK500 CDFN2 Domeflat1 . f i t s ”
gdb1=”mdpdb\$geomap/ m c s d i s t c r r 1 f e b 0 7 n e w . dbs ” gmp1=”m c s d i s t c r r 1 f e b 0 7 n e w . gmp”
# c a l i b r a t i o n data f o r c h i p 2
bpm2=”mdpdb\$bpm/ nlbpm2 FF64r . f i t s ” f l a t 2 =”HK500 CDFN2 Domeflat2 . f i t s ”
gdb2=”mdpdb\$geomap/ m c s d i s t c r r 2 f e b 0 7 n e w . dbs ” gmp2=” m c s d i s t c r r 2 f e b 0 7 n e w . gmp”
coordlist mdpfile
prefix 2 prefix + ” ” + MDP
+ ”˙fits”
bpm[12 ] flat[12 ]
A.2.3 共通処理
mktemplate.sh
AB
A (B )
B(A)
h s e l e c t @ f i l e s . l i s t $ I ”@’DET−ID ’ = 1 & K DITCNT = 1” > f i l e s 1 a . l i s t h s e l e c t @ f i l e s . l i s t $ I ”@’DET−ID ’ = 1 & K DITCNT = 2” > f i l e s 1 b . l i s t h s e l e c t @ f i l e s . l i s t $ I ”@’DET−ID ’ = 2 & K DITCNT = 1” > f i l e s 2 a . l i s t h s e l e c t @ f i l e s . l i s t $ I ”@’DET−ID ’ = 2 & K DITCNT = 2” > f i l e s 2 b . l i s t
mdpproc
u n l e a r n mdpproc
i r a f . mdpproc . bpmask=”mdpdb$bpm/ nlbpm1 FF64r . f i t s ”
i r a f . mdpproc . gdb=”mdpdb$geomap/ m c s d i s t c r r 1 f e b 0 7 n e w . dbs ” i r a f . mdpproc . gmp=”m c s d i s t c r r 1 f e b 0 7 n e w . gmp”
i r a f . mdpproc . f l a t =”HK500 CDFN2 Domeflat1 . f i t s ” mdpproc @obj1a . l s t @obj1b . l s t
Multi Extention FITS (MEF) 0 ab
1 ( )
A.2. 63
A.2.4 データの切り出し
mdptran MDP FITS
−−> mdptran CDFN MASK02 . mdp CDFN MASK02tr1 . mdp gcflabcrMCSA00057147 . f i t s
1 MDP 2 MDP 3
FITS
maskplot ds9 FITS
s e l e c t shape then r e t u r n o r ’ q ’ t o q u i t :
return
c l i c k t h e o b j e c t p o s i t i o n : ds9
HK500 16692˚A
3 4 5
‘q’+return MDP ds9
type q t o w r i t e mdp d a t a b a s e t o f i l e :
‘q’+return MDP q
MDP cutspec
mdpproc
−−> c u t s p e c @gcobj1a . l s t CDFN MASK02tr1 . mdp name=ALE03 205
A.2.5 個別処理
mdpcombine mdptran
ds9 4
64 A MCSMDP
−−> u n l e a r n mdpcombine
−−> i r a f . mdpcombine . c o o r d l i s t =”mdpdb$ohlist / l i s t N S H K 500 ”
−−> i r a f . mdpcombine . ymin=42
−−> i r a f . mdpcombine . ymax=2003
−−> mdpcombine @ g c f l a b 1 . l i s t CDFN MASK02tr1 . mdp HK500 MODS11−0390 MODS11−0390 id mode=manual
1 “id mode=manual” OH
mktemplate sed MDP
mdpcombine
! awk ’{p r i n t $8}’ CDFN MASK02tr1 . mdp | s e d ’ s /\( .∗ \) / mdpcombine\
\@ g c f l a b 1\. l i s t\ ’CDFN MASK02’ t r 1\. mdp\ ’ HK500 ’ \1 \1 / ’
mdpcombine mdpcombine
−−> i r a f . mdpcombine . id mode =’ auto ’
−−> i r a f . mdpcombine . refname =’MODS11−0390 ’
−−> i r a f . mdpcombine . b g i n t e r=no
−−> mdpcombine @ g c f l a b 1 . l i s t CDFN MASK02tr1 . mdp HK500 MODS11−0094 MODS11−0094
A.2.6 標準星解析、フラックス較正
MOS mdpproc
mdpcombine MDP ymin, ymax
−−> u n l e a r n mdpproc
−−> i r a f . mdpproc . bpmask=”mdpdb$bpm/ nlbpm2 FF64r . f i t s ”
−−> i r a f . mdpproc . gdb=”mdpdb$geomap/ m c s d i s t c r r 2 f e b 0 7 n e w . dbs ”
−−> i r a f . mdpproc . gmp=”m c s d i s t c r r 2 f e b 0 7 n e w . gmp”
−−> i r a f . mdpproc . f l a t =”HK500 CDFN2 Domeflat2 . f i t s ”
−−> mdpproc MCSA00057114 . f i t s MCSA00057116 . f i t s
−−> u n l e a r n mdpcombine
−−> i r a f . mdpcombine . c o o r d l i s t =”mdpdb$ohlist / l i s t N S H K 500 ”
−−> mdpcombine gcflabMCSA00057114 . f i t s INDEF HK500 M53735 M53735 id mode=manual c o o r d l i s t =”mdpdb$ohlist / l i s t N S H K 500 ”
c a l f r a m e=”INDEF” ymin=900 ymax=1030
A.2. 65
• idenfity
id section ”line #”
• background ”s”
apall 1 rescurve mdpfcalib 2
apall MEF
imcopy MEF mdpfcalib
−−> imcopy HK500 M53735 . f i t s [ 0 ] HK500 M53735 0 . f i t s
−−> a p a l l HK500 M53735 0 . f i t s
−−> r e s c u r v e HK500 M53735 0 . ms . f i t s 8 . 8 5 6 resc CDFN2 HK500 . f i t s
−−> m d p f c a l i b HK500 MODS11−0390. f i t s resc CDFN2 HK500 . f i t s HK500 MODS11−0390 f l . f i t s
A.2.7 スペクトルの 1 次元化と輝線フィット
matplotlib mcsmdp ipython
$ mcsmdp −−i p y t h o n
mdp2dplot mdpextract 2 mdp2dplot 2
(center) (width) nsum
mdpextract SPECFIT
fittype “Halpha”
([NII]-Hα) “OIII” ([OIII]-Hβ) “emission” SPECFIT
SPECFIT
mdpextract mdp2dplot
In [ 1 ] : mcsmdp
In [ 2 ] : u n l e a r n mdp2dplot
In [ 3 ] : i r a f . mdp2dplot . x l a b e l=r ’ $\lambda {\rm obs} \l e f t [\AA\r i g h t ] $ ’ In [ 4 ] : i r a f . mdp2dplot . y l a b e l=
r ’ $F \lambda\, [\t i m e s 1 0 ˆ{ −18}\, e r g\,cmˆ{−2}s ˆ{−1}\AAˆ{ −1}] $ ’ In [ 5 ] : i r a f . mdp2dplot . b o t t o m p l o t=F a l s e
66 A MCSMDP In [ 6 ] : i r a f . mdp2dplot . y s c a l e =1e18
In [ 7 ] : i r a f . mdp2dplot . l i n e l i s t =” m d p d b $ l i n e l i s t / n e b u l a e t e x . g a i a ” In [ 8 ] : i r a f . mdp2dplot . l i n e w i d t h =1
In [ 9 ] : i r a f . mdp2dplot . l i n e a l p h a =1 In [ 1 0 ] : i r a f . mdp2dplot . l s t y l e =’k−−’ In [ 1 1 ] : i r a f . mdp2dplot . ymin=0
In [ 1 2 ] : mdp2dplot HK500 MODS11−0390 f l 59 21990 width =800 nsum=5 In [ 1 3 ] : mdpextract HK500 MODS11−0390 f l 59 21990
HK500 MODS11−0390 f l H a l p h a width =800 nsum=5 s u b f r a c =0.5 sub fwhm=42 In [ 1 4 ] : mdp2dplot HK500 MODS11−0390 f l H a l p h a 59 21990 p l o t f i t+
ymin=0 width =800
In [ 1 5 ] : pylab . s a v e f i g ( ’ HK500 MODS11−0390 f l H a l p h a . png ’ ) A.1