faulthorintersect.h

Go to the documentation of this file.

#ifndef faulthorintersect_h
#define faulthorintersect_h
                                                                                
/*+
________________________________________________________________________
(C) dGB Beheer B.V.; (LICENSE) http://opendtect.org/OpendTect_license.txt
Author:        Yuancheng Liu
Date:          March 2010
RCS:           $Id$
________________________________________________________________________

-*/

#include "geometrymod.h"
#include "gendefs.h"
#include "paralleltask.h"
#include "executor.h"

class Coord3List;
class Coord3ListImpl;

namespace Geometry
{

class BinIDSurface;
class IndexedShape;
class ExplFaultStickSurface;
class FaultStickSet;


mExpClass(Geometry) FaultBinIDSurfaceIntersector
{
public:
                                FaultBinIDSurfaceIntersector(float horshift,
                                        const BinIDSurface&, 
                                        const ExplFaultStickSurface&,
                                        Coord3List&);
                                ~FaultBinIDSurfaceIntersector() {}

    void                        compute();      
                                
                                //The shape is optional, if not set, we still
                                //compute intersections, stored in crdlist_    
    void                        setShape(const IndexedShape&);
    const IndexedShape*         getShape(bool takeover=true);

protected:

    float                       zshift_;
    Coord3List&                 crdlist_;
    const BinIDSurface&         surf_;
    const IndexedShape*         output_;
    const ExplFaultStickSurface& eshape_;

private:
    void                        sortPointsToLine(TypeSet<Coord3>&,
                                                 TypeSet<Coord3>&);
    const Coord3                findNearestPoint(const Coord3&,
                                                 TypeSet<Coord3>&);
    bool                        findMin(TypeSet<Coord3>&,int&,bool);
    int                         optimizeOrder(TypeSet<Coord3>&);

};

mExpClass(Geometry) BulkFaultBinIDSurfaceIntersector : public Executor
{ mODTextTranslationClass(BulkFaultBinIDSurfaceIntersector)
public:
                                BulkFaultBinIDSurfaceIntersector(
                                    float horshift,
                                    BinIDSurface*,
                                    ObjectSet<FaultStickSet>&,
                                    ObjectSet<Coord3ListImpl>& );

                                ~BulkFaultBinIDSurfaceIntersector() {}

    virtual od_int64            totalNr() const { return totalnr_; }
    virtual od_int64            nrDone() const { return nrdone_; }
    uiString                    uiMessage() const
                                { return tr( "Fault-Horizon "
                                                "Intersector Calculator" ); }
    uiString                    uiNrDoneText() const
                                    { return tr( "Faults handled" ); }

    int                         nextStep();

    const ObjectSet<Coord3ListImpl>     getIntersectionPoints()
                                                { return crdlistset_; }

protected:

    float                                       zshift_;
    ObjectSet<Coord3ListImpl>&                  crdlistset_;
    const BinIDSurface*                         surf_;
    ObjectSet<FaultStickSet>&                   fssset_;
    od_int64                                    nriter_;

    int                                         totalnr_;
    int                                         nrdone_;
};


};


#endif