compareAntennaPositionCorrections
Return to directory of Todd's CASA extensions
This command computes the difference (in mm) between two antenna baseline solutions
on a per antenna basis. The input format is the same as for gencal, and can be used
to compare the output from two different executions of
correctMyAntennaPositions.
It returns dictionaries of the vector differences and magnitudes, keyed by antenna names.
Usage:
au.compareAntennaPositionCorrections(antenna1='', parameter1='', antenna2='', parameter2=''
gencal1='', gencal2='')
Inputs:
- antenna1='DV03,DV04'
- antenna2='DV03,DV04'
- parameter1=[-0.00011,3.965e-08,0.000126,-0.0001939,-2.638e-05,0.0001716]
- parameter2=[-0.00012,3.966e-08,0.000136,-0.0001839,-2.648e-05,0.0002716]
- gencal1: a textfile from which to locate an antpos gencal command and read the antenna and parameter
- gencal2: a textfile from which to locate an antpos gencal command and read the antenna and parameter
Example:
CASA <2>: au.compareAntennaPositionCorrections(\
antenna1 = 'PM01,PM03,PM02,DV11,DV10,DV13,DV12,DA41,DV14,DA43,DV15',
antenna2='PM01,PM03,PM02,DV06,DV10,DV13,DV12,DA41,DV03,DA43,DV14,DV11,DV15',
parameter1=[
-0.000111298679883,3.96572241512e-08,0.000126355064337, # PM01
-0.000193905085325,-2.63825058937e-05,0.000171698629856, # PM03
-6.5671095879e-05,-1.64214945227e-05,8.76620880636e-05, # PM02
0.000164386350662,-0.000342965126038,-0.000527559313923, # DV11
0.00013935379684,-0.00034146849066,-0.000170580577105, # DV10
-8.14814120531e-05,-1.16489827633e-05,0.00020522531122, # DV13
9.39136426995e-05,3.4540914983e-05,0.000139599974595, # DV12
-3.2177194953e-06,-0.000219983048737,8.33319500089e-05, # DA41
-4.16459692004e-05,-1.19805462951e-06,0.000186379124138, # DV14
-1.64100158093e-05,-8.66703933952e-05,3.38814741404e-05, # DA43
1.00582838058e-06,-1.0309740901e-06,2.00234353542e-08],
parameter2=[-0.000108526088297,-0.000314249657094,1.90567225218e-05, # PM01
-0.000369362067431,-8.08034092188e-05,7.96969980001e-05, # PM03
-0.000227167271078,0.000142025761306,0.000128640793264, # PM02
-0.00011939695105,1.67889520526e-05,-9.28249210119e-06, # DV06
0.00013935379684,-0.00034146849066,-0.000170580577105, # DV10
-0.000178528949618,8.29827040434e-05,0.000201276037842, # DV13
9.59411263466e-05,-5.79198822379e-05,0.00025811418891, # DV12
-0.000190026126802,0.000115232542157,0.000217058695853, # DA41
-4.80054877698e-05,0.000225332565606,0.000178911723197, # DV03
.000319970455725,-0.000502289491245,-0.000319080200653, # DA43
-0.000223453178561,0.000216139004078,0.000189656552503, # DV14
0.000164386350662,-0.000342965126038,-0.000527559313923, # DV11
0.000573113560677,-0.000985723920166,-0.000957812648267])
PM01 -0.003 +0.314 +0.107 total = +0.332 mm
PM03 +0.175 +0.054 +0.092 total = +0.205 mm
PM02 +0.161 -0.158 -0.041 total = +0.230 mm
DV11 +0.000 +0.000 +0.000 total = +0.000 mm
DV10 +0.000 +0.000 +0.000 total = +0.000 mm
DV13 +0.097 -0.095 +0.004 total = +0.136 mm
DV12 -0.002 +0.092 -0.119 total = +0.150 mm
DA41 +0.187 -0.335 -0.134 total = +0.406 mm
DV14 +0.182 -0.217 -0.003 total = +0.283 mm
DA43 -0.336 +0.416 +0.353 total = +0.641 mm
DV15 -0.572 +0.985 +0.958 total = +1.488 mm
Standard deviation = +0.399 mm
Out[2]:
({'DA41': array([ 0.18680841, -0.33521559, -0.13372675]),
'DA43': array([-0.33638047, 0.4156191 , 0.35296167]),
'DV10': array([ 0., 0., 0.]),
'DV11': array([ 0., 0., 0.]),
'DV12': array([-0.00202748, 0.0924608 , -0.11851421]),
'DV13': array([ 0.09704754, -0.09463169, 0.00394927]),
'DV14': array([ 0.18180721, -0.21733706, -0.00327743]),
'DV15': array([-0.57210773, 0.98469295, 0.95783267]),
'PM01': array([-0.00277259, 0.31428931, 0.10729834]),
'PM02': array([ 0.16149618, -0.15844726, -0.04097871]),
'PM03': array([ 0.17545698, 0.0544209 , 0.09200163])},
{'DA41': 0.40638616606982575,
'DA43': 0.64068166822540007,
'DV10': 0.0,
'DV11': 0.0,
'DV12': 0.15032873546785858,
'DV13': 0.13560596393782323,
'DV14': 0.28337254630575698,
'DV15': 1.4880762354078025,
'PM01': 0.332112021033518,
'PM02': 0.22992564398834003,
'PM03': 0.20545337075018663})
-- Todd Hunter
--
ToddHunter - 2012-03-27