Changeset 492

Timestamp:

02/21/05 18:03:36 (20 years ago)

Author:

kil064

Message:

add function polarizationLabel

Location:

trunk/src

Files:

• SDPol.cc (modified) (3 diffs)
• SDPol.h (modified) (2 diffs)

trunk/src/SDPol.cc

@@ -328 +328 @@
    if (toStokes) {
       if (val==0) {
-         stokes = Stokes::I;
+          stokes = Stokes::I;
       } else if (val==1) {
          stokes = Stokes::Q;
@@ -361 +361 @@
 }
 
+
+
+String SDPolUtil::polarizationLabel (uInt polIdx, Bool linear, Bool stokes, Bool linPol)
+{   
+   Stokes::StokesTypes type = Stokes::Undefined;
+   if (stokes) {
+      switch (polIdx) {
+         case 0:
+           {
+              type = Stokes::I;
+           }
+           break;
+         case 1:
+           {
+              if (linPol) {
+                 type = Stokes::Plinear;
+              } else {
+                 type = Stokes::Q;
+              }
+           }
+           break;
+         case 2:
+           {
+              if (linPol) {
+                 type = Stokes::Pangle;
+              } else {
+                 type = Stokes::U;
+              }
+           }
+           break;
+         case 3:
+           {
+              type = Stokes::V;
+           }
+           break;
+         default:  
+           {
+               throw(AipsError("Unknown Stokes type"));
+           }
+      }
+   } else {
+      if (linear) {
+         switch (polIdx) {
+            case 0:
+              {
+                 type = Stokes::XX;
+              }
+              break;
+            case 1:
+              {
+                 type = Stokes::YY;
+              }
+              break;
+            case 2:
+              {
+                 type = Stokes::XY;              // Really Real(XY)
+                 return String("Real(XY)");
+              }
+              break;
+            case 3:
+              {
+                 type = Stokes::YX;              // Really Imag(XY)
+                 return String("Imag(XY)");
+              }
+              break;
+            default:  
+              {
+                  throw(AipsError("Unknown linear polarization type"));
+              }
+         }
+      } else {  
+         switch (polIdx) {
+            case 0:
+              {
+                 type = Stokes::RR;
+              }
+              break;
+            case 1:
+              {
+                 type = Stokes::LL;
+              }
+            case 2:
+              {
+                 type = Stokes::RL;               // Really Real(RL)
+                 return String("Real(RL)");
+              }
+              break;
+            case 3:
+              {
+                 type = Stokes::LR;               // Really Imag(RL)
+                 return String("Imag(RL)");
+              }
+              break;
+            default:  
+              {
+                  throw(AipsError("Unknown circular polarization type"));
+              }
+         }
+      }
+   }
+//
+   return SDPolUtil::stokesString(type);
+} 
+
+
+
+
+// private
+
+String SDPolUtil::stokesString (Stokes::StokesTypes type)
+{
+  return Stokes::name (type);
+}
+
+
 Array<casa::uChar> SDPolUtil::andArrays (const Array<casa::uChar>& in1,
                                          const Array<casa::uChar>& in2)
@@ -404 +519 @@
 }
 
+
+
+
+

trunk/src/SDPol.h

@@ -194 +194 @@
    static casa::Stokes::StokesTypes convertStokes(casa::Int val, casa::Bool toStokes, 
                                                   casa::Bool linear);
+
+// Return a label for the given polarization index (0-3).   
+   static casa::String polarizationLabel (casa::uInt polIdx, casa::Bool linear, 
+                                          casa::Bool stokes, casa::Bool linPol);
 
 
@@ -220 +224 @@
 
 private:
+// Return a label for the given StokesType
+   static casa::String stokesString (casa::Stokes::StokesTypes stokes);
+
+// specializations
   static casa::Array<casa::Float> andArrays (const casa::Array<casa::Float>& in1,
                                              const casa::Array<casa::Float>& in2) 
                                              {return (in1+in2)/casa::Float(2.0);}
-
   static casa::Array<casa::Bool> andArrays (const casa::Array<casa::Bool>& in1,
                                             const casa::Array<casa::Bool>& in2)