kriging.h

Go to the documentation of this file.

/*
 * $Id: kriging.h 17611 2014-09-26 20:39:48Z predoehl $
 *
 * Tab size:  4
 */

#ifndef KRIGING_H_INCLUDED_PREDOEHL_UOFARIZONAVISION
#define KRIGING_H_INCLUDED_PREDOEHL_UOFARIZONAVISION 1

#include <l_mt_cpp/l_mt_mutexlock.h>
#include <m_cpp/m_vector.h>
#include <topo_cpp/nedgrid.h>
#include <vector>

namespace kjb
{

class Kriging_interpolator
{
    const double m_char_sf2;  
    char m_utm_zone;          
    std::vector< TopoFusion::pt > m_training;
    Vector m_cov_inv_targ;


    mutable bool kcache_valid;

    mutable TopoFusion::pt m_query_zbad;

    mutable TopoFusion::pt m_query_zgood;

    mutable Vector k_cache;


    void refresh_cache(const TopoFusion::pt&) const;

public:
    Kriging_interpolator(
        double,
        const Ned13_one_degree_grid::IntegralLL&,
        Ned13_grid_cache*,
        Pthread_mutex* = 0
    );

    double operator()(const TopoFusion::pt& utm) const
    {
        refresh_cache(utm);
        return dot( m_cov_inv_targ, k_cache );
    }

    double diff_e(const TopoFusion::pt&) const;

    double diff_n(const TopoFusion::pt&) const;
};


}

#endif /* KRIGING_H_INCLUDED_PREDOEHL_UOFARIZONAVISION */