History of CSV/OpticalPointingAnalysis (No. 4)

• 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 - 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
  • Preparations
    • 1. MATLAB
  • Data reduction
    • 1. MATLAB
    • 2. 器差パラメータのフィット結果
    • 3. Plots
  • Troubleshooting
    • "??? Index exceeds matrix dimensions." during the execution of MATLAB scripts

Preparations†

1. MATLAB†

• rsync MATLAB programs

  asteXX@aste-ms1: cd ~/optpon/matlab
  asteXX@aste-ms1: rsync -av --delete /home/astecal/astepon/matlab/12_ASTE_20XX/ 12_YYYYMMDD/    (YYYYMMDD=current date; e.g. 2016-08-07 -> 20160807)
  ...

• Data copy [#a5845ebf]
• astecntl@aste-1c:/home/astecntl/optponにデータファイルができているのを確認 (200908102207.tabなど、日付と時刻のついた.tabファイル)

  asteXX@aste-ms1: cd optpon/matlab/tabs_i
  asteXX@aste-ms1: pwd
  /home/asteXX/optpon/matlab/12_ASTE_YYYYMMDD/tabs_i
  asteXX@aste-ms1: rsync -e ssh -avz astecntl@aste-1c:./optpon/"*.tab" ./

Data reduction†

1. MATLAB†

12_ASTE_YYYYMMDD/tabs_i に入っている .tab ファイルが順に解析される。

• Run MATLAB

  asteXX@aste-ms1: cd ~/optpon/matlab/12_ASTE_YYYYMMDD/
  asteXX@aste-ms1: /usr/local/matlab/R2009a/bin/matlab &
  matlab >> pwd
  /home/aste16/optpon/matlab/12_ASTE_YYYYMMDD

• Run either of the MATLAB scripts according to the OPT programs

  [iOPT_r17c]
  matlab >> step_02_tab2mat2pm_iOPT_A22c
  [tracking_c2n]
  matlab >> stp_72_tab2mat2pm_v2

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

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

  ans =
  -0.00022  0.00026 -0.00041 -0.00000 -0.00002 -0.00010 -0.00007  0.00007  0.00000 -0.00002  0.00004 -0.00003  0.00010  0.00001 
  ans =
  -0.80449  0.91804 -1.46657 -0.01760 -0.06549 -0.35213 -0.25378  0.25857   0.00461 -0.05692  0.13902 -0.12545  0.34517  0.02084 

3. Plots†

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

• Figure1: Az error vs. El error

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

• Figure2: ?

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

• Figure101: 観測点

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

• Figure201:
• #ref(Figure201_20150406.png);
  • 左上: 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.
• Figure301: Figure201にポインティングモデルを適用した場合のプロット

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

• Figure401

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

• Figure501

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

• Figure601

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

Troubleshooting†

"??? Index exceeds matrix dimensions." during the execution of MATLAB scripts†

??? 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 ファイルの中身がおかしい可能性。中身をチェックして、データが欠損している行を削除。