1 /*****************************************************************************
5 ** Purpose: View & Canvas routines for CTSim program
6 ** Programmer: Kevin Rosenberg
7 ** Date Started: July 2000
9 ** This is part of the CTSim program
10 ** Copyright (C) 1983-2000 Kevin Rosenberg
12 ** $Id: views.cpp,v 1.19 2000/09/02 05:10:39 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
32 // For compilers that support precompilation, includes "wx/wx.h".
33 #include "wx/wxprec.h"
43 #if !wxUSE_DOC_VIEW_ARCHITECTURE
44 #error You must set wxUSE_DOC_VIEW_ARCHITECTURE to 1 in setup.h!
48 #include "wx/progdlg.h"
55 #include "dlgprojections.h"
56 #include "dlgreconstruct.h"
58 #include "backprojectors.h"
59 #include "reconstruct.h"
64 BEGIN_EVENT_TABLE(ImageFileCanvas, wxScrolledWindow)
65 EVT_MOUSE_EVENTS(ImageFileCanvas::OnMouseEvent)
69 ImageFileCanvas::ImageFileCanvas (ImageFileView* v, wxFrame *frame, const wxPoint& pos, const wxSize& size, const long style)
70 : wxScrolledWindow(frame, -1, pos, size, style)
76 ImageFileCanvas::OnDraw(wxDC& dc)
79 m_pView->OnDraw(& dc);
83 ImageFileCanvas::OnMouseEvent(wxMouseEvent& event)
91 wxPoint pt(event.GetLogicalPosition(dc));
93 if (event.LeftIsDown()) {
94 const ImageFile& rIF = m_pView->GetDocument()->getImageFile();
95 ImageFileArrayConst v = rIF.getArray();
99 if (pt.x >= 0 && pt.x < nx && pt.y >= 0 & pt.y < ny) {
101 os << "Image value (" << pt.x << "," << pt.y << ") = " << v[pt.x][ny - 1 - pt.y] << "\n";
102 *theApp->getLog() << os.str().c_str();
104 *theApp->getLog() << "Mouse out of image range (" << pt.x << "," << pt.y << ")\n";
112 IMPLEMENT_DYNAMIC_CLASS(ImageFileView, wxView)
114 BEGIN_EVENT_TABLE(ImageFileView, wxView)
115 EVT_MENU(IFMENU_FILE_PROPERTIES, ImageFileView::OnProperties)
116 EVT_MENU(IFMENU_VIEW_SCALE_MINMAX, ImageFileView::OnScaleMinMax)
117 EVT_MENU(IFMENU_VIEW_SCALE_AUTO, ImageFileView::OnScaleAuto)
120 ImageFileView::ImageFileView(void)
121 : wxView(), m_canvas(NULL), m_frame(NULL), m_bMinSpecified(false), m_bMaxSpecified(false)
125 ImageFileView::~ImageFileView(void)
130 ImageFileView::OnProperties (wxCommandEvent& event)
132 double min, max, mean, mode, median, stddev;
133 const ImageFile& rIF = GetDocument()->getImageFile();
134 if (rIF.nx() == 0 || rIF.ny() == 0)
135 *theApp->getLog() << "Properties: empty imagefile\n";
137 const string& rFilename = rIF.getFilename();
138 rIF.statistics (min, max, mean, mode, median, stddev);
140 os << "file: " << rFilename << "\nmin: "<<min<<"\nmax: "<<max<<"\nmean: "<<mean<<"\nmode: "<<mode<<"\nstddev: "<<stddev << "\n";
141 *theApp->getLog() << os.str().c_str();
142 wxMessageDialog dialogMsg (m_frame, os.str().c_str(), "Imagefile Properties", wxOK | wxICON_INFORMATION);
143 dialogMsg.ShowModal();
148 ImageFileView::OnScaleAuto (wxCommandEvent& event)
150 const ImageFile& rIF = GetDocument()->getImageFile();
151 DialogAutoScaleParameters dialogAutoScale (m_frame, rIF, m_dAutoScaleFactor);
152 int iRetVal = dialogAutoScale.ShowModal();
153 if (iRetVal == wxID_OK) {
154 m_bMinSpecified = true;
155 m_bMaxSpecified = true;
157 dialogAutoScale.getMinMax (&dMin, &dMax);
160 m_dAutoScaleFactor = dialogAutoScale.getAutoScaleFactor();
161 OnUpdate (this, NULL);
166 ImageFileView::OnScaleMinMax (wxCommandEvent& event)
168 const ImageFile& rIF = GetDocument()->getImageFile();
170 if (! m_bMinSpecified && ! m_bMaxSpecified)
171 rIF.getMinMax (min, max);
178 DialogGetImageMinMax dialogMinMax (m_frame, rIF, min, max);
179 int retVal = dialogMinMax.ShowModal();
180 if (retVal == wxID_OK) {
181 m_bMinSpecified = true;
182 m_bMaxSpecified = true;
183 m_dMinPixel = dialogMinMax.getMinimum();
184 m_dMaxPixel = dialogMinMax.getMaximum();
185 OnUpdate (this, NULL);
191 ImageFileView::CreateCanvas (wxView *view, wxFrame *parent)
193 ImageFileCanvas* pCanvas;
195 parent->GetClientSize(&width, &height);
197 pCanvas = new ImageFileCanvas (dynamic_cast<ImageFileView*>(view), parent, wxPoint(0, 0), wxSize(width, height), 0);
199 pCanvas->SetScrollbars(20, 20, 50, 50);
200 pCanvas->SetBackgroundColour(*wxWHITE);
207 ImageFileView::CreateChildFrame(wxDocument *doc, wxView *view)
209 wxDocChildFrame *subframe = new wxDocChildFrame(doc, view, theApp->getMainFrame(), -1, "ImageFile Frame", wxPoint(-1, -1), wxSize(0, 0), wxDEFAULT_FRAME_STYLE);
211 wxMenu *file_menu = new wxMenu;
213 file_menu->Append(MAINMENU_FILE_CREATE_PHANTOM, "Cr&eate Phantom...");
214 file_menu->Append(wxID_OPEN, "&Open...");
215 file_menu->Append(wxID_SAVE, "&Save");
216 file_menu->Append(wxID_SAVEAS, "Save &As...");
217 file_menu->Append(wxID_CLOSE, "&Close");
219 file_menu->AppendSeparator();
220 file_menu->Append(IFMENU_FILE_PROPERTIES, "P&roperties");
222 file_menu->AppendSeparator();
223 file_menu->Append(wxID_PRINT, "&Print...");
224 file_menu->Append(wxID_PRINT_SETUP, "Print &Setup...");
225 file_menu->Append(wxID_PREVIEW, "Print Pre&view");
227 wxMenu *view_menu = new wxMenu;
228 view_menu->Append(IFMENU_VIEW_SCALE_MINMAX, "Display Scale &Set...");
229 view_menu->Append(IFMENU_VIEW_SCALE_AUTO, "Display Scale &Auto...");
231 wxMenu *help_menu = new wxMenu;
232 help_menu->Append(MAINMENU_HELP_ABOUT, "&About");
234 wxMenuBar *menu_bar = new wxMenuBar;
236 menu_bar->Append(file_menu, "&File");
237 menu_bar->Append(view_menu, "&View");
238 menu_bar->Append(help_menu, "&Help");
240 subframe->SetMenuBar(menu_bar);
242 subframe->Centre(wxBOTH);
249 ImageFileView::OnCreate(wxDocument *doc, long WXUNUSED(flags) )
251 m_frame = CreateChildFrame(doc, this);
254 m_bMinSpecified = false;
255 m_bMaxSpecified = false;
256 m_dAutoScaleFactor = 1.;
259 m_frame->GetClientSize(&width, &height);
260 m_frame->SetTitle("ImageFileView");
261 m_canvas = CreateCanvas(this, m_frame);
264 int x, y; // X requires a forced resize
265 m_frame->GetSize(&x, &y);
266 m_frame->SetSize(-1, -1, x, y);
276 ImageFileView::OnDraw (wxDC* dc)
279 dc->DrawBitmap(m_bitmap, 0, 0, false);
284 ImageFileView::OnUpdate(wxView *WXUNUSED(sender), wxObject *WXUNUSED(hint) )
286 const ImageFile& rIF = dynamic_cast<ImageFileDocument*>(GetDocument())->getImageFile();
287 ImageFileArrayConst v = rIF.getArray();
290 if (v != NULL && nx != 0 && ny != 0) {
291 if (! m_bMinSpecified || ! m_bMaxSpecified) {
293 rIF.getMinMax (min, max);
294 if (! m_bMinSpecified)
296 if (! m_bMaxSpecified)
299 double scaleWidth = m_dMaxPixel - m_dMinPixel;
301 unsigned char* imageData = new unsigned char [nx * ny * 3];
302 for (int ix = 0; ix < nx; ix++) {
303 for (int iy = 0; iy < ny; iy++) {
304 double scaleValue = ((v[ix][iy] - m_dMinPixel) / scaleWidth) * 255;
305 int intensity = static_cast<int>(scaleValue + 0.5);
306 intensity = clamp (intensity, 0, 255);
307 int baseAddr = ((ny - 1 - iy) * nx + ix) * 3;
308 imageData[baseAddr] = imageData[baseAddr+1] = imageData[baseAddr+2] = intensity;
311 wxImage image (nx, ny, imageData, true);
312 m_bitmap = image.ConvertToBitmap();
316 xSize = clamp (xSize, 0, 800);
317 ySize = clamp (ySize, 0, 800);
318 m_frame->SetClientSize (xSize, ySize);
319 m_canvas->SetScrollbars(20, 20, nx/20, ny/20);
320 m_canvas->SetBackgroundColour(*wxWHITE);
328 ImageFileView::OnClose (bool deleteWindow)
330 if (!GetDocument()->Close())
334 m_canvas->m_pView = NULL;
336 wxString s(theApp->GetAppName());
338 m_frame->SetTitle(s);
354 PhantomCanvas::PhantomCanvas (PhantomView* v, wxFrame *frame, const wxPoint& pos, const wxSize& size, const long style)
355 : wxScrolledWindow(frame, -1, pos, size, style)
361 PhantomCanvas::OnDraw(wxDC& dc)
364 m_pView->OnDraw(& dc);
370 IMPLEMENT_DYNAMIC_CLASS(PhantomView, wxView)
372 BEGIN_EVENT_TABLE(PhantomView, wxView)
373 EVT_MENU(PHMMENU_FILE_PROPERTIES, PhantomView::OnProperties)
374 EVT_MENU(PHMMENU_PROCESS_RASTERIZE, PhantomView::OnRasterize)
375 EVT_MENU(PHMMENU_PROCESS_PROJECTIONS, PhantomView::OnProjections)
378 PhantomView::PhantomView(void)
379 : wxView(), m_canvas(NULL), m_frame(NULL)
381 m_iDefaultNDet = 367;
382 m_iDefaultNView = 320;
383 m_iDefaultNSample = 2;
384 m_dDefaultRotation = 2;
385 m_dDefaultFocalLength = 2;
386 m_dDefaultFieldOfView = 1;
387 m_iDefaultGeometry = Scanner::GEOMETRY_PARALLEL;
388 m_iDefaultTrace = Trace::TRACE_NONE;
391 PhantomView::~PhantomView(void)
396 PhantomView::OnProperties (wxCommandEvent& event)
398 const int idPhantom = GetDocument()->getPhantomID();
399 const string& namePhantom = GetDocument()->getPhantomName();
401 os << "Phantom " << namePhantom << " (" << idPhantom << ")\n";
402 *theApp->getLog() << os.str().c_str();
404 const Phantom& rPhantom = GetDocument()->getPhantom();
411 PhantomView::OnProjections (wxCommandEvent& event)
413 DialogGetProjectionParameters dialogProjection (m_frame, m_iDefaultNDet, m_iDefaultNView, m_iDefaultNSample, m_dDefaultRotation, m_dDefaultFocalLength, m_dDefaultFieldOfView, m_iDefaultGeometry, m_iDefaultTrace);
414 int retVal = dialogProjection.ShowModal();
415 if (retVal == wxID_OK) {
416 m_iDefaultNDet = dialogProjection.getNDet();
417 m_iDefaultNView = dialogProjection.getNView();
418 m_iDefaultNSample = dialogProjection.getNSamples();
419 m_iDefaultTrace = dialogProjection.getTrace();
420 m_dDefaultRotation = dialogProjection.getRotAngle();
421 m_dDefaultFocalLength = dialogProjection.getFocalLengthRatio();
422 m_dDefaultFieldOfView = dialogProjection.getFieldOfViewRatio();
423 wxString sGeometry = dialogProjection.getGeometry();
424 m_iDefaultGeometry = Scanner::convertGeometryNameToID (sGeometry.c_str());
426 if (m_iDefaultNDet > 0 && m_iDefaultNView > 0 && sGeometry != "") {
427 const Phantom& rPhantom = GetDocument()->getPhantom();
428 ProjectionFileDocument* pProjectionDoc = dynamic_cast<ProjectionFileDocument*>(theApp->getDocManager()->CreateDocument("untitled.pj", wxDOC_SILENT));
429 Projections& rProj = pProjectionDoc->getProjections();
430 Scanner theScanner (rPhantom, sGeometry.c_str(), m_iDefaultNDet, m_iDefaultNView, m_iDefaultNSample, m_dDefaultRotation, m_dDefaultFocalLength, m_dDefaultFieldOfView);
431 if (theScanner.fail()) {
432 *theApp->getLog() << "Failed making scanner: " << theScanner.failMessage().c_str() << "\n";
435 rProj.initFromScanner (theScanner);
436 m_dDefaultRotation /= PI; // convert back to PI units
438 if (m_iDefaultTrace > Trace::TRACE_CONSOLE) {
439 ProjectionsDialog dialogProjections (theScanner, rProj, rPhantom, m_iDefaultTrace, dynamic_cast<wxWindow*>(m_frame));
440 for (int iView = 0; iView < rProj.nView(); iView++) {
443 if (dialogProjections.isCancelled() || ! dialogProjections.projectView (iView)) {
444 pProjectionDoc->DeleteAllViews();
449 while (dialogProjections.isPaused()) {
455 wxProgressDialog dlgProgress (wxString("Projection"), wxString("Projection Progress"), rProj.nView() + 1, m_frame, wxPD_CAN_ABORT);
456 for (int i = 0; i < rProj.nView(); i++) {
457 theScanner.collectProjections (rProj, rPhantom, i, 1, true, m_iDefaultTrace);
458 if (! dlgProgress.Update (i+1)) {
459 pProjectionDoc->DeleteAllViews();
466 os << "Projections for " << rPhantom.name() << ": nDet=" << m_iDefaultNDet << ", nView=" << m_iDefaultNView << ", nSamples=" << m_iDefaultNSample << ", RotAngle=" << m_dDefaultRotation << ", FocalLengthRatio=" << m_dDefaultFocalLength << ", FieldOfViewRatio=" << m_dDefaultFieldOfView << ", Geometry=" << sGeometry.c_str() << "\n";
467 rProj.setRemark (os.str());
468 *theApp->getLog() << os.str().c_str();
471 pProjectionDoc->Modify(true);
472 pProjectionDoc->UpdateAllViews(this);
474 if (wxView* pView = pProjectionDoc->GetFirstView())
475 if (wxFrame* pFrame = pView->GetFrame()) {
486 PhantomView::OnRasterize (wxCommandEvent& event)
488 DialogGetRasterParameters dialogRaster (m_frame, 256, 256, 1);
489 int retVal = dialogRaster.ShowModal();
490 if (retVal == wxID_OK) {
491 int xSize = dialogRaster.getXSize();
492 int ySize = dialogRaster.getYSize();
493 int nSamples = dialogRaster.getNSamples();
496 if (xSize > 0 && ySize > 0) {
497 const Phantom& rPhantom = GetDocument()->getPhantom();
498 ImageFileDocument* pRasterDoc = dynamic_cast<ImageFileDocument*>(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT));
499 ImageFile& imageFile = pRasterDoc->getImageFile();
501 imageFile.setArraySize (xSize, ySize);
502 rPhantom.convertToImagefile (imageFile, nSamples, Trace::TRACE_NONE);
503 pRasterDoc->Modify(true);
504 pRasterDoc->UpdateAllViews(this);
507 os << "Rasterize Phantom " << rPhantom.name() << ": XSize=" << xSize << ", YSize=" << ySize << ", nSamples=" << nSamples << "\n";
508 *theApp->getLog() << os.str().c_str();
515 PhantomView::CreateCanvas (wxView *view, wxFrame *parent)
517 PhantomCanvas* pCanvas;
519 parent->GetClientSize(&width, &height);
521 pCanvas = new PhantomCanvas (dynamic_cast<PhantomView*>(view), parent, wxPoint(0, 0), wxSize(width, height), 0);
523 pCanvas->SetBackgroundColour(*wxWHITE);
530 PhantomView::CreateChildFrame(wxDocument *doc, wxView *view)
532 wxDocChildFrame *subframe = new wxDocChildFrame(doc, view, theApp->getMainFrame(), -1, "Phantom Frame", wxPoint(10, 10), wxSize(256, 256), wxDEFAULT_FRAME_STYLE);
534 wxMenu *file_menu = new wxMenu;
536 file_menu->Append(MAINMENU_FILE_CREATE_PHANTOM, "Cr&eate Phantom...");
537 file_menu->Append(wxID_OPEN, "&Open...");
538 file_menu->Append(wxID_CLOSE, "&Close");
540 file_menu->AppendSeparator();
541 file_menu->Append(PHMMENU_FILE_PROPERTIES, "P&roperties");
543 file_menu->AppendSeparator();
544 file_menu->Append(wxID_PRINT, "&Print...");
545 file_menu->Append(wxID_PRINT_SETUP, "Print &Setup...");
546 file_menu->Append(wxID_PREVIEW, "Print Pre&view");
548 wxMenu *process_menu = new wxMenu;
549 process_menu->Append(PHMMENU_PROCESS_RASTERIZE, "&Rasterize...");
550 process_menu->Append(PHMMENU_PROCESS_PROJECTIONS, "&Projections...");
552 wxMenu *help_menu = new wxMenu;
553 help_menu->Append(MAINMENU_HELP_CONTENTS, "&Contents");
554 help_menu->Append(MAINMENU_HELP_ABOUT, "&About");
556 wxMenuBar *menu_bar = new wxMenuBar;
558 menu_bar->Append(file_menu, "&File");
559 menu_bar->Append(process_menu, "&Process");
560 menu_bar->Append(help_menu, "&Help");
562 subframe->SetMenuBar(menu_bar);
564 subframe->Centre(wxBOTH);
571 PhantomView::OnCreate(wxDocument *doc, long WXUNUSED(flags) )
573 m_frame = CreateChildFrame(doc, this);
577 m_frame->GetClientSize(&width, &height);
578 m_frame->SetTitle("PhantomView");
579 m_canvas = CreateCanvas(this, m_frame);
582 int x, y; // X requires a forced resize
583 m_frame->GetSize(&x, &y);
584 m_frame->SetSize(-1, -1, x, y);
595 PhantomView::OnUpdate(wxView *WXUNUSED(sender), wxObject *WXUNUSED(hint) )
602 PhantomView::OnClose (bool deleteWindow)
604 if (!GetDocument()->Close())
608 m_canvas->m_pView = NULL;
610 wxString s(wxTheApp->GetAppName());
612 m_frame->SetTitle(s);
625 PhantomView::OnDraw (wxDC* dc)
628 m_canvas->GetClientSize (&xsize, &ysize);
629 SGPDriver driver (dc, xsize, ysize);
631 const Phantom& rPhantom = GetDocument()->getPhantom();
632 sgp.setColor (C_RED);
638 ProjectionFileCanvas::ProjectionFileCanvas (ProjectionFileView* v, wxFrame *frame, const wxPoint& pos, const wxSize& size, const long style)
639 : wxScrolledWindow(frame, -1, pos, size, style)
645 ProjectionFileCanvas::OnDraw(wxDC& dc)
648 m_pView->OnDraw(& dc);
651 // ProjectionFileView
653 IMPLEMENT_DYNAMIC_CLASS(ProjectionFileView, wxView)
655 BEGIN_EVENT_TABLE(ProjectionFileView, wxView)
656 EVT_MENU(PJMENU_FILE_PROPERTIES, ProjectionFileView::OnProperties)
657 EVT_MENU(PJMENU_PROCESS_RECONSTRUCT, ProjectionFileView::OnReconstruct)
660 ProjectionFileView::ProjectionFileView(void)
661 : wxView(), m_canvas(NULL), m_frame(NULL)
665 ProjectionFileView::~ProjectionFileView(void)
670 ProjectionFileView::OnProperties (wxCommandEvent& event)
672 const Projections& rProj = GetDocument()->getProjections();
674 rProj.printScanInfo(os);
675 *theApp->getLog() << os.str().c_str();
676 wxMessageDialog dialogMsg (m_frame, os.str().c_str(), "Projection File Properties", wxOK | wxICON_INFORMATION);
677 dialogMsg.ShowModal();
682 ProjectionFileView::OnReconstruct (wxCommandEvent& event)
685 DialogGetReconstructionParameters dialogReconstruction (m_frame, 256, 256, SignalFilter::FILTER_ABS_BANDLIMIT, 1., ProcessSignal::FILTER_METHOD_RFFTW, ProcessSignal::FILTER_GENERATION_INVERSE_FOURIER, 1, Backprojector::INTERP_LINEAR, 1, Backprojector::BPROJ_IDIFF3);
687 DialogGetReconstructionParameters dialogReconstruction (m_frame, 256, 256, SignalFilter::FILTER_ABS_BANDLIMIT, 1., ProcessSignal::FILTER_METHOD_CONVOLUTION, ProcessSignal::FILTER_GENERATION_DIRECT, 1, Backprojector::INTERP_LINEAR, 1, Backprojector::BPROJ_IDIFF3);
689 int retVal = dialogReconstruction.ShowModal();
690 if (retVal == wxID_OK) {
691 int xSize = dialogReconstruction.getXSize();
692 int ySize = dialogReconstruction.getYSize();
693 wxString optFilterName = dialogReconstruction.getFilterName();
694 double optFilterParam = dialogReconstruction.getFilterParam();
695 wxString optFilterMethodName = dialogReconstruction.getFilterMethodName();
696 int optZeropad = dialogReconstruction.getZeropad();
697 wxString optFilterGenerationName = dialogReconstruction.getFilterGenerationName();
698 wxString optInterpName = dialogReconstruction.getInterpName();
699 int optInterpParam = dialogReconstruction.getInterpParam();
700 wxString optBackprojectName = dialogReconstruction.getBackprojectName();
701 int iTrace = dialogReconstruction.getTrace();
702 if (xSize > 0 && ySize > 0) {
703 ImageFileDocument* pReconDoc = dynamic_cast<ImageFileDocument*>(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT));
704 ImageFile& imageFile = pReconDoc->getImageFile();
705 const Projections& rProj = GetDocument()->getProjections();
706 imageFile.setArraySize (xSize, ySize);
709 Reconstructor* pReconstruct = new Reconstructor (rProj, imageFile, optFilterName.c_str(), optFilterParam, optFilterMethodName.c_str(), optZeropad, optFilterGenerationName.c_str(), optInterpName.c_str(), optInterpParam, optBackprojectName.c_str(), iTrace);
710 if (iTrace > Trace::TRACE_CONSOLE) {
711 ReconstructDialog* pDlgReconstruct = new ReconstructDialog (*pReconstruct, rProj, imageFile, iTrace, m_frame);
712 for (int iView = 0; iView < rProj.nView(); iView++) {
715 if (pDlgReconstruct->isCancelled() || ! pDlgReconstruct->reconstructView (iView)) {
716 delete pDlgReconstruct;
718 pReconDoc->DeleteAllViews();
723 while (pDlgReconstruct->isPaused()) {
729 wxProgressDialog dlgProgress (wxString("Reconstruction"), wxString("Reconstruction Progress"), rProj.nView() + 1, m_frame, wxPD_CAN_ABORT);
730 for (int i = 0; i < rProj.nView(); i++) {
731 pReconstruct->reconstructView (i, 1);
732 if (! dlgProgress.Update(i + 1)) {
734 pReconDoc->DeleteAllViews();
739 pReconDoc->Modify(true);
740 pReconDoc->UpdateAllViews(this);
742 os << "Reconstruct " << rProj.getFilename() << ": xSize=" << xSize << ", ySize=" << ySize << ", Filter=" << optFilterName.c_str() << ", FilterParam=" << optFilterParam << ", FilterMethod=" << optFilterMethodName.c_str() << ", FilterGeneration=" << optFilterGenerationName.c_str() << ", Zeropad=" << optZeropad << ", Interpolation=" << optInterpName.c_str() << ", InterpolationParam=" << optInterpParam << ", Backprojection=" << optBackprojectName.c_str() << "\n";
743 *theApp->getLog() << os.str().c_str();
744 imageFile.labelAdd (rProj.getLabel());
745 imageFile.labelAdd (Array2dFileLabel::L_HISTORY, os.str().c_str(), timerRecon.timerEnd());
751 ProjectionFileCanvas*
752 ProjectionFileView::CreateCanvas (wxView *view, wxFrame *parent)
754 ProjectionFileCanvas* pCanvas;
756 parent->GetClientSize(&width, &height);
758 pCanvas = new ProjectionFileCanvas (dynamic_cast<ProjectionFileView*>(view), parent, wxPoint(0, 0), wxSize(width, height), 0);
760 pCanvas->SetScrollbars(20, 20, 50, 50);
761 pCanvas->SetBackgroundColour(*wxWHITE);
768 ProjectionFileView::CreateChildFrame(wxDocument *doc, wxView *view)
770 wxDocChildFrame *subframe = new wxDocChildFrame(doc, view, theApp->getMainFrame(), -1, "Projection Frame", wxPoint(10, 10), wxSize(0, 0), wxDEFAULT_FRAME_STYLE);
772 wxMenu *file_menu = new wxMenu;
774 file_menu->Append(MAINMENU_FILE_CREATE_PHANTOM, "Cr&eate Phantom...");
775 file_menu->Append(wxID_OPEN, "&Open...");
776 file_menu->Append(wxID_SAVE, "&Save");
777 file_menu->Append(wxID_SAVEAS, "Save &As...");
778 file_menu->Append(wxID_CLOSE, "&Close");
780 file_menu->AppendSeparator();
781 file_menu->Append(PJMENU_FILE_PROPERTIES, "P&roperties");
783 file_menu->AppendSeparator();
784 file_menu->Append(wxID_PRINT, "&Print...");
785 file_menu->Append(wxID_PRINT_SETUP, "Print &Setup...");
786 file_menu->Append(wxID_PREVIEW, "Print Pre&view");
788 wxMenu *process_menu = new wxMenu;
789 process_menu->Append(PJMENU_PROCESS_RECONSTRUCT, "R&econstruct...");
791 wxMenu *help_menu = new wxMenu;
792 help_menu->Append(MAINMENU_HELP_CONTENTS, "&Contents");
793 help_menu->AppendSeparator();
794 help_menu->Append(MAINMENU_HELP_ABOUT, "&About");
796 wxMenuBar *menu_bar = new wxMenuBar;
798 menu_bar->Append(file_menu, "&File");
799 menu_bar->Append(process_menu, "&Process");
800 menu_bar->Append(help_menu, "&Help");
802 subframe->SetMenuBar(menu_bar);
804 subframe->Centre(wxBOTH);
811 ProjectionFileView::OnCreate(wxDocument *doc, long WXUNUSED(flags) )
813 m_frame = CreateChildFrame(doc, this);
817 m_frame->GetClientSize(&width, &height);
818 m_frame->SetTitle("ProjectionFileView");
819 m_canvas = CreateCanvas(this, m_frame);
822 int x, y; // X requires a forced resize
823 m_frame->GetSize(&x, &y);
824 m_frame->SetSize(-1, -1, x, y);
834 ProjectionFileView::OnDraw (wxDC* dc)
837 dc->DrawBitmap (m_bitmap, 0, 0, false);
842 ProjectionFileView::OnUpdate(wxView *WXUNUSED(sender), wxObject *WXUNUSED(hint) )
844 const Projections& rProj = GetDocument()->getProjections();
845 const int nDet = rProj.nDet();
846 const int nView = rProj.nView();
847 if (nDet != 0 && nView != 0) {
848 const DetectorArray& detarray = rProj.getDetectorArray(0);
849 const DetectorValue* detval = detarray.detValues();
850 double min = detval[0];
851 double max = detval[0];
852 for (int iy = 0; iy < nView; iy++) {
853 const DetectorArray& detarray = rProj.getDetectorArray(iy);
854 const DetectorValue* detval = detarray.detValues();
855 for (int ix = 0; ix < nDet; ix++) {
856 if (min > detval[ix])
858 else if (max < detval[ix])
863 unsigned char* imageData = new unsigned char [nDet * nView * 3];
864 double scale = (max - min) / 255;
865 for (int iy = 0; iy < nView; iy++) {
866 const DetectorArray& detarray = rProj.getDetectorArray(iy);
867 const DetectorValue* detval = detarray.detValues();
868 for (int ix = 0; ix < nDet; ix++) {
869 int intensity = static_cast<int>(((detval[ix] - min) / scale) + 0.5);
870 intensity = clamp(intensity, 0, 255);
871 int baseAddr = (iy * nDet + ix) * 3;
872 imageData[baseAddr] = imageData[baseAddr+1] = imageData[baseAddr+2] = intensity;
875 wxImage image (nDet, nView, imageData, true);
876 m_bitmap = image.ConvertToBitmap();
880 xSize = clamp (xSize, 0, 800);
881 ySize = clamp (ySize, 0, 800);
882 m_frame->SetClientSize (xSize, ySize);
883 m_canvas->SetScrollbars (20, 20, nDet/20, nView/20);
891 ProjectionFileView::OnClose (bool deleteWindow)
893 if (!GetDocument()->Close())
897 m_canvas->m_pView = NULL;
899 wxString s(wxTheApp->GetAppName());
901 m_frame->SetTitle(s);