Changeset 2186 for trunk/src/MathUtils.cpp

Timestamp:

06/07/11 23:49:53 (13 years ago)

Author:

WataruKawasaki

Message:

New Development: Yes

JIRA Issue: Yes CAS-3149

Ready for Test: Yes

Interface Changes: Yes

What Interface Changed: scantable.*sinusoid_baseline() params

Test Programs:

Put in Release Notes: Yes

Module(s):

Description: (1) Implemented an automated sinusoidal fitting functionality

(2) FFT available with scantable.fft()
(3) fixed a bug of parsing 'edge' param used by linefinder.
(4) a function to show progress status for row-based iterations.

File:

• trunk/src/MathUtils.cpp (modified) (1 diff)

trunk/src/MathUtils.cpp

@@ -248 +248 @@
   }
 }
+
+void mathutil::doZeroOrderInterpolation(casa::Vector<casa::Float>& data, 
+                                        std::vector<bool>& mask) {
+  int fstart = -1;
+  int fend = -1;
+  for (uInt i = 0; i < mask.size(); ++i) {
+    if (!mask[i]) {
+      fstart = i;
+      while (!mask[i] && i < mask.size()) {
+        fend = i;
+        i++;
+      }
+    }
+
+    // execute interpolation as the following criteria:
+    // (1) for a masked region inside the spectrum, replace the spectral 
+    //     values with the mean of those at the two channels just outside 
+    //     the both edges of the masked region. 
+    // (2) for a masked region at the spectral edge, replace the values 
+    //     with the one at the nearest non-masked channel. 
+    //     (ZOH, but bilateral)
+    Float interp = 0.0;
+    if (fstart-1 > 0) {
+      interp = data[fstart-1];
+      if (fend+1 < Int(data.nelements())) {
+        interp = (interp + data[fend+1]) / 2.0;
+      }
+    } else {
+      interp = data[fend+1];
+    }
+    if (fstart > -1 && fend > -1) {
+      for (int j = fstart; j <= fend; ++j) {
+        data[j] = interp;
+      }
+    }
+
+    fstart = -1;
+    fend = -1;
+  }
+}