source: trunk/src/STMath.h@ 1097

Last change on this file since 1097 was 1078, checked in by mar637, 18 years ago

added scan averaging to average_channel

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 5.6 KB
Line 
1//
2// C++ Interface: STMath
3//
4// Description:
5//
6//
7// Author: Malte Marquarding <asap@atnf.csiro.au>, (C) 2006
8//
9// Copyright: See COPYING file that comes with this distribution
10//
11//
12#ifndef ASAPSTMATH_H
13#define ASAPSTMATH_H
14
15#include <string>
16#include <map>
17
18#include <casa/aips.h>
19#include <casa/Utilities/CountedPtr.h>
20#include <casa/BasicSL/String.h>
21#include <casa/Arrays/Vector.h>
22#include <scimath/Mathematics/InterpolateArray1D.h>
23
24#include "Scantable.h"
25#include "STDefs.h"
26#include "STPol.h"
27#include "Logger.h"
28
29namespace asap {
30
31/**
32Mathmatical operations on Scantable objects
33
34@author Malte Marquarding
35*/
36class STMath : private Logger {
37public:
38
39 typedef casa::InterpolateArray1D<casa::Double,
40 casa::Float>::InterpolationMethod imethod;
41
42 typedef std::map<std::string, imethod> imap;
43
44 STMath(bool insitu=true);
45
46 virtual ~STMath();
47
48 /**
49 * set the @attr insitu attribute
50 * @param b
51 */
52 bool insitu() const { return insitu_;};
53 void setInsitu(bool b) { insitu_ = b; };
54
55 casa::CountedPtr<Scantable>
56 average( const std::vector<casa::CountedPtr<Scantable> >& in,
57 const std::vector<bool>& mask = std::vector<bool>(),
58 const std::string& weight = "NONE",
59 const std::string& avmode = "SCAN");
60
61 casa::CountedPtr<Scantable>
62 averageChannel( const casa::CountedPtr<Scantable> & in,
63 const std::string& mode = "MEDIAN",
64 const std::string& avmode = "SCAN");
65
66 casa::CountedPtr< Scantable >
67 averagePolarisations( const casa::CountedPtr< Scantable > & in,
68 const std::vector<bool>& mask,
69 const std::string& weight );
70
71 casa::CountedPtr<Scantable>
72 unaryOperate( const casa::CountedPtr<Scantable>& in, float val,
73 const std::string& mode, bool tsys=false );
74
75 casa::CountedPtr<Scantable> autoQuotient(const casa::CountedPtr<Scantable>& in,
76 const std::string& mode = "NEAREST",
77 bool preserve = true);
78
79 casa::CountedPtr<Scantable> quotient( const casa::CountedPtr<Scantable>& on,
80 const casa::CountedPtr<Scantable>& off,
81 bool preserve = true );
82
83 casa::CountedPtr<Scantable>
84 freqSwitch( const casa::CountedPtr<Scantable>& in );
85
86 std::vector<float> statistic(const casa::CountedPtr<Scantable>& in,
87 const std::vector<bool>& mask,
88 const std::string& which);
89
90 casa::CountedPtr<Scantable> bin( const casa::CountedPtr<Scantable>& in,
91 int width=5);
92 casa::CountedPtr<Scantable>
93 resample(const casa::CountedPtr<Scantable>& in,
94 const std::string& method, float width);
95
96 casa::CountedPtr<Scantable>
97 smooth(const casa::CountedPtr<Scantable>& in, const std::string& kernel,
98 float width);
99
100 casa::CountedPtr<Scantable>
101 gainElevation(const casa::CountedPtr<Scantable>& in,
102 const std::vector<float>& coeff,
103 const std::string& fileName,
104 const std::string& method);
105 casa::CountedPtr<Scantable>
106 convertFlux(const casa::CountedPtr<Scantable>& in, float d,
107 float etaap, float jyperk);
108
109 casa::CountedPtr<Scantable> opacity(const casa::CountedPtr<Scantable>& in,
110 float tau);
111
112 casa::CountedPtr<Scantable>
113 merge(const std::vector<casa::CountedPtr<Scantable> >& in);
114
115 casa::CountedPtr<Scantable>
116 invertPhase( const casa::CountedPtr<Scantable>& in);
117
118 casa::CountedPtr<Scantable>
119 rotateXYPhase( const casa::CountedPtr<Scantable>& in, float phase);
120
121 casa::CountedPtr<Scantable>
122 rotateLinPolPhase( const casa::CountedPtr<Scantable>& in, float phase);
123
124 /// @todo frequency alignment
125
126 casa::CountedPtr<Scantable>
127 swapPolarisations(const casa::CountedPtr<Scantable>& in);
128
129 casa::CountedPtr<Scantable>
130 frequencyAlign( const casa::CountedPtr<Scantable>& in,
131 const std::string& refTime = "",
132 const std::string& method = "cubic" );
133
134 casa::CountedPtr<Scantable>
135 convertPolarisation( const casa::CountedPtr<Scantable>& in,
136 const std::string& newtype);
137
138private:
139 casa::CountedPtr<Scantable> applyToPol( const casa::CountedPtr<Scantable>& in,
140 STPol::polOperation fptr,
141 casa::Float phase);
142
143 static imethod stringToIMethod(const std::string& in);
144 static WeightType stringToWeight(const std::string& in);
145
146 void scaleByVector(casa::Table& in,
147 const casa::Vector<casa::Float>& factor,
148 bool dotsys);
149
150 void scaleFromAsciiTable(casa::Table& in, const std::string& filename,
151 const std::string& method,
152 const casa::Vector<casa::Float>& xout,
153 bool dotsys);
154
155 void scaleFromTable(casa::Table& in, const casa::Table& table,
156 const std::string& method,
157 const casa::Vector<casa::Float>& xout, bool dotsys);
158
159 void convertBrightnessUnits(casa::CountedPtr<Scantable>& in,
160 bool tokelvin, float cfac);
161
162 casa::CountedPtr< Scantable >
163 getScantable(const casa::CountedPtr< Scantable >& in, bool droprows);
164
165 casa::MaskedArray<casa::Float>
166 maskedArray( const casa::Vector<casa::Float>& s,
167 const casa::Vector<casa::uChar>& f );
168 casa::Vector<casa::uChar>
169 flagsFromMA(const casa::MaskedArray<casa::Float>& ma);
170
171 bool insitu_;
172};
173
174}
175#endif
Note: See TracBrowser for help on using the repository browser.