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