VecFilter.cpp

 
/* SPDX-License-Identifier: CC0-1.0 */
 
#include <float.h>
#include <string>
#include <vector>
#include "VecFilter.h"
#include "IException.h"
#include "IString.h"
 
using namespace std;
const int MARKER = -999999;
namespace Isis {
  VecFilter::VecFilter() {}
 
  VecFilter::~VecFilter() {};
 
  vector<double> VecFilter::LowPass(vector<double> invec, int boxsize) {
    vector<double> outvec;
 
    // Clear the output vector
    //outvec.resize(0);
 
    // Boxcar size must be odd and greater than 1
    if((boxsize % 2) == 0) {
      string m = "Boxcar size must be odd and greater than 1 in [VecFilter::LowPass]";
      throw IException(IException::Programmer, m, _FILEINFO_);
    }
 
    // Perform lowpass filter
    int halfwidth = boxsize / 2;
    int vecsize = (int)invec.size() - 1;
    for(int i = 0; i <= vecsize; i++) {
      int i1 = i - halfwidth;
      if(i1 < 0) i1 = 0;
      int i2 = i + halfwidth;
      if(i2 > vecsize) i2 = vecsize;
      int npts = 0;
      double sum = 0.0;
      for(int j = i1; j <= i2; j++) {
        if(invec[j] != 0.0) {
          sum = sum + invec[j];
          npts++;
        }
      }
      if(npts > 0) {
        outvec.push_back(sum / npts);
      }
      else {
        outvec.push_back(0.0);
      }
    }
    return outvec;
  }
 
  vector<double> VecFilter::HighPass(vector<double> invec1, vector<double> invec2) {
    vector<double> outvec;
 
    // Both vectors must be the same size
    if(invec1.size() != invec2.size()) {
      string m = "Both vectors must be the same size in [VecFilter::HighPass]";
      throw IException(IException::Programmer, m, _FILEINFO_);
    }
 
    int vecsize = (int)invec1.size() - 1;
    for(int i = 0; i <= vecsize; i++) {
      if(invec1[i] != 0.0 && invec2[i] != 0.0) {
        outvec.push_back(invec1[i] - invec2[i]);
      }
      else {
        outvec.push_back(0.0);
      }
    }
    return outvec;
  }
 
  vector<double> VecFilter::HighPass(vector<double> pdInVector1, vector<double> pdInVector2,
                 vector<int> piValidPntsVector, int piMaxPoints, const QString & psMode)
  {
    vector<double> dOutVector;
 
    // Both vectors must be the same size
    if(pdInVector1.size() != pdInVector2.size()) {
      string m = "Both vectors must be the same size in [VecFilter::HighPass]";
      throw IException(IException::Programmer, m, _FILEINFO_);
    }
 
    int iSize = (int)pdInVector1.size();
    for(int i = 0; i < iSize; i++) {
      if(pdInVector1[i] != 0.0 && pdInVector2[i] != 0.0 && piValidPntsVector[i]==piMaxPoints) {
        if(psMode == "SUBTRACT") {
          dOutVector.push_back(pdInVector1[i] - pdInVector2[i]);
        }
        else {
          dOutVector.push_back(pdInVector1[i] / pdInVector2[i]);
        }
      }
      else {
        dOutVector.push_back(MARKER);
      }
    }
    return dOutVector;
  }
 
} // end namespace isis