12/15/2005 [[[ 5 ]]]
Slit 5 (ch1)
SCUBA 850.53
SXDS zJ
imarith ../data/MCSA00014965.fits - ../data/MCSA00014967.fits 14965-14967.fits
imarith ../data/MCSA00014966.fits - ../data/MCSA00014968.fits 14966-14968.fits
imarith ../data/MCSA00014576.fits - ../data/MCSA00014578.fits 14576-14578.fits
-----------------------------------
Slit 6 zJ totak= 6600s = 110min
-----------------------------------
600sの引き算画像で写っているので、それに合わせてなるべくぎりぎりのサイズで切り出す。
day1:
imcopy ../data/MCSA00014575.fits[*,705:810] 5o_zJ_14575.fits zJ500 300s
imcopy ../data/MCSA00014577.fits[*,705:810] 5o_zJ_14577.fits
imcopy ../data/MCSA00014579.fits[*,705:810] 5o_zJ_14579.fits 600s
imcopy ../data/MCSA00014581.fits[*,705:810] 5o_zJ_14581.fits
imcopy ../data/MCSA00014583.fits[*,705:810] 5o_zJ_14583.fits
imcopy ../data/MCSA00014585.fits[*,705:810] 5o_zJ_14585.fits
imcopy ../data/MCSA00014587.fits[*,705:810] 5o_zJ_14587.fits 存在しない
imcopy ../data/MCSA00014589.fits[*,705:810] 5o_zJ_14589.fits
imcopy ../data/MCSA00014591.fits[*,705:810] 5o_zJ_14591.fits
day2:
imcopy ../data/MCSA00014965.fits[*,705:810] 5o_zJ_14965.fits zJ500 600sec
imcopy ../data/MCSA00014967.fits[*,705:810] 5o_zJ_14967.fits
imcopy ../data/MCSA00014969.fits[*,705:810] 5o_zJ_14969.fits
imcopy ../data/MCSA00014971.fits[*,705:810] 5o_zJ_14971.fits
---------
差し引き
---------
day1:
imarith 5o_zJ_14575.fits - 5o_zJ_14577.fits 5o_zJ_75-77.fits
imarith 5o_zJ_14579.fits - 5o_zJ_14581.fits 5o_zJ_79-81.fits
imarith 5o_zJ_14583.fits - 5o_zJ_14585.fits 5o_zJ_83-85.fits
imarith 5o_zJ_14589.fits - 5o_zJ_14591.fits 5o_zJ_89-91.fits
day2:
imarith 5o_zJ_14965.fits - 5o_zJ_14967.fits 5o_zJ_65-67.fits
imarith 5o_zJ_14969.fits - 5o_zJ_14971.fits 5o_zJ_69-71.fits
---------
波長較正
---------
for day1:
identify 5o_zJ_14579.fits section="middle line" coordli="OHlines_zJ500.dat" nsum=5 function="cheby"
17 order=4 rms=0.6189
reidentify 5o_zJ_14579.fits 5o_zJ_14579.fits section="middle line" step=10 nsum=5 overrid+ interac+ nlost=3
スペクトルきれい。4次にするとreidも比較的きれい
fitcoord 5o_zJ_14579 function="cheby" xorder=1 yorder=1 interac+
rms=2.116 （yorder=4でよかったが、あえて2でとめた）
X order = 4
Y order = 2
(1, 1) = 19673.76 (2048, 1) = 8437.482
(1, 106) = 19673.99 (2048, 106) = 8369.956
transform 5o_zJ_75-77.fits 5o_zJ_75-77_tr.fits 5o_zJ_14579 interpt="linear"
transform 5o_zJ_79-81.fits 5o_zJ_79-81_tr.fits 5o_zJ_14579 interpt="linear"
transform 5o_zJ_83-85.fits 5o_zJ_83-85_tr.fits 5o_zJ_14579 interpt="linear"
transform 5o_zJ_89-91.fits 5o_zJ_89-91_tr.fits 5o_zJ_14579 interpt="linear"
for day2:
identify 5o_zJ_14965.fits section="middle line" coordli="OHlines_zJ500.dat" nsum=5 function="cheby"
15 order=4 0.5084
短いほうはどうせ波長見ないのでフィットしない。もったいないので絵はカットしない。
reidentify 5o_zJ_14965.fits 5o_zJ_14965.fits section="middle line" step=10 nsum=5 overrid+ interac+ nlost=3
fitcoord 5o_zJ_14965 function="cheby" xorder=1 yorder=1 interac+
2.004
X order = 1
Y order = 1
(1, 1) = 19692.4 (2048, 1) = 8118.201
(1, 106) = 19688.91 (2048, 106) = 8114.715
transform 5o_zJ_65-67.fits 5o_zJ_65-67_tr.fits 5o_zJ_14965 interpt="linear"
transform 5o_zJ_69-71.fits 5o_zJ_69-71_tr.fits 5o_zJ_14965 interpt="linear"
------
FLAT
------
DISPAXISが入ってなかった。objのときは入ってたのに。ヘッダはダークになっている。
画像確認。day1/day2で同じ領域を使える。
for day1:
imcopy ../data/MCSA00014823.fits[*,705:810] 5f_zJ_14823_1.fits
imcopy ../data/MCSA00014825.fits[*,705:810] 5f_zJ_14825_1.fits
imcopy ../data/MCSA00014827.fits[*,705:810] 5f_zJ_14827_1.fits
imcombine 5f_zJ_14823_1.fits,5f_zJ_14825_1.fits,5f_zJ_14827_1.fits 5f_zj_1_p25.fits combine="median" reject="none"
hedit 5f_zj_1_p25.fits DISPAXIS 1 add+ ver-
transform 5f_zj_1_p25.fits 5f_zj_1_p25_tr.fits 5o_zJ_14579 interpt="linear"
blkavg 5f_zj_1_p25_tr.fits 5f_zj_1_p25_trblk.fits 1 105 option="average"
imcopy 5f_zj_1_p25_trblk.fits[*,1] 5f_zj_1_p25_trblk1.fits
imarith 5f_zj_1_p25_tr.fits / 5f_zj_1_p25_trblk1.fits 5f_zj_1_p25_trn.fits
ccdmask 5f_zj_1_p25_trn.fits 5f_zj_1_p25_trn
imcopy 5f_zj_1_p25_trn.fits 5f_zj_1_p25_trnm.fits
fixpix 5f_zj_1_p25_trnm.fits 5f_zj_1_p25_trn.pl << flat for day1&2
--------------
MASK for OBJ
--------------
day1:
ccdmask 5o_zJ_75-77_tr.fits 5o_zJ_75-77_tr
ccdmask 5o_zJ_79-81_tr.fits 5o_zJ_79-81_tr
ccdmask 5o_zJ_83-85_tr.fits 5o_zJ_83-85_tr
ccdmask 5o_zJ_89-91_tr.fits 5o_zJ_89-91_tr
imcopy 5o_zJ_75-77_tr.fits 5o_zJ_75-77_trm.fits
imcopy 5o_zJ_79-81_tr.fits 5o_zJ_79-81_trm.fits
imcopy 5o_zJ_83-85_tr.fits 5o_zJ_83-85_trm.fits
imcopy 5o_zJ_89-91_tr.fits 5o_zJ_89-91_trm.fits
fixpix 5o_zJ_75-77_trm.fits 5o_zJ_75-77_tr.pl
fixpix 5o_zJ_79-81_trm.fits 5o_zJ_79-81_tr.pl
fixpix 5o_zJ_83-85_trm.fits 5o_zJ_83-85_tr.pl
fixpix 5o_zJ_89-91_trm.fits 5o_zJ_89-91_tr.pl
day2:
ccdmask 5o_zJ_65-67_tr.fits 5o_zJ_65-67_tr
ccdmask 5o_zJ_69-71_tr.fits 5o_zJ_69-71_tr
imcopy 5o_zJ_65-67_tr.fits 5o_zJ_65-67_trm.fits
imcopy 5o_zJ_69-71_tr.fits 5o_zJ_69-71_trm.fits
fixpix 5o_zJ_65-67_trm.fits 5o_zJ_65-67_tr.pl
fixpix 5o_zJ_69-71_trm.fits 5o_zJ_69-71_tr.pl
------------
OBJ / FLAT
------------
同じフラット
for day1:
imarith 5o_zJ_75-77_trm.fits / 5f_zj_1_p25_trnm.fits 5o_zJ_75-77_trmf.fits
imarith 5o_zJ_79-81_trm.fits / 5f_zj_1_p25_trnm.fits 5o_zJ_79-81_trmf.fits
imarith 5o_zJ_83-85_trm.fits / 5f_zj_1_p25_trnm.fits 5o_zJ_83-85_trmf.fits
imarith 5o_zJ_89-91_trm.fits / 5f_zj_1_p25_trnm.fits 5o_zJ_89-91_trmf.fits
for day2:
imarith 5o_zJ_65-67_trm.fits / 5f_zj_1_p25_trnm.fits 5o_zJ_65-67_trmf.fits
imarith 5o_zJ_69-71_trm.fits / 5f_zj_1_p25_trnm.fits 5o_zJ_69-71_trmf.fits
-----------
重ね合わせ
-----------
day1:
imarith 5o_zJ_75-77_trmf.fits * -1.0 5o_zJ_75-77_trmf-1.fits
imshift 5o_zJ_75-77_trmf-1.fits 5o_zJ_75-77_trmf-1sh.fits 0 44
imarith 5o_zJ_75-77_trmf.fits + 5o_zJ_75-77_trmf-1sh.fits 5o_zJ_75-77_trmf_comb.fits
imarith 5o_zJ_79-81_trmf.fits * -1.0 5o_zJ_79-81_trmf-1.fits
imshift 5o_zJ_79-81_trmf-1.fits 5o_zJ_79-81_trmf-1sh.fits 0 44
imarith 5o_zJ_79-81_trmf.fits + 5o_zJ_79-81_trmf-1sh.fits 5o_zJ_79-81_trmf_comb.fits
imarith 5o_zJ_83-85_trmf.fits * -1.0 5o_zJ_83-85_trmf-1.fits
imshift 5o_zJ_83-85_trmf-1.fits 5o_zJ_83-85_trmf-1sh.fits 0 44
imarith 5o_zJ_83-85_trmf.fits + 5o_zJ_83-85_trmf-1sh.fits 5o_zJ_83-85_trmf_comb.fits
imarith 5o_zJ_89-91_trmf.fits * -1.0 5o_zJ_89-91_trmf-1.fits
imshift 5o_zJ_89-91_trmf-1.fits 5o_zJ_89-91_trmf-1sh.fits 0 44
imarith 5o_zJ_89-91_trmf.fits + 5o_zJ_89-91_trmf-1sh.fits 5o_zJ_89-91_trmf_comb.fits
imcombine 5o_zJ_75-77_trmf_comb.fits,5o_zJ_79-81_trmf_comb.fits,5o_zJ_83-85_trmf_comb.fits,5o_zJ_89-91_trmf_comb.fits 5o_zJ_day1.fits combine="sum" reject="none"
imcopy 5o_zJ_day1.fits[*,50:105] 5o_zJ_day1s.fits
day2:
imarith 5o_zJ_65-67_trmf.fits * -1.0 5o_zJ_65-67_trmf-1.fits
imshift 5o_zJ_65-67_trmf-1.fits 5o_zJ_65-67_trmf-1sh.fits 0 44
imarith 5o_zJ_65-67_trmf.fits + 5o_zJ_65-67_trmf-1sh.fits 5o_zJ_65-67_trmf_comb.fits
imarith 5o_zJ_69-71_trmf.fits * -1.0 5o_zJ_69-71_trmf-1.fits
imshift 5o_zJ_69-71_trmf-1.fits 5o_zJ_69-71_trmf-1sh.fits 0 44
imarith 5o_zJ_69-71_trmf.fits + 5o_zJ_69-71_trmf-1sh.fits 5o_zJ_69-71_trmf_comb.fits
imcombine 5o_zJ_65-67_trmf_comb.fits,5o_zJ_69-71_trmf_comb.fits 5o_zJ_day2.fits combine="sum" reject="none"
imcopy 5o_zJ_day2.fits[*,50:105] 5o_zJ_day2s.fits
----------------------------------------------------------------------------------------
SXDS HK
imarith ../data/MCSA00014988.fits - ../data/MCSA00014990.fits 14988-14990.fits ch2
imarith ../data/MCSA00014987.fits - ../data/MCSA00014989.fits 14987-14989.fits ch1
-----------
Slit 6 HK total = 5400s = 100min
-----------
imcopy ../data/MCSA00014987.fits[800:2048,700:830] 5o_HK_14987.fits HK500 600sec
imcopy ../data/MCSA00014989.fits[800:2048,700:830] 5o_HK_14989.fits
imcopy ../data/MCSA00014995.fits[800:2048,700:830] 5o_HK_14995.fits
imcopy ../data/MCSA00014997.fits[800:2048,700:830] 5o_HK_14997.fits
imcopy ../data/MCSA00014999.fits[800:2048,700:830] 5o_HK_14999.fits
imcopy ../data/MCSA00015001.fits[800:2048,700:830] 5o_HK_15001.fits
imcopy ../data/MCSA00015003.fits[800:2048,700:830] 5o_HK_15003.fits
imcopy ../data/MCSA00015005.fits[800:2048,700:830] 5o_HK_15005.fits
imcopy ../data/MCSA00014983.fits[800:2048,700:830] 5o_HK_14983.fits HK500 300sec
imcopy ../data/MCSA00014985.fits[800:2048,700:830] 5o_HK_14985.fits
imarith 5o_HK_14987.fits - 5o_HK_14989.fits 5o_HK_87-89.fits
imarith 5o_HK_14995.fits - 5o_HK_14997.fits 5o_HK_95-97.fits
imarith 5o_HK_14999.fits - 5o_HK_15001.fits 5o_HK_99-01.fits
imarith 5o_HK_15003.fits - 5o_HK_15005.fits 5o_HK_03-05.fits
imarith 5o_HK_14983.fits - 5o_HK_14985.fits 5o_HK_83-85.fits 300s
identify 5o_HK_14987.fits section="middle line" coordli="OHlines_HK500.dat" nsum=5 function="cheby"
7+both order=4 0.4814
reidentify 5o_HK_14987.fits 5o_HK_14987.fits section="middle line" step=10 nsum=5 overrid+ interac+ nlost=3
fitcoord 5o_HK_14987 function="cheby" xorder=1 yorder=1 interac+
3.68
X order = 4
Y order = 2
(1, 1) = 24877.34 (1249, 1) = 15218.74
(1, 131) = 24866.2 (1249, 131) = 15200.44
transform 5o_HK_87-89.fits 5o_HK_87-89_tr.fits 5o_HK_14987 interpt="linear"
transform 5o_HK_95-97.fits 5o_HK_95-97_tr.fits 5o_HK_14987 interpt="linear"
transform 5o_HK_99-01.fits 5o_HK_99-01_tr.fits 5o_HK_14987 interpt="linear"
transform 5o_HK_03-05.fits 5o_HK_03-05_tr.fits 5o_HK_14987 interpt="linear"
transform 5o_HK_83-85.fits 5o_HK_83-85_tr.fits 5o_HK_14987 interpt="linear"
------
FLAT
------
imcopy ../data/MCSA00014875.fits[800:2048,700:830] 5f_HK_14875.fits p0mm
imcopy ../data/MCSA00014877.fits[800:2048,700:830] 5f_HK_14877.fits
imcopy ../data/MCSA00014879.fits[800:2048,700:830] 5f_HK_14879.fits
imcombine 5f_HK_14875.fits,5f_HK_14877.fits,5f_HK_14879.fits 5f_HK_p25.fits combine="median" reject="none"
hedit 5f_HK_p25.fits DISPAXIS 1 add+ ver-
transform 5f_HK_p25.fits 5f_HK_p25_tr.fits 5o_HK_14987 interpt="linear"
blkavg 5f_HK_p25_tr.fits 5f_HK_p25_trblk.fits 1 131 option="average"
imcopy 5f_HK_p25_trblk.fits[*,1] 5f_HK_p25_trblk1.fits
imarith 5f_HK_p25_tr.fits / 5f_HK_p25_trblk1.fits 5f_HK_p25_trn.fits
ccdmask 5f_HK_p25_trn.fits 5f_HK_p25_trn
imcopy 5f_HK_p25_trn.fits 5f_HK_p25_trnm.fits
fixpix 5f_HK_p25_trnm.fits 5f_HK_p25_trn.pl << flat
-------------
MASK for OBJ
-------------
ccdmask 5o_HK_87-89_tr.fits 5o_HK_87-89_tr
ccdmask 5o_HK_95-97_tr.fits 5o_HK_95-97_tr
ccdmask 5o_HK_99-01_tr.fits 5o_HK_99-01_tr
ccdmask 5o_HK_03-05_tr.fits 5o_HK_03-05_tr
ccdmask 5o_HK_83-85_tr.fits 5o_HK_83-85_tr
imcopy 5o_HK_87-89_tr.fits 5o_HK_87-89_trm.fits
imcopy 5o_HK_95-97_tr.fits 5o_HK_95-97_trm.fits
imcopy 5o_HK_99-01_tr.fits 5o_HK_99-01_trm.fits
imcopy 5o_HK_03-05_tr.fits 5o_HK_03-05_trm.fits
imcopy 5o_HK_83-85_tr.fits 5o_HK_83-85_trm.fits
fixpix 5o_HK_87-89_trm.fits 5o_HK_87-89_tr.pl
fixpix 5o_HK_95-97_trm.fits 5o_HK_95-97_tr.pl
fixpix 5o_HK_99-01_trm.fits 5o_HK_99-01_tr.pl
fixpix 5o_HK_03-05_trm.fits 5o_HK_03-05_tr.pl
fixpix 5o_HK_83-85_trm.fits 5o_HK_83-85_tr.pl
------------
OBJ / FLAT
------------
imarith 5o_HK_87-89_trm.fits / 5f_HK_p25_trnm.fits 5o_HK_87-89_trmf.fits
imarith 5o_HK_95-97_trm.fits / 5f_HK_p25_trnm.fits 5o_HK_95-97_trmf.fits
imarith 5o_HK_99-01_trm.fits / 5f_HK_p25_trnm.fits 5o_HK_99-01_trmf.fits
imarith 5o_HK_03-05_trm.fits / 5f_HK_p25_trnm.fits 5o_HK_03-05_trmf.fits
imarith 5o_HK_83-85_trm.fits / 5f_HK_p25_trnm.fits 5o_HK_83-85_trmf.fits
imarith 5o_HK_87-89_trmf.fits * -1.0 5o_HK_87-89_trmf-1.fits
imshift 5o_HK_87-89_trmf-1.fits 5o_HK_87-89_trmf-1sh.fits 0 44
imarith 5o_HK_87-89_trmf.fits + 5o_HK_87-89_trmf-1sh.fits 5o_HK_87-89_trmf_comb.fits
imarith 5o_HK_95-97_trmf.fits * -1.0 5o_HK_95-97_trmf-1.fits
imshift 5o_HK_95-97_trmf-1.fits 5o_HK_95-97_trmf-1sh.fits 0 44
imarith 5o_HK_95-97_trmf.fits + 5o_HK_95-97_trmf-1sh.fits 5o_HK_95-97_trmf_comb.fits
imarith 5o_HK_99-01_trmf.fits * -1.0 5o_HK_99-01_trmf-1.fits
imshift 5o_HK_99-01_trmf-1.fits 5o_HK_99-01_trmf-1sh.fits 0 44
imarith 5o_HK_99-01_trmf.fits + 5o_HK_99-01_trmf-1sh.fits 5o_HK_99-01_trmf_comb.fits
imarith 5o_HK_03-05_trmf.fits * -1.0 5o_HK_03-05_trmf-1.fits
imshift 5o_HK_03-05_trmf-1.fits 5o_HK_03-05_trmf-1sh.fits 0 44
imarith 5o_HK_03-05_trmf.fits + 5o_HK_03-05_trmf-1sh.fits 5o_HK_03-05_trmf_comb.fits
imarith 5o_HK_83-85_trmf.fits * -1.0 5o_HK_83-85_trmf-1.fits
imshift 5o_HK_83-85_trmf-1.fits 5o_HK_83-85_trmf-1sh.fits 0 44
imarith 5o_HK_83-85_trmf.fits + 5o_HK_83-85_trmf-1sh.fits 5o_HK_83-85_trmf_comb.fits
imcombine 5o_HK_87-89_trmf_comb.fits,5o_HK_95-97_trmf_comb.fits,5o_HK_99-01_trmf_comb.fits,5o_HK_03-05_trmf_comb.fits,5o_HK_83-85_trmf_comb.fits 5o_HK_trmf_all.fits combine="sum" reject="none"
imcopy 5o_HK_trmf_all.fits[*,55:130] 5o_HK_trmf_alls.fits