Changeset 2012 for trunk/src/Scantable.h

Timestamp:

02/25/11 16:51:50 (13 years ago)

Author:

WataruKawasaki

Message:

New Development: Yes

JIRA Issue: Yes (CAS-2373, CAS-2620)

Ready for Test: Yes

Interface Changes: Yes

What Interface Changed: For Scantable::polyBaseline(), parameters and return value have been changed.

Test Programs:

Put in Release Notes: Yes

Module(s): sdbaseline, sd.linefinder

Description: (1) CAS-2373-related:

(1.1) Modified Scantable::polyBaseline() to have the row-based loop inside.

Now it fits and subtracts baseline for all rows and also output info
about the fitting result to logger and text file, while in the
previous version this method just did baseline fitting/subtraction
for one row only and had to be put inside a row-based loop at the
python side ("poly_baseline()" in scantable.py) and result output had
also to be controlled at the python side. Using a test data from NRO
45m telescope (348,000 rows, 512 channels), the processing time of
scantable.poly_baseline() has reduced from 130 minutes to 5-6 minutes.

(1.2) For accelerating the another polynomial fitting function, namely

scantable.auto_poly_baseline(), added a method
Scantable::autoPolyBaseline(). This basically does the same thing
with Scantable::polyBaseline(), but uses linefinfer also to
automatically flag the line regions.

(1.3) To make linefinder usable in Scantable class, added a method

linefinder.setdata(). This makes it possible to apply linefinder
for a float-list data given as a parameter, without setting scantable,
while it was indispensable to set scantable to use linefinger previously.

(2) CAS-2620-related:

Added Scantable::cubicSplineBaseline() and autoCubicSplineBaseline() for
fit baseline using the cubic spline function. Parameters include npiece
(number of polynomial pieces), clipthresh (clipping threshold), and
clipniter (maximum iteration number).

File:

• trunk/src/Scantable.h (modified) (3 diffs)

trunk/src/Scantable.h

@@ -16 +16 @@
 #include <string>
 #include <vector>
+#include <iostream>
+#include <fstream>
 // AIPS++
 #include <casa/aips.h>
@@ -493 +495 @@
   bool getFlagtraFast(int whichrow);
 
-  void polyBaselineBatch(const std::vector<bool>& mask, int order);
-  STFitEntry polyBaseline(const std::vector<bool>& mask, int order, int rowno);
+  void polyBaseline(const std::vector<bool>& mask, 
+                    int order, 
+                    bool outLogger=false, 
+                    const std::string& blfile="");
+  void autoPolyBaseline(const std::vector<bool>& mask,
+                        int order, 
+                        const std::vector<int>& edge, 
+                        float threshold=3.0, 
+                        int chanAvgLimit=1, 
+                        bool outLogger=false, 
+                        const std::string& blfile="");
+  void cubicSplineBaseline(const std::vector<bool>& mask,
+                           int nPiece,
+                           float thresClip, 
+                           int nIterClip,
+                           bool outLogger=false,
+                           const std::string& blfile="");
+  void autoCubicSplineBaseline(const std::vector<bool>& mask,
+                               int nPiece,
+                               float thresClip, 
+                               int nIterClip,
+                               const std::vector<int>& edge, 
+                               float threshold=3.0, 
+                               int chanAvgLimit=1, 
+                               bool outLogger=false, 
+                               const std::string& blfile="");
+  float getRms(const std::vector<bool>& mask, int whichrow);
+  std::string formatBaselineParams(const std::vector<float>& params, 
+                                   const std::vector<bool>& fixed, 
+                                   float rms, 
+                                   const std::string& masklist, 
+                                   int whichrow, 
+                                   bool verbose=false) const;
+  std::string formatPiecewiseBaselineParams(const std::vector<int>& ranges, 
+                                            const std::vector<float>& params, 
+                                            const std::vector<bool>& fixed, 
+                                            float rms, 
+                                            const std::string& masklist, 
+                                            int whichrow, 
+                                            bool verbose=false) const;
+
 
 private:
@@ -608 +649 @@
                                                       const casa::Array<T2>&);
 
-  void doPolyBaseline(const std::vector<bool>& mask, int order, int rowno, Fitter& fitter);
+  void fitBaseline(const std::vector<bool>& mask, int whichrow, Fitter& fitter);
+  std::vector<float> doCubicSplineFitting(const std::vector<float>& data, 
+                                          const std::vector<bool>& mask,
+                                          std::vector<int>& sectionRanges,
+                                          std::vector<float>& params,
+                                          int nPiece=2,
+                                          float thresClip=3.0, 
+                                          int nIterClip=1);
+  bool hasSameNchanOverIFs();
+  std::string getMaskRangeList(const std::vector<bool>& mask, 
+                                int whichrow, 
+                                const casa::String& coordInfo, 
+                                bool hasSameNchan, 
+                                int firstIF, 
+                                bool silent=false);
+  std::vector<int> getMaskEdgeIndices(const std::vector<bool>& mask, bool getStartIndices=true);
+  std::string formatBaselineParamsHeader(int whichrow, const std::string& masklist, bool verbose) const;
+  std::string formatBaselineParamsFooter(float rms, bool verbose) const;
+  std::vector<bool> getCompositeChanMask(int whichrow, const std::vector<bool>& inMask);
+  //std::vector<bool> getCompositeChanMask(int whichrow, const std::vector<bool>& inMask, const std::vector<int>& edge, const int minEdgeSize, STLineFinder& lineFinder);
 
   void applyChanFlag( casa::uInt whichrow, const std::vector<bool>& msk, casa::uChar flagval);