Sky 引きに関係する Tips¶
HSC pipeline デフォルトの sky 引きを行いたくない場合¶
HSC pipeline では、sky 引きが singleFrameDriver.py と coaddDriver.py の 2 つの解析過程で行われています。 このうち singleFrameDriver.py で引かれた sky は coaddDriver.py で戻すことができます。
# singleFrameDriver.py で引かれた sky を戻して coaddDriver.py を実行
coaddDriver.py $home/HSC --calib=$home/HSC/CALIB --rerun test --id filter=HSC-I tract=0 --selectId visit=124776..124790:2 ccd=0..103 --config makeCoaddTempExp.bgSubtracted=False
detectCoaddSources.detection.reEstimateBackground=False
# オプション
# --config makeCoaddTempExp.bgSubtracted=False:デフォルトは True で、reduceFrames.py の sky 引きを残している
# detectCoaddSources.detection.reEstimateBackground=False : デフォルトは True。True のままだと background の見積もりがうまくいかない。
sky が戻された画像データは ~/HSC/rerun/[rerun]/deepCoadd/[filter]/[tract]/[patch] 以下に warp-[filter]-[tract]-[patch]-[visit].fits という名前で各 visit 毎の warp 画像として確認することができます。 calexp データについては、CCD 間の段差が目立ってしまう場合があります (例えば、各ショットの CCD の隙間がうっすら見える、等)。
Warning
例えば、一度 coaddDriver.py を実行しており、rerun 以下に deepCoadd ディレクトリができている場合、 二度目の coaddDriver.py では既存の deepCoadd 以下にある warp 画像を読み込むため新たに生成してくれません。 そこで、deepCoadd 下にある [filter] ディレクトリ名を変更するとこの問題は回避でき、 sky を戻した warp 画像を生成してくれます。
singleFrameDriver.py の sky estimate の領域のサイズを変更したい場合¶
singleFrameDriver.py の実行時に以下のオプションを追加することで変更できます。
# sky estimate の領域サイズを 256 × 256 ピクセルに変更した場合
singleFrameDriver.py $home/hsc --calib=$home/hsc/CALIB --rerun test --id filter=HSC-I visit=902798..902808:2 --config processCcd.isr.doFringe=False processCcd.detection.background.binSize=256 processCcd.detection.background.binSize=256
# オプション
# --config processCcd.detection.background.binSize=256 processCcd.calibrate.background.binSize=256
skycorrection.pyをしないでcoaddDriver.pyを実行する場合¶
coaddDriver.py の実行時に”–config makeCoaddTempExp.doApplySkyCorr=False”
https://community.lsst.org/t/sky-subtraction/2415
# skycorrection.pyをしないでcoaddDriver.pyを実行する場合
coaddDriver.py $home/hsc --calib=$home/hsc/CALIB --rerun test --id filter=HSC-I visit=902798..902808:2 --config makeCoaddTempExp.doApplySkyCorr=False
# オプション
# --config makeCoaddTempExp.doApplySkyCorr=False