Segmentation.h

#pragma once

using namespace boost::numeric;
/*waehlt auf einfache Weise Saatpkte aus. Kein member, gehoert so nicht dazu.
return: 1: kontur war leer.*/
int  makeSeeds ( std::vector< std::vector< cv::Vec2i >*>& seedsIn, const std::vector<cv::Point2d>   &cont,
                                        const int offsetX, const int offsetY, const cv::Mat& img);

/*benoetigt GMRes, dies benoetigt ublas, verarbeitet grosse Matrizen. 
Diese lassen sich im debug-Modus _nicht verarbeiten_! 
Z.B. kommt er nicht in endlicher Zeit nicht aus der Residuumsberechnung r=Ax-b raus.
Kann mir das nur so erklaeren, dass Matrix-Vektor-Mult fuer sparse unter debug als nicht-sparse implementiert ist.*/

class SegCoeff: public CoeffFunc
{
private:
   //wichtig ist, dass das Ziel der Pointer nicht durch Parallelisierung veraendert wird. 
    //dafuer traegt class Segmentation sorge. 
    ublas::vector<double>* const
            conductivitiesHorizPtr;
    ublas::vector<double>* const
            conductivitiesVertPtr;
    const int cols;
public:
    virtual double operator()( int i, int j) const ;

    SegCoeff(    ublas::vector<double>* const
            conductivitiesHoriz,
            ublas::vector<double>* const
            conductivitiesVert, const int cols);
    ~SegCoeff(void);
 };

class Segmentation
{
private:
    const cv::Mat& img;
    cv::Mat& segmentedImg;
    const std::vector< std::vector< cv::Vec2i >*>&  seeds;

    ublas::vector<double> 
        conductivitiesHoriz;
    ublas::vector<double> 
        conductivitiesVert;

    void    mkConductivities();
public:
    int operator()() ;
    Segmentation(   cv::Mat& segmentedImg, const cv::Mat& img, 
                    const std::vector< std::vector< cv::Vec2i >*>&  seeds);
    ~Segmentation(void);
};