|| || || || || || || ||

---

ColorGaussResponses.h

Go to the documentation of this file.

#ifndef Impala_Core_Array_ColorGaussResponses_h
#define Impala_Core_Array_ColorGaussResponses_h

#include <vector>
#include "Util/Math.h"
#include "Core/Array/ColorSegmentationInvariant.h"
#include "Core/Array/ColorSegmentationAlgorithm.h"
#include "Core/Array/ArrayPool.h"
#include "Core/Array/ColorSpace.h"
#include "Core/Array/ProjectRange.h"
#include "Core/Array/Set.h"
#include "Core/Array/RecGauss.h"
#include "Core/Array/Norm2.h"

namespace Impala
{
namespace Core
{
namespace Array
{


inline std::vector<Array2dScalarReal64*>
ColorGaussResponses(Array2dVec3Real64* im)
{
    const int nrScales = 4;
    double scales[nrScales];

    scales[0] = 1.0 ;
    scales[1] = 2.0 ;
    scales[2] = 3.5 ;
    scales[3] = 7.0 ;

    // convert to opponent color space
    ColorSpace(im, im, Element::RGB, Element::OOO);

    int width = im->CW();
    int height = im->CH();

    // declare and initialize temporary images
    std::vector<Array2dScalarReal64*> imList;
    int borderSize = 14; // make sure this is large enough!
    ArrayPool imOOO(width, height, 0, 0);
    imOOO.AllocSR64(3, borderSize, borderSize);
    imOOO.AllocC64(3, borderSize, borderSize);
    for (int ii=0 ; ii<3 ; ii++)
    {
        Array2dScalarReal64* tmpSR64 = imOOO.GetSR64(ii);
        ProjectRange(tmpSR64, im, ii+1);
        Array2dComplex64* tmpC64 = imOOO.GetC64(ii);
        Set(tmpC64, tmpSR64);
    }

    for (int si=0; si<nrScales-1; si++)
    {
        for (int derx=0 ; derx<=1 ; derx++)
        {
            for (int dery=0 ; dery<=derx ; dery++)
            {
                for (int pixDim=1 ; pixDim<=3 ; pixDim++)
                {
                    Array2dScalarReal64* res = 0;
                    RecGauss(res, imOOO.GetSR64(pixDim-1), scales[si],
                             scales[si], derx, dery, 3);
                    Norm2(res, res);
                    RecGauss(res, res, 1.95*scales[si], 1.95*scales[si],
                             0, 0, 3);
                    imList.push_back(res);
std::cout << "imList.size = " << imList.size() << std::endl;
                }
            }
        }
    }

    return imList;
}

} // namespace Array
} // namespace Core
} // namespace Impala

#endif

---