SplineFill.h

#ifndef SplineFill_h
#define SplineFill_h
 
/* SPDX-License-Identifier: CC0-1.0 */                                                                  
 
#include <string>
#include <vector>
 
#include "IString.h"
#include "Module.h"
#include "NumericalApproximation.h"
#include "SpecialPixel.h"
#include "IException.h"
 
namespace Isis {
 
  class SplineFill : public Module {
 
    public:
      //  Constructors and Destructor
      SplineFill() : Module("SplineFill"), _filled(0) { }
 
      SplineFill(const Module &c) :
                     Module("SplineFill", c), _filled(0) {
        fill(c.ref());
        _history.add(formHistory());
      }
 
      SplineFill(const HiVector &v) :
                     Module("SplineFill"), _filled(0) {
        fill(v);
        _history.add(formHistory());
      }
      SplineFill(const HiVector &v, const HiHistory &h) :
                     Module("SplineFill", h), _filled(0) {
        fill(v);
        _history.add(formHistory());
      }
 
      virtual ~SplineFill() { }
 
      void Process(const HiVector &v) {
        fill(v);
        _history.clear();
        _history.add(formHistory());
      }
 
      inline int Filled() const { return (_filled); }
 
    private:
      int   _filled;         
 
      QString formHistory() {
        QString cfilled(toString(_filled));
        return (QString("SplineFill(Cubic,Filled[" + cfilled + "])"));
      }
 
      void fill(const HiVector &v) {
        NumericalApproximation spline(NumericalApproximation::CubicNatural);
        for (int i = 0 ; i < v.dim() ; i++) {
          if (!IsSpecial(v[i])) {
            spline.AddData(i, v[i]);
          }
        }
 
        //  Compute the spline and fill missing data
        HiVector vout(v.dim());
        _filled = 0;
        for (int j = 0 ; j < v.dim() ; j++) {
          if (IsSpecial(v[j])) {
            vout[j] = spline.Evaluate(j,NumericalApproximation::NearestEndpoint);
            _filled++;
          }
          else {
            vout[j] = v[j];
          }
        }
 
        _data = vout;
        return;
      }
  };
 
}     // namespace Isis
#endif