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.31 2000/12/19 21:37:51 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 ******************************************************************************/
28 // For compilers that support precompilation, includes "wx/wx.h".
29 #include "wx/wxprec.h"
39 #if !wxUSE_DOC_VIEW_ARCHITECTURE
40 #error You must set wxUSE_DOC_VIEW_ARCHITECTURE to 1 in setup.h!
44 #include "wx/progdlg.h"
51 #include "dlgprojections.h"
52 #include "dlgreconstruct.h"
53 #include "backprojectors.h"
54 #include "reconstruct.h"
57 #if defined(MSVC) || HAVE_SSTREAM
\r
60 #include <sstream_subst>
\r
66 BEGIN_EVENT_TABLE(ImageFileCanvas, wxScrolledWindow)
67 EVT_MOUSE_EVENTS(ImageFileCanvas::OnMouseEvent)
71 ImageFileCanvas::ImageFileCanvas (ImageFileView* v, wxFrame *frame, const wxPoint& pos, const wxSize& size, const long style)
72 : wxScrolledWindow(frame, -1, pos, size, style)
80 ImageFileCanvas::OnDraw(wxDC& dc)
83 m_pView->OnDraw(& dc);
87 ImageFileCanvas::DrawRubberBandCursor (wxDC& dc, int x, int y)
\r
89 const ImageFile& rIF = m_pView->GetDocument()->getImageFile();
\r
90 ImageFileArrayConst v = rIF.getArray();
\r
94 dc.SetLogicalFunction (wxINVERT);
\r
95 dc.SetPen (*wxGREEN_PEN);
\r
96 dc.DrawLine (0, y, nx, y);
\r
97 dc.DrawLine (x, 0, x, ny);
\r
98 dc.SetLogicalFunction (wxCOPY);
\r
102 ImageFileCanvas::GetCurrentCursor (int& x, int& y)
\r
107 if (m_xCursor >= 0 && m_yCursor >= 0)
\r
114 ImageFileCanvas::OnMouseEvent(wxMouseEvent& event)
122 wxPoint pt(event.GetLogicalPosition(dc));
124 if (event.RightIsDown()) {
\r
125 const ImageFile& rIF = m_pView->GetDocument()->getImageFile();
\r
126 ImageFileArrayConst v = rIF.getArray();
\r
130 if (pt.x >= 0 && pt.x < nx && pt.y >= 0 && pt.y < ny) {
\r
131 std::ostringstream os;
\r
132 os << "Image value (" << pt.x << "," << pt.y << ") = " << v[pt.x][ny - 1 - pt.y] << "\n";
\r
133 *theApp->getLog() << os.str().c_str();
\r
135 *theApp->getLog() << "Mouse out of image range (" << pt.x << "," << pt.y << ")\n";
\r
137 else if (event.LeftIsDown()) {
\r
138 const ImageFile& rIF = m_pView->GetDocument()->getImageFile();
\r
139 ImageFileArrayConst v = rIF.getArray();
\r
143 if (pt.x >= 0 && pt.x < nx && pt.y >= 0 && pt.y < ny) {
\r
144 if (m_xCursor >= 0 && m_yCursor >= 0) {
\r
145 DrawRubberBandCursor (dc, m_xCursor, m_yCursor);
\r
147 DrawRubberBandCursor (dc, pt.x, pt.y);
\r
150 std::ostringstream os;
\r
151 os << "Selected column" << pt.x << " and row" << pt.y << "\n";
\r
152 *theApp->getLog() << os.str().c_str();
\r
154 *theApp->getLog() << "Mouse out of image range (" << pt.x << "," << pt.y << ")\n";
\r
160 IMPLEMENT_DYNAMIC_CLASS(ImageFileView, wxView)
162 BEGIN_EVENT_TABLE(ImageFileView, wxView)
163 EVT_MENU(IFMENU_FILE_PROPERTIES, ImageFileView::OnProperties)
164 EVT_MENU(IFMENU_VIEW_SCALE_MINMAX, ImageFileView::OnScaleMinMax)
165 EVT_MENU(IFMENU_VIEW_SCALE_AUTO, ImageFileView::OnScaleAuto)
\r
166 EVT_MENU(IFMENU_PLOT_ROW, ImageFileView::OnPlotRow)
167 EVT_MENU(IFMENU_PLOT_COL, ImageFileView::OnPlotCol)
\r
170 ImageFileView::ImageFileView(void)
171 : wxView(), m_canvas(NULL), m_frame(NULL), m_bMinSpecified(false), m_bMaxSpecified(false)
175 ImageFileView::~ImageFileView(void)
180 ImageFileView::OnProperties (wxCommandEvent& event)
182 double min, max, mean, mode, median, stddev;
183 const ImageFile& rIF = GetDocument()->getImageFile();
184 if (rIF.nx() == 0 || rIF.ny() == 0)
185 *theApp->getLog() << "Properties: empty imagefile\n";
187 const std::string& rFilename = rIF.getFilename();
188 rIF.statistics (min, max, mean, mode, median, stddev);
189 std::ostringstream os;
190 os << "file: " << rFilename << "\nmin: "<<min<<"\nmax: "<<max<<"\nmean: "<<mean<<"\nmedian: "<<median<<"\nmode: "<<mode<<"\nstddev: "<<stddev << "\n";
191 *theApp->getLog() << os.str().c_str();
192 wxMessageDialog dialogMsg (m_frame, os.str().c_str(), "Imagefile Properties", wxOK | wxICON_INFORMATION);
193 dialogMsg.ShowModal();
198 ImageFileView::OnScaleAuto (wxCommandEvent& event)
200 const ImageFile& rIF = GetDocument()->getImageFile();
201 DialogAutoScaleParameters dialogAutoScale (m_frame, rIF, m_dAutoScaleFactor);
202 int iRetVal = dialogAutoScale.ShowModal();
203 if (iRetVal == wxID_OK) {
204 m_bMinSpecified = true;
205 m_bMaxSpecified = true;
207 dialogAutoScale.getMinMax (&dMin, &dMax);
210 m_dAutoScaleFactor = dialogAutoScale.getAutoScaleFactor();
211 OnUpdate (this, NULL);
216 ImageFileView::OnScaleMinMax (wxCommandEvent& event)
218 const ImageFile& rIF = GetDocument()->getImageFile();
220 if (! m_bMinSpecified && ! m_bMaxSpecified)
221 rIF.getMinMax (min, max);
228 DialogGetImageMinMax dialogMinMax (m_frame, rIF, min, max);
229 int retVal = dialogMinMax.ShowModal();
230 if (retVal == wxID_OK) {
231 m_bMinSpecified = true;
232 m_bMaxSpecified = true;
233 m_dMinPixel = dialogMinMax.getMinimum();
234 m_dMaxPixel = dialogMinMax.getMaximum();
235 OnUpdate (this, NULL);
241 ImageFileView::CreateCanvas (wxView *view, wxFrame *parent)
243 ImageFileCanvas* pCanvas;
245 parent->GetClientSize(&width, &height);
247 pCanvas = new ImageFileCanvas (dynamic_cast<ImageFileView*>(view), parent, wxPoint(0, 0), wxSize(width, height), 0);
249 pCanvas->SetScrollbars(20, 20, 50, 50);
250 pCanvas->SetBackgroundColour(*wxWHITE);
257 ImageFileView::CreateChildFrame(wxDocument *doc, wxView *view)
259 wxDocChildFrame *subframe = new wxDocChildFrame(doc, view, theApp->getMainFrame(), -1, "ImageFile Frame", wxPoint(-1, -1), wxSize(0, 0), wxDEFAULT_FRAME_STYLE);
261 wxMenu *file_menu = new wxMenu;
263 file_menu->Append(MAINMENU_FILE_CREATE_PHANTOM, "Cr&eate Phantom...");
264 file_menu->Append(wxID_OPEN, "&Open...");
265 file_menu->Append(wxID_SAVE, "&Save");
266 file_menu->Append(wxID_SAVEAS, "Save &As...");
267 file_menu->Append(wxID_CLOSE, "&Close");
269 file_menu->AppendSeparator();
270 file_menu->Append(IFMENU_FILE_PROPERTIES, "P&roperties");
272 file_menu->AppendSeparator();
273 file_menu->Append(wxID_PRINT, "&Print...");
274 file_menu->Append(wxID_PRINT_SETUP, "Print &Setup...");
275 file_menu->Append(wxID_PREVIEW, "Print Pre&view");
277 wxMenu *view_menu = new wxMenu;
278 view_menu->Append(IFMENU_VIEW_SCALE_MINMAX, "Display Scale &Set...");
279 view_menu->Append(IFMENU_VIEW_SCALE_AUTO, "Display Scale &Auto...");
281 wxMenu *plot_menu = new wxMenu;
\r
282 plot_menu->Append (IFMENU_PLOT_ROW, "Plot &Row");
\r
283 plot_menu->Append (IFMENU_PLOT_COL, "Plot &Column");
\r
285 wxMenu *help_menu = new wxMenu;
286 help_menu->Append(MAINMENU_HELP_ABOUT, "&About");
288 wxMenuBar *menu_bar = new wxMenuBar;
290 menu_bar->Append(file_menu, "&File");
291 menu_bar->Append(view_menu, "&View");
\r
292 menu_bar->Append(plot_menu, "&Plot");
293 menu_bar->Append(help_menu, "&Help");
295 subframe->SetMenuBar(menu_bar);
297 subframe->Centre(wxBOTH);
304 ImageFileView::OnCreate(wxDocument *doc, long WXUNUSED(flags) )
306 m_frame = CreateChildFrame(doc, this);
309 m_bMinSpecified = false;
310 m_bMaxSpecified = false;
311 m_dAutoScaleFactor = 1.;
314 m_frame->GetClientSize (&width, &height);
315 m_frame->SetTitle("ImageFileView");
316 m_canvas = CreateCanvas (this, m_frame);
318 int x, y; // X requires a forced resize
319 m_frame->GetSize(&x, &y);
320 m_frame->SetSize(-1, -1, x, y);
321 m_frame->SetFocus();
\r
329 ImageFileView::OnDraw (wxDC* dc)
332 dc->DrawBitmap(m_bitmap, 0, 0, false);
334 int xCursor, yCursor;
\r
335 if (m_canvas->GetCurrentCursor (xCursor, yCursor))
\r
336 m_canvas->DrawRubberBandCursor (*dc, xCursor, yCursor);
\r
341 ImageFileView::OnUpdate (wxView *WXUNUSED(sender), wxObject *WXUNUSED(hint) )
343 const ImageFile& rIF = dynamic_cast<ImageFileDocument*>(GetDocument())->getImageFile();
344 ImageFileArrayConst v = rIF.getArray();
347 if (v != NULL && nx != 0 && ny != 0) {
348 if (! m_bMinSpecified || ! m_bMaxSpecified) {
350 rIF.getMinMax (min, max);
351 if (! m_bMinSpecified)
353 if (! m_bMaxSpecified)
356 double scaleWidth = m_dMaxPixel - m_dMinPixel;
358 unsigned char* imageData = new unsigned char [nx * ny * 3];
359 for (int ix = 0; ix < nx; ix++) {
360 for (int iy = 0; iy < ny; iy++) {
361 double scaleValue = ((v[ix][iy] - m_dMinPixel) / scaleWidth) * 255;
362 int intensity = static_cast<int>(scaleValue + 0.5);
363 intensity = clamp (intensity, 0, 255);
364 int baseAddr = ((ny - 1 - iy) * nx + ix) * 3;
365 imageData[baseAddr] = imageData[baseAddr+1] = imageData[baseAddr+2] = intensity;
368 wxImage image (nx, ny, imageData, true);
369 m_bitmap = image.ConvertToBitmap();
373 xSize = clamp (xSize, 0, 800);
374 ySize = clamp (ySize, 0, 800);
375 m_frame->SetClientSize (xSize, ySize);
376 m_canvas->SetScrollbars(20, 20, nx/20, ny/20);
377 m_canvas->SetBackgroundColour(*wxWHITE);
\r
385 ImageFileView::OnClose (bool deleteWindow)
387 if (!GetDocument()->Close())
391 m_canvas->m_pView = NULL;
393 wxString s(theApp->GetAppName());
395 m_frame->SetTitle(s);
407 #include "wx/plot.h"
\r
409 class VectorPlotCurve : wxPlotCurve {
\r
412 std::vector<double> m_points;
\r
415 VectorPlotCurve (const double* const points, int n, double minVal, double maxVal)
\r
416 : wxPlotCurve (minVal, minVal, maxVal)
\r
418 for (int i = 0; i < n; i++)
\r
419 m_points.push_back(points[i]);
\r
422 wxInt32 GetStartX ()
\r
426 { return m_nPoints - 1; }
\r
428 double GetY (wxInt32 x)
\r
429 { return m_points[x]; }
\r
434 ImageFileView::OnPlotRow (wxCommandEvent& event)
\r
436 int xCursor, yCursor;
\r
437 if (! m_canvas->GetCurrentCursor (xCursor, yCursor)) {
\r
438 // os << "No row selected. Please use left mouse button on image to select row\n";
\r
442 const ImageFile& rIF = dynamic_cast<ImageFileDocument*>(GetDocument())->getImageFile();
\r
443 ImageFileArrayConst v = rIF.getArray();
\r
447 if (v != NULL && yCursor < ny) {
\r
448 double* pVec = new double [nx];
\r
449 double minVal = v[0][yCursor];
\r
450 double maxVal = minVal;
\r
451 for (int i = 0; i < nx; i++) {
\r
452 pVec[i] = v[i][yCursor];
\r
453 if (pVec[i] < minVal)
\r
455 else if (pVec[i] > maxVal)
\r
458 VectorPlotCurve* plotCurve = new VectorPlotCurve (pVec, nx, minVal, maxVal);
\r
459 wxPlotWindow* plotWindow = new wxPlotWindow (m_canvas, NULL, wxPoint(0,0), wxSize(500,300));
\r
460 plotWindow->Add (reinterpret_cast<wxPlotCurve*>(plotCurve));
\r
461 plotWindow->Enable(TRUE);
\r
462 plotWindow->Show();
\r
468 ImageFileView::OnPlotCol (wxCommandEvent& event)
\r
470 int xCursor, yCursor;
\r
471 if (! m_canvas->GetCurrentCursor (xCursor, yCursor)) {
\r
472 // os << "No column selected. Please use left mouse button on image to select column\n";
\r
476 const ImageFile& rIF = dynamic_cast<ImageFileDocument*>(GetDocument())->getImageFile();
\r
477 ImageFileArrayConst v = rIF.getArray();
\r
481 if (v != NULL && xCursor < nx) {
\r
482 double* pVec = new double [ny];
\r
483 double minVal = v[xCursor][0];
\r
484 double maxVal = minVal;
\r
485 for (int i = 0; i < ny; i++) {
\r
486 if (minVal < v[xCursor][i])
\r
487 minVal = v[xCursor][i];
\r
488 else if (maxVal > v[xCursor][i])
\r
489 maxVal = v[xCursor][i];
\r
490 pVec[i] = v[xCursor][i];
\r
493 VectorPlotCurve* plotCurve = new VectorPlotCurve (pVec, ny, minVal, maxVal);
\r
494 wxPlotWindow plotWindow (m_canvas, NULL, wxPoint(0,0), wxSize(500,300));
\r
495 plotWindow.Add (reinterpret_cast<wxPlotCurve*>(plotCurve));
\r
502 PhantomCanvas::PhantomCanvas (PhantomView* v, wxFrame *frame, const wxPoint& pos, const wxSize& size, const long style)
503 : wxScrolledWindow(frame, -1, pos, size, style)
509 PhantomCanvas::OnDraw(wxDC& dc)
512 m_pView->OnDraw(& dc);
518 IMPLEMENT_DYNAMIC_CLASS(PhantomView, wxView)
520 BEGIN_EVENT_TABLE(PhantomView, wxView)
521 EVT_MENU(PHMMENU_FILE_PROPERTIES, PhantomView::OnProperties)
522 EVT_MENU(PHMMENU_PROCESS_RASTERIZE, PhantomView::OnRasterize)
523 EVT_MENU(PHMMENU_PROCESS_PROJECTIONS, PhantomView::OnProjections)
526 PhantomView::PhantomView(void)
527 : wxView(), m_canvas(NULL), m_frame(NULL)
529 m_iDefaultNDet = 367;
530 m_iDefaultNView = 320;
531 m_iDefaultNSample = 2;
532 m_dDefaultRotation = 2;
533 m_dDefaultFocalLength = 2;
534 m_dDefaultFieldOfView = 1;
535 m_iDefaultGeometry = Scanner::GEOMETRY_PARALLEL;
536 m_iDefaultTrace = Trace::TRACE_NONE;
539 PhantomView::~PhantomView(void)
544 PhantomView::OnProperties (wxCommandEvent& event)
546 const int idPhantom = GetDocument()->getPhantomID();
547 const wxString& namePhantom = GetDocument()->getPhantomName();
548 std::ostringstream os;
549 os << "Phantom " << namePhantom.c_str() << " (" << idPhantom << ")\n";
550 *theApp->getLog() << os.str().c_str();
552 const Phantom& rPhantom = GetDocument()->getPhantom();
559 PhantomView::OnProjections (wxCommandEvent& event)
561 DialogGetProjectionParameters dialogProjection (m_frame, m_iDefaultNDet, m_iDefaultNView, m_iDefaultNSample, m_dDefaultRotation, m_dDefaultFocalLength, m_dDefaultFieldOfView, m_iDefaultGeometry, m_iDefaultTrace);
562 int retVal = dialogProjection.ShowModal();
563 if (retVal == wxID_OK) {
564 m_iDefaultNDet = dialogProjection.getNDet();
565 m_iDefaultNView = dialogProjection.getNView();
566 m_iDefaultNSample = dialogProjection.getNSamples();
567 m_iDefaultTrace = dialogProjection.getTrace();
568 m_dDefaultRotation = dialogProjection.getRotAngle();
569 m_dDefaultFocalLength = dialogProjection.getFocalLengthRatio();
570 m_dDefaultFieldOfView = dialogProjection.getFieldOfViewRatio();
571 wxString sGeometry = dialogProjection.getGeometry();
572 m_iDefaultGeometry = Scanner::convertGeometryNameToID (sGeometry.c_str());
574 if (m_iDefaultNDet > 0 && m_iDefaultNView > 0 && sGeometry != "") {
575 const Phantom& rPhantom = GetDocument()->getPhantom();
576 ProjectionFileDocument* pProjectionDoc = dynamic_cast<ProjectionFileDocument*>(theApp->getDocManager()->CreateDocument("untitled.pj", wxDOC_SILENT));
577 Projections& rProj = pProjectionDoc->getProjections();
578 Scanner theScanner (rPhantom, sGeometry.c_str(), m_iDefaultNDet, m_iDefaultNView, m_iDefaultNSample, m_dDefaultRotation, m_dDefaultFocalLength, m_dDefaultFieldOfView);
579 if (theScanner.fail()) {
580 *theApp->getLog() << "Failed making scanner: " << theScanner.failMessage().c_str() << "\n";
583 rProj.initFromScanner (theScanner);
584 m_dDefaultRotation /= PI; // convert back to PI units
586 if (m_iDefaultTrace > Trace::TRACE_CONSOLE) {
587 ProjectionsDialog dialogProjections (theScanner, rProj, rPhantom, m_iDefaultTrace, dynamic_cast<wxWindow*>(m_frame));
588 for (int iView = 0; iView < rProj.nView(); iView++) {
591 if (dialogProjections.isCancelled() || ! dialogProjections.projectView (iView)) {
592 pProjectionDoc->DeleteAllViews();
597 while (dialogProjections.isPaused()) {
603 wxProgressDialog dlgProgress (wxString("Projection"), wxString("Projection Progress"), rProj.nView() + 1, m_frame, wxPD_CAN_ABORT);
604 for (int i = 0; i < rProj.nView(); i++) {
605 theScanner.collectProjections (rProj, rPhantom, i, 1, true, m_iDefaultTrace);
606 if (! dlgProgress.Update (i+1)) {
607 pProjectionDoc->DeleteAllViews();
613 std::ostringstream os;
614 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();
615 rProj.setRemark (os.str());
616 *theApp->getLog() << os.str().c_str() << "\n";
620 ProjectionFileView* projView = dynamic_cast<ProjectionFileView*>(pProjectionDoc->GetFirstView());
\r
622 projView->getFrame()->SetFocus();
\r
623 projView->OnUpdate (projView, NULL);
\r
625 if (wxView* pView = pProjectionDoc->GetFirstView()) {
626 if (wxFrame* pFrame = pView->GetFrame()) {
630 theApp->getDocManager()->ActivateView (pView, true, false);
633 pProjectionDoc->Modify(true);
634 pProjectionDoc->UpdateAllViews(this);
641 PhantomView::OnRasterize (wxCommandEvent& event)
643 DialogGetRasterParameters dialogRaster (m_frame, 256, 256, 1);
644 int retVal = dialogRaster.ShowModal();
645 if (retVal == wxID_OK) {
646 int xSize = dialogRaster.getXSize();
647 int ySize = dialogRaster.getYSize();
648 int nSamples = dialogRaster.getNSamples();
651 if (xSize > 0 && ySize > 0) {
652 const Phantom& rPhantom = GetDocument()->getPhantom();
653 ImageFileDocument* pRasterDoc = dynamic_cast<ImageFileDocument*>(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT));
654 ImageFile& imageFile = pRasterDoc->getImageFile();
656 imageFile.setArraySize (xSize, ySize);
657 wxProgressDialog dlgProgress (wxString("Rasterize"), wxString("Rasterization Progress"), imageFile.nx() + 1, m_frame, wxPD_CAN_ABORT);
658 for (unsigned int i = 0; i < imageFile.nx(); i++) {
659 rPhantom.convertToImagefile (imageFile, nSamples, Trace::TRACE_NONE, i, 1, true);
660 if (! dlgProgress.Update(i+1)) {
661 pRasterDoc->DeleteAllViews();
665 pRasterDoc->Modify(true);
666 pRasterDoc->UpdateAllViews(this);
667 ImageFileView* rasterView = dynamic_cast<ImageFileView*>(pRasterDoc->GetFirstView());
\r
669 rasterView->getFrame()->SetFocus();
\r
670 rasterView->OnUpdate (rasterView, NULL);
\r
673 std::ostringstream os;
674 os << "Rasterize Phantom " << rPhantom.name() << ": XSize=" << xSize << ", YSize=" << ySize << ", nSamples=" << nSamples << "\n";
675 *theApp->getLog() << os.str().c_str();
682 PhantomView::CreateCanvas (wxView *view, wxFrame *parent)
684 PhantomCanvas* pCanvas;
686 parent->GetClientSize(&width, &height);
688 pCanvas = new PhantomCanvas (dynamic_cast<PhantomView*>(view), parent, wxPoint(0, 0), wxSize(width, height), 0);
690 pCanvas->SetBackgroundColour(*wxWHITE);
697 PhantomView::CreateChildFrame(wxDocument *doc, wxView *view)
699 wxDocChildFrame *subframe = new wxDocChildFrame(doc, view, theApp->getMainFrame(), -1, "Phantom Frame", wxPoint(10, 10), wxSize(256, 256), wxDEFAULT_FRAME_STYLE);
701 wxMenu *file_menu = new wxMenu;
703 file_menu->Append(MAINMENU_FILE_CREATE_PHANTOM, "Cr&eate Phantom...");
704 file_menu->Append(wxID_OPEN, "&Open...");
705 file_menu->Append(wxID_CLOSE, "&Close");
707 file_menu->AppendSeparator();
708 file_menu->Append(PHMMENU_FILE_PROPERTIES, "P&roperties");
710 file_menu->AppendSeparator();
711 file_menu->Append(wxID_PRINT, "&Print...");
712 file_menu->Append(wxID_PRINT_SETUP, "Print &Setup...");
713 file_menu->Append(wxID_PREVIEW, "Print Pre&view");
715 wxMenu *process_menu = new wxMenu;
716 process_menu->Append(PHMMENU_PROCESS_RASTERIZE, "&Rasterize...");
717 process_menu->Append(PHMMENU_PROCESS_PROJECTIONS, "&Projections...");
719 wxMenu *help_menu = new wxMenu;
720 help_menu->Append(MAINMENU_HELP_CONTENTS, "&Contents");
721 help_menu->Append(MAINMENU_HELP_ABOUT, "&About");
723 wxMenuBar *menu_bar = new wxMenuBar;
725 menu_bar->Append(file_menu, "&File");
726 menu_bar->Append(process_menu, "&Process");
727 menu_bar->Append(help_menu, "&Help");
729 subframe->SetMenuBar(menu_bar);
731 subframe->Centre(wxBOTH);
738 PhantomView::OnCreate(wxDocument *doc, long WXUNUSED(flags) )
740 m_frame = CreateChildFrame(doc, this);
744 m_frame->GetClientSize(&width, &height);
745 m_frame->SetTitle("PhantomView");
746 m_canvas = CreateCanvas(this, m_frame);
749 int x, y; // X requires a forced resize
750 m_frame->GetSize(&x, &y);
751 m_frame->SetSize(-1, -1, x, y);
762 PhantomView::OnUpdate(wxView *WXUNUSED(sender), wxObject *WXUNUSED(hint) )
769 PhantomView::OnClose (bool deleteWindow)
771 if (!GetDocument()->Close())
775 m_canvas->m_pView = NULL;
777 wxString s(wxTheApp->GetAppName());
779 m_frame->SetTitle(s);
792 PhantomView::OnDraw (wxDC* dc)
795 m_canvas->GetClientSize (&xsize, &ysize);
796 SGPDriver driver (dc, xsize, ysize);
798 const Phantom& rPhantom = GetDocument()->getPhantom();
799 sgp.setColor (C_RED);
805 ProjectionFileCanvas::ProjectionFileCanvas (ProjectionFileView* v, wxFrame *frame, const wxPoint& pos, const wxSize& size, const long style)
806 : wxScrolledWindow(frame, -1, pos, size, style)
812 ProjectionFileCanvas::OnDraw(wxDC& dc)
815 m_pView->OnDraw(& dc);
818 // ProjectionFileView
820 IMPLEMENT_DYNAMIC_CLASS(ProjectionFileView, wxView)
822 BEGIN_EVENT_TABLE(ProjectionFileView, wxView)
823 EVT_MENU(PJMENU_FILE_PROPERTIES, ProjectionFileView::OnProperties)
824 EVT_MENU(PJMENU_PROCESS_RECONSTRUCT, ProjectionFileView::OnReconstruct)
827 ProjectionFileView::ProjectionFileView(void)
828 : wxView(), m_canvas(NULL), m_frame(NULL)
832 m_iDefaultFilter = SignalFilter::FILTER_ABS_BANDLIMIT;
833 m_dDefaultFilterParam = 1.;
835 m_iDefaultFilterMethod = ProcessSignal::FILTER_METHOD_RFFTW;
836 m_iDefaultFilterGeneration = ProcessSignal::FILTER_GENERATION_INVERSE_FOURIER;
838 m_iDefaultFilterMethod = ProcessSignal::FILTER_METHOD_CONVOLUTION;
839 m_iDefaultFilterGeneration = ProcessSignal::FILTER_GENERATION_DIRECT;
841 m_iDefaultZeropad = 1;
842 m_iDefaultBackprojector = Backprojector::BPROJ_IDIFF3;
843 m_iDefaultInterpolation = Backprojector::INTERP_LINEAR;
844 m_iDefaultInterpParam = 1;
845 m_iDefaultTrace = Trace::TRACE_NONE;
848 ProjectionFileView::~ProjectionFileView(void)
853 ProjectionFileView::OnProperties (wxCommandEvent& event)
855 const Projections& rProj = GetDocument()->getProjections();
856 std::ostringstream os;
857 rProj.printScanInfo(os);
858 *theApp->getLog() << os.str().c_str();
859 wxMessageDialog dialogMsg (m_frame, os.str().c_str(), "Projection File Properties", wxOK | wxICON_INFORMATION);
860 dialogMsg.ShowModal();
865 ProjectionFileView::OnReconstruct (wxCommandEvent& event)
867 DialogGetReconstructionParameters dialogReconstruction (m_frame, m_iDefaultNX, m_iDefaultNY, m_iDefaultFilter, m_dDefaultFilterParam, m_iDefaultFilterMethod, m_iDefaultFilterGeneration, m_iDefaultZeropad, m_iDefaultInterpolation, m_iDefaultInterpParam, m_iDefaultBackprojector, m_iDefaultTrace);
869 int retVal = dialogReconstruction.ShowModal();
870 if (retVal == wxID_OK) {
871 m_iDefaultNX = dialogReconstruction.getXSize();
872 m_iDefaultNY = dialogReconstruction.getYSize();
873 wxString optFilterName = dialogReconstruction.getFilterName();
874 m_iDefaultFilter = SignalFilter::convertFilterNameToID (optFilterName.c_str());
875 m_dDefaultFilterParam = dialogReconstruction.getFilterParam();
876 wxString optFilterMethodName = dialogReconstruction.getFilterMethodName();
877 m_iDefaultFilterMethod = ProcessSignal::convertFilterMethodNameToID(optFilterMethodName.c_str());
878 m_iDefaultZeropad = dialogReconstruction.getZeropad();
879 wxString optFilterGenerationName = dialogReconstruction.getFilterGenerationName();
880 m_iDefaultFilterGeneration = ProcessSignal::convertFilterGenerationNameToID (optFilterGenerationName.c_str());
881 wxString optInterpName = dialogReconstruction.getInterpName();
882 m_iDefaultInterpolation = Backprojector::convertInterpNameToID (optInterpName.c_str());
883 m_iDefaultInterpParam = dialogReconstruction.getInterpParam();
884 wxString optBackprojectName = dialogReconstruction.getBackprojectName();
885 m_iDefaultBackprojector = Backprojector::convertBackprojectNameToID (optBackprojectName.c_str());
886 m_iDefaultTrace = dialogReconstruction.getTrace();
887 if (m_iDefaultNX > 0 && m_iDefaultNY > 0) {
888 ImageFileDocument* pReconDoc = dynamic_cast<ImageFileDocument*>(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT));
889 ImageFile& imageFile = pReconDoc->getImageFile();
890 const Projections& rProj = GetDocument()->getProjections();
891 imageFile.setArraySize (m_iDefaultNX, m_iDefaultNY);
893 if (m_iDefaultFilterMethod != ProcessSignal::FILTER_METHOD_CONVOLUTION && m_iDefaultFilterGeneration == ProcessSignal::FILTER_GENERATION_DIRECT && rProj.geometry() != Scanner::GEOMETRY_PARALLEL) {
894 wxMessageBox ("Sorry!\nCurrently, frequency-based filtering with direct filter generation is not support for geometries other than parallel.\nAborting command.", "Not Supported", wxOK | wxICON_WARNING, m_frame);
898 SGPDriver* pSGPDriver = NULL;
900 wxMemoryDC* pDCPlot = NULL;
902 if (m_iDefaultTrace >= Trace::TRACE_PLOT) {
903 bitmap.Create (500, 500);
904 pDCPlot = new wxMemoryDC;
905 pDCPlot->SelectObject (bitmap);
906 pSGPDriver = new SGPDriver (dynamic_cast<wxDC*>pDCPlot, 500, 500);
907 pSGP = new SGP (*pSGPDriver);
909 Reconstructor* pReconstruct = new Reconstructor (rProj, imageFile, optFilterName.c_str(), m_dDefaultFilterParam, optFilterMethodName.c_str(), m_iDefaultZeropad, optFilterGenerationName.c_str(), optInterpName.c_str(), m_iDefaultInterpParam, optBackprojectName.c_str(), m_iDefaultTrace, pSGP);
912 Reconstructor* pReconstruct = new Reconstructor (rProj, imageFile, optFilterName.c_str(), m_dDefaultFilterParam, optFilterMethodName.c_str(), m_iDefaultZeropad, optFilterGenerationName.c_str(), optInterpName.c_str(), m_iDefaultInterpParam, optBackprojectName.c_str(), m_iDefaultTrace);
916 if (m_iDefaultTrace > Trace::TRACE_CONSOLE) {
917 ReconstructDialog* pDlgReconstruct = new ReconstructDialog (*pReconstruct, rProj, imageFile, m_iDefaultTrace, m_frame);
918 for (int iView = 0; iView < rProj.nView(); iView++) {
921 if (pDlgReconstruct->isCancelled() || ! pDlgReconstruct->reconstructView (iView)) {
922 delete pDlgReconstruct;
924 pReconDoc->DeleteAllViews();
929 while (pDlgReconstruct->isPaused()) {
934 delete pDlgReconstruct;
936 wxProgressDialog dlgProgress (wxString("Reconstruction"), wxString("Reconstruction Progress"), rProj.nView() + 1, m_frame, wxPD_CAN_ABORT);
937 for (int i = 0; i < rProj.nView(); i++) {
938 pReconstruct->reconstructView (i, 1);
939 if (! dlgProgress.Update(i + 1)) {
941 pReconDoc->DeleteAllViews();
947 pReconDoc->Modify(true);
948 pReconDoc->UpdateAllViews(this);
949 ImageFileView* rasterView = dynamic_cast<ImageFileView*>(pReconDoc->GetFirstView());
\r
951 rasterView->getFrame()->SetFocus();
\r
952 rasterView->OnUpdate (rasterView, NULL);
\r
954 std::ostringstream os;
955 os << "Reconstruct " << rProj.getFilename() << ": xSize=" << m_iDefaultNX << ", ySize=" << m_iDefaultNY << ", Filter=" << optFilterName.c_str() << ", FilterParam=" << m_dDefaultFilterParam << ", FilterMethod=" << optFilterMethodName.c_str() << ", FilterGeneration=" << optFilterGenerationName.c_str() << ", Zeropad=" << m_iDefaultZeropad << ", Interpolation=" << optInterpName.c_str() << ", InterpolationParam=" << m_iDefaultInterpParam << ", Backprojection=" << optBackprojectName.c_str();
956 *theApp->getLog() << os.str().c_str() << "\n";
957 imageFile.labelAdd (rProj.getLabel());
958 imageFile.labelAdd (Array2dFileLabel::L_HISTORY, os.str().c_str(), timerRecon.timerEnd());
964 ProjectionFileCanvas*
965 ProjectionFileView::CreateCanvas (wxView *view, wxFrame *parent)
967 ProjectionFileCanvas* pCanvas;
969 parent->GetClientSize(&width, &height);
971 pCanvas = new ProjectionFileCanvas (dynamic_cast<ProjectionFileView*>(view), parent, wxPoint(0, 0), wxSize(width, height), 0);
973 pCanvas->SetScrollbars(20, 20, 50, 50);
974 pCanvas->SetBackgroundColour(*wxWHITE);
981 ProjectionFileView::CreateChildFrame(wxDocument *doc, wxView *view)
983 wxDocChildFrame *subframe = new wxDocChildFrame(doc, view, theApp->getMainFrame(), -1, "Projection Frame", wxPoint(10, 10), wxSize(0, 0), wxDEFAULT_FRAME_STYLE);
985 wxMenu *file_menu = new wxMenu;
987 file_menu->Append(MAINMENU_FILE_CREATE_PHANTOM, "Cr&eate Phantom...");
988 file_menu->Append(wxID_OPEN, "&Open...");
989 file_menu->Append(wxID_SAVE, "&Save");
990 file_menu->Append(wxID_SAVEAS, "Save &As...");
991 file_menu->Append(wxID_CLOSE, "&Close");
993 file_menu->AppendSeparator();
994 file_menu->Append(PJMENU_FILE_PROPERTIES, "P&roperties");
996 file_menu->AppendSeparator();
997 file_menu->Append(wxID_PRINT, "&Print...");
998 file_menu->Append(wxID_PRINT_SETUP, "Print &Setup...");
999 file_menu->Append(wxID_PREVIEW, "Print Pre&view");
1001 wxMenu *process_menu = new wxMenu;
1002 process_menu->Append(PJMENU_PROCESS_RECONSTRUCT, "R&econstruct...");
1004 wxMenu *help_menu = new wxMenu;
1005 help_menu->Append(MAINMENU_HELP_CONTENTS, "&Contents");
1006 help_menu->AppendSeparator();
1007 help_menu->Append(MAINMENU_HELP_ABOUT, "&About");
1009 wxMenuBar *menu_bar = new wxMenuBar;
1011 menu_bar->Append(file_menu, "&File");
1012 menu_bar->Append(process_menu, "&Process");
1013 menu_bar->Append(help_menu, "&Help");
1015 subframe->SetMenuBar(menu_bar);
1017 subframe->Centre(wxBOTH);
1024 ProjectionFileView::OnCreate(wxDocument *doc, long WXUNUSED(flags) )
1026 m_frame = CreateChildFrame(doc, this);
1030 m_frame->GetClientSize(&width, &height);
1031 m_frame->SetTitle("ProjectionFileView");
1032 m_canvas = CreateCanvas(this, m_frame);
1035 int x, y; // X requires a forced resize
1036 m_frame->GetSize(&x, &y);
1037 m_frame->SetSize(-1, -1, x, y);
1040 m_frame->Show(true);
1047 ProjectionFileView::OnDraw (wxDC* dc)
1050 dc->DrawBitmap (m_bitmap, 0, 0, false);
1055 ProjectionFileView::OnUpdate(wxView *WXUNUSED(sender), wxObject *WXUNUSED(hint) )
1057 const Projections& rProj = GetDocument()->getProjections();
1058 const int nDet = rProj.nDet();
1059 const int nView = rProj.nView();
1060 if (nDet != 0 && nView != 0) {
1061 const DetectorArray& detarray = rProj.getDetectorArray(0);
1062 const DetectorValue* detval = detarray.detValues();
1063 double min = detval[0];
1064 double max = detval[0];
1065 for (int iy = 0; iy < nView; iy++) {
1066 const DetectorArray& detarray = rProj.getDetectorArray(iy);
1067 const DetectorValue* detval = detarray.detValues();
1068 for (int ix = 0; ix < nDet; ix++) {
1069 if (min > detval[ix])
1071 else if (max < detval[ix])
1076 unsigned char* imageData = new unsigned char [nDet * nView * 3];
1077 double scale = (max - min) / 255;
1078 for (int iy2 = 0; iy2 < nView; iy2++) {
1079 const DetectorArray& detarray = rProj.getDetectorArray (iy2);
1080 const DetectorValue* detval = detarray.detValues();
1081 for (int ix = 0; ix < nDet; ix++) {
1082 int intensity = static_cast<int>(((detval[ix] - min) / scale) + 0.5);
1083 intensity = clamp(intensity, 0, 255);
1084 int baseAddr = (iy2 * nDet + ix) * 3;
1085 imageData[baseAddr] = imageData[baseAddr+1] = imageData[baseAddr+2] = intensity;
1088 wxImage image (nDet, nView, imageData, true);
1089 m_bitmap = image.ConvertToBitmap();
1093 xSize = clamp (xSize, 0, 800);
1094 ySize = clamp (ySize, 0, 800);
1095 m_frame->SetClientSize (xSize, ySize);
1096 m_canvas->SetScrollbars (20, 20, nDet/20, nView/20);
1100 m_canvas->Refresh();
1104 ProjectionFileView::OnClose (bool deleteWindow)
1106 if (!GetDocument()->Close())
1110 m_canvas->m_pView = NULL;
1112 wxString s(wxTheApp->GetAppName());
1114 m_frame->SetTitle(s);
1130 PlotFileCanvas::PlotFileCanvas (PlotFileView* v, wxFrame *frame, const wxPoint& pos, const wxSize& size, const long style)
1131 : wxScrolledWindow(frame, -1, pos, size, style)
1137 PlotFileCanvas::OnDraw(wxDC& dc)
1140 m_pView->OnDraw(& dc);
1145 IMPLEMENT_DYNAMIC_CLASS(PlotFileView, wxView)
1147 BEGIN_EVENT_TABLE(PlotFileView, wxView)
1148 EVT_MENU(PJMENU_FILE_PROPERTIES, PlotFileView::OnProperties)
1151 PlotFileView::PlotFileView(void)
1152 : wxView(), m_canvas(NULL), m_frame(NULL)
1156 PlotFileView::~PlotFileView(void)
1161 PlotFileView::OnProperties (wxCommandEvent& event)
1163 const PlotFile& rProj = GetDocument()->getPlot();
1164 std::ostringstream os;
\r
1165 os << "Columns: " << rProj.getNumColumns() << ", Records: " << rProj.getNumRecords() << "\n";
1166 *theApp->getLog() << os.str().c_str();
1167 wxMessageDialog dialogMsg (m_frame, os.str().c_str(), "Plot File Properties", wxOK | wxICON_INFORMATION);
1168 dialogMsg.ShowModal();
1173 PlotFileView::CreateCanvas (wxView *view, wxFrame *parent)
1175 PlotFileCanvas* pCanvas;
1177 parent->GetClientSize(&width, &height);
1179 pCanvas = new PlotFileCanvas (dynamic_cast<PlotFileView*>(view), parent, wxPoint(0, 0), wxSize(width, height), 0);
1181 pCanvas->SetBackgroundColour(*wxWHITE);
1188 PlotFileView::CreateChildFrame(wxDocument *doc, wxView *view)
1190 wxDocChildFrame *subframe = new wxDocChildFrame(doc, view, theApp->getMainFrame(), -1, "Plot Frame", wxPoint(10, 10), wxSize(500, 300), wxDEFAULT_FRAME_STYLE);
1192 wxMenu *file_menu = new wxMenu;
1194 file_menu->Append(MAINMENU_FILE_CREATE_PHANTOM, "Cr&eate Phantom...");
1195 file_menu->Append(wxID_OPEN, "&Open...");
1196 file_menu->Append(wxID_SAVE, "&Save");
1197 file_menu->Append(wxID_SAVEAS, "Save &As...");
1198 file_menu->Append(wxID_CLOSE, "&Close");
1200 file_menu->AppendSeparator();
1201 file_menu->Append(PJMENU_FILE_PROPERTIES, "P&roperties");
1203 file_menu->AppendSeparator();
1204 file_menu->Append(wxID_PRINT, "&Print...");
1205 file_menu->Append(wxID_PRINT_SETUP, "Print &Setup...");
1206 file_menu->Append(wxID_PREVIEW, "Print Pre&view");
1208 wxMenu *help_menu = new wxMenu;
1209 help_menu->Append(MAINMENU_HELP_CONTENTS, "&Contents");
1210 help_menu->AppendSeparator();
1211 help_menu->Append(MAINMENU_HELP_ABOUT, "&About");
1213 wxMenuBar *menu_bar = new wxMenuBar;
1215 menu_bar->Append(file_menu, "&File");
1216 menu_bar->Append(help_menu, "&Help");
1218 subframe->SetMenuBar(menu_bar);
1220 subframe->Centre(wxBOTH);
1227 PlotFileView::OnCreate(wxDocument *doc, long WXUNUSED(flags) )
1229 m_frame = CreateChildFrame(doc, this);
1233 m_frame->GetClientSize(&width, &height);
1234 m_frame->SetTitle ("Plot File");
1235 m_canvas = CreateCanvas (this, m_frame);
1238 int x, y; // X requires a forced resize
1239 m_frame->GetSize(&x, &y);
1240 m_frame->SetSize(-1, -1, x, y);
1243 m_frame->Show(true);
1250 PlotFileView::OnDraw (wxDC* dc)
1253 dc->DrawBitmap (m_bitmap, 0, 0, false);
1258 PlotFileView::OnUpdate(wxView *WXUNUSED(sender), wxObject *WXUNUSED(hint) )
1260 const PlotFile& rPlot = GetDocument()->getPlotFile();
1261 const int iNColumns = rPlot.getNumColumns();
\r
1262 const int iNRecords = rPlot.getNumRecords();
1263 if (iNColumns > 0 && iNRecords > 0) {
1264 // int xSize, ySize;
\r
1265 // m_frame->SetClientSize (xSize, ySize);
1266 // m_canvas->SetScrollbars (20, 20, nDet/20, nView/20);
1270 m_canvas->Refresh();
1274 PlotFileView::OnClose (bool deleteWindow)
1276 if (!GetDocument()->Close())
1280 m_canvas->m_pView = NULL;
1282 wxString s(wxTheApp->GetAppName());
1284 m_frame->SetTitle(s);