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

---

MakeGaussian2d.h

Go to the documentation of this file.

#ifndef Impala_Core_Array_MakeGaussian2d_h
#define Impala_Core_Array_MakeGaussian2d_h

#include "Util/Math.h"
#include "Core/Array/Arrays.h"

namespace Impala
{
namespace Core
{
namespace Array
{


#ifndef M_PI
#define M_PI    3.14159265358979323846
#endif


//#define MAX(a,b) ((a) > (b) ? (a) : (b))
#define ABS(a)   ((a) < 0.0 ? -(a) : (a))


inline Array2dScalarReal64*
MakeGaussian2d(double st, double sr, double phi, int dert, double n)
{
    int         xx, yy, width, height;
    double      l, divido;

    phi = phi*M_PI/180.0;
    if (phi == 0.0) {
        width  = (int)(n*st);
        height = (int)(n*sr);
    } else {
        l = sqrt(sr*sr/(cos(phi)*cos(phi)) + st*st/(sin(phi)*sin(phi)));
        width = (int)(n*ABS(st*st/tan(phi) + sr*sr*tan(phi))/l + 0.5);
        l = sqrt(st*st/(cos(phi)*cos(phi)) + sr*sr/(sin(phi)*sin(phi)));
        height = (int)(n*ABS(sr*sr/tan(phi) + st*st*tan(phi))/l + 0.5);
    }
    divido = 1./(2*M_PI*sr*st);

    double* fData = new double[(2*width+1) * (2*height+1)];

    int i=0;
    for (yy=0; yy<2*height+1; yy++) {
        for (xx=0; xx<2*width+1; xx++) {
            if (dert==2) {
                divido = (pow((cos(phi)*(xx-width)+sin(phi)*(yy-height))/st,
                                        2.0)-1.)/(2*M_PI*sr*st*st);
            }
            fData[i] = (double)(divido * exp(-0.5 *
                (pow(( cos(phi)*(xx-width)+sin(phi)*(yy-height))/st,2.0)+
                pow((-sin(phi)*(xx-width)+cos(phi)*(yy-height))/sr,2.0))));
            i++;
        }
    }
    return ArrayCreate<Array2dScalarReal64>(2*width+1, 2*height+1, 0, 0, fData);
}

#undef ABS

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

#endif

---