source: GTP/trunk/Lib/Vis/Preprocessing/src/CombinedPreprocessor.cpp @ 2066

Revision 2066, 4.4 KB checked in by mattausch, 18 years ago (diff)

worked on integration manual

  • Property svn:executable set to *
Line 
1#include "SceneGraph.h"
2#include "KdTree.h"
3#include "X3dExporter.h"
4#include "Environment.h"
5#include "Polygon3.h"
6#include "VssRay.h"
7#include "RssTree.h"
8#include "ViewCellsManager.h"
9#include "RenderSimulator.h"
10#include "GlRenderer.h"
11#include "SamplingStrategy.h"
12#include "PreprocessorThread.h"
13#include "CombinedPreprocessor.h"
14#include "RayCaster.h"
15
16
17
18namespace GtpVisibilityPreprocessor {
19
20
21const bool pruneInvalidRays = true;
22
23 
24  CombinedPreprocessor::CombinedPreprocessor():
25        Preprocessor()
26  {
27        // this should increase coherence of the samples
28//      Environment::GetSingleton()->GetIntValue("RssPreprocessor.vssSamples", mTotalSamples);
29//      Environment::GetSingleton()->GetIntValue("RssPreprocessor.samplesPerPass", mSamplesPerPass);
30//      Environment::GetSingleton()->GetIntValue("RssPreprocessor.vssSamplesPerPass",
31//                                                                                       mSamplesPerEvaluation);
32       
33       
34  }
35
36
37
38bool
39CombinedPreprocessor::ComputeVisibility()
40{
41 
42  Debug << "Preprocessor type: combined" << endl;
43  cout<<"Combined Preprocessor started\n"<<flush;
44 
45  // if not already loaded, construct view cells from file
46  if (!mLoadViewCells) {
47        // construct view cells using it's own set of samples
48        mViewCellsManager->Construct(this);
49       
50        //-- several visualizations and statistics
51        Debug << "view cells construction finished: " << endl;
52        mViewCellsManager->PrintStatistics(Debug);
53  }
54 
55 
56  // now decode distribution string
57  char buff[1024];
58  Environment::GetSingleton()->GetStringValue("RssPreprocessor.distributions",
59                                                                                          buff);
60 
61 
62  mMixtureDistribution = new MixtureDistribution(*this);
63  mMixtureDistribution->Construct(buff);
64 
65  SimpleRayContainer rays;
66  VssRayContainer vssRays;
67
68  int lastEvaluation = 0;
69  int totalVssRays = 0;
70  mPass = 0;
71 
72  vector<VssRayContainer> rayBuffer;
73  rays.reserve(mSamplesPerPass);
74  rayBuffer.reserve(2*mSamplesPerPass);
75 
76  long startTime = GetTime();
77
78  for (int i=0; i < mTotalSamples; i += mSamplesPerPass, mPass++) {
79        mRayCaster->InitPass();
80       
81        cout<<"Progress :"<<i/1e6f<<"M rays, "<<(100.0f*i)/mTotalSamples<<"%"<<endl;
82       
83        rays.clear();
84        vssRays.clear();
85       
86       
87        mMixtureDistribution->GenerateSamples(mSamplesPerPass, rays);
88       
89        bool doubleRays = true;
90        CastRays(rays, vssRays, doubleRays, pruneInvalidRays);
91        totalVssRays += (int)vssRays.size();
92
93        mMixtureDistribution->ComputeContributions(vssRays);
94
95
96        if (mExportRays) {
97          Debug<<"Exporting rays..."<<endl<<flush;
98          char filename[256];
99
100          VssRayContainer contributingRays;
101          vssRays.GetContributingRays(contributingRays, 0);
102          sprintf(filename, "rss-crays-%04d.x3d", mPass);
103          ExportRays(filename, contributingRays, mExportNumRays);
104
105          sprintf(filename, "rss-rays-i%04d.x3d", mPass);
106          ExportRays(filename, vssRays, mExportNumRays);
107
108
109          if (mExportAnimation) {
110                rayBuffer.push_back(VssRayContainer());
111                vssRays.SelectRays(mExportNumRays, rayBuffer[rayBuffer.size()-1], true);
112          }
113
114          Debug<<"done."<<endl<<flush;
115        }
116
117
118        mMixtureDistribution->UpdateDistributions(vssRays);
119       
120        if (i - lastEvaluation >= mSamplesPerEvaluation) {
121          Real time = TimeDiff(startTime, GetTime());
122          mViewCellsManager->PrintPvsStatistics(mStats);
123          mStats <<
124                "#Pass\n" <<mPass<<endl<<
125                "#Time\n" << time <<endl<<
126                "#TotalSamples\n" <<i<<endl<<
127                "#RssSamples\n" <<totalVssRays<<endl;
128
129          float last = 0.0f;
130          for (int k=0; k < 6; k++) {
131                float ratio = 0.0f;
132                if (k < mMixtureDistribution->mDistributions.size()) {
133                  ratio = mMixtureDistribution->mDistributions[k]->mRatio-last;
134                  last = mMixtureDistribution->mDistributions[k]->mRatio;;
135                }
136                mStats <<
137                  "#Distribution"<<k<<endl<<
138                  ratio<<endl;
139          }
140         
141          lastEvaluation = i;
142         
143          if (renderer) {
144                ComputeRenderError();
145          }
146        }
147       
148       
149        Real time = TimeDiff(startTime, GetTime());
150        if (mTotalTime!= -1 && time/1000 > mTotalTime)
151          break;
152
153#if 0
154        cerr<<"deleting rays"<<endl;
155        //clear all unreferenced rays
156        for (int j=0; j < vssRays.size(); j++) {
157          if (vssRays[j]->RefCount() == 0)
158                delete vssRays[j];
159        }
160        cerr<<"done."<<endl;
161#endif 
162       
163        //      if (!mMixtureDistribution->RequiresRays())
164        //        CLEAR_CONTAINER(vssRays);
165
166       
167  }
168 
169  EvalViewCellHistogram();
170
171  if (mExportAnimation && mExportRays) {
172        char filename[64];
173        sprintf(filename, "rss-rays-i.x3d");
174        //rayBuffer.resize(mPass);
175        ExportRayAnimation(filename, rayBuffer);
176  }
177
178  return true;
179}
180
181}
Note: See TracBrowser for help on using the repository browser.