emor.h

Go to the documentation of this file.

// -*- c-basic-offset: 4 -*-
#ifndef VIGRA_EXT_EMOR_H
#define VIGRA_EXT_EMOR_H
 
#include <assert.h>
 
#include <vector>
 
#include <vigra/numerictraits.hxx>
#include <vigra_ext/utils.h>
 
//#define DEBUG_WRITE_FILES
 
namespace vigra_ext{
 
namespace EMoR
{
    extern IMPEX double f0[1024];
    extern IMPEX double h[25][1024];
 
    template <class VECTOR>
    inline void createEMoRLUT(const std::vector<float> & params, VECTOR & lut)
    {
        typedef typename VECTOR::value_type VT;
 
        VT s = (VT) vigra_ext::LUTTraits<VT>::max();
 
// lookup tables
        size_t nDim = params.size();
        assert(nDim < 26);
        lut.resize(1024);
        for (int i=0; i<1024; ++i) {
            double t = vigra_ext::EMoR::f0[i];
            for (size_t j=0; j < nDim; j++) {
                t += params[j] * vigra_ext::EMoR::h[j][i];
            }
            lut[i] = vigra::NumericTraits<VT>::fromRealPromote(t*s);
        }
    }
}
 
} // namespace
 
#endif // VIGRA_EXT_VIGNETTING_CORRECTION_H