NIC Web
staff/地球照偏光観測
をテンプレートにして作成
開始行:
[[staff]]
NICでの地球照偏光観測についてのページ(PI: 高橋)
* 観測手順 [#ic7284ad]
** 事前準備 [#k10dabec]
- なゆたファイル駆動用のトラックファイルを作り、UCCの
/export/home/observer/NHAO/FILEMOVE
に保存しておく
*** 詳細 [#q78d3d37]
-- 個人 MacBookPro (inyanga) でターミナルを開き
$ cd /Users/junon/Sync/Documents/20YYdo/obsprep/Earthsh...
-- lunar2radec.awk をエディタで開き、観測対象の月座標 (lm...
-- 2021年6月:ターゲット位置を変更した。Montanes-Rodrigıue...
--- 上弦期: Grimaldi side の patch #3: 月座標 (lat=-5.0...
--- 下弦期: Crisium side の patch #3: 月座標 (lat=0.00,...
-- トラックファイルを作り、%%gamera%%(<--ftpがない) OBS1...
$ ./prepLunarObs.py [start_date(yy-mm-dd) start_time(h...
(例)./prepLunarObs.py 2021-05-26 10:00:00
--- 引数で指定した開始時刻から5時間分(15秒間隔)のトラック...
-- OBS1b に ssh ログイン
-- ucc に ftp ログインして、上記のファイル置き場に移動、...
** 導入まで [#w53c41fe]
- ucc_run: 近赤外線観測モード、視野回転=(モード= AZ/EL同...
- calfocus でざっくり焦点合わせ
- shutterはopenにしておく
- HWP挿入
$ pol in
- ステラナビで月に近い、6-8等くらいの星に望遠鏡を向ける
- 以下のコマンドで、星を偏光視野に入れる。そのときのオフ...
$TL 3 2; finpol
$xyoffset2 (左側の偏光窓に星を移動)
$TL 3 2
- 偏光素子を挿入
$darkshutter.py -p
- focus値を +0.2 して、ざっくり焦点を合わせる
- autofocus で焦点合わせ
- ucc_run ファイル駆動モードで、月の地球照側に pointing
- さきほどのオフセット量を適用
- テスト撮影しながら、視野を調整。月の端が視野の真ん中に...
-- オフセット適用後は、望遠鏡は月の内側に向きすぎているこ...
--- 2021年5月以降、月座標でオフセットする
obs1b% lunarOffset.sh 月座標で北へのオフセット量 月座標...
オフセット量は秒角。月座標での東西と、赤道座標での東西(R...
- 地球照観測
$ PLo 3 (exp) Earthshine (repeat) ;fin
- 30分ごとを目安にスカイ観測
-- スカイ領域にpointing。上弦期は del_RA>0 (東に)、下弦期...
$ nayuta -do (del_RA) 0
-- スカイ観測
$ PLo 3 (exp) Earthshine_sky 1 ;fin
- 天候が安定していたら、nicpc:~/takahashi/obsScript/ESpol...
vi ESpol.sh
で変数を編集(exptime, dRAなど)。(地球照 x $series回 + ス...
* status [#o57a13ba]
- 第1期観測
| 日付 | 上弦/下弦側 (月齢)| 露出時間(s)xセット数 | 天...
| 2019.05.29 | 下 (24d) | 120s x 3 | 晴れ | bad | 再済 |...
| 2019.06.10 | 上 (6d) | 120s x 2 | 曇り | NG | 済み | ...
| 2019.06.25 | 下(22d) | (45s, 30s, 20s) x many | 薄曇...
| 2019.07.07 | 上 (4d) | 180s x 2 | 薄曇り | たぶんNG |...
| 2019.08.25 | 下 (24d) | 120s x 12 | 快晴 | good | 済...
| 2019.08.26 | 下 (25d) | 60s x 8 | 薄曇り | so so | 済...
| 2019.10.06 | 上 (7d) | 60s x 3 | 薄曇り->晴れ | bad |...
| 2019.10.21 | 下 (23d) | 60s x8, 120s x7 | 曇り->晴れ ...
| 2019.10.22 | 下 (24d) | 60s x many | 曇り->晴れ | bad...
| 2019.11.01 | 上 (4d) | 60s x 3 | 曇り->晴れ | so so |...
| 2019.11.05 | 上 (8d) | 60s x 13 | 晴れ | good | 再済 ...
| 2019.11.06 | 上 (9d) | 60s x 16 | 晴れ->少し雲? | goo...
| 2019.11.20 | 下 (24d) | 120s x 3 | 薄曇り | so so | ...
| 2019.11.21 | 下 (25d) | 180s x 10 | 晴れ | good | 再...
| 2019.12.04 | 上 (7d) | 60s x 7 | 薄曇り | so so | 済 ...
| 2019.12.06 | 上 (10d) | 5s x many | 曇り | NG | - | -...
| 2019.12.07 | 上 (11d) | 30s x 8 | 曇り->晴れ | so so...
| 2019.12.08 | 上 (12d) | 25s x 17 | 晴れ | so so | 済...
| 2019.12.09 | 上 (13d) | 25s x 25 | 晴れ | so so | 済...
| 2019.12.10 | 上 (14d) | 2s x | 薄曇り | NG | - | -...
| 2019.12.14 | 下 (18d) | 20s x | 曇り | NG | - | - ...
| 2019.12.15 | 下 (19d) | 2s x many | 薄曇り | NG | - ...
| 2019.12.18 | 下 (22d) | 120s, 90s x (3h分) | 晴れ->...
| 2019.12.19 | 下 (23d) | 90s x many (2.5h分) | 晴れ->...
| 2019.12.20 | 下 (24d) | 15s x many | 薄曇り | たぶ...
| 2019.12.22 | 下 (26d) | 100s x 2 | 晴れ->薄曇り | b...
| 2019.12.28 | 上(2d) | 20s x 7 | 薄曇り<->晴れ | so ...
| 2020.01.02 | 上(7d) | 90s x 13 | 薄曇り | bad | ...
| 2020.01.03 | 上(8d) | 60s x (4hくらい) | 晴れ | good ...
| 2020.01.04 | 上(9d) | 60s x 2 | 曇り | bad | 済 | 済...
| 2020.01.05 | 上(10d) | 40s x 11 | 曇り<->晴れ | so so ...
| 2020.01.15 | 下(20d) | 45s x 約5時間分 | 薄曇り | so s...
| 2020.01.17 | 下(23d) | 30s x | 薄曇り | bad - so so ...
| 2020.01.18 | 下(24d) | 90s x | 快晴 | good | 済 | ...
| 2020.01.30 | 上(5d) | 60s x 6 | 曇り | bad | 済 | 済 ...
| 2020.02.01 | 上(8d) | 60s x 7 | 晴れ/曇り | so so| 済...
| 2020.02.02 | 上(8d) | 60s x 4 | 曇り | bad| 済 | しな...
| 2020.02.03 | 上(9d) | 60s x 7 | 曇り | bad| 済 | 済 |...
| 2020.02.04 | 上(10d) | 30s x 8 | 晴れ | so so| 済 | ...
| 2020.02.05 | 上(11d) | 30s x 8 | 晴れ、時々雲 | so so|...
| 2020.02.06 | 上(12d) | 15s x 7, 25s x 3 | 薄曇り | ba...
| 2020.02.13 | 下(20d) | 30s x 10? | 晴れ | bad| 済 | ...
| 2020.02.26 | 上(3d) | 30s x | 晴れ | good| 済 | 済 ...
| 2020.03.02 | 上(8d) | 60s x (2.5h分) | 晴れ | good| ...
| 2020.03.14 | 下(20d) | 30s x 6 | 薄曇り | so so | 済 ...
| 2020.04.02 | 上(9d) | 30s x 3 | 薄曇り | so so | 済 |...
| 2020.04.03 | 上(10d) | 20s x (5.5h分) | 薄曇り | so s...
| 2020.04.29 | 上(6d) | 60s x (4h分) | 晴れ(ごく薄い雲...
- 第2期観測
| 日付 | 上弦/下弦側 (月齢)| 露出時間(s)xセット数 | 天...
| 2021.05.22 | 上(10d) | 60s x (h分) | 晴れ->曇 | | -...
| 2021.07.19 | 上(9d) | 30s x (h分) | 晴れ | | - |- ...
| 2021.09.15 | 上(9d) | 60s x (2h分) | 晴れ | | - |-...
| 2021.12.13 | 上(9d) | 30s x (7h分) | 晴れ | | - |-...
| 2022.01.10 | 上(7d) | 50s x (3h分) | 薄曇り/晴れ | | ...
| 2022.01.11 | 上(8d) | 50s x (3h分) | 薄曇り/晴れ | | -...
| 2022.01.24 | 下(22d) | 90s x (1h分) | 晴れ->曇り | | -...
| 2022.02.24 | 下(24d) | 120s x (2h分) | 快晴 | | - |- ...
| 2022.02.25 | 下(25d) | 180s x (1h分) | 快晴 | | - |- ...
| 2022.05.07 | 上(7d) | 75s x (2h分) | 晴 | | - |- | ...
| 2022.6.07 | 上(8d) | 70s x (実質1h分) | 晴 | | - |- ...
| 2024.2.14 | 上(4d) | 60s x (a few sets)| 曇り | | - |...
* 長時間データ [#hcc54d4f]
| date | delta_t | UT range | wax/wan | P | F_ocean| com...
| 191121 | 2.2 h | 18.5-21| wan | 変化なし | 80->95 % |...
| 191218 | 3.2 h | 16-18.5 | wan | 増加? | 75->95 % | -...
| 191219 | h | | wan | 増加? | | - |
| 200103 | 4.5 h | 10-14.5 | wax | 減少 | 70->40 %| - |
| 200115 | 4.9 h | 15-20 | wan | ? | 75->100| obs data ...
| 200302 | 2.7 h | 9.5-12.5 | wax | 減少 | 60->40| - |
| 200403 | 5.4 h | 10.5-15.5 | wax | 減少 | 48?->30->35...
| 200429 | 5.4 h | 10.5-14.5 | wax | ほとんど変化なし |...
* 解析 [#f8696597]
- 画像処理(観測中にやる)
@gamera
$ runnicred.sh -g[o] tonight
- 偏光導出
@jabulani2
$cd /mnt/hdd1/research/Earthshine/nic/YYMMDD
$../../task/calnicpolES.py --sky upper/lower [-r] *.ps...
-- calnicpolES.py
usage: ./calnicpolES.py fits1 fits2 ... [-s upper/lower]...
positional arguments:
fits FITS files
optional arguments:
-h, --help show this help message and exit
-s {upper,lower}, --sky {upper,lower}
sky position: "upper" (larger Y) ...
-r, --no_rename do not rename final result files.
- 全結果グラフ化
@jabulani2
$cd /mnt/hdd1/research/Earthshine/nic/results
$../../task/plotResults.py [-r] [-t]
-- plotResults.py
usage: ./plotResults.py [-r] [-t] [-h]
optional arguments:
-h, --help show this help message and exit
-r, --no_rename do not rename final result files.
-t, --two_dim use 2D qu results
* To Dos [#hfce6288]
- calnicpolES.py 改良
-- 2D q, u を計算するオプション (済)
-- q, u の時刻を独立に (19/12/6 済)
-- skyの正しい位置(上 or 下)をコマンドライン引数で与える ...
-- NG データ自動除外 (19/12/6 済)
-- 露出平均時刻を出して、位相角や散乱面方位角を出す(19/12...
-- 0D と 2D を一緒にプロットする (19/12/7 済)
-- python3化(20/02/04 済)
-- 結果の採否基準を明確に。例えば
--- |theta-N| > 30d
--- P > 2% かつ sigmaP/P > 1
- Plot PD vs phase angle, PA vs normal to scat. plane (19...
- plotResults.py
-- 結果まとめ表を作る (19/12/15 済)
-- JHK別に PD vs phase ang (済)
- sampling region を変えてみる(下弦期「月の海」を避ける)
- MODIS から地球の scene を抽出する
終了行:
[[staff]]
NICでの地球照偏光観測についてのページ(PI: 高橋)
* 観測手順 [#ic7284ad]
** 事前準備 [#k10dabec]
- なゆたファイル駆動用のトラックファイルを作り、UCCの
/export/home/observer/NHAO/FILEMOVE
に保存しておく
*** 詳細 [#q78d3d37]
-- 個人 MacBookPro (inyanga) でターミナルを開き
$ cd /Users/junon/Sync/Documents/20YYdo/obsprep/Earthsh...
-- lunar2radec.awk をエディタで開き、観測対象の月座標 (lm...
-- 2021年6月:ターゲット位置を変更した。Montanes-Rodrigıue...
--- 上弦期: Grimaldi side の patch #3: 月座標 (lat=-5.0...
--- 下弦期: Crisium side の patch #3: 月座標 (lat=0.00,...
-- トラックファイルを作り、%%gamera%%(<--ftpがない) OBS1...
$ ./prepLunarObs.py [start_date(yy-mm-dd) start_time(h...
(例)./prepLunarObs.py 2021-05-26 10:00:00
--- 引数で指定した開始時刻から5時間分(15秒間隔)のトラック...
-- OBS1b に ssh ログイン
-- ucc に ftp ログインして、上記のファイル置き場に移動、...
** 導入まで [#w53c41fe]
- ucc_run: 近赤外線観測モード、視野回転=(モード= AZ/EL同...
- calfocus でざっくり焦点合わせ
- shutterはopenにしておく
- HWP挿入
$ pol in
- ステラナビで月に近い、6-8等くらいの星に望遠鏡を向ける
- 以下のコマンドで、星を偏光視野に入れる。そのときのオフ...
$TL 3 2; finpol
$xyoffset2 (左側の偏光窓に星を移動)
$TL 3 2
- 偏光素子を挿入
$darkshutter.py -p
- focus値を +0.2 して、ざっくり焦点を合わせる
- autofocus で焦点合わせ
- ucc_run ファイル駆動モードで、月の地球照側に pointing
- さきほどのオフセット量を適用
- テスト撮影しながら、視野を調整。月の端が視野の真ん中に...
-- オフセット適用後は、望遠鏡は月の内側に向きすぎているこ...
--- 2021年5月以降、月座標でオフセットする
obs1b% lunarOffset.sh 月座標で北へのオフセット量 月座標...
オフセット量は秒角。月座標での東西と、赤道座標での東西(R...
- 地球照観測
$ PLo 3 (exp) Earthshine (repeat) ;fin
- 30分ごとを目安にスカイ観測
-- スカイ領域にpointing。上弦期は del_RA>0 (東に)、下弦期...
$ nayuta -do (del_RA) 0
-- スカイ観測
$ PLo 3 (exp) Earthshine_sky 1 ;fin
- 天候が安定していたら、nicpc:~/takahashi/obsScript/ESpol...
vi ESpol.sh
で変数を編集(exptime, dRAなど)。(地球照 x $series回 + ス...
* status [#o57a13ba]
- 第1期観測
| 日付 | 上弦/下弦側 (月齢)| 露出時間(s)xセット数 | 天...
| 2019.05.29 | 下 (24d) | 120s x 3 | 晴れ | bad | 再済 |...
| 2019.06.10 | 上 (6d) | 120s x 2 | 曇り | NG | 済み | ...
| 2019.06.25 | 下(22d) | (45s, 30s, 20s) x many | 薄曇...
| 2019.07.07 | 上 (4d) | 180s x 2 | 薄曇り | たぶんNG |...
| 2019.08.25 | 下 (24d) | 120s x 12 | 快晴 | good | 済...
| 2019.08.26 | 下 (25d) | 60s x 8 | 薄曇り | so so | 済...
| 2019.10.06 | 上 (7d) | 60s x 3 | 薄曇り->晴れ | bad |...
| 2019.10.21 | 下 (23d) | 60s x8, 120s x7 | 曇り->晴れ ...
| 2019.10.22 | 下 (24d) | 60s x many | 曇り->晴れ | bad...
| 2019.11.01 | 上 (4d) | 60s x 3 | 曇り->晴れ | so so |...
| 2019.11.05 | 上 (8d) | 60s x 13 | 晴れ | good | 再済 ...
| 2019.11.06 | 上 (9d) | 60s x 16 | 晴れ->少し雲? | goo...
| 2019.11.20 | 下 (24d) | 120s x 3 | 薄曇り | so so | ...
| 2019.11.21 | 下 (25d) | 180s x 10 | 晴れ | good | 再...
| 2019.12.04 | 上 (7d) | 60s x 7 | 薄曇り | so so | 済 ...
| 2019.12.06 | 上 (10d) | 5s x many | 曇り | NG | - | -...
| 2019.12.07 | 上 (11d) | 30s x 8 | 曇り->晴れ | so so...
| 2019.12.08 | 上 (12d) | 25s x 17 | 晴れ | so so | 済...
| 2019.12.09 | 上 (13d) | 25s x 25 | 晴れ | so so | 済...
| 2019.12.10 | 上 (14d) | 2s x | 薄曇り | NG | - | -...
| 2019.12.14 | 下 (18d) | 20s x | 曇り | NG | - | - ...
| 2019.12.15 | 下 (19d) | 2s x many | 薄曇り | NG | - ...
| 2019.12.18 | 下 (22d) | 120s, 90s x (3h分) | 晴れ->...
| 2019.12.19 | 下 (23d) | 90s x many (2.5h分) | 晴れ->...
| 2019.12.20 | 下 (24d) | 15s x many | 薄曇り | たぶ...
| 2019.12.22 | 下 (26d) | 100s x 2 | 晴れ->薄曇り | b...
| 2019.12.28 | 上(2d) | 20s x 7 | 薄曇り<->晴れ | so ...
| 2020.01.02 | 上(7d) | 90s x 13 | 薄曇り | bad | ...
| 2020.01.03 | 上(8d) | 60s x (4hくらい) | 晴れ | good ...
| 2020.01.04 | 上(9d) | 60s x 2 | 曇り | bad | 済 | 済...
| 2020.01.05 | 上(10d) | 40s x 11 | 曇り<->晴れ | so so ...
| 2020.01.15 | 下(20d) | 45s x 約5時間分 | 薄曇り | so s...
| 2020.01.17 | 下(23d) | 30s x | 薄曇り | bad - so so ...
| 2020.01.18 | 下(24d) | 90s x | 快晴 | good | 済 | ...
| 2020.01.30 | 上(5d) | 60s x 6 | 曇り | bad | 済 | 済 ...
| 2020.02.01 | 上(8d) | 60s x 7 | 晴れ/曇り | so so| 済...
| 2020.02.02 | 上(8d) | 60s x 4 | 曇り | bad| 済 | しな...
| 2020.02.03 | 上(9d) | 60s x 7 | 曇り | bad| 済 | 済 |...
| 2020.02.04 | 上(10d) | 30s x 8 | 晴れ | so so| 済 | ...
| 2020.02.05 | 上(11d) | 30s x 8 | 晴れ、時々雲 | so so|...
| 2020.02.06 | 上(12d) | 15s x 7, 25s x 3 | 薄曇り | ba...
| 2020.02.13 | 下(20d) | 30s x 10? | 晴れ | bad| 済 | ...
| 2020.02.26 | 上(3d) | 30s x | 晴れ | good| 済 | 済 ...
| 2020.03.02 | 上(8d) | 60s x (2.5h分) | 晴れ | good| ...
| 2020.03.14 | 下(20d) | 30s x 6 | 薄曇り | so so | 済 ...
| 2020.04.02 | 上(9d) | 30s x 3 | 薄曇り | so so | 済 |...
| 2020.04.03 | 上(10d) | 20s x (5.5h分) | 薄曇り | so s...
| 2020.04.29 | 上(6d) | 60s x (4h分) | 晴れ(ごく薄い雲...
- 第2期観測
| 日付 | 上弦/下弦側 (月齢)| 露出時間(s)xセット数 | 天...
| 2021.05.22 | 上(10d) | 60s x (h分) | 晴れ->曇 | | -...
| 2021.07.19 | 上(9d) | 30s x (h分) | 晴れ | | - |- ...
| 2021.09.15 | 上(9d) | 60s x (2h分) | 晴れ | | - |-...
| 2021.12.13 | 上(9d) | 30s x (7h分) | 晴れ | | - |-...
| 2022.01.10 | 上(7d) | 50s x (3h分) | 薄曇り/晴れ | | ...
| 2022.01.11 | 上(8d) | 50s x (3h分) | 薄曇り/晴れ | | -...
| 2022.01.24 | 下(22d) | 90s x (1h分) | 晴れ->曇り | | -...
| 2022.02.24 | 下(24d) | 120s x (2h分) | 快晴 | | - |- ...
| 2022.02.25 | 下(25d) | 180s x (1h分) | 快晴 | | - |- ...
| 2022.05.07 | 上(7d) | 75s x (2h分) | 晴 | | - |- | ...
| 2022.6.07 | 上(8d) | 70s x (実質1h分) | 晴 | | - |- ...
| 2024.2.14 | 上(4d) | 60s x (a few sets)| 曇り | | - |...
* 長時間データ [#hcc54d4f]
| date | delta_t | UT range | wax/wan | P | F_ocean| com...
| 191121 | 2.2 h | 18.5-21| wan | 変化なし | 80->95 % |...
| 191218 | 3.2 h | 16-18.5 | wan | 増加? | 75->95 % | -...
| 191219 | h | | wan | 増加? | | - |
| 200103 | 4.5 h | 10-14.5 | wax | 減少 | 70->40 %| - |
| 200115 | 4.9 h | 15-20 | wan | ? | 75->100| obs data ...
| 200302 | 2.7 h | 9.5-12.5 | wax | 減少 | 60->40| - |
| 200403 | 5.4 h | 10.5-15.5 | wax | 減少 | 48?->30->35...
| 200429 | 5.4 h | 10.5-14.5 | wax | ほとんど変化なし |...
* 解析 [#f8696597]
- 画像処理(観測中にやる)
@gamera
$ runnicred.sh -g[o] tonight
- 偏光導出
@jabulani2
$cd /mnt/hdd1/research/Earthshine/nic/YYMMDD
$../../task/calnicpolES.py --sky upper/lower [-r] *.ps...
-- calnicpolES.py
usage: ./calnicpolES.py fits1 fits2 ... [-s upper/lower]...
positional arguments:
fits FITS files
optional arguments:
-h, --help show this help message and exit
-s {upper,lower}, --sky {upper,lower}
sky position: "upper" (larger Y) ...
-r, --no_rename do not rename final result files.
- 全結果グラフ化
@jabulani2
$cd /mnt/hdd1/research/Earthshine/nic/results
$../../task/plotResults.py [-r] [-t]
-- plotResults.py
usage: ./plotResults.py [-r] [-t] [-h]
optional arguments:
-h, --help show this help message and exit
-r, --no_rename do not rename final result files.
-t, --two_dim use 2D qu results
* To Dos [#hfce6288]
- calnicpolES.py 改良
-- 2D q, u を計算するオプション (済)
-- q, u の時刻を独立に (19/12/6 済)
-- skyの正しい位置(上 or 下)をコマンドライン引数で与える ...
-- NG データ自動除外 (19/12/6 済)
-- 露出平均時刻を出して、位相角や散乱面方位角を出す(19/12...
-- 0D と 2D を一緒にプロットする (19/12/7 済)
-- python3化(20/02/04 済)
-- 結果の採否基準を明確に。例えば
--- |theta-N| > 30d
--- P > 2% かつ sigmaP/P > 1
- Plot PD vs phase angle, PA vs normal to scat. plane (19...
- plotResults.py
-- 結果まとめ表を作る (19/12/15 済)
-- JHK別に PD vs phase ang (済)
- sampling region を変えてみる(下弦期「月の海」を避ける)
- MODIS から地球の scene を抽出する
ページ名: