source: branches/Release2.1.2/src/STMath.h@ 2357

Last change on this file since 2357 was 1310, checked in by mar637, 18 years ago

merge from trunk, to get binary operator changes

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 7.8 KB
RevLine 
[805]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
[2]14
[38]15#include <string>
[805]16#include <map>
17
[130]18#include <casa/aips.h>
[81]19#include <casa/Utilities/CountedPtr.h>
[805]20#include <casa/BasicSL/String.h>
21#include <casa/Arrays/Vector.h>
22#include <scimath/Mathematics/InterpolateArray1D.h>
[294]23
[805]24#include "Scantable.h"
[834]25#include "STDefs.h"
[896]26#include "STPol.h"
[894]27#include "Logger.h"
[2]28
[83]29namespace asap {
[2]30
[805]31/**
[1106]32 * Mathmatical operations on Scantable objects
33 * @author Malte Marquarding
[805]34*/
[890]35class STMath : private Logger {
[716]36public:
[1106]37 // typedef for long method name
[805]38 typedef casa::InterpolateArray1D<casa::Double,
39 casa::Float>::InterpolationMethod imethod;
[299]40
[1106]41 // typedef for std::map
[805]42 typedef std::map<std::string, imethod> imap;
[177]43
[1106]44/**
45 * whether to operate on the given Scantable or return a new one
46 * @param insitu the toggle for this behaviour
47 */
[805]48 STMath(bool insitu=true);
[221]49
[995]50 virtual ~STMath();
[227]51
[805]52 /**
[1295]53 * get the currnt @attr inistu_ state
[805]54 */
55 bool insitu() const { return insitu_;};
[1143]56
[1106]57 /**
58 * set the currnt @attr inistu state
59 * @param b the new state
60 */
[805]61 void setInsitu(bool b) { insitu_ = b; };
[262]62
[1106]63
[1140]64 /**
65 * average a vector of Scantables
66 * @param in the vector of Scantables to average
[1295]67 * @param mask an optional mask to apply on specific weights
[1140]68 * @param weight weighting scheme
69 * @param avmode the mode ov averaging. Per "SCAN" or "ALL".
70 * @return a casa::CountedPtr<Scantable> which either holds a new Scantable
71 * or returns the imput pointer.
72 */
[805]73 casa::CountedPtr<Scantable>
74 average( const std::vector<casa::CountedPtr<Scantable> >& in,
[858]75 const std::vector<bool>& mask = std::vector<bool>(),
76 const std::string& weight = "NONE",
[977]77 const std::string& avmode = "SCAN");
[1066]78
[1140]79 /**
80 * median average a vector of Scantables. See also STMath::average
81 * @param in the Scantable to average
82 * @param mode the averaging mode. Currently only "MEDIAN"
83 * @param avmode the mode ov averaging. Per "SCAN" or "ALL".
84 * @return a casa::CountedPtr<Scantable> which either holds a new Scantable
85 * or returns the imput pointer.
86 */
[1069]87 casa::CountedPtr<Scantable>
88 averageChannel( const casa::CountedPtr<Scantable> & in,
[1078]89 const std::string& mode = "MEDIAN",
90 const std::string& avmode = "SCAN");
[1069]91
[1140]92 /**
[1145]93 * Average polarisations together. really only useful if only linears are
94 * available.
[1140]95 * @param in the input Scantable
96 * @param mask an optional mask if weight allows one
97 * @param weight weighting scheme
[1143]98 * @return
[1140]99 */
[940]100 casa::CountedPtr< Scantable >
101 averagePolarisations( const casa::CountedPtr< Scantable > & in,
[1140]102 const std::vector<bool>& mask,
103 const std::string& weight );
[234]104
[1145]105 /**
106 * Average beams together.
107 * @param in the input Scantable
108 * @param mask an optional mask if weight allows one
109 * @param weight weighting scheme
110 * @return
111 */
112 casa::CountedPtr< Scantable >
113 averageBeams( const casa::CountedPtr< Scantable > & in,
114 const std::vector<bool>& mask,
115 const std::string& weight );
116
[805]117 casa::CountedPtr<Scantable>
118 unaryOperate( const casa::CountedPtr<Scantable>& in, float val,
119 const std::string& mode, bool tsys=false );
[169]120
[1310]121 casa::CountedPtr<Scantable>
122 binaryOperate( const casa::CountedPtr<Scantable>& left,
123 const casa::CountedPtr<Scantable>& right,
124 const std::string& mode);
125
[1066]126 casa::CountedPtr<Scantable> autoQuotient(const casa::CountedPtr<Scantable>& in,
127 const std::string& mode = "NEAREST",
128 bool preserve = true);
129
130 casa::CountedPtr<Scantable> quotient( const casa::CountedPtr<Scantable>& on,
131 const casa::CountedPtr<Scantable>& off,
[805]132 bool preserve = true );
[169]133
[805]134 casa::CountedPtr<Scantable>
135 freqSwitch( const casa::CountedPtr<Scantable>& in );
[716]136
[805]137 std::vector<float> statistic(const casa::CountedPtr<Scantable>& in,
138 const std::vector<bool>& mask,
139 const std::string& which);
[457]140
[805]141 casa::CountedPtr<Scantable> bin( const casa::CountedPtr<Scantable>& in,
142 int width=5);
143 casa::CountedPtr<Scantable>
144 resample(const casa::CountedPtr<Scantable>& in,
145 const std::string& method, float width);
[457]146
[805]147 casa::CountedPtr<Scantable>
148 smooth(const casa::CountedPtr<Scantable>& in, const std::string& kernel,
[995]149 float width);
[169]150
[805]151 casa::CountedPtr<Scantable>
152 gainElevation(const casa::CountedPtr<Scantable>& in,
[867]153 const std::vector<float>& coeff,
[805]154 const std::string& fileName,
[995]155 const std::string& method);
[805]156 casa::CountedPtr<Scantable>
157 convertFlux(const casa::CountedPtr<Scantable>& in, float d,
158 float etaap, float jyperk);
[152]159
[805]160 casa::CountedPtr<Scantable> opacity(const casa::CountedPtr<Scantable>& in,
161 float tau);
[144]162
[841]163 casa::CountedPtr<Scantable>
164 merge(const std::vector<casa::CountedPtr<Scantable> >& in);
165
[896]166 casa::CountedPtr<Scantable>
167 invertPhase( const casa::CountedPtr<Scantable>& in);
[912]168
[896]169 casa::CountedPtr<Scantable>
170 rotateXYPhase( const casa::CountedPtr<Scantable>& in, float phase);
[912]171
[896]172 casa::CountedPtr<Scantable>
173 rotateLinPolPhase( const casa::CountedPtr<Scantable>& in, float phase);
174
175 casa::CountedPtr<Scantable>
176 swapPolarisations(const casa::CountedPtr<Scantable>& in);
177
[917]178 casa::CountedPtr<Scantable>
179 frequencyAlign( const casa::CountedPtr<Scantable>& in,
[927]180 const std::string& refTime = "",
181 const std::string& method = "cubic" );
[917]182
[992]183 casa::CountedPtr<Scantable>
184 convertPolarisation( const casa::CountedPtr<Scantable>& in,
185 const std::string& newtype);
[1143]186
[1140]187 casa::CountedPtr<Scantable>
[1143]188 mxExtract( const casa::CountedPtr<Scantable>& in,
189 const std::string& srctype = "on");
[992]190
[1192]191 /**
192 * "hard" flag the data, this flags everything selected in setSelection()
193 * @param frequency the frequency to remove
194 * @param width the number of lags to flag left to the side of the frequency
195 */
196 casa::CountedPtr<Scantable>
197 lagFlag( const casa::CountedPtr<Scantable>& in, double frequency,
[1200]198 double width);
[1192]199
[805]200private:
[896]201 casa::CountedPtr<Scantable> applyToPol( const casa::CountedPtr<Scantable>& in,
202 STPol::polOperation fptr,
203 casa::Float phase);
204
[805]205 static imethod stringToIMethod(const std::string& in);
206 static WeightType stringToWeight(const std::string& in);
[146]207
[805]208 void scaleByVector(casa::Table& in,
209 const casa::Vector<casa::Float>& factor,
210 bool dotsys);
[152]211
[805]212 void scaleFromAsciiTable(casa::Table& in, const std::string& filename,
213 const std::string& method,
214 const casa::Vector<casa::Float>& xout,
215 bool dotsys);
[162]216
[805]217 void scaleFromTable(casa::Table& in, const casa::Table& table,
218 const std::string& method,
219 const casa::Vector<casa::Float>& xout, bool dotsys);
[227]220
[805]221 void convertBrightnessUnits(casa::CountedPtr<Scantable>& in,
222 bool tokelvin, float cfac);
[227]223
[805]224 casa::CountedPtr< Scantable >
225 getScantable(const casa::CountedPtr< Scantable >& in, bool droprows);
[230]226
[805]227 casa::MaskedArray<casa::Float>
228 maskedArray( const casa::Vector<casa::Float>& s,
229 const casa::Vector<casa::uChar>& f );
230 casa::Vector<casa::uChar>
231 flagsFromMA(const casa::MaskedArray<casa::Float>& ma);
[230]232
[805]233 bool insitu_;
234};
[234]235
[805]236}
[165]237#endif
Note: See TracBrowser for help on using the repository browser.