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

---

Resample.h

Go to the documentation of this file.

#ifndef Impala_Core_Array_Resample_h
#define Impala_Core_Array_Resample_h

#include <vector>
#include "Util/Math.h"
#include "Core/Array/ProjectRange.h"
//#include "Core/Array/GaussDerivative.h"
#include "Core/Array/RecGauss.h"
#include "Core/Array/Set.h"
#include "Core/Array/Scale.h"
#include "Core/Array/MakeFrom3Images.h"
#include "Core/Array/ArrayListDelete.h"

namespace Impala
{
namespace Core
{
namespace Array
{


//template<class DstArrayT, class SrcArrayT>
inline void
//Resample(DstArrayT*& dst, SrcArrayT* src, int newWidth, int newHeight,
Resample(Array2dVec3UInt8*& dst, Array2dVec3UInt8* src, int newWidth,
         int newHeight)
{
    double scaleX = (double) newWidth / src->CW();
    double scaleY = (double) newHeight / src->CH();

    // reduce by factor of 2 == sigma 0.75
    //double sigmaX = (1.0 / scaleX) * 0.375;
    //double sigmaY = (1.0 / scaleY) * 0.375;
    double sigmaX = 0.75;
    double sigmaY = sigmaX;

    std::vector<Array2dScalarReal64*> resList;
    for (int d=1 ; d<=src->ElemSize() ; d++)
    {
        Array2dScalarUInt8* comp = 0;
        ProjectRange(comp, src, d);
        double precision = 3.0;
        Array2dScalarReal64* compReal = 0;
        //Set(compReal, comp);
        Scale(compReal, comp, scaleX, scaleY, Geometry::LINEAR, true);
        Array2dScalarReal64* p = 0;
        RecGauss(p, compReal, sigmaX, sigmaY, 0, 0, precision);
        resList.push_back(p);
        delete comp;
        delete compReal;
    }
    MakeFrom3Images(dst, resList[0], resList[1], resList[2]);
    ArrayListDelete(&resList);
}

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

#endif

---