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

---

ComputeWeibull.h

Go to the documentation of this file.

#ifndef Impala_Core_Histogram_ComputeWeibull_h
#define Impala_Core_Histogram_ComputeWeibull_h

#include "Core/Geometry/RectangleSet.h"
#include "Core/Histogram/MakeHistogram1dSet.h"
#include "Core/Histogram/FitWeibullMarginal.h"

namespace Impala
{
namespace Core
{
namespace Histogram
{


template <class VecElemT, class HistElemT>
void
ComputeWeibull(Vector::VectorTem<VecElemT>& v, Vector::VectorTem<HistElemT> hist,
               HistElemT histLow, HistElemT histHigh,
               bool storeMu, bool storeA)
{
    FitWeibullMarginal<VecElemT> fit(hist.GetData(), hist.Size(), histLow,
                                     histHigh);
    v[0] = fit.Beta();
    v[1] = fit.Gamma();
    int idx = 2;
    if (storeMu)
        v[idx++] = fit.Mu();
    if (storeA)
        v[idx] = fit.A2();
}

template <class VecElemT, class HistElemT>
void
ComputeWeibull(Vector::VectorTem<VecElemT>& v, Histogram1dTem<HistElemT> hist,
               bool storeMu, bool storeA)
{
    FitWeibullMarginal<VecElemT> fit(hist.GetData(), hist.GetBinCount(),
                                     hist.GetLow(), hist.GetHigh());
    v[0] = fit.Beta();
    v[1] = fit.Gamma();
    int idx = 2;
    if (storeMu)
        v[idx++] = fit.Mu();
    if (storeA)
        v[idx] = fit.A2();
}

template <class VecElemT, class HistArrayT>
void
ComputeWeibull(Vector::VectorTem<VecElemT>& v, Histogram1dSet<HistArrayT>* hSet,
               bool storeMu, bool storeA)
{
    typedef typename Histogram1dSet<HistArrayT>::HistT HistType;

    int nrPar = 2 + ((storeMu) ? 1 : 0) + ((storeA) ? 1 : 0);
    for (int i=0 ; i<hSet->Size() ; i++)
    {
        HistType hist = hSet->Hist(i);
        FitWeibullMarginal<VecElemT> fit(hist.GetData(), hist.GetBinCount(),
                                         hist.GetLow(), hist.GetHigh());
        int idx = i*nrPar;
        v[idx++] = fit.Beta();
        v[idx++] = fit.Gamma();
        if (storeMu)
            v[idx++] = fit.Mu();
        if (storeA)
            v[idx++] = fit.A2();
    }
}

template <class VecArrayT, class ImArrayT, class HistArrayT>
void
ComputeWeibullSet(Vector::VectorSet<VecArrayT>* vSet, Array::ArraySet<ImArrayT>& imSet,
                  const Geometry::RectangleSet& rects, Histogram1dSet<HistArrayT>* hSet,
                  bool storeMu, bool storeA)
{
    typedef typename Vector::VectorSet<VecArrayT>::ElemT VecElemT;
    typedef typename Vector::VectorSet<VecArrayT>::VectorT VectorT;

    for (int i=0 ; i<rects.size() ; i++)
    {
        MakeHistogram1dSet<HistArrayT,ImArrayT>(hSet, imSet, rects[i]);
        VectorT v = vSet->GetVector(i, true);
        ComputeWeibull<VecElemT,HistArrayT>(v, hSet, storeMu, storeA);
    }
    vSet->SetSize(rects.size());
}

} // namespace Histogram
} // namespace Core
} // namespace Impala

#endif

---