source: trunk/src/STMathWrapper.h@ 1317

Last change on this file since 1317 was 1308, checked in by mar637, 18 years ago

re-enabled the scantable binary operators

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 6.0 KB
Line 
1//
2// C++ Interface: STMathWrapper
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#ifndef ASAPSTMATHWRAPPER_H
13#define ASAPSTMATHWRAPPER_H
14
15#include <vector>
16#include <string>
17
18#include <casa/Utilities/CountedPtr.h>
19
20#include "STMath.h"
21#include "Scantable.h"
22#include "ScantableWrapper.h"
23
24namespace asap {
25
26/**
27Wrapper class to handle ScantableWrapper
28
29@author Malte Marquarding
30*/
31class STMathWrapper : public STMath {
32public:
33 STMathWrapper() {;}
34 STMathWrapper(bool insitu) : STMath(insitu) {;}
35
36 virtual ~STMathWrapper() {;}
37
38 ScantableWrapper
39 average( const std::vector<ScantableWrapper>& in,
40 const std::vector<bool>& mask,
41 const std::string& weight,
42 const std::string& avmode )
43 {
44 std::vector<casa::CountedPtr<Scantable> > sts;
45 for (unsigned int i=0; i<in.size(); ++i) sts.push_back(in[i].getCP());
46 return ScantableWrapper(STMath::average(sts, mask, weight, avmode));
47 }
48
49 ScantableWrapper
50 averageChannel( const ScantableWrapper& in,
51 const std::string& mode = "MEDIAN",
52 const std::string& avmode = "NONE")
53 {
54 return ScantableWrapper(STMath::averageChannel(in.getCP(), mode, avmode));
55 }
56
57 ScantableWrapper
58 averagePolarisations( const ScantableWrapper& in,
59 const std::vector<bool>& mask,
60 const std::string& weight)
61 { return ScantableWrapper(STMath::averagePolarisations(in.getCP(),mask, weight));}
62
63 ScantableWrapper
64 averageBeams( const ScantableWrapper& in,
65 const std::vector<bool>& mask,
66 const std::string& weight)
67
68 { return ScantableWrapper(STMath::averageBeams(in.getCP(),mask, weight));}
69
70 ScantableWrapper
71 unaryOperate( const ScantableWrapper& in, float val,
72 const std::string& mode, bool tsys=false )
73 { return ScantableWrapper(STMath::unaryOperate(in.getCP(), val, mode, tsys)); }
74
75 ScantableWrapper binaryOperate( const ScantableWrapper& left,
76 const ScantableWrapper& right,
77 const std::string& mode)
78 { return ScantableWrapper( STMath::binaryOperate( left.getCP(), right.getCP(),
79 mode ) ); }
80
81
82 ScantableWrapper autoQuotient( const ScantableWrapper& in,
83 const std::string& mode = "NEAREST",
84 bool preserve = true )
85 { return ScantableWrapper(STMath::autoQuotient(in.getCP(), mode, preserve)); }
86
87 ScantableWrapper quotient( const ScantableWrapper& on,
88 const ScantableWrapper& off,
89 bool preserve = true )
90 { return ScantableWrapper( STMath::quotient( on.getCP(), off.getCP(),
91 preserve ) ); }
92
93 ScantableWrapper
94 freqSwitch( const ScantableWrapper& in )
95 { return ScantableWrapper(STMath::freqSwitch(in.getCP())); }
96
97 std::vector<float> statistic(const ScantableWrapper& in,
98 const std::vector<bool>& mask,
99 const std::string& which)
100 { return STMath::statistic(in.getCP(), mask, which); }
101
102 ScantableWrapper bin( const ScantableWrapper& in, int width=5)
103 { return ScantableWrapper(STMath::bin(in.getCP(), width)); }
104
105 ScantableWrapper
106 resample(const ScantableWrapper& in,
107 const std::string& method, float width)
108 { return ScantableWrapper(STMath::resample(in.getCP(), method, width)); }
109
110 ScantableWrapper
111 smooth(const ScantableWrapper& in, const std::string& kernel, float width)
112 { return ScantableWrapper(STMath::smooth(in.getCP(), kernel, width)); }
113
114 ScantableWrapper
115 gainElevation(const ScantableWrapper& in,
116 const std::vector<float>& coeff,
117 const std::string& filename,
118 const std::string& method)
119
120 { return
121 ScantableWrapper(STMath::gainElevation(in.getCP(), coeff, filename, method)); }
122
123 ScantableWrapper
124 convertFlux(const ScantableWrapper& in, float d,
125 float etaap, float jyperk)
126 { return ScantableWrapper(STMath::convertFlux(in.getCP(), d, etaap, jyperk)); }
127
128 ScantableWrapper opacity(const ScantableWrapper& in,
129 float tau)
130 { return ScantableWrapper(STMath::opacity(in.getCP(), tau)); }
131
132 ScantableWrapper
133 merge(const std::vector<ScantableWrapper >& in)
134
135 {
136 std::vector<casa::CountedPtr<Scantable> > sts;
137 for (unsigned int i=0; i<in.size(); ++i) sts.push_back(in[i].getCP());
138 return ScantableWrapper(STMath::merge(sts)); }
139
140 ScantableWrapper rotateXYPhase( const ScantableWrapper& in, float angle)
141 { return ScantableWrapper(STMath::rotateXYPhase(in.getCP(), angle)); }
142
143 ScantableWrapper rotateLinPolPhase( const ScantableWrapper& in, float angle)
144 { return ScantableWrapper(STMath::rotateLinPolPhase(in.getCP(), angle)); }
145
146 ScantableWrapper invertPhase( const ScantableWrapper& in )
147 { return ScantableWrapper(STMath::invertPhase(in.getCP())); }
148
149 ScantableWrapper swapPolarisations( const ScantableWrapper& in )
150 { return ScantableWrapper(STMath::swapPolarisations(in.getCP())); }
151
152 ScantableWrapper frequencyAlign( const ScantableWrapper& in,
153 const std::string& refTime,
154 const std::string& method )
155 { return ScantableWrapper(STMath::frequencyAlign(in.getCP())); }
156
157 ScantableWrapper convertPolarisation( const ScantableWrapper& in,
158 const std::string& newtype )
159 { return ScantableWrapper(STMath::convertPolarisation(in.getCP(),newtype)); }
160
161 ScantableWrapper mxExtract( const ScantableWrapper& in,
162 const std::string& scantype="on" )
163 { return ScantableWrapper(STMath::mxExtract(in.getCP(),scantype)); }
164
165 ScantableWrapper lagFlag( const ScantableWrapper& in,
166 double frequency, double width )
167 { return ScantableWrapper(STMath::lagFlag(in.getCP(), frequency, width)); }
168
169};
170
171}
172
173#endif
Note: See TracBrowser for help on using the repository browser.