source: trunk/src/STMath.h @ 1068

Last change on this file since 1068 was 1066, checked in by mar637, 18 years ago

Enhancement Ticket #50; arbitrary quotients.

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 5.4 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/**
32Mathmatical operations on Scantable objects
[2]33
[805]34@author Malte Marquarding
35*/
[890]36class STMath : private Logger {
[716]37public:
[171]38
[805]39  typedef casa::InterpolateArray1D<casa::Double,
40                                   casa::Float>::InterpolationMethod imethod;
[299]41
[805]42  typedef std::map<std::string, imethod> imap;
[177]43
[805]44  STMath(bool insitu=true);
[221]45
[995]46  virtual ~STMath();
[227]47
[805]48  /**
49   * set the @attr insitu attribute
50   * @param b
51   */
52  bool insitu() const { return insitu_;};
53  void setInsitu(bool b) { insitu_ = b; };
[262]54
[805]55  casa::CountedPtr<Scantable>
56    average( const std::vector<casa::CountedPtr<Scantable> >& in,
[858]57             const std::vector<bool>& mask = std::vector<bool>(),
58             const std::string& weight = "NONE",
[977]59             const std::string& avmode = "SCAN");
[1066]60
[940]61  casa::CountedPtr< Scantable >
62    averagePolarisations( const casa::CountedPtr< Scantable > & in,
63                                const std::vector<bool>& mask,
64                                const std::string& weight );
[234]65
[805]66  casa::CountedPtr<Scantable>
67    unaryOperate( const casa::CountedPtr<Scantable>& in, float val,
68                  const std::string& mode, bool tsys=false );
[169]69
[1066]70  casa::CountedPtr<Scantable> autoQuotient(const casa::CountedPtr<Scantable>& in,
71                                           const std::string& mode = "NEAREST",
72                                           bool preserve = true);
73
74  casa::CountedPtr<Scantable> quotient( const casa::CountedPtr<Scantable>& on,
75                                        const casa::CountedPtr<Scantable>& off,
[805]76                                        bool preserve = true );
[169]77
[805]78  casa::CountedPtr<Scantable>
79    freqSwitch( const casa::CountedPtr<Scantable>& in );
[716]80
[805]81  std::vector<float> statistic(const casa::CountedPtr<Scantable>& in,
82                               const std::vector<bool>& mask,
83                               const std::string& which);
[457]84
[805]85  casa::CountedPtr<Scantable> bin( const casa::CountedPtr<Scantable>& in,
86                                   int width=5);
87  casa::CountedPtr<Scantable>
88    resample(const casa::CountedPtr<Scantable>& in,
89             const std::string& method, float width);
[457]90
[805]91  casa::CountedPtr<Scantable>
92    smooth(const casa::CountedPtr<Scantable>& in, const std::string& kernel,
[995]93                      float width);
[169]94
[805]95  casa::CountedPtr<Scantable>
96    gainElevation(const casa::CountedPtr<Scantable>& in,
[867]97                  const std::vector<float>& coeff,
[805]98                  const std::string& fileName,
[995]99                  const std::string& method);
[805]100  casa::CountedPtr<Scantable>
101    convertFlux(const casa::CountedPtr<Scantable>& in, float d,
102                float etaap, float jyperk);
[152]103
[805]104  casa::CountedPtr<Scantable> opacity(const casa::CountedPtr<Scantable>& in,
105                                      float tau);
[144]106
[841]107  casa::CountedPtr<Scantable>
108    merge(const std::vector<casa::CountedPtr<Scantable> >& in);
109
[896]110  casa::CountedPtr<Scantable>
111    invertPhase( const casa::CountedPtr<Scantable>& in);
[912]112
[896]113  casa::CountedPtr<Scantable>
114    rotateXYPhase( const casa::CountedPtr<Scantable>& in, float phase);
[912]115
[896]116  casa::CountedPtr<Scantable>
117    rotateLinPolPhase( const casa::CountedPtr<Scantable>& in, float phase);
118
[805]119  /// @todo frequency alignment
[152]120
[896]121  casa::CountedPtr<Scantable>
122    swapPolarisations(const casa::CountedPtr<Scantable>& in);
123
[917]124  casa::CountedPtr<Scantable>
125    frequencyAlign( const casa::CountedPtr<Scantable>& in,
[927]126                    const std::string& refTime = "",
127                    const std::string& method = "cubic" );
[917]128
[992]129  casa::CountedPtr<Scantable>
130    convertPolarisation( const casa::CountedPtr<Scantable>& in,
131                         const std::string& newtype);
132
[805]133private:
[896]134  casa::CountedPtr<Scantable>  applyToPol( const casa::CountedPtr<Scantable>& in,
135                                           STPol::polOperation fptr,
136                                           casa::Float phase);
137
[805]138  static imethod stringToIMethod(const std::string& in);
139  static WeightType stringToWeight(const std::string& in);
[146]140
[805]141  void scaleByVector(casa::Table& in,
142                     const casa::Vector<casa::Float>& factor,
143                     bool dotsys);
[152]144
[805]145  void scaleFromAsciiTable(casa::Table& in, const std::string& filename,
146                           const std::string& method,
147                           const casa::Vector<casa::Float>& xout,
148                           bool dotsys);
[162]149
[805]150  void scaleFromTable(casa::Table& in, const casa::Table& table,
151                      const std::string& method,
152                      const casa::Vector<casa::Float>& xout, bool dotsys);
[227]153
[805]154  void convertBrightnessUnits(casa::CountedPtr<Scantable>& in,
155                              bool tokelvin, float cfac);
[227]156
[805]157  casa::CountedPtr< Scantable >
158    getScantable(const casa::CountedPtr< Scantable >& in, bool droprows);
[230]159
[805]160  casa::MaskedArray<casa::Float>
161    maskedArray( const casa::Vector<casa::Float>& s,
162                 const casa::Vector<casa::uChar>& f );
163  casa::Vector<casa::uChar>
164    flagsFromMA(const casa::MaskedArray<casa::Float>& ma);
[230]165
[805]166  bool insitu_;
167};
[234]168
[805]169}
[165]170#endif
Note: See TracBrowser for help on using the repository browser.