casa6

#ifndef POINTINGINTERPOLATIONTVI_H_

#define POINTINGINTERPOLATIONTVI_H_

​

​

​

#include <casacore/casa/aips.h>

​

#include <casacore/casa/Containers/Record.h>

​

#include <msvis/MSVis/VisibilityIterator2.h>

#include <msvis/MSVis/ViImplementation2.h>

#include <msvis/MSVis/TransformingVi2.h>

​

#include <casacore/measures/Measures/MDirection.h>

​

​

namespace casa {

​

namespace vi {

​

using casacore::Vector;

using casacore::Double;

​

​

​

class PointingInterpolationTVI : public TransformingVi2 {

​

public:

    class Interpolator {

    public:

        using PointingTimes = Vector<Double>;      // pointingsCount[antId]

        using PointingDir   = Vector<Double>;      // 2

        using PointingDirs  = Vector<PointingDir>; // pointingsCount[antId]

        Interpolator();

​

        enum class InterpMethod {

            NEAREST,

            CUBIC_SPLINE

        };

​

        void setInterpMethod(InterpMethod method);

        InterpMethod getInterpMethod() const;

​

        void setData(const Vector<PointingTimes> &antsTimes, // nAnt

                     const Vector<PointingDirs> &antsDirs,   // nAnt

                     const Vector<bool> &antSelected);       // nAnt

​

        // Interpolation

        Vector<Double> pointingDir(int ant,double time) const;

        double nearestPointingTimeStamp() const;

    private:

        void clear();

        void init(const Vector<bool> &antSelected = Vector<bool>());

        using SplineSegmentCoeffs = Vector<Double>;           // 4

        using DirSegmentCoeffs = Vector<SplineSegmentCoeffs>; // 2

        using SplineCoeffs = Vector<DirSegmentCoeffs>;        // pointingsCount[antId] - 1

        void computeSplineCoeffs(const PointingTimes& timeStamps,

                             const PointingDirs& dirs,

                             SplineCoeffs& coeffs);

        Vector<PointingTimes> antsTimes_;        // nAnt

        Vector<bool> isSelected_;                // nAnt

        Vector<bool> isInterpolated_;            // nAnt

        Vector<SplineCoeffs> antsSplinesCoeffs_; // nAnt