GroundDistances.h

Go to the documentation of this file.

/******************************************************************************
 *
 *  Copyright (C) 2008 Max Hermann
 *
 *  This file is part of the "CEDRIC Event Display" application.
 *
 *****************************************************************************/

#ifndef REDUKT_GROUNDDISTANCES
#define REDUKT_GROUNDDISTANCES

#include "EarthMoversDistance.h" // for class GroundDistance
#include "Feature.h"             // for class AbstractFeature
#include <cmath>     // for fabs(), pow(), sqrt()
#include <algorithm> // for min()

template<int P>
struct GroundDistance_Lp: public GroundDistance<AbstractFeature>
{
        float distance( const AbstractFeature& f1, const AbstractFeature& f2 ) const
        {
                int d = std::min( f1.dim(), f2.dim() );         
                float dist = 0.f;
                
                for( int i=0; i < d; ++i )
                        dist += std::pow( fabs(f1[i]-f2[i]), P );
                
                return std::pow( dist, 1.0/P );
        }
};

template<>
float GroundDistance_Lp<1>::distance( const AbstractFeature& f1, const AbstractFeature& f2 ) const
{
        int d = std::min( f1.dim(), f2.dim() );         
        float dist = 0.f;

        for( int i=0; i < d; ++i )
                dist += std::fabs( f1[i] - f2[i] );
        
        return dist;
}

template<>
float GroundDistance_Lp<2>::distance( const AbstractFeature& f1, const AbstractFeature& f2 ) const
{
        int d = std::min( f1.dim(), f2.dim() );         
        float dist = 0.f;
        
        for( int i=0; i < d; ++i )
                dist += std::fabs( f1[i]*f1[i] - f2[i]*f2[i] );
        
        return std::sqrt(dist);
}

#endif // REDUKT_GROUNDDISTANCES