filters.h

Go to the documentation of this file.

/****************************************************************************************
 * SCYLLARUS : C++ HYPERSPECTRAL PROCESSING LIBRARY
 * filters.h - Image processing routines for use with HS images.
 *
 * Copyright 2013 NICTA (Jeremy Oorloff, Keith Grochow)
 *
 ***************************************************************************************/

#ifndef FILTERS_H_
#define FILTERS_H_


#include <armadillo>
#include <iso646.h>
#include "logger.h"

namespace scyl
{
    template<class T> inline void clip_lower(T &a, float f_lim, float f_set)
    {
        for (arma::uword i = 0; i < a.n_elem; i++)
        {
            if (a[i] <= f_lim)
            {
                a[i] = f_set;
            }
        }
    }

    template<class T> inline void clip_lower(T& a, float f_lim)
    {
        clip_lower(a, f_lim, f_lim);
    }

    arma::fmat conv(const arma::fmat & a,
                          arma::fvec   kernel);

    arma::fmat conv2(const arma::fmat & a,
                           arma::fmat   kernel);

    arma::fmat wiener2(const arma::fmat & a,
                             arma::fvec   nhood);

    arma::fmat median(const arma::fmat & a,
                             arma::u32   iSize);

    void savitzky_golay(arma::fcube & in,
                                int   window=5);

    void moving_average(arma::fcube & in,
                                int   window=3);
}

#endif /* FILTERS_H_ */