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

---

Laplacian.h

Go to the documentation of this file.

#ifndef Impala_Core_Array_Laplacian_h
#define Impala_Core_Array_Laplacian_h

#include "Core/Array/GaussDerivative.h"
#include "Core/Array/RecGauss.h"
#include "Core/Array/MulVal.h"
#include "Core/Array/Abs.h"
#include "Core/Array/Add.h"
#include "Core/Array/Set.h"

namespace Impala
{
namespace Core
{
namespace Array
{

void
Laplacian(Array2dScalarReal64*& dst, Array2dScalarReal64* src, Real64 sigma, 
          int gamma, bool useRecGauss, double precision)
{
    /* the Laplacian is multiplied by sigma^gamma */

    Array2dScalarReal64* Lxx = 0;
    Array2dScalarReal64* Lyy = 0;

    if(useRecGauss)
    {
        RecGauss(Lxx, src, sigma, sigma, 2, 0, precision);
        RecGauss(Lyy, src, sigma, sigma, 0, 2, precision);
    } 
    else 
    {
        GaussDerivative(Lxx, src, sigma, 2, 0, precision);
        GaussDerivative(Lyy, src, sigma, 0, 2, precision);
    }
    
    Add(dst, Lxx, Lyy);
    delete Lxx;
    delete Lyy;
    Real64 sigma2 = 1;
    for(int i = 0; i < gamma; i++)
        sigma2 = sigma2 * sigma;
    MulVal(dst, dst, sigma2);
}


// multi-channel version
template<class SrcArrayT>
void
Laplacian(Array2dScalarReal64*& dst, SrcArrayT* src, Real64 sigma, int gamma, 
          bool useRecGauss, double precision)
{
    for (int d=1 ; d<=src->ElemSize() ; d++) 
    {
        Array2dScalarReal64* imIntensity = 0;
        ProjectRange(imIntensity, src, d);
        Array2dScalarReal64* temp = 0;
        Laplacian(temp, imIntensity, sigma, gamma, useRecGauss, precision);
        if(d == 1) 
        {
            Set(dst, temp);
        }
        else 
        {
            Add(dst, dst, temp);
        }
        delete imIntensity;
        delete temp;
    }
}

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

#endif

---