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

---

Weibull.h

Go to the documentation of this file.

#ifndef Impala_Core_Feature_Weibull_h
#define Impala_Core_Feature_Weibull_h

#include "Core/Vector/VectorSet.h"
#include "Core/Histogram/ComputeWeibull.h"
#include "Core/Array/InvWiccest.h"
#include "Core/Array/RGB2Intensity.h"
#include "Core/Array/DivVal.h"
#include "Core/Histogram/MakeHistogram1dSet.h"

namespace Impala
{
namespace Core
{
namespace Feature
{


Vector::VectorSet<Core::Array::Array2dScalarReal64>*
Weibull(Core::Array::Array2dVec3UInt8* im, Real64 sigma)
{
    typedef Array::Array2dScalarReal64 Array2dScalarReal64;
    typedef Vector::VectorTem<Real64> VectorReal64;
    typedef Array::ArraySet<Array2dScalarReal64> InvSetType;
    typedef Histogram::Histogram1dSet<Array2dScalarReal64> HistSetType;
    typedef HistSetType::HistT HistType;
    typedef Vector::VectorSet<Array2dScalarReal64> VecSetType;

    int histBinCount = 1001;
    Real64 histLow = -1;
    Real64 histHigh = 1;
    bool doRot = false;
    bool doC = false;
    int borderSize = 15;

    Geometry::Rectangle rect(15, 15, im->CW() - 15, im->CH() - 15);
    Array::ArraySet<Array2dScalarReal64> invSet;
    Array::InvWiccest(invSet, im, doRot, doC, sigma, 3, true);
    //Array::WriteRawList(invSet, "outW.raw", true);

    HistSetType hSet(histLow, histHigh, histBinCount, invSet.Size());
    Histogram::MakeHistogram1dSet(&hSet, invSet, rect, true);
    VecSetType* res = new VecSetType(true, 2, invSet.Size());
    for (int h=0 ; h<hSet.Size() ; h++)
    {
        VectorReal64 v(2);
        Histogram::ComputeWeibull(v, hSet.Hist(h), false, false);
        res->AddVector(v);
    }

    invSet.Delete();
    return res;
}


Vector::VectorSet<Core::Array::Array2dScalarReal64>*
WeibullIRGB(Core::Array::Array2dVec3UInt8* im, Real64 sigma, bool doMu = false,
            bool doAnderson = false, bool doArjan = false)
{
    typedef Array::Array2dScalarReal64 Array2dScalarReal64;
    typedef Vector::VectorTem<Real64> VectorReal64;
    typedef Array::ArraySet<Array2dScalarReal64> InvSetType;
    typedef Histogram::Histogram1dSet<Array2dScalarReal64> HistSetType;
    typedef HistSetType::HistT HistType;
    typedef Vector::VectorSet<Array2dScalarReal64> VecSetType;

    int histBinCount = 1001;
    Real64 histLow = -1;
    Real64 histHigh = 1;
    int borderSize = 15;

    int vecSize = (doMu) ? 3 : 2;
    if (doAnderson)
        vecSize++;
    Geometry::Rectangle rect(15, 15, im->CW() - 15, im->CH() - 15);
    if (doArjan)
        rect = Geometry::Rectangle(5, 5, 240, im->CH() - 5);
    Array::ArraySet<Array2dScalarReal64> invSet;
    Array::Array2dVec3Real64* imR = 0;
    Array::Set(imR, im);
    Array::DivVal(imR, imR, Array::Element::Vec3Real64(255, 255, 255));
    Array::Array2dScalarReal64* l = 0;
    Array::RGB2Intensity(l, imR);
    RecGauss(invSet.Array("Lx"), l, sigma, sigma, 1, 0, 3.0);
    RecGauss(invSet.Array("Ly"), l, sigma, sigma, 0, 1, 3.0);
    std::vector<Array2dScalarReal64*> rgbList;
    rgbList = InvCompDer(imR, sigma, 3.0, true);
    Array::Set(invSet.Array("Rx"), rgbList[1]);
    Array::Set(invSet.Array("Ry"), rgbList[2]);
    Array::Set(invSet.Array("Gx"), rgbList[4]);
    Array::Set(invSet.Array("Gy"), rgbList[5]);
    Array::Set(invSet.Array("Bx"), rgbList[7]);
    Array::Set(invSet.Array("By"), rgbList[8]);
    Array::ArrayListDelete(&rgbList);
    //Array::WriteRawList(invSet, "outIRGB.raw", true);

    HistSetType hSet(histLow, histHigh, histBinCount, invSet.Size());
    Histogram::MakeHistogram1dSet(&hSet, invSet, rect, true);
    VecSetType* res = new VecSetType(true, vecSize, invSet.Size());
    for (int h=0 ; h<hSet.Size() ; h++)
    {
        VectorReal64 v(vecSize);
        Histogram::ComputeWeibull(v, hSet.Hist(h), doMu, doAnderson);
        res->AddVector(v);
    }

    delete l;
    delete imR;
    invSet.Delete();
    return res;
}


} // namespace Feature
} // namespace Core
} // namespace Impala

#endif

---