Ignore:
Timestamp:
12/13/06 23:12:27 (18 years ago)
Author:
bittner
Message:

combined preprocessor

Location:
GTP/trunk/Lib/Vis/Preprocessing/src
Files:
2 added
9 edited

Legend:

Unmodified
Added
Removed
  • GTP/trunk/Lib/Vis/Preprocessing/src/Makefile

    r1884 r1891  
    11############################################################################# 
    22# Makefile for building: preprocessor 
    3 # Generated by qmake (2.00a) (Qt 4.1.2) on: st 13. XII 00:50:24 2006 
     3# Generated by qmake (2.00a) (Qt 4.1.2) on: st 13. XII 23:09:08 2006 
    44# Project:  preprocessor.pro 
    55# Template: app 
  • GTP/trunk/Lib/Vis/Preprocessing/src/PreprocessorFactory.cpp

    r1486 r1891  
    66#include "RenderSampler.h" 
    77#include "GvsPreprocessor.h" 
     8#include "CombinedPreprocessor.h" 
    89 
    910 
     
    1718                return new VssPreprocessor(); 
    1819        } 
    19          
     20 
     21        if (preprocessorType == "combined") 
     22        { 
     23                return new CombinedPreprocessor(); 
     24        } 
     25 
    2026        if (preprocessorType == "rss") 
    2127        { 
  • GTP/trunk/Lib/Vis/Preprocessing/src/RssPreprocessor.cpp

    r1884 r1891  
    435435  Randomize(0); 
    436436 
    437   if (1) { 
     437  if (0) { 
    438438        Exporter *exporter = Exporter::GetExporter("samples.wrl"); 
    439439    exporter->SetFilled(); 
     
    501501  char *curr = buff; 
    502502  mUseRssTree = false; 
     503 
    503504  while (1) { 
    504505        char *e = strchr(curr,'+'); 
     
    542543  mMixtureDistribution->Init(); 
    543544  bool oldGenerator = false; 
     545 
     546  const int batch = 100000; 
    544547   
    545548  if (mLoadInitialSamples) { 
     
    579582           
    580583        } else { 
    581           const int batch = 1000; 
    582584          for (int i=0; i < mInitialSamples; i+=batch) { 
    583585                rays.clear(); 
  • GTP/trunk/Lib/Vis/Preprocessing/src/RssTree.cpp

    r1884 r1891  
    2929#include "SceneGraph.h" 
    3030#include "SamplingStrategy.h" 
     31#include "ViewCellsManager.h" 
    3132 
    3233 
     
    27332734// if small very high importance of the last sample 
    27342735// if 1.0f then weighs = 1 1/2 1/3 1/4 
    2735 float passSampleWeightDecay = 1.0f; 
     2736float passSampleWeightDecay = 2.0f; 
    27362737//float passSampleWeightDecay = 1000.0f; 
    27372738 
     
    27402741{ 
    27412742  int passDiff = mCurrentPass - pass; 
     2743  if (passDiff < 0) { 
     2744        cerr<<"Error: passdif < 0"<<passDiff<<endl; 
     2745        exit(1); 
     2746  } 
    27422747  float weight; 
    2743   if (0)  
     2748  if (1)  
    27442749        weight = 1.0f/(passDiff + passSampleWeightDecay); 
    27452750  else  
     
    27882793RssTree::ComputeImportance(RssTreeLeaf *leaf)  
    27892794{ 
    2790 #if 1 
     2795#if 0 
    27912796  if (leaf->mTotalRays) 
    27922797        leaf->mImportance = leaf->mTotalContribution / leaf->mTotalRays; 
     
    28562861         
    28572862        if (sumWeights != 0.0f) { 
     2863#if 1 
    28582864          leaf->mImportance = 
    2859                         sqr(((weightAbsContributions*sumContributions + 
    2860                                                 (1.0f - weightAbsContributions)*sumRelContributions)/sumWeights)); 
    2861            
    2862           //      leaf->mImportance = 
    2863           //            (weightAbsContributions*maxContribution + 
    2864           //             (1.0f - weightAbsContributions)*maxRelContribution)/sumWeights; 
    2865            
     2865                pow(((weightAbsContributions*sumContributions + 
     2866                          (1.0f - weightAbsContributions)*sumRelContributions)/sumWeights), 3); 
     2867#else  
     2868          leaf->mImportance = 
     2869                (weightAbsContributions*maxContribution + 
     2870                 (1.0f - weightAbsContributions)*maxRelContribution)/sumWeights; 
     2871#endif 
    28662872        } else 
    28672873          leaf->mImportance = 0.0f; 
     
    28772883RssTreeNode::GetImportance()  const 
    28782884{ 
    2879   //  return sqr(mImportance); 
     2885  // return mImportance*mImportance; 
    28802886  return mImportance; 
    28812887} 
     
    31203126  float r[5]; 
    31213127 
    3122   if (mRssTree == NULL) 
    3123         return false; 
    3124    
    31253128  mHalton.GetNext(5, r); 
    31263129 
    3127   SimpleRayContainer rays; 
    3128   mRssTree->GenerateRay(r, rays); 
    3129   ray = rays[0]; 
    3130    
    3131   return false; 
     3130  if (mRssTree == NULL) { 
     3131        // use direction based distribution 
     3132        Vector3 origin, direction; 
     3133        mPreprocessor.mViewCellsManager->GetViewPoint(origin, Vector3(r[0], r[1], r[2])); 
     3134         
     3135        direction = UniformRandomVector(r[3], r[4]); 
     3136 
     3137        const float pdf = 1.0f; 
     3138        ray = SimpleRay(origin, direction, RSS_BASED_DISTRIBUTION, pdf); 
     3139         
     3140  } else { 
     3141         
     3142        SimpleRayContainer rays; 
     3143        mRssTree->GenerateRay(r, rays); 
     3144        ray = rays[0]; 
     3145  } 
     3146   
     3147  return true; 
    31323148} 
    31333149 
     
    31413157        /// compute view cell contribution of rays if view cells manager already constructed 
    31423158        //  mViewCellsManager->ComputeSampleContributions(mVssRays, true, false); 
    3143          
    31443159        mRssTree->Construct(mPreprocessor.mObjects, vssRays); 
    31453160        mRssTree->stat.Print(mPreprocessor.mStats); 
     
    31473162         
    31483163  } else { 
     3164        mRssTree->SetPass(mPass); 
    31493165        Debug<<"Adding rays...\n"<<flush; 
    31503166        mRssTree->AddRays(vssRays); 
  • GTP/trunk/Lib/Vis/Preprocessing/src/SamplingStrategy.cpp

    r1888 r1891  
    55#include "ViewCellsManager.h" 
    66#include "AxisAlignedBox3.h" 
    7  
     7#include "RssTree.h" 
    88 
    99namespace GtpVisibilityPreprocessor { 
     
    475475  for (int i=0; i < mDistributions.size(); i++) { 
    476476        // small non-zero value 
    477         mDistributions[i]->mRays = 1; 
     477        mDistributions[i]->mTotalRays = 0; 
    478478        // unit contribution per ray 
    479         mDistributions[i]->mContribution = 1.0f; 
     479        mDistributions[i]->mTotalContribution = 0.0f; 
    480480  } 
    481481  UpdateRatios(); 
     
    490490 
    491491  // pickup a distribution 
    492   for (int i=0; i < mDistributions.size(); i++) 
     492  for (int i=0; i < mDistributions.size()-1; i++) 
    493493        if (r < mDistributions[i]->mRatio) 
    494494          break; 
    495    
     495 
    496496  return mDistributions[i]->GenerateSample(ray); 
    497497} 
     
    505505  VssRayContainer::iterator it = vssRays.begin(); 
    506506 
     507  for (i=0; i < mDistributions.size(); i++) { 
     508        mDistributions[i]->mContribution = 0; 
     509        mDistributions[i]->mRays = 0; 
     510  } 
     511 
    507512  for(; it != vssRays.end(); ++it) { 
    508513        VssRay *ray = *it; 
    509         for (i=0; i < mDistributions.size()-1; i++) 
     514        for (i=0; i < mDistributions.size()-1; i++) { 
    510515          if (mDistributions[i]->mType == ray->mDistribution) 
    511516                break; 
     517        } 
     518  
    512519        float contribution = 
    513           mPreprocessor.mViewCellsManager->ComputeSampleContributions(vssRays, true, false); 
     520          mPreprocessor.mViewCellsManager->ComputeSampleContribution(*ray, true, false); 
    514521        mDistributions[i]->mContribution += contribution; 
    515522        mDistributions[i]->mRays ++; 
     523 
     524        mDistributions[i]->mTotalContribution += contribution; 
     525        mDistributions[i]->mTotalRays ++; 
    516526  } 
    517527 
    518528  // now update the distributions with all the rays 
    519   for (i=0; i < mDistributions.size(); i++) 
     529  for (i=0; i < mDistributions.size(); i++) { 
    520530        mDistributions[i]->Update(vssRays); 
     531  } 
    521532   
    522533   
     
    531542  int i; 
    532543  for (i=0; i < mDistributions.size(); i++) { 
    533         float importance = mDistributions[i]->mContribution/mDistributions[i]->mRays; 
     544        cout<<i<<": c="<<mDistributions[i]->mContribution<< 
     545          " rays="<<mDistributions[i]->mRays<<endl; 
     546        float importance = 0.0f; 
     547        if (mDistributions[i]->mRays != 0) { 
     548          importance = pow(mDistributions[i]->mContribution/mDistributions[i]->mRays, 3); 
     549        } 
     550        if (importance < Limits::Small) 
     551          importance = Limits::Small; 
    534552        mDistributions[i]->mRatio = importance; 
    535553        sum += importance; 
    536554  } 
    537    
     555 
     556  const float minratio = 0.05f; 
     557  float threshold = minratio*sum; 
     558  for (i=0; i < mDistributions.size(); i++) { 
     559        if (mDistributions[i]->mRatio < threshold) { 
     560          sum += threshold - mDistributions[i]->mRatio; 
     561          mDistributions[i]->mRatio = threshold; 
     562        } 
     563  } 
    538564   
    539565  mDistributions[0]->mRatio /= sum; 
     566         
    540567  for (i=1; i < mDistributions.size(); i++) { 
    541568        float r = mDistributions[i]->mRatio / sum; 
     
    548575  cout<<endl; 
    549576} 
    550 } 
    551  
    552  
     577 
     578 
     579 
     580bool 
     581MixtureDistribution::Construct(char *str) 
     582{ 
     583  char *curr = str; 
     584 
     585  while (1) { 
     586        char *e = strchr(curr,'+'); 
     587        if (e!=NULL) { 
     588          *e=0; 
     589        } 
     590         
     591        if (strcmp(curr, "rss")==0) { 
     592          mDistributions.push_back(new RssBasedDistribution(mPreprocessor)); 
     593        } else 
     594          if (strcmp(curr, "object")==0) { 
     595                mDistributions.push_back(new ObjectBasedDistribution(mPreprocessor)); 
     596          } else 
     597                if (strcmp(curr, "spatial")==0) { 
     598                  mDistributions.push_back(new SpatialBoxBasedDistribution(mPreprocessor)); 
     599                } else 
     600                  if (strcmp(curr, "global")==0) { 
     601                        mDistributions.push_back(new GlobalLinesDistribution(mPreprocessor)); 
     602                  } else 
     603                        if (strcmp(curr, "direction")==0) { 
     604                          mDistributions.push_back(new DirectionBasedDistribution(mPreprocessor)); 
     605                        } else 
     606                          if (strcmp(curr, "object_direction")==0) { 
     607                                mDistributions.push_back(new ObjectDirectionBasedDistribution(mPreprocessor)); 
     608                          } else 
     609                                if (strcmp(curr, "reverse_object")==0) { 
     610                                  mDistributions.push_back(new ReverseObjectBasedDistribution(mPreprocessor)); 
     611                                } else 
     612                                  if (strcmp(curr, "reverse_viewspace_border")==0) { 
     613                                        mDistributions.push_back(new ReverseViewSpaceBorderBasedDistribution(mPreprocessor)); 
     614                                  }  
     615         
     616        if (e==NULL) 
     617          break; 
     618        curr = e+1; 
     619  } 
     620 
     621  Init(); 
     622  return true; 
     623} 
     624 
     625} 
     626 
     627 
  • GTP/trunk/Lib/Vis/Preprocessing/src/SamplingStrategy.h

    r1888 r1891  
    7373        int mRays; 
    7474        float mContribution; 
    75         float mTime; 
    76         float mRatio; 
    7775 
    7876  int mTotalRays; 
    7977  float mTotalContribution; 
     78 
     79  float mTime; 
     80        float mRatio; 
     81 
    8082protected: 
    8183 
     
    159161  SpatialBoxBasedDistribution(Preprocessor &preprocessor): 
    160162        SamplingStrategy(preprocessor){ 
    161         mType = DIRECTION_BOX_BASED_DISTRIBUTION; 
     163        mType = SPATIAL_BOX_BASED_DISTRIBUTION; 
    162164  } 
    163165   
     
    172174         ViewSpaceBorderBasedDistribution(Preprocessor &preprocessor): 
    173175           SamplingStrategy(preprocessor){ 
    174            mType = SPATIAL_BOX_BASED_DISTRIBUTION; 
     176           mType = VIEWSPACE_BORDER_BASED_DISTRIBUTION; 
    175177         } 
    176178          
     
    185187         ReverseViewSpaceBorderBasedDistribution(Preprocessor &preprocessor): 
    186188           SamplingStrategy(preprocessor){ 
    187            mType = SPATIAL_BOX_BASED_DISTRIBUTION; 
     189           mType = REVERSE_VIEWSPACE_BORDER_BASED_DISTRIBUTION; 
    188190         } 
    189191          
     
    199201 public: 
    200202  ViewCellBorderBasedDistribution(Preprocessor &preprocessor): 
    201         SamplingStrategy(preprocessor) {} 
     203        SamplingStrategy(preprocessor) { 
     204        mType = VIEWCELL_BORDER_BASED_DISTRIBUTION; 
     205 
     206  } 
    202207   
    203208  virtual bool GenerateSample(SimpleRay &ray); 
     
    261266  UpdateRatios(); 
    262267 
     268  // construct distribution mixyure from a string describing the required distributions 
     269  bool 
     270  Construct(char *str); 
     271 
    263272   
    264273}; 
  • GTP/trunk/Lib/Vis/Preprocessing/src/ViewCellsManager.cpp

    r1889 r1891  
    3030 
    3131#define KD_PVS_AREA (1e-5f) 
    32  
    33  
    3432 
    3533namespace GtpVisibilityPreprocessor { 
  • GTP/trunk/Lib/Vis/Preprocessing/src/default.env

    r1884 r1891  
    4848#       type sampling 
    4949#       type vss 
    50         type rss 
     50#       type rss 
     51        type combined 
    5152#       type render 
    5253        detectEmptyViewSpace true 
     
    8889} 
    8990 
     91 
    9092RssPreprocessor { 
    9193 
    9294#       distributions rss+spatial+object_direction 
    93         distributions object_direction+spatial 
     95        distributions rss+object_direction+spatial+object+direction 
    9496#       distributions object_direction 
    9597 
    96         samplesPerPass 1000 
     98        samplesPerPass 200000 
    9799        initialSamples 2000000 
    98100        vssSamples 50000000 
    99         vssSamplesPerPass 2000000 
     101        vssSamplesPerPass 500000 
    100102        useImportanceSampling true 
    101103 
     
    133135        maxCostRatio 1.0 
    134136        maxRayContribution 1.0 
    135         maxRays         1000000 
     137        maxRays         10000000 
    136138        maxTotalMemory  400 
    137139        maxStaticMemory 200 
  • GTP/trunk/Lib/Vis/Preprocessing/src/preprocessor.pro

    r1883 r1891  
    109109InternalRayCaster.cpp IntelRayCaster.cpp \ 
    110110RayCaster.cpp PreprocessorFactory.cpp GvsPreprocessor.cpp \ 
    111 Trackball.cpp ObjExporter.cpp SubdivisionCandidate.cpp QtGlViewer.cpp Mailable.cpp 
     111Trackball.cpp ObjExporter.cpp SubdivisionCandidate.cpp QtGlViewer.cpp \ 
     112Mailable.cpp \ 
     113CombinedPreprocessor.cpp 
    112114 
    113115 
Note: See TracChangeset for help on using the changeset viewer.