junctiondetector.h

#ifndef JUNCTIONDETECTOR_H
#define JUNCTIONDETECTOR_H

#include <QtGui>
#include <vector>
#include <cmath>

extern "C" {
    #include "emd.h"
}

#define MAX_RECURSIVE_STEPS 1
#define MIN_EIGENVALUE_VALUE 0.1f

using namespace std;

class JunctionDetector
{
public:

    static void detectJunctions(const QImage &image,
        const QImage &cannyEdgeImage,
        const QImage &gradientImage,
        QImage &junctionImage,
        QImage &outputJunctionImage);

    static void detectXJunctions(
        QImage &image,
        QImage &junctionImage,
        QImage &outputXJunctionImage,
        std::vector<std::vector<float>> &luminances);

    struct ColorSignature {
        float N;
        QVector3D pixel;
    };


private:
    static void calculateColorsignature(
        const vector<QVector3D> &imagePixels,
        signature_t &signature);

    static void recursiveSplit(
        const vector<QVector3D> &imagePixels,
        const QMatrix3x3 &R, const QVector3D &m,
        vector<ColorSignature> &colorSignatures);

    static QMatrix3x3 multiplyVectors(
        const QVector3D *vec,
        const QVector3D *vecT);
    
    static void calcCubicEquation(
        const QMatrix3x3 &A,
        float root[3]);
    
    static void eigenvectorSolver(
        const QMatrix3x3 &A,
        float evalue[3],
        QVector3D res[3]);
    
    static int computeRank(QMatrix3x3 M);
    
    static void getComplement1 (
        const QVector3D &U,
        const QVector3D &V,
        QVector3D &W);
    
    static void getComplement2 (
        QVector3D U,
        QVector3D &V,
        QVector3D &W);

    static bool colorSignatureSortFunction(
        const ColorSignature &sig1,
        const ColorSignature &sig2);

    static void growRegion(
        const int &seedX, const int &seedY,
        QImage &image,
        const QColor &initialColor,
        QImage &neighbours);

    static float colorDiff(
        const QColor &color1,
        const QColor &color2);

    static unsigned int recursiveDepth;
    static const QVector2D normalizedVector;
};

#endif