History of CSV/OpticalPointingAnalysis(No. 1)

History
View the source.
CSV/OpticalPointingAnalysis has been deleted.
- 1 (2016-08-05 (Fri) 17:33:12)
- 2 (2016-08-06 (Sat) 23:48:55)
- 3 (2016-08-07 (Sun) 12:42:28)
- 4 (2016-08-07 (Sun) 14:13:20)
- 5 (2016-08-07 (Sun) 19:46:37)
- 6 (2016-08-10 (Wed) 16:07:18)
- 7 (2016-08-11 (Thu) 12:15:59)

Commissioning and Science Verification

Optical Pointing > Analysis†

What's this page?: This page describes the way to analyze optical pointing data.

Revision history

2010-04-13 - modified by B. Hatsukade (astepon2010_red.txt)
2010-05-17 - BH - modified for wiki by B. Hatsukade
2014-04-23 - updated by BH and YW
2016-08-01 - modified for Linux COSMOS3 by Kamazaki,T.

Table of contents

Optical Pointing > Analysis
Troubleshoot
- matlab スクリプトを走らせた際、以下のようなエラー

↑

0. MATLAB準備†

解析ファイル用ディレクトリ作成

mkdir ~/matlab/12_ASTE_2014PM
mkdir ~/matlab/12_ASTE_2014PM/tabs_i

出力ファイル保存用ディレクトリ作成

cd ~/matlab/12_ASTE_2014PM
mkdir Figd
mkdir Figd/figs_i
mkdir tab2mat_i

解析スクリプトをコピー（以前の解析ディレクトリから）
```
cp -p ~/matlab/12_ASTE_2013PM/*.m .
```

↑

1. データコピー†

アカウント: astepon
aste-2c の /home/astepon/optpon にデータファイルができているのを確認 (200908102207.tabなど、日付と時刻のついた.tabファイル)
```
% ssh astepon@aste-2c
aste-2c{astepon}41: ttt
aste-2c{astepon}41: cd /home/astepon/optpon
```

aste-mx2.mtk.nao.ac.jp の /home/astepon/matlab/ に 12_ASTE_2014PM/tabs_i ディレクトリを作成し、データファイルをコピー
```
cd ~/matlab/12_ASTE_2014PM/tabs_i
scp -p astepon@aste-2c:/home/astepon/optpon/201009090238.tab .
```

↑

2. 解析†

aste-mx2.mtk.nao.ac.jp:/home/astepon/matlab/12_ASTE_2014PM/tabs_i に入っている .tab ファイルが順に解析される。

MATLABを起動
```
cd ~/matlab/12_ASTE_2014PM/
matlab & 
```

↑

方法1†

MATLABツールバーの File -> Open から、解析ツール step_02_tab2mat2pm_iOPT_A22c.m (浮田さん作成) を選択

F5キーでプログラムを実行

※ 測定時に使用したプログラムによって出力形式が異なり、解析ツールも異なる：

--------------------------------------------------
測定プログラム     解析プログラム
iOPT_r17c      ==> step_02_tab2mat2pm_iOPT_A22c.m
tracking_c2n   ==> stp_72_tab2mat2pm_v2.m
--------------------------------------------------

↑

方法2†

MATLAB terminalで

>> cd /home/astepon/matlab/12_ASTE_2010PM
>> step_02_tab2mat2pm_iOPT_A22c  [enter] (.m 拡張子はいらない)

↑

3. 結果†

↑

器差パラメータのフィット結果†

ポインティングモデルの16パラメータのうち、P9, P10を除いた14パラメータに対する補正値がmatlab上に出力される。
複数のファイルを解析した場合は、まず個々の結果が表示され、最後にそれらを平均した値が表示される。
上段はdegree, 下段はarcsec。器差更新には上段の値を現在使用されている値に足し合わせることになる。

ans =
0.00041  0.00017  0.00081 -0.00006  0.00001 -0.00011  0.00005 -0.00003 -0.00002 -0.00002 -0.00016  0.00012  0.00003  0.00005  
ans =
1.48944  0.59760  2.92666 -0.21296  0.04278 -0.40678  0.16991 -0.09039 -0.08448 -0.05949 -0.57215  0.43172  0.09684  0.18882

↑

プロット†

プロットは、/home/astepon/matlab/12_ASTE_2010PM/Figd/figs_i に保存される。
複数のファイルを解析した場合は、個々のファイルの結果が順に表示されます。一つのファイルのプロットが表示された後、次のファイルに移ると、以前のプロットが上書きされます。個々の結果を保存したい場合は、ディレクトリ（/home/astepon/matlab/12_ASTE_2010PM/tabs_i）にあるファイルを一つにしておく必要があります。
複数のファイルを解析した場合でも、プロットに関してはそれらを平均した結果は表示されません。

Figure1: Az error vs. El error

#ref(): File not found: "Figure1.png" at page "CSV/OpticalPointingAnalysis"

Figure2: ?

#ref(): File not found: "Figure2.png" at page "CSV/OpticalPointingAnalysis"

Figure101: 観測点

#ref(): File not found: "Figure101.png" at page "CSV/OpticalPointingAnalysis"

Figure201:

左上: Az vs. Az error. filename.
右上: El vs. Az error. Az error rms.
左下: Az vs. EL error. number of targets, Az El error rms.
右下: El vs. EL error. EL error rms.
#ref(): File not found: "Figure201.png" at page "CSV/OpticalPointingAnalysis"

Figure301: Figure201にポインティングモデルを適用した場合のプロット

#ref(): File not found: "Figure301.png" at page "CSV/OpticalPointingAnalysis"

Figure401

#ref(): File not found: "Figure401.png" at page "CSV/OpticalPointingAnalysis"

Figure501

#ref(): File not found: "Figure501.png" at page "CSV/OpticalPointingAnalysis"

Figure601

#ref(): File not found: "Figure601.png" at page "CSV/OpticalPointingAnalysis"

↑

5．器差ファイル更新†

いくつかの測定データに対するフィット結果を平均し元の器差パラメータに加える。(P9,P10の2パラメータを固定し、後の14パラメータでフィットしている。）
「主鏡器差ファイル更新手順」にしたがって器差ファイルを入れ替える。

↑

Troubleshoot†

↑

matlab スクリプトを走らせた際、以下のようなエラー†

??? Index exceeds matrix dimensions.

Error in ==> step_02_tab2mat2pm_iOPT_A22c at 26
        az = str2double(nextline(34:41));    dat(icount,1)=az;

??? Error using ==> feof
Invalid file identifier.  Use fopen to generate a valid file identifier.

Error in ==> step_02_tab2mat2pm_iOPT_A22c at 18
    while feof(fid) == 0

.tab ファイルの中身がおかしい可能性。中身をチェックして、データが欠損している行を削除。

History of CSV/OpticalPointingAnalysis (No. 1)