Diff of CSV/OpticalPointingAnalysis - ASTE Observation Wiki

The added line is THIS COLOR.
The deleted line is THIS COLOR.
CSV/OpticalPointingAnalysis has been deleted.
Deleting diff of CSV/OpticalPointingAnalysis
[[Commissioning and Science Verification]]

* Optical Pointing > Analysis [#g0858f28]

: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|
#contents

----

** Preparations [#a6f00c07]

+ Login to the ASTE analysis computer (aste-ms1) with asteXX (e.g. aste16)
 > ssh -Y asteXX@aste-ms1
-- Don't forget to add '-Y' option to forward windows.
+ Create "matlab" directory and rsync MATLAB programs
 asteXX@aste-ms1: mkdir -p ~/optpon/matlab
 asteXX@aste-ms1: cd ~/optpon/matlab
 asteXX@aste-ms1: pwd
 /home/asteXX/optpon/matlab
 asteXX@aste-ms1: rsync -av --delete ~astecal/astepon/matlab/12_ASTE_20XX/ ./12_ASTE_YYYYMMDD/    (YYYYMMDD=current date; e.g. 2016-08-07 -> 20160807)
 ...

** Data reduction [#jc6a2663]

+ Rsync data files from ASTE to Mitaka
 asteXX@aste-ms1: cd ~/optpon/matlab/12_ASTE_YYYYMMDD
 asteXX@aste-ms1: pwd
 /home/asteXX/optpon/matlab/12_ASTE_YYYYMMDD
 asteXX@aste-ms1: rsync -e ssh -avz astecntl@aste-1c:./optpon/"*.tab" ./
 ...
-- There will be data files, whose names are "YYYYMMDDHHMM.tab" (year+month+day+hour+minute) like "201608111234.tab", under aste-c1c:~astecntl/optpon.
+ Concatenate all .tab files taken on the day
 asteXX@aste-ms1: cat *.tab ./tabs_i/YYYYMMDD0000.tab
-- All .tab files under the "tabs_i" directory are individually reduced by MATLAB and then, it outputs individual results of the files and their averaged parameters. ''It is necessary to put all .tab files taken on the day into a .tab file in order to do fitting for the whole data of the files.''
+ Run MATLAB and execute either of the MATLAB scripts according to the OPT programs
 asteXX@aste-ms1: cd ~/optpon/matlab/12_ASTE_YYYYMMDD/
 asteXX@aste-ms1: /usr/local/matlab/R2009a/bin/matlab &
 matlab >> pwd
 /home/asteXX/optpon/matlab/12_ASTE_YYYYMMDD
+
 [iOPT_r17c]
 matlab >> step_02_tab2mat2pm_iOPT_A22c
+
 [tracking_c2n] (This program is NOT used since at least 2015.)
 matlab >> stp_72_tab2mat2pm_v2

*** 4. 器差パラメータのフィット結果 [#a9d29924]
- ポインティングモデルの16パラメータのうち、P9, P10を除いた14パラメータに対する補正値がmatlab上に出力される。
- 複数のファイルを解析した場合は、まず個々の結果が表示され、最後にそれらを平均した値が表示される。
- 上段はdegree, 下段はarcsec。器差更新には上段の値を現在使用されている値に足し合わせることになる。
   P1       P2       P3       P4       P5       P6       P7       P8       P11      P12      P13      P14      P15      P16
 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 

*** 5. Plots [#h3d6a78f]
- プロットは、~/optpon/matlab/12_ASTE_YYYYMMDD/Figd/figs_iに保存される。
- 複数のファイルを解析した場合は、個々のファイルの結果が順に表示されます。一つのファイルのプロットが表示された後、次のファイルに移ると、以前のプロットが上書きされます。個々の結果を保存したい場合は、ディレクトリ（/home/astepon/matlab/12_ASTE_2010PM/tabs_i）にあるファイルを一つにしておく必要があります。
- 複数のファイルを解析した場合でも、プロットに関してはそれらを平均した結果は表示されません。
- 全てのデータをまとめた結果をプロットしたい場合は、catでファイルを1つに結合して処理して下さい。
- tabs_i内にあるファイル・ディレクトリを全て処理しようとしますので、処理したくないファイル・ディレクトリをtabs_iに置かないようにして下さい。

- Figure1: Az error vs. El error
#ref(Figure1_20150406.png);
- Figure2: ?
#ref(Figure2_20150406.png);
- Figure101: 観測点
#ref(Figure101_20150406.png);
- Figure201: 
#ref(Figure201_20150406.png);
-- Upper-Left: Az vs. Az error. filename. 
-- Upper-Right: El vs. Az error. Az error rms. 
-- Lower-Left: Az vs. EL error. number of targets, Az El error rms. 
-- Lower-Right: El vs. EL error. EL error rms. 
- Figure301: Figure201にポインティングモデルを適用した場合のプロット
#ref(Figure301_20150406.png);
- Figure401
#ref(Figure401_20150406.png);
- Figure501
#ref(Figure501_20150406.png);
- Figure601
-- This plot will be displayed, when total number of observed stars exceeds 118.
-
#ref(Figure601_20150406.png);

** Troubleshooting [#h28e6671]
- "??? Index exceeds matrix dimensions." during the execution of the MATLAB script
 ??? 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
-- There could be invalid lines in .tab files. Search such lines and delete them.

// End of the page