r480: no message
[ctsim.git] / src / graph3dview.cpp
index eb5cbf4a8afd8a72dfa5c9c989e3e158b7fae250..c9b28f70e0bf6e4fd4ac571d682a94fce2c38f30 100644 (file)
@@ -9,7 +9,7 @@
 **  This is part of the CTSim program
 **  Copyright (c) 1983-2001 Kevin Rosenberg
 **
-**  $Id: graph3dview.cpp,v 1.4 2001/01/31 01:01:22 kevin Exp $
+**  $Id: graph3dview.cpp,v 1.5 2001/02/02 00:46:38 kevin Exp $
 **
 **  This program is free software; you can redistribute it and/or modify
 **  it under the terms of the GNU General Public License (version 2) as
 #include <sstream_subst>
 #endif
 
+// Rainbow: Purple->Blue->Cyan->Green->Yellow->Red = (1,0,1)-(0,0,1)-(0,1,1)-(0,1,0)-(1,1,0)-(1,0,0)
+static void 
+intensityToColor (double dIntensity, float vecColor[3])
+{
+  double dRange = dIntensity * 5;
+  int iRange = static_cast<int>(floor (dRange));
+  double dFrac = dRange - iRange;
+  
+  switch (iRange) {
+  case 0:
+    vecColor[0] = 1 - dFrac; vecColor[1] = 0; vecColor[2] = 1;
+    break;
+  case 1:
+    vecColor[0] = 0; vecColor[1] = dFrac; vecColor[2] = 1;
+    break;
+  case 2:
+    vecColor[0] = 0; vecColor[1] = 1; vecColor[2] = 1 - dFrac;
+    break;
+  case 3:
+    vecColor[0] = dFrac; vecColor[1] = 1; vecColor[2] = 0;
+    break;
+  case 4:
+    vecColor[0] = 0; vecColor[1] = 1 - dFrac; vecColor[2] = 0;
+    break;
+  case 5:
+    vecColor[0] = 1; vecColor[1] = 0; vecColor[2] = 0;
+    break;
+  }
+}
 
 //***********************************************************************
 // Function: CalculateVectorNormal
@@ -110,17 +139,25 @@ IMPLEMENT_DYNAMIC_CLASS(Graph3dFileView, wxView)
 
 BEGIN_EVENT_TABLE(Graph3dFileView, wxView)
 EVT_MENU(IFMENU_FILE_PROPERTIES, Graph3dFileView::OnProperties)
+EVT_MENU(GRAPH3D_VIEW_LIGHTING, Graph3dFileView::OnLighting)
+EVT_MENU(GRAPH3D_VIEW_COLOR, Graph3dFileView::OnColor)
+EVT_MENU(GRAPH3D_VIEW_SMOOTH, Graph3dFileView::OnSmooth)
+EVT_MENU(GRAPH3D_VIEW_SURFACE, Graph3dFileView::OnSurface)
 END_EVENT_TABLE()
 
 Graph3dFileView::Graph3dFileView ()
-: m_pFileMenu(NULL)
+: m_pFileMenu(NULL), m_pViewMenu(NULL)
 {
-  m_bUseVertexArrays = GL_FALSE;
-  m_bDoubleBuffer = GL_TRUE;
-  m_bSmooth = GL_TRUE;
-  m_bLighting = GL_TRUE;
+  m_bUseVertexArrays = false;
+  m_bDoubleBuffer = true;
+  m_bSmooth = true;
+  m_bLighting = true;
+  m_bSurface = true;
+  m_bLighting = true;
+  m_bColor = true;
   m_dXRotate = 0;
   m_dYRotate = 0;
+  m_dZRotate = 0;
 }
 
 Graph3dFileView::~Graph3dFileView()
@@ -152,6 +189,10 @@ Graph3dFileView::OnCreate (wxDocument *doc, long WXUNUSED(flags) )
   m_pFrame->Show(true);
   Activate(true);
   
+  m_pViewMenu->Check (GRAPH3D_VIEW_COLOR, m_bColor);
+  m_pViewMenu->Check (GRAPH3D_VIEW_LIGHTING, m_bLighting);
+  m_pViewMenu->Check (GRAPH3D_VIEW_SMOOTH, m_bSmooth);
+  m_pViewMenu->Check (GRAPH3D_VIEW_SURFACE, m_bSurface);
   return true;
 } 
 
@@ -188,67 +229,112 @@ Graph3dFileView::CreateCanvas (wxFrame* parent)
 void
 Graph3dFileView::DrawSurface()
 {
-  int nVertices = GetDocument()->m_nVertices;
-  glTripleFloat* pVertices = GetDocument()->m_pVertices;
-  glTripleFloat* pNormals = GetDocument()->m_pNormals;
-
 #ifdef GL_EXT_vertex_array
   if (m_bUseVertexArrays) {
-    glDrawArraysEXT( GL_TRIANGLE_STRIP, 0, nVertices );
+    //    glDrawArraysEXT( GL_TRIANGLE_STRIP, 0, nVertices );
   }
   else {
 #endif
-       double edge = 1.;       
-       unsigned int nx = GetDocument()->m_nx;
-  unsigned int ny = GetDocument()->m_ny;
-  const ImageFileArray v = GetDocument()->m_array;
-  if (nx == 0 || ny == 0 || ! v)
-    return;
-
-  double dMin = v[0][0];
-  double dMax = dMin;
-  unsigned int ix;
-  for (ix = 0; ix < nx; ix++)
-    for (unsigned int iy = 0; iy < ny; iy++)
-      if (v[ix][iy] < dMin)
-        dMin = v[ix][iy];
-      else if (v[ix][iy] > dMax)
-        dMax = v[ix][iy];
-
-  double actOffset = dMin;
-  double actScale = 0.3 * sqrt(nx*nx+ny*ny) / (dMax - dMin);
-
-  glRotatef( m_dYRotate, 0.0, 0.0, 1.0 );
-  glRotatef( m_dXRotate, 1.0, 0.0, 0.0 );
-  glTranslatef (-static_cast<double>(nx) / 2, -static_cast<double>(ny) / 2, 0);
-
-//     glNewList(opnListNum++,GL_COMPILE);             
-               for (ix = 0; ix < nx-1; ix++) {                 
-                       for (unsigned int iy = 0; iy < ny-1; iy++) {                    
-                               
-         float p1[3], p2[3], p3[3], p4[3];
-             float n1[3], n2[3], n3[3], n4[3];
-                               glBegin(GL_LINE_LOOP);
-                                       
-                               p1[0] = ix;  p1[1] = actScale * (v[ix][iy] + actOffset); p1[2] = iy;
-                               p2[0] = ix+1; p2[1] = actScale * (v[ix+1][iy] + actOffset); p2[2] = iy; 
-                               p3[0] = ix+1; p3[1] = actScale * (v[ix+1][iy+1] + actOffset); p3[2] = iy;
-                               p4[0] = ix;  p4[1] = actScale * (v[ix][iy+1] + actOffset); p4[2] = iy;
-                                                                                       
-                               n1[0] = -(p2[1] - p1[1])*(p3[2] - p1[2]) + (p2[2] - p1[2])*(p3[1] - p2[1]);
-                               n1[1] = -(p2[2] - p1[2])*(p3[0] - p2[0]) + (p2[0] - p1[0])*(p3[2] - p2[2]);
-                               n1[2] = -(p2[0] - p1[0])*(p3[1] - p2[1]) + (p2[1] - p1[1])*(p3[0] - p2[0]);
-                                       
-                               glVertex3fv(p1); glNormal3fv(n1);                                       
-                               glVertex3fv(p2); glNormal3fv(n1);                                       
-                               glVertex3fv(p3); glNormal3fv(n1);                                       
-                               glVertex3fv(p4); glNormal3fv(n1);                                                                                                                                               
-        glEnd();                               
-                       }                       
-                       
-               }
-       glEndList();
-               
+    double edge = 1.;  
+    unsigned int nx = GetDocument()->nx();
+    unsigned int ny = GetDocument()->ny();
+    const ImageFileArrayConst v = GetDocument()->getArray();
+    if (nx == 0 || ny == 0 || ! v)
+      return;
+    
+        glRotatef( m_dXRotate, 1.0, 0.0, 0.0 );
+        glRotatef( m_dZRotate, 0.0, 1.0, 0.0 );
+        glRotatef( m_dYRotate, 0.0, 0.0, 1.0 );
+        glTranslatef (-static_cast<double>(nx) / 2., 0, -static_cast<double>(ny) / 2.);
+        
+    InitMaterials(); 
+    
+    if (m_bSmooth) {
+      glShadeModel (GL_SMOOTH);
+    } else {
+      glShadeModel (GL_FLAT);
+    }
+    
+    if (m_bLighting) {
+      glEnable (GL_LIGHTING);
+    } else {
+      glDisable (GL_LIGHTING);
+    }
+    
+    double dMin = v[0][0];
+    double dMax = dMin;
+    unsigned int ix;
+    for (ix = 0; ix < nx; ix++)
+      for (unsigned int iy = 0; iy < ny; iy++)
+        if (v[ix][iy] < dMin)
+          dMin = v[ix][iy];
+        else if (v[ix][iy] > dMax)
+          dMax = v[ix][iy];
+        
+        double dIntensityScale = dMax - dMin;
+        double actOffset = dMin;
+        double actScale = 0.3 * sqrt(nx*nx+ny*ny) / (dMax - dMin);
+        
+        //     glNewList(opnListNum++,GL_COMPILE);             
+        if (! m_bColor)
+          glColor3f (1.0, 1.0, 1.0);
+        
+        glDisable (GL_CULL_FACE);
+        for (ix = 0; ix < nx-1; ix++) {                        
+          for (unsigned int iy = 0; iy < ny-1; iy++) {                 
+            
+            float p1[3], p2[3], p3[3], p4[3];
+            float n1[3], n2[3], n3[3], n4[3];
+            if (m_bSurface)
+              glBegin(GL_QUADS);
+            else
+              glBegin(GL_LINE_LOOP);
+            
+            p1[0] = ix;  p1[1] = actScale * (v[ix][iy] + actOffset); p1[2] = iy;
+            p2[0] = ix+1; p2[1] = actScale * (v[ix+1][iy] + actOffset); p2[2] = iy; 
+            p3[0] = ix+1; p3[1] = actScale * (v[ix+1][iy+1] + actOffset); p3[2] = iy+1;
+            p4[0] = ix;  p4[1] = actScale * (v[ix][iy+1] + actOffset); p4[2] = iy+1;
+            
+            //                         n1[0] = -(p2[1] - p1[1])*(p3[2] - p1[2]) + (p2[2] - p1[2])*(p3[1] - p2[1]);
+            //                         n1[1] = -(p2[2] - p1[2])*(p3[0] - p2[0]) + (p2[0] - p1[0])*(p3[2] - p2[2]);
+            //                         n1[2] = -(p2[0] - p1[0])*(p3[1] - p2[1]) + (p2[1] - p1[1])*(p3[0] - p2[0]);
+            CalculateVectorNormal (p1, p2, p4, &n1[0], &n1[1], &n1[2]);
+            //CalculateVectorNormal (p2, p1, p3, &n2[0], &n2[1], &n2[2])
+            //CalculateVectorNormal (p3, p2, p4, &n1[0], &n1[1], &n1[2])
+            double dIntensity1, dIntensity2, dIntensity3, dIntensity4;
+            if (m_bColor) {
+              dIntensity1 =    (v[ix][iy] - dMin) / dIntensityScale;
+              dIntensity2 =    (v[ix+1][iy] - dMin) / dIntensityScale;
+              dIntensity3 =    (v[ix+1][iy+1] - dMin) / dIntensityScale;
+              dIntensity4 =    (v[ix][iy+1] - dMin) / dIntensityScale;
+            }
+            float vecColor[3];
+            if (m_bColor) {
+              intensityToColor (dIntensity1, vecColor);
+              glColor3fv (vecColor);
+            }
+            glVertex3fv (p1); glNormal3fv (n1);                                        
+            if (m_bColor) {
+              intensityToColor (dIntensity2, vecColor);
+              glColor3fv (vecColor);
+            }
+            glVertex3fv (p2); glNormal3fv (n1);                                        
+            if (m_bColor) {
+              intensityToColor (dIntensity3, vecColor);
+              glColor3fv (vecColor);
+            }
+            glVertex3fv (p3); glNormal3fv (n1);                                        
+            if (m_bColor) {
+              intensityToColor (dIntensity4, vecColor);
+              glColor3fv (vecColor);
+            }
+            glVertex3fv (p4); glNormal3fv (n1);                                                                                                                                                
+            glEnd();                           
+          }                    
+          
+        }
+        glEndList();
+        
 #ifdef GL_EXT_vertex_array
   }
 #endif
@@ -258,10 +344,43 @@ Graph3dFileView::DrawSurface()
 void
 Graph3dFileView::OnProperties (wxCommandEvent& event)
 {
-    std::ostringstream os;
-    *theApp->getLog() << ">>>>\n" << os.str().c_str() << "<<<<\n";
-    wxMessageDialog dialogMsg (getFrameForChild(), os.str().c_str(), "Imagefile Properties", wxOK | wxICON_INFORMATION);
-    dialogMsg.ShowModal();
+  std::ostringstream os;
+  *theApp->getLog() << ">>>>\n" << os.str().c_str() << "<<<<\n";
+  wxMessageDialog dialogMsg (getFrameForChild(), os.str().c_str(), "Imagefile Properties", wxOK | wxICON_INFORMATION);
+  dialogMsg.ShowModal();
+}
+
+void
+Graph3dFileView::OnLighting (wxCommandEvent& event)
+{
+  m_bLighting = ! m_bLighting;
+  m_pViewMenu->Check (GRAPH3D_VIEW_LIGHTING, m_bLighting);
+  
+  m_pCanvas->Refresh();
+}
+
+void
+Graph3dFileView::OnSurface (wxCommandEvent& event)
+{
+  m_bSurface = ! m_bSurface;
+  m_pViewMenu->Check (GRAPH3D_VIEW_SURFACE, m_bSurface);
+  m_pCanvas->Refresh();
+}
+
+void
+Graph3dFileView::OnColor (wxCommandEvent& event)
+{
+  m_bColor = ! m_bColor;
+  m_pViewMenu->Check (GRAPH3D_VIEW_COLOR, m_bColor);
+  m_pCanvas->Refresh();
+}
+
+void
+Graph3dFileView::OnSmooth (wxCommandEvent& event)
+{
+  m_bSmooth = ! m_bSmooth;
+  m_pViewMenu->Check (GRAPH3D_VIEW_SMOOTH, m_bSmooth);
+  m_pCanvas->Refresh();
 }
 
 
@@ -274,6 +393,9 @@ Graph3dFileView::OnDraw (wxDC* dc)
 #endif
   
   Draw();
+  std::ostringstream os;
+  os << "Xangle=" << m_dXRotate << ", Yangle=" << m_dYRotate << ", Zangle=" << m_dZRotate;
+  m_statusBar.SetStatusText (os.str().c_str());
   m_pCanvas->SwapBuffers();
 }
 
@@ -293,38 +415,134 @@ Graph3dFileView::Draw ()
 void
 Graph3dFileView::InitMaterials()
 {
+  if (! GetDocument())
+    return;
+  int nx = GetDocument()->nx();
+  int ny = GetDocument()->ny();
+  
+#if 1
   static float ambient[] = {0.1, 0.1, 0.1, 1.0};
-  static float diffuse[] = {0.5, 1.0, 1.0, 1.0};
-  static float position0[] = {0.0, 0.0, 20.0, 0.0};
-  static float position1[] = {0.0, 0.0, -20.0, 0.0};
-  static float front_mat_shininess[] = {60.0};
-  static float front_mat_specular[] = {0.2, 0.2, 0.2, 1.0};
-  static float front_mat_diffuse[] = {0.5, 0.28, 0.38, 1.0};
+  static float diffuse[] = {1.0, 1.0, 1.0, 1.0};
+  static float position0[] = {0, 0, -nx/4, 0, 0.0};
+  static float position1[] = {nx/2, ny/2, nx/4, 0.0};
+  //  static float position0[] = {0.0, 0.0, 20.0, 0.0};
+  //  static float position1[] = {0.0, 0.0, -20.0, 0.0};
+  static float front_mat_shininess[] = {5.0};
+  static float front_mat_specular[] = {0.1, 0.1, 0.1, 1.0};
+  static float front_mat_diffuse[] = {0.3, 0.3, 0.3, 1.0};
   /*
   static float back_mat_shininess[] = {60.0};
-  static float back_mat_specular[] = {0.5, 0.5, 0.2, 1.0};
-  static float back_mat_diffuse[] = {1.0, 1.0, 0.2, 1.0};
+  static float back_mat_specular[] = {0.2, 0.2, 0.2, 1.0};
+  static float back_mat_diffuse[] = {1.0, 1.0, 1.0, 1.0};
   */
   static float lmodel_ambient[] = {1.0, 1.0, 1.0, 1.0};
   static float lmodel_twoside[] = {GL_FALSE};
   
-  glLightfv(GL_LIGHT0, GL_AMBIENT, ambient);
-  glLightfv(GL_LIGHT0, GL_DIFFUSE, diffuse);
-  glLightfv(GL_LIGHT0, GL_POSITION, position0);
-  glEnable(GL_LIGHT0);
+  glLightfv (GL_LIGHT0, GL_AMBIENT, ambient);
+  glLightfv (GL_LIGHT0, GL_DIFFUSE, diffuse);
+  glLightfv (GL_LIGHT0, GL_POSITION, position0);
+  glEnable (GL_LIGHT0);
   
-  glLightfv(GL_LIGHT1, GL_AMBIENT, ambient);
-  glLightfv(GL_LIGHT1, GL_DIFFUSE, diffuse);
-  glLightfv(GL_LIGHT1, GL_POSITION, position1);
-  glEnable(GL_LIGHT1);
+  glLightfv (GL_LIGHT1, GL_AMBIENT, ambient);
+  glLightfv (GL_LIGHT1, GL_DIFFUSE, diffuse);
+  glLightfv (GL_LIGHT1, GL_POSITION, position1);
+  glEnable (GL_LIGHT1);
+  
+  glLightModelfv (GL_LIGHT_MODEL_AMBIENT, lmodel_ambient);
+  glLightModelfv (GL_LIGHT_MODEL_TWO_SIDE, lmodel_twoside);
+  
+  glMaterialfv (GL_FRONT_AND_BACK, GL_SHININESS, front_mat_shininess);
+  glMaterialfv (GL_FRONT_AND_BACK, GL_SPECULAR, front_mat_specular);
+  glMaterialfv (GL_FRONT_AND_BACK, GL_DIFFUSE, front_mat_diffuse);
+  
+  glColorMaterial (GL_FRONT_AND_BACK, GL_DIFFUSE);
+  //  glColorMaterial (GL_FRONT_AND_BACK, GL_SPECULAR);
+  glEnable(GL_COLOR_MATERIAL);
+#else
+  GLfloat impLPos[]  = {1.0, 1.0, 1.0, 0.0};
+  
+  GLfloat defaultLightAmb   [] = {.2, .2, .2, 1.0};
+  GLfloat defaultLightDiff  [] = {.2, .2, .2, 1.0};
+  GLfloat defaultLightSpec  [] = { .3, .3, .3, 1.0};
+  
+  GLfloat defaultGlobalAmb [] = {.3, .3, .3, 1.0};
+  GLfloat defaultGlobalDiff[] = {.3, .3, .3, 1.0};
+  
+  GLfloat defaultMatShine[] = {  30.0 };
+  GLfloat defaultMatSpec[]  = { .4, .4, .4, 1.0};
+  GLfloat defaultMatAmb[]   = { .3, .3, .3, 1.0};
+  GLfloat defaultMatDiff[]  = { .5, .5, .5, 1.0};
+  
+  GLfloat brassMatAmb[]   = { .33, .22, .03, 1.0};
+  GLfloat brassMatDiff[]  = { .78, .57, .11, 1.0};
+  GLfloat brassMatSpec[]  = { .99, .91, .81, 1.0};
+  GLfloat brassMatShine[] = {  27.8 };
+  
+  GLfloat emeraldMatAmb[]   = { .021, .1745 , .021, 1.0};
+  GLfloat emeraldMatDiff[]  = { .075, .6142 , .075, 1.0 };
+  GLfloat emeraldMatSpec[]  = { .633, .7278 , .633, 1.0 };
+  GLfloat emeraldMatShine[] = {  76.8 };
+  
+  GLfloat slateMatAmb[]   = { .02, .02 , .02, 1.0 };
+  GLfloat slateMatDiff[]  = { .02, .01 , .01, 1.0 };
+  GLfloat slateMatSpec[]  = { .4,  .4 ,  .4 , 1.0 };
+  GLfloat slateMatShine[] = { .768 };
+  
+  //       double opnX = nx, opnY = ny, opnZ = z;
+  //       eyeX = 1; eyeY = 1, eyeZ = 1;
+  
+  impLPos[0] = nx/2.; impLPos[1]= ny/2.; impLPos[2] = 0.;
+  //opnListNum = 1;
+  //impGraphicsFlag = IMP__3D;
+  
+  //       glutInitDisplayMode (GLUT_DOUBLE| GLUT_RGB | GLUT_DEPTH | GLUT_ACCUM);
+  //       glutInitWindowSize (IMP_WIN_X, IMP_WIN_Y);
+  //  glutInitWindowPosition (100, 100);
+  //       glutCreateWindow ("- imp3D graphics -" );
+  
+  glClearColor(0.0, 0.0, 0.0, 0.0);
+  
+  glShadeModel (GL_SMOOTH);
+  glBlendFunc(GL_SRC_ALPHA,GL_ONE_MINUS_SRC_ALPHA);
+  glHint(GL_LINE_SMOOTH, GL_DONT_CARE);
+  glEnable(GL_NORMALIZE);
+  
+  
+  glEnable(GL_DEPTH_TEST);
+  
+  glLightfv(GL_LIGHT0, GL_AMBIENT, defaultLightAmb);
+  glLightfv(GL_LIGHT0, GL_DIFFUSE, defaultLightDiff);
+  glLightfv(GL_LIGHT0, GL_SPECULAR,defaultLightSpec);
+  
+  glLightfv(GL_LIGHT1, GL_AMBIENT, defaultLightAmb);
+  glLightfv(GL_LIGHT1, GL_DIFFUSE, defaultLightDiff);
+  glLightfv(GL_LIGHT1, GL_SPECULAR,defaultLightSpec);
+  
+  glLightfv(GL_LIGHT2, GL_AMBIENT , defaultLightAmb);
+  glLightfv(GL_LIGHT2, GL_DIFFUSE , defaultLightDiff);
+  glLightfv(GL_LIGHT2, GL_SPECULAR, defaultLightSpec);
+  
+  glLightfv(GL_LIGHT0, GL_POSITION,impLPos);
+  glLightfv(GL_LIGHT1, GL_POSITION,impLPos);
+  glLightfv(GL_LIGHT2, GL_POSITION,impLPos);
+  
+  glMaterialfv(GL_FRONT_AND_BACK, GL_AMBIENT  , defaultMatAmb);
+  glMaterialfv(GL_FRONT_AND_BACK, GL_DIFFUSE  , defaultMatDiff);
+  glMaterialfv(GL_FRONT_AND_BACK, GL_SPECULAR , defaultMatSpec);
+  glMaterialfv(GL_FRONT_AND_BACK, GL_SHININESS, defaultMatShine);
+  
+  glLightModelfv(GL_LIGHT_MODEL_AMBIENT, defaultGlobalAmb);
+  
+  glColorMaterial(GL_FRONT_AND_BACK, GL_DIFFUSE);
+  
+  glEnable(GL_COLOR_MATERIAL);
   
-  glLightModelfv(GL_LIGHT_MODEL_AMBIENT, lmodel_ambient);
-  glLightModelfv(GL_LIGHT_MODEL_TWO_SIDE, lmodel_twoside);
   glEnable(GL_LIGHTING);
+  glEnable(GL_LIGHT1);
+  glEnable(GL_LIGHT2);
+  glEnable(GL_LIGHT0);
+#endif
   
-  glMaterialfv(GL_FRONT_AND_BACK, GL_SHININESS, front_mat_shininess);
-  glMaterialfv(GL_FRONT_AND_BACK, GL_SPECULAR, front_mat_specular);
-  glMaterialfv(GL_FRONT_AND_BACK, GL_DIFFUSE, front_mat_diffuse);
 }
 
 
@@ -333,17 +551,8 @@ Graph3dFileView::InitGL ()
 {
   glClearColor(0.0, 0.0, 0.0, 0.0);
   
-  glShadeModel(GL_SMOOTH);
   glEnable(GL_DEPTH_TEST);
   
-  InitMaterials();
-  
-  glMatrixMode(GL_PROJECTION);
-  glLoadIdentity();
-  glOrtho (-300, 300, -300, 300, 200, -200);
-  glMatrixMode(GL_MODELVIEW);
-  glLoadIdentity();
 }
 
 void 
@@ -352,34 +561,19 @@ Graph3dFileView::OnUpdate (wxView *WXUNUSED(sender), wxObject *WXUNUSED(hint) )
   int nVertices = GetDocument()->m_nVertices;
   glTripleFloat* pVertices = GetDocument()->m_pVertices;
   glTripleFloat* pNormals = GetDocument()->m_pNormals;
-
-#if 0    
-  const ImageFile& rIF = GetDocument()->getImageFile();
-  ImageFileArrayConst v = rIF.getArray();
-  int nx = rIF.nx();
-  int ny = rIF.ny();
-  if (v != NULL && nx != 0 && ny != 0) {
-    unsigned char* imageData = new unsigned char [nx * ny * 3];
-    for (int ix = 0; ix < nx; ix++) {
-      for (int iy = 0; iy < ny; iy++) {
-        double scaleValue = ((v[ix][iy] - m_dMinPixel) / scaleWidth) * 255;
-        int intensity = static_cast<int>(scaleValue + 0.5);
-        intensity = clamp (intensity, 0, 255);
-        int baseAddr = ((ny - 1 - iy) * nx + ix) * 3;
-        imageData[baseAddr] = imageData[baseAddr+1] = imageData[baseAddr+2] = intensity;
-      }
-    }
-    wxImage image (nx, ny, imageData, true);
-    m_bitmap = image.ConvertToBitmap();
-    delete imageData;
-    int xSize = nx;
-    int ySize = ny;
-    ySize = clamp (ySize, 0, 800);
-    m_pFrame->SetClientSize (xSize, ySize);
-    m_pCanvas->SetScrollbars(20, 20, nx/20, ny/20);
-    m_pCanvas->SetBackgroundColour(*wxWHITE);
-  } 
-#endif
+  
+  glMatrixMode(GL_PROJECTION);
+  glLoadIdentity();
+  if (! GetDocument())
+    return;
+  int nx = GetDocument()->nx();
+  int ny = GetDocument()->ny();
+  int maxDim = maxValue<int> (nx, ny);
+  
+  glOrtho (-maxDim * 0.71, maxDim * 0.71, -maxDim * 0.71, maxDim * 0.71, maxDim * 0.71, -maxDim * 0.71);
+  glMatrixMode(GL_MODELVIEW);
+  glLoadIdentity();
+  
   
 #ifdef GL_EXT_vertex_array
   if (m_bUseVertexArrays) {
@@ -389,6 +583,7 @@ Graph3dFileView::OnUpdate (wxView *WXUNUSED(sender), wxObject *WXUNUSED(hint) )
     glEnable( GL_NORMAL_ARRAY_EXT );
   }
 #endif
+  
   if (m_pCanvas)
     m_pCanvas->Refresh();
 }
@@ -434,13 +629,14 @@ Graph3dFileView::CreateChildFrame (wxDocument *doc, wxView *view)
 #endif
   theApp->setIconForFrame (subframe);
   
+  m_statusBar.Create (subframe, -1);
+  subframe->SetStatusBar (&m_statusBar);
+  
   m_pFileMenu = new wxMenu;
   
   m_pFileMenu->Append(MAINMENU_FILE_CREATE_PHANTOM, "Cr&eate Phantom...\tCtrl-P");
   m_pFileMenu->Append(MAINMENU_FILE_CREATE_FILTER, "Create &Filter...\tCtrl-F");
   m_pFileMenu->Append(wxID_OPEN, "&Open...\tCtrl-O");
-  m_pFileMenu->Append(wxID_SAVE, "&Save\tCtrl-S");
-  m_pFileMenu->Append(wxID_SAVEAS, "Save &As...");
   m_pFileMenu->Append(wxID_CLOSE, "&Close\tCtrl-W");
   
   m_pFileMenu->AppendSeparator();
@@ -457,6 +653,12 @@ Graph3dFileView::CreateChildFrame (wxDocument *doc, wxView *view)
   GetDocumentManager()->FileHistoryAddFilesToMenu(m_pFileMenu);
   GetDocumentManager()->FileHistoryUseMenu(m_pFileMenu);
   
+  m_pViewMenu = new wxMenu;
+  m_pViewMenu->Append(GRAPH3D_VIEW_SURFACE, "&Surface\tCtrl-U", "", true);
+  m_pViewMenu->Append(GRAPH3D_VIEW_SMOOTH, "&Smooth\tCtrl-M", "", true);
+  m_pViewMenu->Append(GRAPH3D_VIEW_COLOR, "&Color\tCtrl-R", "", true);
+  m_pViewMenu->Append(GRAPH3D_VIEW_LIGHTING, "&Lighting\tCtrl-L", "", true);
+  
   wxMenu *help_menu = new wxMenu;
   help_menu->Append(MAINMENU_HELP_CONTENTS, "&Contents\tF1");
   help_menu->Append(MAINMENU_HELP_TOPICS, "&Topics\tCtrl-H");
@@ -465,6 +667,7 @@ Graph3dFileView::CreateChildFrame (wxDocument *doc, wxView *view)
   wxMenuBar *menu_bar = new wxMenuBar;
   
   menu_bar->Append(m_pFileMenu, "&File");
+  menu_bar->Append(m_pViewMenu, "&View");
   menu_bar->Append(help_menu, "&Help");
   
   subframe->SetMenuBar(menu_bar);
@@ -473,16 +676,15 @@ Graph3dFileView::CreateChildFrame (wxDocument *doc, wxView *view)
   
   wxAcceleratorEntry accelEntries[10];
   accelEntries[0].Set (wxACCEL_CTRL, static_cast<int>('O'), wxID_OPEN);
-  accelEntries[1].Set (wxACCEL_CTRL, static_cast<int>('S'), wxID_SAVE);
-  accelEntries[2].Set (wxACCEL_CTRL, static_cast<int>('W'), wxID_CLOSE);
-  accelEntries[3].Set (wxACCEL_CTRL, static_cast<int>('H'), MAINMENU_HELP_TOPICS);
-  accelEntries[4].Set (wxACCEL_CTRL, static_cast<int>('P'), MAINMENU_FILE_CREATE_PHANTOM);
-  accelEntries[5].Set (wxACCEL_CTRL, static_cast<int>('F'), MAINMENU_FILE_CREATE_FILTER);
-  accelEntries[6].Set (wxACCEL_NORMAL, WXK_F1, MAINMENU_HELP_CONTENTS);
-  accelEntries[7].Set (wxACCEL_CTRL, static_cast<int>('A'), IFMENU_VIEW_SCALE_AUTO);
-  accelEntries[8].Set (wxACCEL_CTRL, static_cast<int>('U'), IFMENU_VIEW_SCALE_FULL);
-  accelEntries[9].Set (wxACCEL_CTRL, static_cast<int>('E'), IFMENU_VIEW_SCALE_MINMAX);
-  wxAcceleratorTable accelTable (10, accelEntries);
+  accelEntries[1].Set (wxACCEL_CTRL, static_cast<int>('H'), MAINMENU_HELP_TOPICS);
+  accelEntries[2].Set (wxACCEL_CTRL, static_cast<int>('P'), MAINMENU_FILE_CREATE_PHANTOM);
+  accelEntries[3].Set (wxACCEL_CTRL, static_cast<int>('F'), MAINMENU_FILE_CREATE_FILTER);
+  accelEntries[4].Set (wxACCEL_NORMAL, WXK_F1, MAINMENU_HELP_CONTENTS);
+  accelEntries[5].Set (wxACCEL_CTRL, static_cast<int>('U'), GRAPH3D_VIEW_SURFACE);
+  accelEntries[6].Set (wxACCEL_CTRL, static_cast<int>('R'), GRAPH3D_VIEW_COLOR);
+  accelEntries[7].Set (wxACCEL_CTRL, static_cast<int>('L'), GRAPH3D_VIEW_LIGHTING);
+  accelEntries[8].Set (wxACCEL_CTRL, static_cast<int>('M'), GRAPH3D_VIEW_SMOOTH);
+  wxAcceleratorTable accelTable (9, accelEntries);
   subframe->SetAcceleratorTable (accelTable);
   
   return subframe;
@@ -546,7 +748,8 @@ Graph3dFileCanvas::OnChar(wxKeyEvent& event)
   if (! m_pView)
     return;
   
-  switch(event.KeyCode()) {
+  wxCommandEvent dummyEvent;
+  switch (event.KeyCode()) {
   case WXK_LEFT:
        m_pView->m_dYRotate -= 15.0;
     break;
@@ -559,21 +762,20 @@ Graph3dFileCanvas::OnChar(wxKeyEvent& event)
   case WXK_DOWN:
     m_pView->m_dXRotate -= 15.0;
     break;
+  case 'z': case 'Z':
+    m_pView->m_dZRotate += 15.0;
+    break;
+  case 'x': case 'X':
+    m_pView->m_dZRotate -= 15.0;
+    break;
   case 's': case 'S':
-    m_pView->m_bSmooth = !m_pView->m_bSmooth;
-    if (m_pView->m_bSmooth) {
-      glShadeModel(GL_SMOOTH);
-    } else {
-      glShadeModel(GL_FLAT);
-    }
+    m_pView->OnSmooth (dummyEvent);
     break;
   case 'l': case 'L':
-    m_pView->m_bLighting = !m_pView->m_bLighting;
-    if (m_pView->m_bLighting) {
-      glEnable(GL_LIGHTING);
-    } else {
-      glDisable(GL_LIGHTING);
-    }
+    m_pView->OnLighting (dummyEvent);
+    break;
+  case 'c': case 'C':
+    m_pView->OnColor (dummyEvent);
     break;
   default:
     {
@@ -598,8 +800,11 @@ Graph3dFileCanvas::OnMouseEvent(wxMouseEvent& event)
   static int dragging = 0;
   static float last_x, last_y;
   
+  if (! m_pView)
+    return;
+  
   if(event.LeftIsDown()) {
-    if(!dragging) {
+    if(! dragging) {
       dragging = 1;
     } else {
       m_pView->m_dXRotate += (event.GetX() - last_x)*1.0;