1 | //
|
---|
2 | // C++ Implementation: STPolStokes
|
---|
3 | //
|
---|
4 | // Description:
|
---|
5 | //
|
---|
6 | //
|
---|
7 | // Author: Malte Marquarding <Malte.Marquarding@csiro.au>, (C) 2006
|
---|
8 | //
|
---|
9 | // Copyright: See COPYING file that comes with this distribution
|
---|
10 | //
|
---|
11 | //
|
---|
12 |
|
---|
13 | #include <casa/Arrays/ArrayMath.h>
|
---|
14 | #include <casa/BasicMath/Math.h>
|
---|
15 | #include <casa/Exceptions/Error.h>
|
---|
16 | #include <casa/BasicSL/Constants.h>
|
---|
17 | #include "STPolStokes.h"
|
---|
18 |
|
---|
19 | using namespace casa;
|
---|
20 |
|
---|
21 | namespace asap {
|
---|
22 |
|
---|
23 | Factory<STPol,STPolStokes> STPolStokes::myFactory;
|
---|
24 |
|
---|
25 | STPolStokes::~STPolStokes()
|
---|
26 | {
|
---|
27 | }
|
---|
28 |
|
---|
29 |
|
---|
30 | Vector<Float> asap::STPolStokes::getStokes( uint index )
|
---|
31 | {
|
---|
32 | return getSpectrum(index);
|
---|
33 | }
|
---|
34 |
|
---|
35 | Vector<Float> asap::STPolStokes::getLinPol( uInt index )
|
---|
36 | {
|
---|
37 | if ( index < 0 || index >4 ) throw(AipsError("LinPol index out of range"));
|
---|
38 | Vector<Float> out;
|
---|
39 | if ( nspec() == 4) {
|
---|
40 | switch(index) {
|
---|
41 | case 1:
|
---|
42 | out = Vector<Float>(sqrt(pow(getSpectrum(1),Float(2.0))
|
---|
43 | +pow(getSpectrum(2), Float(2.0))));
|
---|
44 | break;
|
---|
45 | case 2:
|
---|
46 | out = Vector<Float>(Float(180.0/C::pi/2.0)
|
---|
47 | * atan2(getSpectrum(2),getSpectrum(1)));
|
---|
48 | break;
|
---|
49 | default:
|
---|
50 | out = getSpectrum(index);
|
---|
51 | }
|
---|
52 | }
|
---|
53 | return out;
|
---|
54 | }
|
---|
55 |
|
---|
56 | Vector<Float> asap::STPolStokes::getLinear(uInt index )
|
---|
57 | {
|
---|
58 | Vector<Float> out;
|
---|
59 | switch(index) {
|
---|
60 | case 0:
|
---|
61 | out = (getSpectrum(0)+getSpectrum(1))/Float(2.0);
|
---|
62 | case 1:
|
---|
63 | out = (getSpectrum(0)-getSpectrum(1))/Float(2.0);
|
---|
64 | break;
|
---|
65 | default:
|
---|
66 | out = Vector<Float>();
|
---|
67 | }
|
---|
68 | return out;
|
---|
69 |
|
---|
70 | }
|
---|
71 |
|
---|
72 | Vector<Float> asap::STPolStokes::getCircular(uInt index )
|
---|
73 | {
|
---|
74 | // convert to stokes I/ V first
|
---|
75 | //
|
---|
76 | // We use the convention
|
---|
77 | // I = (RR+LL) // definition changed
|
---|
78 |
|
---|
79 | if ( index == 2 || index ==3 ) throw(AipsError("Re/Imag RL not implemented"));
|
---|
80 | Vector<Float> out;
|
---|
81 | switch(index) {
|
---|
82 | case 0:
|
---|
83 | out = (getSpectrum(0) + getSpectrum(3));
|
---|
84 | break;
|
---|
85 | case 1:
|
---|
86 | out = (getSpectrum(0) - getSpectrum(3));
|
---|
87 | break;
|
---|
88 | default:
|
---|
89 | out = Vector<Float>();
|
---|
90 | }
|
---|
91 | return out;
|
---|
92 | }
|
---|
93 |
|
---|
94 | }
|
---|