1 /*****************************************************************************
4 ** Name: dlgreconstruct.cpp
5 ** Purpose: Reconstruction Animation Dialog
6 ** Programmer: Kevin Rosenberg
7 ** Date Started: August 2000
9 ** This is part of the CTSim program
10 ** Copyright (C) 1983-2000 Kevin Rosenberg
12 ** $Id: dlgreconstruct.cpp,v 1.11 2000/12/18 06:32:13 kevin Exp $
14 ** This program is free software; you can redistribute it and/or modify
15 ** it under the terms of the GNU General Public License (version 2) as
16 ** published by the Free Software Foundation.
18 ** This program is distributed in the hope that it will be useful,
19 ** but WITHOUT ANY WARRANTY; without even the implied warranty of
20 ** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
21 ** GNU General Public License for more details.
23 ** You should have received a copy of the GNU General Public License
24 ** along with this program; if not, write to the Free Software
25 ** Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
26 ******************************************************************************/
29 #pragma implementation "dlgreconstruct.h"
32 // For compilers that support precompilation, includes "wx.h".
33 #include "wx/wxprec.h"
41 #include "wx/button.h"
42 #include "wx/stattext.h"
43 #include "wx/layout.h"
46 #include "wx/settings.h"
47 #include "wx/dcclient.h"
51 #include "dlgreconstruct.h"
56 static const int LAYOUT_X_MARGIN = 4;
57 static const int LAYOUT_Y_MARGIN = 4;
60 const int ReconstructDialog::ID_BTN_PAUSE = 19998;
\r
61 const int ReconstructDialog::ID_BTN_STEP = 19999;
\r
62 const int ReconstructDialog::MAX_IMAGE_X = 400;
\r
63 const int ReconstructDialog::MAX_IMAGE_Y = 400;
\r
66 BEGIN_EVENT_TABLE(ReconstructDialog, wxDialog)
67 EVT_BUTTON(wxID_CANCEL, ReconstructDialog::OnCancel)
68 EVT_BUTTON(ID_BTN_PAUSE, ReconstructDialog::OnPause)
69 EVT_BUTTON(ID_BTN_STEP, ReconstructDialog::OnStep)
70 EVT_CLOSE(ReconstructDialog::OnClose)
71 EVT_PAINT(ReconstructDialog::OnPaint)
74 IMPLEMENT_CLASS(ReconstructDialog, wxDialog)
77 ReconstructDialog::ReconstructDialog (Reconstructor& rReconstruct, const Projections& rProj, ImageFile& rIF, const int iTrace, wxWindow *parent)
78 : 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)
83 while ( m_parentTop && m_parentTop->GetParent() )
84 m_parentTop = m_parentTop->GetParent();
86 m_btnAbort = new wxButton(this, wxID_CANCEL, _("Cancel"));
87 wxLayoutConstraints* c = new wxLayoutConstraints;
88 c->right.SameAs(this, wxRight, 2*LAYOUT_X_MARGIN);
89 c->bottom.SameAs(this, wxBottom, 2*LAYOUT_Y_MARGIN);
91 wxSize sizeBtn = wxButton::GetDefaultSize();
92 c->width.Absolute(sizeBtn.x);
93 c->height.Absolute(sizeBtn.y);
95 m_btnAbort->SetConstraints(c);
97 m_btnPause = new wxButton (this, ID_BTN_PAUSE, wxString("Pause"));
98 wxLayoutConstraints* cPause = new wxLayoutConstraints;
99 cPause->right.SameAs(this, wxRight, 3*LAYOUT_X_MARGIN + sizeBtn.x);
100 cPause->bottom.SameAs(this, wxBottom, 2*LAYOUT_Y_MARGIN);
101 cPause->width.Absolute(sizeBtn.x);
102 cPause->height.Absolute(sizeBtn.y);
103 m_btnPause->SetConstraints(cPause);
105 m_btnStep = new wxButton (this, ID_BTN_STEP, wxString("Step"));
106 wxLayoutConstraints* cStep = new wxLayoutConstraints;
107 cStep->right.SameAs(this, wxRight, 5*LAYOUT_X_MARGIN + sizeBtn.x * 2);
108 cStep->bottom.SameAs(this, wxBottom, 2*LAYOUT_Y_MARGIN);
109 cStep->width.Absolute(sizeBtn.x);
110 cStep->height.Absolute(sizeBtn.y);
111 m_btnStep->SetConstraints(cStep);
118 wxSize sizeDlg (m_nxGraph, m_nyGraph);
119 m_nxImage = m_rImageFile.nx();
120 if (m_nxImage > MAX_IMAGE_X)
121 m_nxImage = MAX_IMAGE_X;
122 m_nyImage = m_rImageFile.ny();
123 if (m_nyImage > MAX_IMAGE_Y)
124 m_nyImage = MAX_IMAGE_Y;
126 sizeDlg.x += m_nxImage;
127 sizeDlg.y = max (sizeDlg.y, m_nyImage);
129 m_iClientX = sizeDlg.x;
130 m_iClientY = sizeDlg.y;
131 SetClientSize (sizeDlg);
133 Centre(wxCENTER_FRAME | wxBOTH);
\r
136 m_parentTop->Enable(FALSE);
\r
139 Enable(TRUE); // enable this window
\r
141 m_bitmap.Create (m_iClientX, m_iClientY); // save a copy of screen
142 m_pDC = dynamic_cast<wxDC*> (new wxClientDC (this));
\r
144 this->GetClientSize(&x, &y);
\r
145 m_pSGPDriver = new SGPDriver (m_pDC, x, y);
\r
146 m_pSGP = new SGP (*m_pSGPDriver);
\r
148 wxYield(); // Update the display
150 m_pSGPDriver->idWX()->SetFont(*wxSWISS_FONT);
154 ReconstructDialog::showView (int iViewNumber, bool bBackprojectView)
156 if ( iViewNumber < m_rProjections.nView() ) {
157 m_iLastView = iViewNumber;
158 ::wxYield(); // update the display
\r
159 m_pSGP->eraseWindow();
\r
162 char szProgress [256];
163 snprintf (szProgress, sizeof(szProgress), "Reconstructing View %d (%.1f%%)", iViewNumber, 100 * iViewNumber / static_cast<double>(m_rProjections.nView()));
164 m_pSGP->setViewport (0, 0, 1, 1);
165 m_pSGP->setWindow (0, 0, 1, 1);
166 m_pSGP->setTextColor (C_LTRED, -1);
167 double dCharHeight = m_pSGP->getCharHeight();
168 m_pSGP->setTextSize (dCharHeight * 2);
169 m_pSGP->moveAbs(0., m_pSGP->getCharHeight());
170 m_pSGP->drawText (szProgress);
171 m_pSGP->setTextSize (dCharHeight);
173 m_pSGP->setViewport (0.0, 0.1, 0.66, 1.);
174 m_rReconstructor.reconstructView (iViewNumber, 1, m_pSGP, bBackprojectView);
176 ImageFileArrayConst v = m_rImageFile.getArray();
177 int xBase = m_nxGraph;
179 if (m_nyGraph > m_nyImage)
180 yBase = (m_nyGraph - m_nyImage) / 2;
181 double minValue = v[0][0];
182 double maxValue = v[0][0];
183 for (int ix = 0; ix < m_nxImage; ix++) {
184 for (int iy = 0; iy < m_nyImage; iy++) {
185 double dPixel = v[ix][iy];
186 if (dPixel < minValue)
188 else if (dPixel > maxValue)
192 unsigned char* imageData = new unsigned char [m_nxImage * m_nyImage * 3];
193 double dScale = 255 / (maxValue - minValue);
194 for (int ix2 = 0; ix2 < m_nxImage; ix2++) {
195 for (int iy = 0; iy < m_nyImage; iy++) {
196 double dPixel = v[ix2][iy];
197 dPixel = (dPixel - minValue) * dScale;
198 int intensity = nearest<int>(dPixel);
199 intensity = clamp (intensity, 0, 255);
200 int baseAddr = ((m_nyImage - 1 - iy) * m_nxImage + ix2) * 3;
201 imageData[baseAddr] = imageData[baseAddr+1] = imageData[baseAddr+2] = intensity;
204 wxImage image (m_nxImage, m_nyImage, imageData, true);
205 wxBitmap bitmap = image.ConvertToBitmap();
206 m_pSGP->getDriver().idWX()->DrawBitmap(bitmap, xBase, yBase, false);
212 ReconstructDialog::reconstructView (int iViewNumber, bool bBackproject)
214 if (iViewNumber <= m_iLastView) // have already done this view
217 if (iViewNumber < m_rProjections.nView()) {
218 ::wxYield(); // update the display
219 showView (iViewNumber, bBackproject);
220 ::wxYield(); // update the display
221 if (m_iTrace >= Trace::TRACE_PLOT) {
225 m_state = Finished; // so that we return TRUE below and
226 } // that [Cancel] handler knew what to do
228 ::wxYield(); // update the display
229 return m_state != Cancelled;
235 void ReconstructDialog::OnCancel (wxCommandEvent& event)
237 if ( m_state == Finished ) {
238 // this means that the count down is already finished and we're being
239 // shown as a modal dialog - so just let the default handler do the job
242 // request to cancel was received, the next time Update() is called we
246 // update the button state immediately so that the user knows that the
247 // request has been noticed
248 m_btnAbort->Disable();
254 ReconstructDialog::OnPause (wxCommandEvent& event)
256 if ( m_state == Finished ) {
257 // this means that the count down is already finished and we're being
258 // shown as a modal dialog - so just let the default handler do the job
260 } else if (m_state == Continue) {
\r
261 m_memoryDC.SelectObject (m_bitmap); // in memoryDC
\r
262 m_pSGP->setDC (&m_memoryDC);
\r
263 m_memoryDC.SetFont (*wxSWISS_FONT);
\r
264 showView (m_iLastView, false);
\r
266 m_btnPause->SetLabel (wxString("Resume"));
267 m_pSGP->setDC (m_pDC);
\r
268 m_memoryDC.SelectObject(wxNullBitmap);
\r
269 } else if (m_state == Paused) {
271 m_btnPause->SetLabel (wxString("Pause"));
277 ReconstructDialog::OnStep (wxCommandEvent& event)
279 if ( m_state == Finished ) {
281 } else if (m_state == Continue) {
\r
282 m_memoryDC.SelectObject (m_bitmap); // in memoryDC
\r
283 m_pSGP->setDC (&m_memoryDC);
\r
284 m_memoryDC.SetFont (*wxSWISS_FONT);
\r
285 showView (m_iLastView, false);
\r
287 m_btnPause->SetLabel (wxString("Resume"));
288 m_pSGP->setDC (m_pDC);
\r
289 m_memoryDC.SelectObject(wxNullBitmap);
\r
291 } else if (m_state == Paused) {
292 m_memoryDC.SelectObject (m_bitmap); // in memoryDC
\r
293 m_pSGP->setDC (&m_memoryDC);
\r
294 m_memoryDC.SetFont (*wxSWISS_FONT);
\r
295 reconstructView (m_iLastView + 1);
296 m_pSGP->setDC (m_pDC);
\r
297 m_memoryDC.SelectObject(wxNullBitmap);
\r
302 void ReconstructDialog::OnClose(wxCloseEvent& event)
304 if ( m_state == Uncancellable )
305 event.Veto(TRUE); // can't close this dialog
306 else if ( m_state == Finished )
307 event.Skip(); // let the default handler close the window as we already terminated
309 m_state = Cancelled; // next Update() will notice it
313 ReconstructDialog::OnPaint (wxPaintEvent& event)
315 wxPaintDC paintDC (this);
\r
316 if (m_state == Paused) {
317 paintDC.DrawBitmap (m_bitmap, 0, 0, false);
\r
322 /////////////////////////////////////////////////////
325 ReconstructDialog::~ReconstructDialog()
328 m_parentTop->Enable(TRUE);