光学ポインティング自動測定の整約

last update 2006/12/6 Handa

1. 実行は、obs@cygで行なう。
2. ポインティングログファイルを取得する。 データは/obsdata/log/pt????????.datという名称で保存されている。 ここで、???の部分は観測開始の年月日である。 以下では、ptlogファイルと呼ぶ。
3. opintprmコマンドの引数は、引数を１つも指定しないと表示されるが、以下の通り。

   第１引数　読み込むべきptlogファイル

   第２引数　使用するlevel10の下限値（指定なき場合は0.0）

   第３引数　使用するlevel50の下限値（指定なき場合は0.0）

   第４引数　残差ファイル名（指定がない場合は作成しない）

   最終引数　--azelcmdの場合にはフィッティングに計算値のAzElの代わりにプログラム指定値を用いる。このオプションは推奨しない。

4. ptlogファイルに以下のコマンドを適用する。(opintprmは~/UltraASTROS/bin/に収録。)
   % opintprm ptlogファイル
   (e.g. % opintprm pn20061201.dat)
5. フィッティングは画面に表示されるとともに、自動的に結果ファイルに書かれる。 結果ファイルの名称は、kisadiff.????????.pfitとなる。 ここで、???の部分は整約実行の年月日である。
6. 画面表示または上記の結果ファイルを見て、器差ファイルを書き換える。 パラメータA1,A2,A3,B1,B2,B3をsctの~/UltraASTROS/dev/kisa.datに差分として 加えればよい。 また、このファイルは電波ポインティングの整約で hoseiAやrevisekisaなどが使用する。
7. 残差を見るためには、以下のコマンドを用いればよい。
   % cat logファイル | awk '{print($3" "$4" "$5" "$6);}' | estres A1 A2 A3 B1 B2 B3 > 残差ファイル
   (e.g. % cat pn130405.t | awk '{print($3" "$4" "$5" "$6);}' | estres A1 A2 A3 B1 B2 B3 > pn130405.pfit.res)
8. gnuplotで以下のコマンドを実行すると残差グラフが描ける。
   > plot '残差ファイル' using 5:6 with po
   (e.g. > plot 'pn130405.pfit.res' using 5:6 with po)