r460: no message
[ctsim.git] / src / dlgreconstruct.cpp
index e9831843568d87a90af89727ac465c0ae928efc2..d0be191674e6315011cdd02a9058bc22c91c39c5 100644 (file)
@@ -7,9 +7,9 @@
 **   Date Started:  August 2000
 **
 **  This is part of the CTSim program
-**  Copyright (C) 1983-2000 Kevin Rosenberg
+**  Copyright (c) 1983-2001 Kevin Rosenberg
 **
-**  $Id: dlgreconstruct.cpp,v 1.5 2000/09/07 14:29:05 kevin Exp $
+**  $Id: dlgreconstruct.cpp,v 1.15 2001/01/28 19:10:18 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
 ******************************************************************************/
 
 #ifdef __GNUG__
-  #pragma implementation "dlgreconstruct.h"
+#pragma implementation "dlgreconstruct.h"
 #endif
 
 // For compilers that support precompilation, includes "wx.h".
 #include "wx/wxprec.h"
 
 #ifdef __BORLANDC__
-    #pragma hdrstop
+#pragma hdrstop
 #endif
 
-#ifndef WX_PRECOMP
-    #include "wx/utils.h"
-    #include "wx/frame.h"
-    #include "wx/button.h"
-    #include "wx/stattext.h"
-    #include "wx/layout.h"
-    #include "wx/event.h"
-    #include "wx/intl.h"
-    #include "wx/settings.h"
-    #include "wx/dcclient.h"
-    #include "wx/timer.h"
-    #include "wx/image.h"
-#endif
+#include "wx/utils.h"
+#include "wx/frame.h"
+#include "wx/button.h"
+#include "wx/stattext.h"
+#include "wx/layout.h"
+#include "wx/event.h"
+#include "wx/intl.h"
+#include "wx/settings.h"
+#include "wx/dcclient.h"
+#include "wx/timer.h"
+#include "wx/image.h"
 
 #include "dlgreconstruct.h"
 #include <algorithm>
 static const int LAYOUT_X_MARGIN = 4;
 static const int LAYOUT_Y_MARGIN = 4;
 
+
+const int ReconstructDialog::ID_BTN_PAUSE = 19998;
+const int ReconstructDialog::ID_BTN_STEP = 19999;
+const int ReconstructDialog::MAX_IMAGE_X = 400;
+const int ReconstructDialog::MAX_IMAGE_Y = 400;
+
+
 BEGIN_EVENT_TABLE(ReconstructDialog, wxDialog)
-   EVT_BUTTON(wxID_CANCEL, ReconstructDialog::OnCancel)
-   EVT_BUTTON(ID_BTN_PAUSE, ReconstructDialog::OnPause)
-   EVT_BUTTON(ID_BTN_STEP, ReconstructDialog::OnStep)
-   EVT_CLOSE(ReconstructDialog::OnClose)
-   EVT_PAINT(ReconstructDialog::OnPaint)
+EVT_BUTTON(wxID_CANCEL, ReconstructDialog::OnCancel)
+EVT_BUTTON(ID_BTN_PAUSE, ReconstructDialog::OnPause)
+EVT_BUTTON(ID_BTN_STEP, ReconstructDialog::OnStep)
+EVT_CLOSE(ReconstructDialog::OnClose)
+EVT_PAINT(ReconstructDialog::OnPaint)
 END_EVENT_TABLE()
 
 IMPLEMENT_CLASS(ReconstructDialog, wxDialog)
 
 
 ReconstructDialog::ReconstructDialog (Reconstructor& rReconstruct, const Projections& rProj, ImageFile& rIF, const int iTrace, wxWindow *parent)
-    : wxDialog(parent, -1, "Reconstruction", wxDefaultPosition), m_rReconstructor(rReconstruct), m_rProjections(rProj), m_rImageFile(rIF), m_pSGPDriver(NULL), m_pSGP(NULL), m_iTrace(iTrace), m_pDC(NULL), m_btnAbort(0), m_btnPause(0), m_btnStep(0)
+: wxDialog(parent, -1, "Reconstruction", wxDefaultPosition), m_rReconstructor(rReconstruct), m_rProjections(rProj), m_rImageFile(rIF), m_pSGPDriver(NULL), m_pSGP(NULL), m_iTrace(iTrace), m_pDC(NULL), m_btnAbort(0), m_btnPause(0), m_btnStep(0)
 {
     m_state = Continue;
     m_iLastView = -1;
@@ -82,13 +87,13 @@ ReconstructDialog::ReconstructDialog (Reconstructor& rReconstruct, const Project
     wxLayoutConstraints* c = new wxLayoutConstraints;
     c->right.SameAs(this, wxRight, 2*LAYOUT_X_MARGIN);
     c->bottom.SameAs(this, wxBottom, 2*LAYOUT_Y_MARGIN);
-
+       
     wxSize sizeBtn = wxButton::GetDefaultSize();
     c->width.Absolute(sizeBtn.x);
     c->height.Absolute(sizeBtn.y);
-
+       
     m_btnAbort->SetConstraints(c);
-
+       
     m_btnPause = new wxButton (this, ID_BTN_PAUSE, wxString("Pause"));
     wxLayoutConstraints* cPause = new wxLayoutConstraints;
     cPause->right.SameAs(this, wxRight, 3*LAYOUT_X_MARGIN + sizeBtn.x);
@@ -96,7 +101,7 @@ ReconstructDialog::ReconstructDialog (Reconstructor& rReconstruct, const Project
     cPause->width.Absolute(sizeBtn.x);
     cPause->height.Absolute(sizeBtn.y);
     m_btnPause->SetConstraints(cPause);
-
+       
     m_btnStep = new wxButton (this, ID_BTN_STEP, wxString("Step"));
     wxLayoutConstraints* cStep = new wxLayoutConstraints;
     cStep->right.SameAs(this, wxRight, 5*LAYOUT_X_MARGIN + sizeBtn.x * 2);
@@ -104,133 +109,126 @@ ReconstructDialog::ReconstructDialog (Reconstructor& rReconstruct, const Project
     cStep->width.Absolute(sizeBtn.x);
     cStep->height.Absolute(sizeBtn.y);
     m_btnStep->SetConstraints(cStep);
-
+       
     SetAutoLayout(TRUE);
     Layout();
-
+       
     m_nxGraph = 500;
     m_nyGraph = 500;
     wxSize sizeDlg (m_nxGraph, m_nyGraph);
     m_nxImage = m_rImageFile.nx();
     if (m_nxImage > MAX_IMAGE_X)
-      m_nxImage = MAX_IMAGE_X;
+               m_nxImage = MAX_IMAGE_X;
     m_nyImage = m_rImageFile.ny();
     if (m_nyImage > MAX_IMAGE_Y)
-      m_nyImage = MAX_IMAGE_Y;
-
+               m_nyImage = MAX_IMAGE_Y;
+       
     sizeDlg.x += m_nxImage;
     sizeDlg.y = max (sizeDlg.y, m_nyImage);
-
+       
     m_iClientX = sizeDlg.x;
     m_iClientY = sizeDlg.y;
     SetClientSize (sizeDlg);
-
-    m_bitmap.Create (m_iClientX, m_iClientY); // save a copy of screen
-    m_pDC = new wxMemoryDC;
-    m_pDC->SelectObject (m_bitmap);       // in memoryDC
-    m_pDC->SetFont (*wxSWISS_FONT);
-    int x, y;
-    this->GetClientSize(&x, &y);
-    m_pSGPDriver = new SGPDriver (dynamic_cast<wxDC*>(m_pDC), x, y);
-    m_pSGP = new SGP (*m_pSGPDriver);
-
+       
     Centre(wxCENTER_FRAME | wxBOTH);
-
+       
     if ( m_parentTop )
-      m_parentTop->Enable(FALSE);
-
+               m_parentTop->Enable(FALSE);
+       
     Show(TRUE);
     Enable(TRUE); // enable this window
-
+       
+    m_bitmap.Create (m_iClientX, m_iClientY); // save a copy of screen
+    m_pDC = dynamic_cast<wxDC*> (new wxClientDC (this));
+    int x, y;
+    this->GetClientSize(&x, &y);
+    m_pSGPDriver = new SGPDriver (m_pDC, x, y);
+    m_pSGP = new SGP (*m_pSGPDriver);
+       
     wxYield();     // Update the display
-
+       
     m_pSGPDriver->idWX()->SetFont(*wxSWISS_FONT);
-#ifdef __WXMAC__
-    MacUpdateImmediately();
-#endif
 }
 
 void
-ReconstructDialog::showView (int iViewNumber)
+ReconstructDialog::showView (int iViewNumber, bool bBackprojectView)
 {
-  if ( iViewNumber < m_rProjections.nView() ) {
-    m_iLastView = iViewNumber;
-    m_pSGP->eraseWindow();
-
-    char szProgress [256];
-    snprintf (szProgress, sizeof(szProgress), "Reconstructing View %d (%.1f%%)", iViewNumber, 100 * iViewNumber / static_cast<double>(m_rProjections.nView()));
-    m_pSGP->setViewport (0, 0, 1, 1);
-    m_pSGP->setWindow (0, 0, 1, 1);
-    m_pSGP->setTextColor (C_LTRED, -1);
-    double dCharHeight = m_pSGP->getCharHeight();
-    m_pSGP->setTextSize (dCharHeight * 2);
-    m_pSGP->moveAbs(0., m_pSGP->getCharHeight());
-    m_pSGP->drawText (szProgress);
-    m_pSGP->setTextSize (dCharHeight);
-
-    m_pSGP->setViewport (0.0, 0.1, 0.66, 1.);
-    m_rReconstructor.reconstructView (iViewNumber, 1, m_pSGP);
-
-    ImageFileArrayConst v = m_rImageFile.getArray();
-    int xBase = m_nxGraph;
-    int yBase = 0;
-    if (m_nyGraph > m_nyImage)
-       yBase = (m_nyGraph - m_nyImage) / 2;
-    double minValue = v[0][0];
-    double maxValue = v[0][0];
-    for (int ix = 0; ix < m_nxImage; ix++) {
-       for (int iy = 0; iy < m_nyImage; iy++) {
-           double dPixel = v[ix][iy];
-           if (dPixel < minValue)
-               minValue = dPixel;
-           else if (dPixel > maxValue)
-               maxValue = dPixel;
+       if ( iViewNumber < m_rProjections.nView() ) {
+               m_iLastView = iViewNumber;
+               ::wxYield();        // update the display
+               m_pSGP->eraseWindow();
+               m_btnPause->Refresh();
+               m_btnStep->Refresh();
+               m_btnAbort->Refresh();
+
+               char szProgress [256];
+               snprintf (szProgress, sizeof(szProgress), "Reconstructing View %d (%.1f%%)", iViewNumber, 100 * iViewNumber / static_cast<double>(m_rProjections.nView()));
+               m_pSGP->setViewport (0, 0, 1, 1);
+               m_pSGP->setWindow (0, 0, 1, 1);
+               m_pSGP->setTextColor (C_LTRED, -1);
+               double dCharHeight = m_pSGP->getCharHeight();
+               m_pSGP->setTextSize (dCharHeight * 2);
+               m_pSGP->moveAbs(0., m_pSGP->getCharHeight());
+               m_pSGP->drawText (szProgress);
+               m_pSGP->setTextSize (dCharHeight);
+               
+               // m_pSGP->setViewport (0.0, 0.1, 0.66, 1.); // viewport equiv now done in subroutine
+               m_rReconstructor.reconstructView (iViewNumber, 1, m_pSGP, bBackprojectView);
+               
+               ImageFileArrayConst v = m_rImageFile.getArray();
+               int xBase = m_nxGraph;
+               int yBase = 0;
+               if (m_nyGraph > m_nyImage)
+                       yBase = (m_nyGraph - m_nyImage) / 2;
+               double minValue = v[0][0];
+               double maxValue = v[0][0];
+               for (int ix = 0; ix < m_nxImage; ix++) {
+                       for (int iy = 0; iy < m_nyImage; iy++) {
+                               double dPixel = v[ix][iy];
+                               if (dPixel < minValue)
+                                       minValue = dPixel;
+                               else if (dPixel > maxValue)
+                                       maxValue = dPixel;
+                       }
+               }
+               unsigned char* imageData = new unsigned char [m_nxImage * m_nyImage * 3];
+               double dScale = 255 / (maxValue - minValue);
+               for (int ix2 = 0; ix2 < m_nxImage; ix2++) {
+                       for (int iy = 0; iy < m_nyImage; iy++) {
+                               double dPixel = v[ix2][iy];
+                               dPixel = (dPixel - minValue) * dScale;
+                               int intensity = nearest<int>(dPixel);
+                               intensity = clamp (intensity, 0, 255);
+                               int baseAddr = ((m_nyImage - 1 - iy) * m_nxImage + ix2) * 3;
+                               imageData[baseAddr] = imageData[baseAddr+1] = imageData[baseAddr+2] = intensity;
+                       }
+               }
+               wxImage image (m_nxImage, m_nyImage, imageData, true);
+               wxBitmap bitmap = image.ConvertToBitmap();
+               m_pSGP->getDriver().idWX()->DrawBitmap(bitmap, xBase, yBase, false);
+               delete imageData;
        }
-    }
-    unsigned char* imageData = new unsigned char [m_nxImage * m_nyImage * 3];
-    double dScale = 255 / (maxValue - minValue);
-    for (int ix = 0; ix < m_nxImage; ix++) {
-       for (int iy = 0; iy < m_nyImage; iy++) {
-           double dPixel = v[ix][iy];
-           dPixel = (dPixel - minValue) * dScale;
-           int intensity = nearest<int>(dPixel);
-           intensity = clamp (intensity, 0, 255);
-           int baseAddr = ((m_nyImage - 1 - iy) * m_nxImage + ix) * 3;
-           imageData[baseAddr] = imageData[baseAddr+1] = imageData[baseAddr+2] = intensity;
-       }
-    }
-    wxImage image (m_nxImage, m_nyImage, imageData, true);
-    wxBitmap bitmap = image.ConvertToBitmap();
-    m_pSGP->getDriver().idWX()->DrawBitmap(bitmap, xBase, yBase, false);
-    delete imageData;
-
-    Refresh();
-  }
 }
 
 bool
-ReconstructDialog::reconstructView (int iViewNumber)
+ReconstructDialog::reconstructView (int iViewNumber, bool bBackproject)
 {
-  if (iViewNumber <= m_iLastView)  // have already done this view
-    return true;
-
-  if (iViewNumber < m_rProjections.nView()) {
-    ::wxYield();        // update the display
-    showView (iViewNumber);
-    ::wxYield();        // update the display
-    if (m_iTrace >= Trace::TRACE_PLOT) {
-      ::wxUsleep(250);
-    }
-  } else {
-    m_state = Finished;    // so that we return TRUE below and 
-  }                        // that [Cancel] handler knew what to do
-  
-#ifdef __WXMAC__
-  MacUpdateImmediately();
-#endif
-
-  ::wxYield();        // update the display
-  return m_state != Cancelled;
+       if (iViewNumber <= m_iLastView)  // have already done this view
+               return true;
+       
+       if (iViewNumber < m_rProjections.nView()) {
+               ::wxYield();        // update the display
+               showView (iViewNumber, bBackproject);
+               ::wxYield();        // update the display
+               if (m_iTrace >= Trace::TRACE_PLOT) {
+                       ::wxUsleep(250);
+               }
+       } else {
+               m_state = Finished;    // so that we return TRUE below and 
+       }                        // that [Cancel] handler knew what to do
+       
+       ::wxYield();        // update the display
+       return m_state != Cancelled;
 }
 
 
@@ -238,70 +236,88 @@ ReconstructDialog::reconstructView (int iViewNumber)
 
 void ReconstructDialog::OnCancel (wxCommandEvent& event)
 {
-  if ( m_state == Finished ) {
-    // this means that the count down is already finished and we're being
-    // shown as a modal dialog - so just let the default handler do the job
-    event.Skip();
-  } else {
-    // request to cancel was received, the next time Update() is called we
-    // will handle it
-    m_state = Cancelled;
-
-    // update the button state immediately so that the user knows that the
-    // request has been noticed
-    m_btnAbort->Disable();
-  }
+       if ( m_state == Finished ) {
+               // this means that the count down is already finished and we're being
+               // shown as a modal dialog - so just let the default handler do the job
+               event.Skip();
+       } else {
+               // request to cancel was received, the next time Update() is called we
+               // will handle it
+               m_state = Cancelled;
+               
+               // update the button state immediately so that the user knows that the
+               // request has been noticed
+               m_btnAbort->Disable();
+       }
 }
 
 
 void 
 ReconstructDialog::OnPause (wxCommandEvent& event)
 {
-  if ( m_state == Finished ) {
-    // this means that the count down is already finished and we're being
-    // shown as a modal dialog - so just let the default handler do the job
-    event.Skip();
-  } else {
-      if (m_state == Continue) {
-       m_state = Paused;
-       m_btnPause->SetLabel (wxString("Resume"));
-      } else if (m_state == Paused) {
-       m_state = Continue;
-       m_btnPause->SetLabel (wxString("Pause"));
-      }
-  }
+       if ( m_state == Finished ) {
+               // this means that the count down is already finished and we're being
+               // shown as a modal dialog - so just let the default handler do the job
+               event.Skip();
+       } else if (m_state == Continue) {
+               m_memoryDC.SelectObject (m_bitmap);       // in memoryDC
+               m_pSGP->setDC (&m_memoryDC);
+               m_memoryDC.SetFont (*wxSWISS_FONT);
+               showView (m_iLastView, false);
+               m_state = Paused;
+               m_btnPause->SetLabel (wxString("Resume"));
+               m_pSGP->setDC (m_pDC);
+               m_memoryDC.SelectObject(wxNullBitmap);
+       } else if (m_state == Paused) {
+               m_state = Continue;
+               m_btnPause->SetLabel (wxString("Pause"));
+       }
 }
 
+
 void 
 ReconstructDialog::OnStep (wxCommandEvent& event)
 {
-  if ( m_state == Finished ) {
-    event.Skip();
-  } else {
-    if (m_state == Continue) {
-      m_state = Paused;
-      m_btnPause->SetLabel (wxString("Resume"));
-    } else if (m_state == Paused) {
-      reconstructView (m_iLastView + 1);
-    }
-  }
+       if ( m_state == Finished ) {
+               event.Skip();
+       } else if (m_state == Continue) {
+               m_memoryDC.SelectObject (m_bitmap);       // in memoryDC
+               m_pSGP->setDC (&m_memoryDC);
+               m_memoryDC.SetFont (*wxSWISS_FONT);
+               showView (m_iLastView, false);
+               m_state = Paused;
+               m_btnPause->SetLabel (wxString("Resume"));
+               m_pSGP->setDC (m_pDC);
+               m_memoryDC.SelectObject(wxNullBitmap);
+               Refresh();
+       } else if (m_state == Paused) {
+               m_memoryDC.SelectObject (m_bitmap);       // in memoryDC
+               m_pSGP->setDC (&m_memoryDC);
+               m_memoryDC.SetFont (*wxSWISS_FONT);
+               reconstructView (m_iLastView + 1);
+               m_pSGP->setDC (m_pDC);
+               m_memoryDC.SelectObject(wxNullBitmap);
+               Refresh();
+       }
 }
 
 void ReconstructDialog::OnClose(wxCloseEvent& event)
 {
     if ( m_state == Uncancellable )
-      event.Veto(TRUE);    // can't close this dialog
+               event.Veto(TRUE);    // can't close this dialog
     else if ( m_state == Finished )
-      event.Skip(); // let the default handler close the window as we already terminated
+               event.Skip(); // let the default handler close the window as we already terminated
     else
-      m_state = Cancelled;          // next Update() will notice it
+               m_state = Cancelled;          // next Update() will notice it
 }
 
 void
 ReconstructDialog::OnPaint (wxPaintEvent& event)
 {
-  wxPaintDC paintDC (this);
-  paintDC.DrawBitmap(m_bitmap, 0, 0, false);
+       wxPaintDC paintDC (this);
+       if (m_state == Paused) {
+               paintDC.DrawBitmap (m_bitmap, 0, 0, false);
+       }
 }
 
 
@@ -310,11 +326,11 @@ ReconstructDialog::OnPaint (wxPaintEvent& event)
 
 ReconstructDialog::~ReconstructDialog()
 {
-  if ( m_parentTop )
-    m_parentTop->Enable(TRUE);
-
-  delete m_pSGP;
-  delete m_pSGPDriver;
-  delete m_pDC;
+       if ( m_parentTop )
+               m_parentTop->Enable(TRUE);
+       
+       delete m_pSGP;
+       delete m_pSGPDriver;
+       delete m_pDC;
 }