source: trunk/src/STMathWrapper.h@ 1079

Last change on this file since 1079 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.1 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 unaryOperate( const ScantableWrapper& in, float val,
65 const std::string& mode, bool tsys=false )
66 { return ScantableWrapper(STMath::unaryOperate(in.getCP(), val, mode, tsys)); }
67
68 ScantableWrapper autoQuotient( const ScantableWrapper& in,
69 const std::string& mode = "NEAREST",
70 bool preserve = true )
71 { return ScantableWrapper(STMath::autoQuotient(in.getCP(), mode, preserve)); }
72
73 ScantableWrapper quotient( const ScantableWrapper& on,
74 const ScantableWrapper& off,
75 bool preserve = true )
76 { return ScantableWrapper( STMath::quotient( on.getCP(), off.getCP(),
77 preserve ) ); }
78
79 ScantableWrapper
80 freqSwitch( const ScantableWrapper& in )
81 { return ScantableWrapper(STMath::freqSwitch(in.getCP())); }
82
83 std::vector<float> statistic(const ScantableWrapper& in,
84 const std::vector<bool>& mask,
85 const std::string& which)
86 { return STMath::statistic(in.getCP(), mask, which); }
87
88 ScantableWrapper bin( const ScantableWrapper& in, int width=5)
89 { return ScantableWrapper(STMath::bin(in.getCP(), width)); }
90
91 ScantableWrapper
92 resample(const ScantableWrapper& in,
93 const std::string& method, float width)
94 { return ScantableWrapper(STMath::resample(in.getCP(), method, width)); }
95
96 ScantableWrapper
97 smooth(const ScantableWrapper& in, const std::string& kernel, float width)
98 { return ScantableWrapper(STMath::smooth(in.getCP(), kernel, width)); }
99
100 ScantableWrapper
101 gainElevation(const ScantableWrapper& in,
102 const std::vector<float>& coeff,
103 const std::string& filename,
104 const std::string& method)
105
106 { return
107 ScantableWrapper(STMath::gainElevation(in.getCP(), coeff, filename, method)); }
108
109 ScantableWrapper
110 convertFlux(const ScantableWrapper& in, float d,
111 float etaap, float jyperk)
112 { return ScantableWrapper(STMath::convertFlux(in.getCP(), d, etaap, jyperk)); }
113
114 ScantableWrapper opacity(const ScantableWrapper& in,
115 float tau)
116 { return ScantableWrapper(STMath::opacity(in.getCP(), tau)); }
117
118 ScantableWrapper
119 merge(const std::vector<ScantableWrapper >& in)
120
121 {
122 std::vector<casa::CountedPtr<Scantable> > sts;
123 for (unsigned int i=0; i<in.size(); ++i) sts.push_back(in[i].getCP());
124 return ScantableWrapper(STMath::merge(sts)); }
125
126 ScantableWrapper rotateXYPhase( const ScantableWrapper& in, float angle)
127 { return ScantableWrapper(STMath::rotateXYPhase(in.getCP(), angle)); }
128
129 ScantableWrapper rotateLinPolPhase( const ScantableWrapper& in, float angle)
130 { return ScantableWrapper(STMath::rotateLinPolPhase(in.getCP(), angle)); }
131
132 ScantableWrapper invertPhase( const ScantableWrapper& in )
133 { return ScantableWrapper(STMath::invertPhase(in.getCP())); }
134
135 ScantableWrapper swapPolarisations( const ScantableWrapper& in )
136 { return ScantableWrapper(STMath::swapPolarisations(in.getCP())); }
137
138 ScantableWrapper frequencyAlign( const ScantableWrapper& in,
139 const std::string& refTime,
140 const std::string& method )
141 { return ScantableWrapper(STMath::frequencyAlign(in.getCP())); }
142
143 ScantableWrapper convertPolarisation( const ScantableWrapper& in,
144 const std::string& newtype )
145 { return ScantableWrapper(STMath::convertPolarisation(in.getCP(),newtype)); }
146
147};
148
149}
150
151#endif
Note: See TracBrowser for help on using the repository browser.