houghtransform.h

Go to the documentation of this file.

#pragma once
/*
________________________________________________________________________

 (C) dGB Beheer B.V.; (LICENSE) http://opendtect.org/OpendTect_license.txt
 Author:        N. Fredman
 Date:          18-12-2002
________________________________________________________________________

*/


#include "algomod.h"
#include "threadlock.h"
#include "arrayndimpl.h"
#include "ranges.h"

class Array3DInfo;
template <class T> class Array3D;
class SequentialTask;
template <class T, class B> class TopList;
class PlaneFrom3DSpaceHoughTransformTask;
class Plane3;

mExpClass(Algo) PlaneFrom3DSpaceHoughTransform
{
public:
                                PlaneFrom3DSpaceHoughTransform();
    virtual                     ~PlaneFrom3DSpaceHoughTransform();

    void                        setResolution( double dangle,
                                               int distsize );
    int                         getParamSpaceSize() const;
    int                         getNrDistVals() const;

    void                        setClipRate( float );
    float                       clipRate() const;

    void                        setData( const Array3D<float>* );
    ObjectSet<SequentialTask>*  createCalculators();

    TopList<unsigned int, unsigned int>*
                                sortParamSpace(int) const;
    Plane3                      getPlane( int plane ) const;
    int                         getNrPointsAfterClip() const; 
protected:
    void                        incParamPos( int normal, double dist );
    float                       cliprate_;
    TypeSet<unsigned int>       calcpositions_;
    Array3DInfo*                datainfo_;

    Array2D<unsigned int>*      paramspace_;
    double                      deltadist_;
    TypeSet<Coord3>*            normals_;
    friend                      class ::PlaneFrom3DSpaceHoughTransformTask;

    Threads::Lock               paramspacelock_;
};


mExpClass(Algo) LineFrom2DSpaceHoughTransform
{
public:
                                LineFrom2DSpaceHoughTransform(
                                        const Array2D<float>&);
                                ~LineFrom2DSpaceHoughTransform();

    bool                        compute();

    void                        setThreshold(float val,bool above_val);
                                /*<Make hough array based on above/below val.*/
    void                        setLineAngleRange(Interval<float>rg);
                                /*<rg is between 0-PI in radian.*/
    void                        setTopList(int nr)      { toplistnr_ = nr; }
                                /*<nr is the number of top confident lines.*/
    Array2D<unsigned char>*     getResult() const       { return result_; }
                                /*<0/1 array, 1 means on a line. */

//protected:

    bool                        setLineFlag(float radius,float theta);

    const Array2D<float>&       input_;
    Array2D<int>*               origcnt_;       
    Array2D<int>*               hougharr_;
    Array2D<unsigned char>*     result_;

    float                       threshold_;
    bool                        abovethreshold_;
    int                         toplistnr_; 
    Interval<float>             anglerg_;
};