source: GTP/trunk/Lib/Vis/Preprocessing/src/SamplingStrategy.h @ 1898

#ifndef _SamplingStategy_H__
#define _SamplingStategy_H__

#include <vector>
using namespace std;

#include "Halton.h"

namespace GtpVisibilityPreprocessor {


class Preprocessor;
struct SimpleRay;
class SimpleRayContainer;
struct VssRayContainer;

/** This class generates a specific sampling strategy.
*/
class SamplingStrategy
{
public:

        /** Sample rays of particular type
        */
        enum 
        {
          DUMMY_DISTRIBUTION = 0,
          DIRECTION_BASED_DISTRIBUTION,
                OBJECT_BASED_DISTRIBUTION,
                DIRECTION_BOX_BASED_DISTRIBUTION,
                SPATIAL_BOX_BASED_DISTRIBUTION,
                VSS_BASED_DISTRIBUTION,
                RSS_BASED_DISTRIBUTION,
                RSS_SILHOUETTE_BASED_DISTRIBUTION,
                OBJECT_DIRECTION_BASED_DISTRIBUTION,
                OBJECTS_INTERIOR_DISTRIBUTION,
                REVERSE_OBJECT_BASED_DISTRIBUTION,
                VIEWCELL_BORDER_BASED_DISTRIBUTION,
                VIEWSPACE_BORDER_BASED_DISTRIBUTION,
                REVERSE_VIEWSPACE_BORDER_BASED_DISTRIBUTION,
          GLOBAL_LINES_DISTRIBUTION,
          GVS
        };

        /** Default constructor
        */
        SamplingStrategy(Preprocessor &preprocessor);

        virtual ~SamplingStrategy();

        /** Each strategy has to implement this function.
        @returns true if generated valid sample.
        */

        virtual int GenerateSamples(const int number, SimpleRayContainer &rays);


        virtual bool GenerateSample(SimpleRay &ray) = 0;
  


  virtual void Update(VssRayContainer &vssRays) {}



  friend bool LowerRatio(const SamplingStrategy *a, const SamplingStrategy *b) {
        return a->mRatio < b->mRatio;
  }

public:

        /// variables usefull for mixed distribution sampling
        int mType;
        int mRays;
        float mContribution;

        int mTotalRays;
        float mTotalContribution;

        float mTime;
        float mRatio;

protected:

        Preprocessor &mPreprocessor;
        static HaltonSequence sHalton;
};


class ObjectBasedDistribution: public SamplingStrategy
{
 public:
  //HaltonSequence mHalton;
  
  ObjectBasedDistribution(Preprocessor &preprocessor):
        SamplingStrategy(preprocessor) {
        mType = OBJECT_BASED_DISTRIBUTION;
  }

private:
  virtual bool GenerateSample(SimpleRay &ray);
};


class ReverseObjectBasedDistribution: public SamplingStrategy
{
 public:
        
         ReverseObjectBasedDistribution(Preprocessor &preprocessor):
           SamplingStrategy(preprocessor) {
           mType = REVERSE_OBJECT_BASED_DISTRIBUTION;
         }

private:
  virtual bool GenerateSample(SimpleRay &ray);
};


class ObjectDirectionBasedDistribution: public SamplingStrategy
{
  //HaltonSequence mHalton;
  
public:
         ObjectDirectionBasedDistribution(Preprocessor &preprocessor):
         SamplingStrategy(preprocessor) {
           mType = OBJECT_DIRECTION_BASED_DISTRIBUTION;
         }
private:
        virtual bool GenerateSample(SimpleRay &ray);
};


class DirectionBasedDistribution: public SamplingStrategy
{
 public:
  DirectionBasedDistribution(Preprocessor &preprocessor):
        SamplingStrategy(preprocessor){
        mType = DIRECTION_BASED_DISTRIBUTION;
  }
private:
  virtual bool GenerateSample(SimpleRay &ray);
};


class DirectionBoxBasedDistribution: public SamplingStrategy
{
 public:
         DirectionBoxBasedDistribution(Preprocessor &preprocessor):
           SamplingStrategy(preprocessor){
           mType = DIRECTION_BOX_BASED_DISTRIBUTION;
         }
           
private:
           virtual bool GenerateSample(SimpleRay &ray);
};


class SpatialBoxBasedDistribution: public SamplingStrategy
{
 public:
  //HaltonSequence mHalton;
  SpatialBoxBasedDistribution(Preprocessor &preprocessor):
        SamplingStrategy(preprocessor){
        mType = SPATIAL_BOX_BASED_DISTRIBUTION;
  }
  
private:
           virtual bool GenerateSample(SimpleRay &ray);
};


class ViewSpaceBorderBasedDistribution: public SamplingStrategy
{
 public:
         ViewSpaceBorderBasedDistribution(Preprocessor &preprocessor):
           SamplingStrategy(preprocessor){
           mType = VIEWSPACE_BORDER_BASED_DISTRIBUTION;
         }
         
private:
         virtual bool GenerateSample(SimpleRay &ray);
};


class ReverseViewSpaceBorderBasedDistribution: public SamplingStrategy
{
 public:
         ReverseViewSpaceBorderBasedDistribution(Preprocessor &preprocessor):
           SamplingStrategy(preprocessor){
           mType = REVERSE_VIEWSPACE_BORDER_BASED_DISTRIBUTION;
         }
         
private:
         virtual bool GenerateSample(SimpleRay &ray);
};




class ViewCellBorderBasedDistribution: public SamplingStrategy
{
 public:
  ViewCellBorderBasedDistribution(Preprocessor &preprocessor):
        SamplingStrategy(preprocessor) {
        mType = VIEWCELL_BORDER_BASED_DISTRIBUTION;

  }
  
  virtual bool GenerateSample(SimpleRay &ray);
};

class GlobalLinesDistribution: public SamplingStrategy
{
public:
  //Halton<4> halton;
  //HaltonSequence mHalton;

  GlobalLinesDistribution(Preprocessor &preprocessor):
        SamplingStrategy(preprocessor) {
        mType = GLOBAL_LINES_DISTRIBUTION;
  }
  
  virtual bool GenerateSample(SimpleRay &ray);
};


/** This strategy generates samples inside of the objects, e.g.,
        for sampling the inside of a colon.
*/
/*class ObjectsInteriorDistribution: public SamplingStrategy
{
 public:
         ObjectsInteriorDistribution(Preprocessor &preprocessor):
         SamplingStrategy(preprocessor) {}
         
         virtual bool GenerateSample(SimpleRay &ray);
};
*/

class MixtureDistribution : public SamplingStrategy 
{
public:
  // halton sequence generator for deciding between distributions
  //HaltonSequence mHalton;
  
  // container for the distributions
  vector<SamplingStrategy *> mDistributions;

  MixtureDistribution(Preprocessor &preprocessor):
        SamplingStrategy(preprocessor)
  {
  }

  // has to called before first usage
  void Init();

  // equalize distribution contributions
  void
  Reset();
  
  // Generate a new sample according to a mixture distribution
  virtual bool GenerateSample(SimpleRay &ray);

  // add contributions of the sample to the strategies
  void ComputeContributions(VssRayContainer &vssRays);

  void
  UpdateRatios();

  // construct distribution mixyure from a string describing the required distributions
  bool
  Construct(char *str);

  
};

};

#endif