HFRCMSpaceInterpolator.h

Go to the documentation of this file.

#ifndef HFRCMSPACEINTERPOLATOR_H_
#define HFRCMSPACEINTERPOLATOR_H_

#include "component/SyncComponent.h"
#include "data/Mtx.h"
#include "HFRadarCompactModelForecaster.h" // for MAX_EOFS, etc.

class UniversalDataReader;

class HFRCMSpaceInterpolator : public SyncEstimationComponent
{
public:
    HFRCMSpaceInterpolator( const Module* module );
    virtual ~HFRCMSpaceInterpolator();

    void initialize( void );

    void run();

    static HFRCMSpaceInterpolator* Instance()
    {
        return Instance_;
    }

    bool lookupEmpiricalOrthogonalFunctions( Mtx& eastEOFs, Mtx& northEOFs, float& eastEM,  float& northEM, float& eastSF,  float& northSF, const float latDeg, const float lonDeg ); // TODO: Is there a reason that I should not make this a static function so that user components to not need to store the pointer to the singleton in order to call the function?

protected:
    int verbosity_;
    const Mtx eofs_;
    const Mtx ensMean_;
    const Mtx gridIdxRev_;
    const Mtx gridLon_;
    const Mtx gridLat_;
    const Mtx scalingFactors_;
    Mtx gridIdx_;
    const int numLocations_, numModes_;

    // readers & writers
    UniversalDataReader* latitudeReader_;
    UniversalDataReader* longitudeReader_;
    ConfigReader* verbosityLevelConfigReader_;

    bool getLowerLeftCorner( int &latitude_index, int &longitude_index, const float latitude_degrees, const float longitude_degrees );

    float interpolateRaveledTwoDimensionalGrid( const Mtx z, const float latitude_degrees, const float longitude_degrees, const int latitude_index, const int longitude_index, const int starting_offset = 0, const int layer = 0 );

    static HFRCMSpaceInterpolator* Instance_;

private:
    // Note that the copy constructor below is private and not given a body.
    // Any attempt to call it will return a compiler error.
    HFRCMSpaceInterpolator( const HFRCMSpaceInterpolator& old ); // disallow copy constructor

};

#endif /* HFRCMSPACEINTERPOLATOR_H_ */