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

---

MakeCooccurenceMatrix.h

Go to the documentation of this file.

#ifndef Impala_Core_Histogram_MakeCooccurenceMatrix_h
#define Impala_Core_Histogram_MakeCooccurenceMatrix_h

#include "Core/Histogram/Histogram2dTem.h"
#include "Core/Array/MakeRoi.h"
#include "Core/Histogram/NgbOutCooccurence.h"
#include "Core/Array/Pattern/PatNgbOp2d.h"
#include "Core/Array/Pattern/PatSetBorder.h"

namespace Impala
{
namespace Core
{
namespace Histogram
{


/*
call function something like this:

    Array::Array2dScalarReal64* im1 = MakeFromValue<Array::Array2dScalarReal64>(0.1, 10, 10);
    Histogram::Histogram2dTem<double> hist(0, 1, 4, 0, 1, 4);
    Histogram::MakeCooccurenceMatrix(&hist, im1, Geometry::Rectangle(0,0,9,9));
    hist.Dump(true);

*/
template<class HistT, class ArrayT>
void
MakeCooccurenceMatrix(HistT* hist, ArrayT* im, Geometry::Rectangle rect, typename ArrayT::ArithType safeBorderValue = -1)
{
    const int borderSize = 1;
    ArrayT* borderedIm = CheckBorderSize(im, borderSize, borderSize);
    typename ArrayT::ArithType dummy = 0;
    Array::Pattern::PatSetBorder(borderedIm, borderSize, borderSize, 
                                 Array::Pattern::BORDERCONSTANT, safeBorderValue,
                                 dummy, dummy, dummy, dummy);
    Geometry::Rectangle borderedRect = GrowBorder(rect, borderSize);
    ArrayT* roi = Array::MakeRoi(borderedIm, borderedRect);
    if(borderedIm != im)
        delete borderedIm;
    roi->GrowBorder(borderSize);

    hist->Clear();
    Histogram::NgbOutCooccurence cooccurence(hist);
    Array::Pattern::PatNgbOp2d(roi, roi, cooccurence, true);
    delete roi;
}

template<class HistT, class ArrayT>
void
MakeCooccurenceMatrix(HistT* hist, ArrayT* im)
{
    MakeCooccurenceMatrix<HistT, ArrayT>(hist, im, Geometry::Rectangle(0,0,im->CW()-1, im->CH()-1));
}

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

#endif

---