Changeset 2939 for GTP/trunk/App/Demos/Vis/FriendlyCulling/src/ShadowMapping.cpp

Timestamp:

09/12/08 22:02:41 (16 years ago)

Author:

mattausch

Message:

lispsm working for standard case (light dir perpendicular to scene), but not even uniform working for all other cases

File:

• GTP/trunk/App/Demos/Vis/FriendlyCulling/src/ShadowMapping.cpp (modified) (14 diffs)

GTP/trunk/App/Demos/Vis/FriendlyCulling/src/ShadowMapping.cpp

@@ -60 +60 @@
         }
 
+        // first create near plane because of precision issues
+        clipPlanes.push_back(frustum.mClipPlanes[4]);
+
         clipPlanes.push_back(frustum.mClipPlanes[0]);
         clipPlanes.push_back(frustum.mClipPlanes[1]);
@@ -65 +68 @@
         clipPlanes.push_back(frustum.mClipPlanes[3]);
         clipPlanes.push_back(frustum.mClipPlanes[5]);
-
-        clipPlanes.push_back(frustum.mClipPlanes[4]);
 
         return Polyhedron::CreatePolyhedron(clipPlanes, sceneBox);
@@ -226 +227 @@
                                                                                  )
 {
-        //return IdentityMatrix();
+        return IdentityMatrix();
 
         AxisAlignedBox3 bounds_ls = GetExtremalPoints(lightSpace, body);
@@ -235 +236 @@
 
         //const float n = 1e2f;
-        const float n = 1e2f;
-        //const float n = ComputeN(bounds_ls) * 100;
+        const float n = 1e6f;
+        //const float n = ComputeN(bounds_ls);
 
         cout << "n: " << n << endl;
@@ -245 +246 @@
         const Vector3 lsNear = lightSpace * nearPt;
 
-        //c start has the x and y coordinate of e, the z coord of the near plane of the light volume
+        // the start point has the x and y coordinate of e, the z coord of the near plane of the light volume
         const Vector3 startPt = Vector3(lsNear.x, lsNear.y, bounds_ls.Max().z);
         //const Vector3 schas = lightSpace * mCamera->GetPosition();
@@ -278 +279 @@
         Matrix4x4 matLispSM;
         
-        matLispSM = GetFrustum(-1.0, 1.0, -1.0, 1.0, n, n + d * 1e10f);
+        matLispSM = GetFrustum(-1.0, 1.0, -1.0, 1.0, n, n + d);
 
         //cout << "lispsm\n" << matLispSM << endl;
@@ -405 +406 @@
 
         VertexArray frustumPoints;
-        IncludeLightVolume(*polyhedron, frustumPoints, mShadowCam->GetDirection(), mSceneBox);
+        //Vector3 lightDir = Vector3(0, 0, -1);
+        IncludeLightVolume(*polyhedron, frustumPoints, mLight->GetDirection(), mSceneBox);
 
 
@@ -445 +447 @@
         lightProj *= frame;
 
-        cout << "here9\n" << lightProj << endl;
-
-        //Matrix4x4 refl2 = ScaleMatrix(-1, -1, 1); 
-        //lightProj = refl2 * lightProj;
-
         const Matrix4x4 matLispSM = 
                 CalcLispSMTransform(lightView * lightProj, extremalPoints, frustumPoints);
@@ -464 +461 @@
 
         lightProj *= transformToGL;
-        //cout << "here4 \n" << lightProj << endl;
 
         AxisAlignedBox3 lightPts = GetExtremalPoints(lightView * lightProj, frustumPoints);
 
-        //cout << "max2: " << lightPts.Max() << endl;
-        //cout << "min2: " << lightPts.Min() << endl;
-
-        AxisAlignedBox3 dummyBox = GetExtremalPoints(IdentityMatrix(), frustumPoints);
-
+        cout << "max: " << lightPts.Max() << " min: " << lightPts.Min() << endl;
         // focus projection matrix on the extremal points => scale to unit cube
-        //Matrix4x4 scaleTranslate = GetFittingProjectionMatrix(extremalPoints);
         Matrix4x4 scaleTranslate = GetFittingProjectionMatrix(lightPts);
-        lightProj = scaleTranslate * lightProj;
-
-        cout << "max: " << lightView * lightProj * dummyBox.Max() << endl;
-        cout << "min: " << lightView * lightProj * dummyBox.Min() << endl;
+
+        cout << "scaleTranslate:\n" << scaleTranslate << endl;
+
+        lightProj = lightProj * scaleTranslate;
+
+        Matrix4x4 mymat = lightView * lightProj;
+
+        AxisAlignedBox3 lightPtsNew = GetExtremalPoints(mymat, frustumPoints);
+
+        cout << "newmax: " << lightPtsNew.Max() << endl;
+        cout << "newmin: " << lightPtsNew.Min() << endl;
 
         // we have to flip the signs in order to tranform to opengl right handed system
@@ -554 +552 @@
                 frustumPoints.push_back(v);
                 
-                // hack: get point surely outside of box
+                // hack: start at point which is guaranteed to be outside of box
                 v -= Magnitude(mSceneBox.Diagonal()) * lightDir;
 
@@ -575 +573 @@
         const float ylen = Magnitude(mSceneBox.Diagonal() * 0.5f);
         
-        //const Vector3 dir = mLight->GetDirection();
-        const Vector3 dir(0, 0, -1);
+        const Vector3 dir = mLight->GetDirection();
+        //const Vector3 dir(0, 0, -1);
 
         // set position so that we can see the whole scene
@@ -627 +625 @@
 
         mLightProjView = lightView * lightProj;
-;
+
         cout << "here3" << endl;
         DEL_PTR(lightPoly);
@@ -690 +688 @@
 void ShadowMap::RenderShadowView(RenderTraverser *renderer, const Matrix4x4 &projView)
 {
-        //const Vector3 dir = mLight->GetDirection();
-        const Vector3 dir(0, 0, -1);
-
+        const Vector3 dir = mLight->GetDirection();
+        //const Vector3 dir(0, 0, -1);
+
+        cout << "dir: " << dir << endl;
         mShadowCam->SetDirection(dir);
 
@@ -735 +734 @@
         lightPoly = CreatePolyhedron(mLightProjView, mSceneBox);
 
+        Frustum frustum(mLightProjView);
+        cout << "near: " << Normalize(frustum.mClipPlanes[4].mNormal) << endl;
+        cout << "far: " << Normalize(frustum.mClipPlanes[5].mNormal) << endl;
+
         glMatrixMode(GL_MODELVIEW);
         glPushMatrix();