From 7aa9fa8966a8ff45125fd1f5e10fbdbcd3a7f9cb Mon Sep 17 00:00:00 2001 From: "Kevin M. Rosenberg" Date: Tue, 2 Jan 2001 13:32:51 +0000 Subject: [PATCH] r341: *** empty log message *** --- src/views.cpp | 2124 ++++++++++++++++++++++++------------------------- 1 file changed, 1062 insertions(+), 1062 deletions(-) diff --git a/src/views.cpp b/src/views.cpp index 454b7d2..1ba30cd 100644 --- a/src/views.cpp +++ b/src/views.cpp @@ -9,7 +9,7 @@ ** This is part of the CTSim program ** Copyright (C) 1983-2000 Kevin Rosenberg ** -** $Id: views.cpp,v 1.46 2001/01/02 08:18:33 kevin Exp $ +** $Id: views.cpp,v 1.47 2001/01/02 13:32:51 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 @@ -53,13 +53,13 @@ #include "backprojectors.h" #include "reconstruct.h" #include "timer.h" - -#if defined(MSVC) || HAVE_SSTREAM -#include -#else -#include -#endif - + +#if defined(MSVC) || HAVE_SSTREAM +#include +#else +#include +#endif + // ImageFileCanvas @@ -71,8 +71,8 @@ END_EVENT_TABLE() ImageFileCanvas::ImageFileCanvas (ImageFileView* v, wxFrame *frame, const wxPoint& pos, const wxSize& size, const long style) : wxScrolledWindow(frame, -1, pos, size, style) { - m_pView = v; - m_xCursor = -1; + m_pView = v; + m_xCursor = -1; m_yCursor = -1; } @@ -81,35 +81,35 @@ ImageFileCanvas::OnDraw(wxDC& dc) { if (m_pView) m_pView->OnDraw(& dc); -} - -void -ImageFileCanvas::DrawRubberBandCursor (wxDC& dc, int x, int y) -{ - const ImageFile& rIF = m_pView->GetDocument()->getImageFile(); - int nx = rIF.nx(); - int ny = rIF.ny(); - - int yPt = ny - y - 1; - dc.SetLogicalFunction (wxINVERT); - dc.SetPen (*wxGREEN_PEN); - dc.DrawLine (0, yPt, nx, yPt); - dc.DrawLine (x, 0, x, ny); - dc.SetLogicalFunction (wxCOPY); -} - -bool -ImageFileCanvas::GetCurrentCursor (int& x, int& y) -{ - x = m_xCursor; - y = m_yCursor; - - if (m_xCursor >= 0 && m_yCursor >= 0) - return true; - else - return false; -} - +} + +void +ImageFileCanvas::DrawRubberBandCursor (wxDC& dc, int x, int y) +{ + const ImageFile& rIF = m_pView->GetDocument()->getImageFile(); + int nx = rIF.nx(); + int ny = rIF.ny(); + + int yPt = ny - y - 1; + dc.SetLogicalFunction (wxINVERT); + dc.SetPen (*wxGREEN_PEN); + dc.DrawLine (0, yPt, nx, yPt); + dc.DrawLine (x, 0, x, ny); + dc.SetLogicalFunction (wxCOPY); +} + +bool +ImageFileCanvas::GetCurrentCursor (int& x, int& y) +{ + x = m_xCursor; + y = m_yCursor; + + if (m_xCursor >= 0 && m_yCursor >= 0) + return true; + else + return false; +} + void ImageFileCanvas::OnMouseEvent(wxMouseEvent& event) { @@ -120,87 +120,87 @@ ImageFileCanvas::OnMouseEvent(wxMouseEvent& event) PrepareDC(dc); wxPoint pt(event.GetLogicalPosition(dc)); - - const ImageFile& rIF = m_pView->GetDocument()->getImageFile(); - ImageFileArrayConst v = rIF.getArray(); - int nx = rIF.nx(); - int ny = rIF.ny(); - const int yPt = ny - 1 - pt.y; - if (event.RightIsDown()) { - if (pt.x >= 0 && pt.x < nx && pt.y >= 0 && pt.y < ny) { - std::ostringstream os; - os << "Image value (" << pt.x << "," << yPt << ") = " << v[pt.x][yPt]; - if (rIF.isComplex()) { - double dImag = rIF.getImaginaryArray()[pt.x][yPt]; - if (dImag < 0) - os << " - " << -dImag; - else - os << " + " << dImag; - os << "i\n"; - } else - os << "\n"; - *theApp->getLog() << os.str().c_str(); - } else - *theApp->getLog() << "Mouse out of image range (" << pt.x << "," << yPt << ")\n"; - } - else if (event.LeftIsDown() || event.LeftUp() || event.RightUp()) { - if (pt.x >= 0 && pt.x < nx && pt.y >= 0 && pt.y < ny) { - if (m_xCursor >= 0 && m_yCursor >= 0) { - DrawRubberBandCursor (dc, m_xCursor, m_yCursor); - } - DrawRubberBandCursor (dc, pt.x, yPt); - m_xCursor = pt.x; - m_yCursor = yPt; - } else - *theApp->getLog() << "Mouse out of image range (" << pt.x << "," << yPt << ")\n"; - } - if (event.LeftUp()) { - std::ostringstream os; - os << "Selected column " << pt.x << " , row " << yPt << "\n"; - *theApp->getLog() << os.str().c_str(); - } + + const ImageFile& rIF = m_pView->GetDocument()->getImageFile(); + ImageFileArrayConst v = rIF.getArray(); + int nx = rIF.nx(); + int ny = rIF.ny(); + const int yPt = ny - 1 - pt.y; + if (event.RightIsDown()) { + if (pt.x >= 0 && pt.x < nx && pt.y >= 0 && pt.y < ny) { + std::ostringstream os; + os << "Image value (" << pt.x << "," << yPt << ") = " << v[pt.x][yPt]; + if (rIF.isComplex()) { + double dImag = rIF.getImaginaryArray()[pt.x][yPt]; + if (dImag < 0) + os << " - " << -dImag; + else + os << " + " << dImag; + os << "i\n"; + } else + os << "\n"; + *theApp->getLog() << os.str().c_str(); + } else + *theApp->getLog() << "Mouse out of image range (" << pt.x << "," << yPt << ")\n"; + } + else if (event.LeftIsDown() || event.LeftUp() || event.RightUp()) { + if (pt.x >= 0 && pt.x < nx && pt.y >= 0 && pt.y < ny) { + if (m_xCursor >= 0 && m_yCursor >= 0) { + DrawRubberBandCursor (dc, m_xCursor, m_yCursor); + } + DrawRubberBandCursor (dc, pt.x, yPt); + m_xCursor = pt.x; + m_yCursor = yPt; + } else + *theApp->getLog() << "Mouse out of image range (" << pt.x << "," << yPt << ")\n"; + } + if (event.LeftUp()) { + std::ostringstream os; + os << "Selected column " << pt.x << " , row " << yPt << "\n"; + *theApp->getLog() << os.str().c_str(); + } } // ImageFileView IMPLEMENT_DYNAMIC_CLASS(ImageFileView, wxView) -BEGIN_EVENT_TABLE(ImageFileView, wxView) +BEGIN_EVENT_TABLE(ImageFileView, wxView) EVT_MENU(IFMENU_FILE_EXPORT, ImageFileView::OnExport) EVT_MENU(IFMENU_FILE_PROPERTIES, ImageFileView::OnProperties) EVT_MENU(IFMENU_VIEW_SCALE_MINMAX, ImageFileView::OnScaleMinMax) -EVT_MENU(IFMENU_VIEW_SCALE_AUTO, ImageFileView::OnScaleAuto) -EVT_MENU(IFMENU_COMPARE_IMAGES, ImageFileView::OnCompare) -EVT_MENU(IFMENU_COMPARE_ROW, ImageFileView::OnCompareRow) -EVT_MENU(IFMENU_COMPARE_COL, ImageFileView::OnCompareCol) -EVT_MENU(IFMENU_FILTER_INVERTVALUES, ImageFileView::OnInvertValues) -EVT_MENU(IFMENU_FILTER_SQUARE, ImageFileView::OnSquare) -EVT_MENU(IFMENU_FILTER_SQRT, ImageFileView::OnSquareRoot) -EVT_MENU(IFMENU_FILTER_LOG, ImageFileView::OnLog) -EVT_MENU(IFMENU_FILTER_EXP, ImageFileView::OnExp) -EVT_MENU(IFMENU_FILTER_FOURIER, ImageFileView::OnFourier) -EVT_MENU(IFMENU_FILTER_INVERSE_FOURIER, ImageFileView::OnInverseFourier) -EVT_MENU(IFMENU_FILTER_SHUFFLEFOURIERTONATURALORDER, ImageFileView::OnShuffleFourierToNaturalOrder) -EVT_MENU(IFMENU_FILTER_SHUFFLENATURALTOFOURIERORDER, ImageFileView::OnShuffleNaturalToFourierOrder) -EVT_MENU(IFMENU_IMAGE_ADD, ImageFileView::OnAdd) -EVT_MENU(IFMENU_IMAGE_SUBTRACT, ImageFileView::OnSubtract) -EVT_MENU(IFMENU_IMAGE_MULTIPLY, ImageFileView::OnMultiply) -EVT_MENU(IFMENU_IMAGE_DIVIDE, ImageFileView::OnDivide) -EVT_MENU(IFMENU_IMAGE_SCALESIZE, ImageFileView::OnScaleSize) -#ifdef HAVE_FFTW -EVT_MENU(IFMENU_FILTER_FFT, ImageFileView::OnFFT) -EVT_MENU(IFMENU_FILTER_IFFT, ImageFileView::OnIFFT) -#endif -EVT_MENU(IFMENU_FILTER_MAGNITUDE, ImageFileView::OnMagnitude) -EVT_MENU(IFMENU_FILTER_PHASE, ImageFileView::OnPhase) -EVT_MENU(IFMENU_PLOT_ROW, ImageFileView::OnPlotRow) -EVT_MENU(IFMENU_PLOT_COL, ImageFileView::OnPlotCol) +EVT_MENU(IFMENU_VIEW_SCALE_AUTO, ImageFileView::OnScaleAuto) +EVT_MENU(IFMENU_COMPARE_IMAGES, ImageFileView::OnCompare) +EVT_MENU(IFMENU_COMPARE_ROW, ImageFileView::OnCompareRow) +EVT_MENU(IFMENU_COMPARE_COL, ImageFileView::OnCompareCol) +EVT_MENU(IFMENU_FILTER_INVERTVALUES, ImageFileView::OnInvertValues) +EVT_MENU(IFMENU_FILTER_SQUARE, ImageFileView::OnSquare) +EVT_MENU(IFMENU_FILTER_SQRT, ImageFileView::OnSquareRoot) +EVT_MENU(IFMENU_FILTER_LOG, ImageFileView::OnLog) +EVT_MENU(IFMENU_FILTER_EXP, ImageFileView::OnExp) +EVT_MENU(IFMENU_FILTER_FOURIER, ImageFileView::OnFourier) +EVT_MENU(IFMENU_FILTER_INVERSE_FOURIER, ImageFileView::OnInverseFourier) +EVT_MENU(IFMENU_FILTER_SHUFFLEFOURIERTONATURALORDER, ImageFileView::OnShuffleFourierToNaturalOrder) +EVT_MENU(IFMENU_FILTER_SHUFFLENATURALTOFOURIERORDER, ImageFileView::OnShuffleNaturalToFourierOrder) +EVT_MENU(IFMENU_IMAGE_ADD, ImageFileView::OnAdd) +EVT_MENU(IFMENU_IMAGE_SUBTRACT, ImageFileView::OnSubtract) +EVT_MENU(IFMENU_IMAGE_MULTIPLY, ImageFileView::OnMultiply) +EVT_MENU(IFMENU_IMAGE_DIVIDE, ImageFileView::OnDivide) +EVT_MENU(IFMENU_IMAGE_SCALESIZE, ImageFileView::OnScaleSize) +#ifdef HAVE_FFTW +EVT_MENU(IFMENU_FILTER_FFT, ImageFileView::OnFFT) +EVT_MENU(IFMENU_FILTER_IFFT, ImageFileView::OnIFFT) +#endif +EVT_MENU(IFMENU_FILTER_MAGNITUDE, ImageFileView::OnMagnitude) +EVT_MENU(IFMENU_FILTER_PHASE, ImageFileView::OnPhase) +EVT_MENU(IFMENU_PLOT_ROW, ImageFileView::OnPlotRow) +EVT_MENU(IFMENU_PLOT_COL, ImageFileView::OnPlotCol) END_EVENT_TABLE() ImageFileView::ImageFileView(void) : wxView(), m_canvas(NULL), m_frame(NULL), m_bMinSpecified(false), m_bMaxSpecified(false) -{ - m_iDefaultExportFormatID = ImageFile::FORMAT_PNG; +{ + m_iDefaultExportFormatID = ImageFile::FORMAT_PNG; } ImageFileView::~ImageFileView(void) @@ -216,24 +216,24 @@ ImageFileView::OnProperties (wxCommandEvent& event) else { const std::string& rFilename = rIF.getFilename(); std::ostringstream os; - double min, max, mean, mode, median, stddev; - rIF.statistics (rIF.getArray(), min, max, mean, mode, median, stddev); - os << "Filename: " << rFilename << "\n"; - os << "Size: (" << rIF.nx() << "," << rIF.ny() << ")\n"; - os << "Data type: "; - if (rIF.isComplex()) - os << "Complex\n"; - else - os << "Real\n"; + double min, max, mean, mode, median, stddev; + rIF.statistics (rIF.getArray(), min, max, mean, mode, median, stddev); + os << "Filename: " << rFilename << "\n"; + os << "Size: (" << rIF.nx() << "," << rIF.ny() << ")\n"; + os << "Data type: "; + if (rIF.isComplex()) + os << "Complex\n"; + else + os << "Real\n"; os << "\nMinimum: "< 0) { - os << "\n"; - rIF.printLabelsBrief (os); - } + if (rIF.isComplex()) { + rIF.statistics (rIF.getImaginaryArray(), min, max, mean, mode, median, stddev); + os << "\nImaginary: min: "< 0) { + os << "\n"; + rIF.printLabelsBrief (os); + } *theApp->getLog() << os.str().c_str(); wxMessageDialog dialogMsg (m_frame, os.str().c_str(), "Imagefile Properties", wxOK | wxICON_INFORMATION); dialogMsg.ShowModal(); @@ -244,19 +244,19 @@ void ImageFileView::OnScaleAuto (wxCommandEvent& event) { const ImageFile& rIF = GetDocument()->getImageFile(); - double min, max, mean, mode, median, stddev; - rIF.statistics(min, max, mean, mode, median, stddev); + double min, max, mean, mode, median, stddev; + rIF.statistics(min, max, mean, mode, median, stddev); DialogAutoScaleParameters dialogAutoScale (m_frame, mean, mode, median, stddev, m_dAutoScaleFactor); int iRetVal = dialogAutoScale.ShowModal(); if (iRetVal == wxID_OK) { m_bMinSpecified = true; m_bMaxSpecified = true; - double dMin, dMax; + double dMin, dMax; if (dialogAutoScale.getMinMax (&dMin, &dMax)) { m_dMinPixel = dMin; m_dMaxPixel = dMax; m_dAutoScaleFactor = dialogAutoScale.getAutoScaleFactor(); - OnUpdate (this, NULL); + OnUpdate (this, NULL); } } } @@ -284,396 +284,396 @@ ImageFileView::OnScaleMinMax (wxCommandEvent& event) OnUpdate (this, NULL); } } - -void -ImageFileView::OnCompare (wxCommandEvent& event) -{ - std::vector vecIF; - theApp->getCompatibleImages (GetDocument(), vecIF); - - if (vecIF.size() == 0) { - wxMessageBox("There are no compatible image files open for comparision", "No comparison images"); - } else { - DialogGetComparisonImage dialogGetCompare(m_frame, "Get Comparison Image", vecIF, true); - - if (dialogGetCompare.ShowModal() == wxID_OK) { - const ImageFile& rIF = GetDocument()->getImageFile(); - ImageFileDocument* pCompareDoc = dialogGetCompare.getImageFileDocument(); - const ImageFile& rCompareIF = pCompareDoc->getImageFile(); - std::ostringstream os; - double min, max, mean, mode, median, stddev; - rIF.statistics (min, max, mean, mode, median, stddev); - os << rIF.getFilename() << ": minimum=" << min << ", maximum=" << max << ", mean=" << mean << ", mode=" << mode << ", median=" << median << ", stddev=" << stddev << "\n"; - rCompareIF.statistics (min, max, mean, mode, median, stddev); - os << pCompareDoc->GetFirstView()->GetFrame()->GetTitle().c_str() << ": minimum=" << min << ", maximum=" << max << ", mean=" << mean << ", mode=" << mode << ", median=" << median << ", stddev=" << stddev << "\n"; - os << "\n"; - double d, r, e; - rIF.comparativeStatistics (rCompareIF, d, r, e); - os << "Comparative Statistics: d=" << d << ", r=" << r << ", e=" << e << "\n"; - *theApp->getLog() << os.str().c_str(); - if (dialogGetCompare.getMakeDifferenceImage()) { - ImageFileDocument* pDifferenceDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT)); - if (! pDifferenceDoc) { - sys_error (ERR_SEVERE, "Unable to create image file"); - return; - } - ImageFile& differenceImage = pDifferenceDoc->getImageFile(); - - differenceImage.setArraySize (rIF.nx(), rIF.ny()); - if (! rIF.subtractImages (rCompareIF, differenceImage)) { - pDifferenceDoc->DeleteAllViews(); - return; - } - - wxString s = GetDocument()->GetFirstView()->GetFrame()->GetTitle() + ": "; - differenceImage.labelsCopy (rIF, s.c_str()); - s = pCompareDoc->GetFirstView()->GetFrame()->GetTitle() + ": "; - differenceImage.labelsCopy (rCompareIF, s.c_str()); - std::ostringstream osLabel; - osLabel << "Compare image " << GetDocument()->GetFirstView()->GetFrame()->GetTitle().c_str() - << " and " << pCompareDoc->GetFirstView()->GetFrame()->GetTitle().c_str() << ": " - << os.str().c_str(); - differenceImage.labelAdd (os.str().c_str()); - if (theApp->getSetModifyNewDocs()) - pDifferenceDoc->Modify(true); - pDifferenceDoc->UpdateAllViews(this); - pDifferenceDoc->GetFirstView()->OnUpdate (this, NULL); - } - wxMessageBox(os.str().c_str(), "Image Comparison"); - } - } -} - -void -ImageFileView::OnInvertValues (wxCommandEvent& event) -{ - ImageFile& rIF = GetDocument()->getImageFile(); - rIF.invertPixelValues (rIF); - rIF.labelAdd ("Invert Pixel Values"); - if (theApp->getSetModifyNewDocs()) - GetDocument()->Modify(TRUE); - GetDocument()->UpdateAllViews(this); -} - -void -ImageFileView::OnSquare (wxCommandEvent& event) -{ - ImageFile& rIF = GetDocument()->getImageFile(); - rIF.square (rIF); - rIF.labelAdd ("Square Pixel Values"); - if (theApp->getSetModifyNewDocs()) - GetDocument()->Modify(TRUE); - GetDocument()->UpdateAllViews(this); -} - -void -ImageFileView::OnSquareRoot (wxCommandEvent& event) -{ - ImageFile& rIF = GetDocument()->getImageFile(); - rIF.sqrt (rIF); - rIF.labelAdd ("Square-root Pixel Values"); - if (theApp->getSetModifyNewDocs()) - GetDocument()->Modify(TRUE); - GetDocument()->UpdateAllViews(this); -} - -void -ImageFileView::OnLog (wxCommandEvent& event) -{ - ImageFile& rIF = GetDocument()->getImageFile(); - rIF.log (rIF); - rIF.labelAdd ("Logrithm base-e Pixel Values"); - if (theApp->getSetModifyNewDocs()) - GetDocument()->Modify(TRUE); - GetDocument()->UpdateAllViews(this); -} - -void -ImageFileView::OnExp (wxCommandEvent& event) -{ - ImageFile& rIF = GetDocument()->getImageFile(); - rIF.exp (rIF); - rIF.labelAdd ("Exponent base-e Pixel Values"); - if (theApp->getSetModifyNewDocs()) - GetDocument()->Modify(TRUE); - GetDocument()->UpdateAllViews(this); -} - -void -ImageFileView::OnAdd (wxCommandEvent& event) -{ - std::vector vecIF; - theApp->getCompatibleImages (GetDocument(), vecIF); - - if (vecIF.size() == 0) { - wxMessageBox ("There are no compatible image files open for comparision", "No comparison images"); - } else { - DialogGetComparisonImage dialogGetCompare (m_frame, "Get Image to Add", vecIF, false); - - if (dialogGetCompare.ShowModal() == wxID_OK) { - ImageFile& rIF = GetDocument()->getImageFile(); - ImageFileDocument* pRHSDoc = dialogGetCompare.getImageFileDocument(); - const ImageFile& rRHSIF = pRHSDoc->getImageFile(); - ImageFileDocument* pNewDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT)); - if (! pNewDoc) { - sys_error (ERR_SEVERE, "Unable to create image file"); - return; - } - ImageFile& newImage = pNewDoc->getImageFile(); - newImage.setArraySize (rIF.nx(), rIF.ny()); - rIF.addImages (rRHSIF, newImage); - std::ostringstream os; - os << "Add image " << GetDocument()->GetFirstView()->GetFrame()->GetTitle().c_str() << " and " - << pRHSDoc->GetFirstView()->GetFrame()->GetTitle().c_str(); - wxString s = GetDocument()->GetFirstView()->GetFrame()->GetTitle() + ": "; - newImage.labelsCopy (rIF, s.c_str()); - s = pRHSDoc->GetFirstView()->GetFrame()->GetTitle() + ": "; - newImage.labelsCopy (rRHSIF, s.c_str()); - newImage.labelAdd (os.str().c_str()); - *theApp->getLog() << os.str().c_str() << "\n"; - if (theApp->getSetModifyNewDocs()) - pNewDoc->Modify(TRUE); - pNewDoc->UpdateAllViews(this); - pNewDoc->GetFirstView()->OnUpdate (this, NULL); - } - } -} - -void -ImageFileView::OnSubtract (wxCommandEvent& event) -{ - std::vector vecIF; - theApp->getCompatibleImages (GetDocument(), vecIF); - - if (vecIF.size() == 0) { - wxMessageBox ("There are no compatible image files open for comparision", "No comparison images"); - } else { - DialogGetComparisonImage dialogGetCompare (m_frame, "Get Image to Subtract", vecIF, false); - - if (dialogGetCompare.ShowModal() == wxID_OK) { - ImageFile& rIF = GetDocument()->getImageFile(); - ImageFileDocument* pRHSDoc = dialogGetCompare.getImageFileDocument(); - const ImageFile& rRHSIF = pRHSDoc->getImageFile(); - ImageFileDocument* pNewDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT)); - if (! pNewDoc) { - sys_error (ERR_SEVERE, "Unable to create image file"); - return; - } - ImageFile& newImage = pNewDoc->getImageFile(); - newImage.setArraySize (rIF.nx(), rIF.ny()); - rIF.subtractImages (rRHSIF, newImage); - std::ostringstream os; - os << "Subtract image " << GetDocument()->GetFirstView()->GetFrame()->GetTitle().c_str() << " and " - << pRHSDoc->GetFirstView()->GetFrame()->GetTitle().c_str(); - wxString s = GetDocument()->GetFirstView()->GetFrame()->GetTitle() + ": "; - newImage.labelsCopy (rIF, s.c_str()); - s = pRHSDoc->GetFirstView()->GetFrame()->GetTitle() + ": "; - newImage.labelsCopy (rRHSIF, s.c_str()); - newImage.labelAdd (os.str().c_str()); - *theApp->getLog() << os.str().c_str() << "\n"; - if (theApp->getSetModifyNewDocs()) - pNewDoc->Modify(TRUE); - pNewDoc->UpdateAllViews(this); - pNewDoc->GetFirstView()->OnUpdate (this, NULL); - } - } -} - -void -ImageFileView::OnMultiply (wxCommandEvent& event) -{ - std::vector vecIF; - theApp->getCompatibleImages (GetDocument(), vecIF); - - if (vecIF.size() == 0) { - wxMessageBox ("There are no compatible image files open for comparision", "No comparison images"); - } else { - DialogGetComparisonImage dialogGetCompare (m_frame, "Get Image to Multiply", vecIF, false); - - if (dialogGetCompare.ShowModal() == wxID_OK) { - ImageFile& rIF = GetDocument()->getImageFile(); - ImageFileDocument* pRHSDoc = dialogGetCompare.getImageFileDocument(); - const ImageFile& rRHSIF = pRHSDoc->getImageFile(); - ImageFileDocument* pNewDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT)); - if (! pNewDoc) { - sys_error (ERR_SEVERE, "Unable to create image file"); - return; - } - ImageFile& newImage = pNewDoc->getImageFile(); - newImage.setArraySize (rIF.nx(), rIF.ny()); - rIF.multiplyImages (rRHSIF, newImage); - std::ostringstream os; - os << "Multiply image " << GetDocument()->GetFirstView()->GetFrame()->GetTitle().c_str() << " and " - << pRHSDoc->GetFirstView()->GetFrame()->GetTitle().c_str(); - wxString s = GetDocument()->GetFirstView()->GetFrame()->GetTitle() + ": "; - newImage.labelsCopy (rIF, s.c_str()); - s = pRHSDoc->GetFirstView()->GetFrame()->GetTitle() + ": "; - newImage.labelsCopy (rRHSIF, s.c_str()); - newImage.labelAdd (os.str().c_str()); - *theApp->getLog() << os.str().c_str() << "\n"; - if (theApp->getSetModifyNewDocs()) - pNewDoc->Modify(TRUE); - pNewDoc->UpdateAllViews(this); - pNewDoc->GetFirstView()->OnUpdate (this, NULL); - } - } -} - -void -ImageFileView::OnDivide (wxCommandEvent& event) -{ - std::vector vecIF; - theApp->getCompatibleImages (GetDocument(), vecIF); - - if (vecIF.size() == 0) { - wxMessageBox ("There are no compatible image files open for comparision", "No comparison images"); - } else { - DialogGetComparisonImage dialogGetCompare (m_frame, "Get Image to Divide", vecIF, false); - - if (dialogGetCompare.ShowModal() == wxID_OK) { - ImageFile& rIF = GetDocument()->getImageFile(); - ImageFileDocument* pRHSDoc = dialogGetCompare.getImageFileDocument(); - const ImageFile& rRHSIF = pRHSDoc->getImageFile(); - ImageFileDocument* pNewDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT)); - if (! pNewDoc) { - sys_error (ERR_SEVERE, "Unable to create image file"); - return; - } - ImageFile& newImage = pNewDoc->getImageFile(); - newImage.setArraySize (rIF.nx(), rIF.ny()); - rIF.divideImages (rRHSIF, newImage); - std::ostringstream os; - os << "Divide image " << GetDocument()->GetFirstView()->GetFrame()->GetTitle().c_str() << " by " - << pRHSDoc->GetFirstView()->GetFrame()->GetTitle().c_str(); - wxString s = GetDocument()->GetFirstView()->GetFrame()->GetTitle() + ": "; - newImage.labelsCopy (rIF, s.c_str()); - s = pRHSDoc->GetFirstView()->GetFrame()->GetTitle() + ": "; - newImage.labelsCopy (rRHSIF, s.c_str()); - newImage.labelAdd (os.str().c_str()); - *theApp->getLog() << os.str().c_str() << "\n"; - if (theApp->getSetModifyNewDocs()) - pNewDoc->Modify(TRUE); - pNewDoc->UpdateAllViews(this); - pNewDoc->GetFirstView()->OnUpdate (this, NULL); - } - } -} - - -#ifdef HAVE_FFTW -void -ImageFileView::OnFFT (wxCommandEvent& event) -{ - ImageFile& rIF = GetDocument()->getImageFile(); - wxProgressDialog dlgProgress (wxString("FFT"), wxString("FFT Progress"), 1, m_frame, wxPD_APP_MODAL); - rIF.fft (rIF); - rIF.labelAdd ("FFT Image"); - m_bMinSpecified = false; - m_bMaxSpecified = false; - if (theApp->getSetModifyNewDocs()) - GetDocument()->Modify(TRUE); - GetDocument()->UpdateAllViews(this); -} - -void -ImageFileView::OnIFFT (wxCommandEvent& event) -{ - ImageFile& rIF = GetDocument()->getImageFile(); - wxProgressDialog dlgProgress (wxString("IFFT"), wxString("IFFT Progress"), 1, m_frame, wxPD_APP_MODAL); - rIF.ifft (rIF); - rIF.labelAdd ("IFFT Image"); - m_bMinSpecified = false; - m_bMaxSpecified = false; - if (theApp->getSetModifyNewDocs()) - GetDocument()->Modify(TRUE); - GetDocument()->UpdateAllViews(this); -} -#endif - -void -ImageFileView::OnFourier (wxCommandEvent& event) -{ - ImageFile& rIF = GetDocument()->getImageFile(); - wxProgressDialog dlgProgress (wxString("Fourier"), wxString("Fourier Progress"), 1, m_frame, wxPD_APP_MODAL); - rIF.fourier (rIF); - rIF.labelAdd ("Fourier Image"); - m_bMinSpecified = false; - m_bMaxSpecified = false; - if (theApp->getSetModifyNewDocs()) - GetDocument()->Modify(TRUE); - GetDocument()->UpdateAllViews(this); -} -void -ImageFileView::OnInverseFourier (wxCommandEvent& event) -{ - ImageFile& rIF = GetDocument()->getImageFile(); - wxProgressDialog dlgProgress (wxString("Inverse Fourier"), wxString("Inverse Fourier Progress"), 1, m_frame, wxPD_APP_MODAL); - rIF.inverseFourier (rIF); - rIF.labelAdd ("Inverse Fourier Image"); - m_bMinSpecified = false; - m_bMaxSpecified = false; - if (theApp->getSetModifyNewDocs()) - GetDocument()->Modify(TRUE); - GetDocument()->UpdateAllViews(this); -} - -void -ImageFileView::OnShuffleNaturalToFourierOrder (wxCommandEvent& event) -{ - ImageFile& rIF = GetDocument()->getImageFile(); - Fourier::shuffleNaturalToFourierOrder (rIF); - rIF.labelAdd ("Shuffle Natural To Fourier Order"); - m_bMinSpecified = false; - m_bMaxSpecified = false; - if (theApp->getSetModifyNewDocs()) - GetDocument()->Modify(TRUE); - GetDocument()->UpdateAllViews(this); -} - -void -ImageFileView::OnShuffleFourierToNaturalOrder (wxCommandEvent& event) -{ - ImageFile& rIF = GetDocument()->getImageFile(); - Fourier::shuffleFourierToNaturalOrder (rIF); - rIF.labelAdd ("Shuffle Fourier To Natural Order"); - m_bMinSpecified = false; - m_bMaxSpecified = false; - if (theApp->getSetModifyNewDocs()) - GetDocument()->Modify(TRUE); - GetDocument()->UpdateAllViews(this); -} - -void -ImageFileView::OnMagnitude (wxCommandEvent& event) -{ - ImageFile& rIF = GetDocument()->getImageFile(); - if (rIF.isComplex()) { - rIF.magnitude (rIF); - rIF.labelAdd ("Magnitude of complex-image"); - m_bMinSpecified = false; - m_bMaxSpecified = false; - if (theApp->getSetModifyNewDocs()) - GetDocument()->Modify(TRUE); - GetDocument()->UpdateAllViews(this); - } -} - -void -ImageFileView::OnPhase (wxCommandEvent& event) -{ - ImageFile& rIF = GetDocument()->getImageFile(); - if (rIF.isComplex()) { - rIF.phase (rIF); - rIF.labelAdd ("Phase of complex-image"); - m_bMinSpecified = false; - m_bMaxSpecified = false; - if (theApp->getSetModifyNewDocs()) - GetDocument()->Modify(TRUE); - GetDocument()->UpdateAllViews(this); - } -} - - + +void +ImageFileView::OnCompare (wxCommandEvent& event) +{ + std::vector vecIF; + theApp->getCompatibleImages (GetDocument(), vecIF); + + if (vecIF.size() == 0) { + wxMessageBox("There are no compatible image files open for comparision", "No comparison images"); + } else { + DialogGetComparisonImage dialogGetCompare(m_frame, "Get Comparison Image", vecIF, true); + + if (dialogGetCompare.ShowModal() == wxID_OK) { + const ImageFile& rIF = GetDocument()->getImageFile(); + ImageFileDocument* pCompareDoc = dialogGetCompare.getImageFileDocument(); + const ImageFile& rCompareIF = pCompareDoc->getImageFile(); + std::ostringstream os; + double min, max, mean, mode, median, stddev; + rIF.statistics (min, max, mean, mode, median, stddev); + os << rIF.getFilename() << ": minimum=" << min << ", maximum=" << max << ", mean=" << mean << ", mode=" << mode << ", median=" << median << ", stddev=" << stddev << "\n"; + rCompareIF.statistics (min, max, mean, mode, median, stddev); + os << pCompareDoc->GetFirstView()->GetFrame()->GetTitle().c_str() << ": minimum=" << min << ", maximum=" << max << ", mean=" << mean << ", mode=" << mode << ", median=" << median << ", stddev=" << stddev << "\n"; + os << "\n"; + double d, r, e; + rIF.comparativeStatistics (rCompareIF, d, r, e); + os << "Comparative Statistics: d=" << d << ", r=" << r << ", e=" << e << "\n"; + *theApp->getLog() << os.str().c_str(); + if (dialogGetCompare.getMakeDifferenceImage()) { + ImageFileDocument* pDifferenceDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT)); + if (! pDifferenceDoc) { + sys_error (ERR_SEVERE, "Unable to create image file"); + return; + } + ImageFile& differenceImage = pDifferenceDoc->getImageFile(); + + differenceImage.setArraySize (rIF.nx(), rIF.ny()); + if (! rIF.subtractImages (rCompareIF, differenceImage)) { + pDifferenceDoc->DeleteAllViews(); + return; + } + + wxString s = GetDocument()->GetFirstView()->GetFrame()->GetTitle() + ": "; + differenceImage.labelsCopy (rIF, s.c_str()); + s = pCompareDoc->GetFirstView()->GetFrame()->GetTitle() + ": "; + differenceImage.labelsCopy (rCompareIF, s.c_str()); + std::ostringstream osLabel; + osLabel << "Compare image " << GetDocument()->GetFirstView()->GetFrame()->GetTitle().c_str() + << " and " << pCompareDoc->GetFirstView()->GetFrame()->GetTitle().c_str() << ": " + << os.str().c_str(); + differenceImage.labelAdd (os.str().c_str()); + if (theApp->getSetModifyNewDocs()) + pDifferenceDoc->Modify(true); + pDifferenceDoc->UpdateAllViews(this); + pDifferenceDoc->GetFirstView()->OnUpdate (this, NULL); + } + wxMessageBox(os.str().c_str(), "Image Comparison"); + } + } +} + +void +ImageFileView::OnInvertValues (wxCommandEvent& event) +{ + ImageFile& rIF = GetDocument()->getImageFile(); + rIF.invertPixelValues (rIF); + rIF.labelAdd ("Invert Pixel Values"); + if (theApp->getSetModifyNewDocs()) + GetDocument()->Modify(TRUE); + GetDocument()->UpdateAllViews(this); +} + +void +ImageFileView::OnSquare (wxCommandEvent& event) +{ + ImageFile& rIF = GetDocument()->getImageFile(); + rIF.square (rIF); + rIF.labelAdd ("Square Pixel Values"); + if (theApp->getSetModifyNewDocs()) + GetDocument()->Modify(TRUE); + GetDocument()->UpdateAllViews(this); +} + +void +ImageFileView::OnSquareRoot (wxCommandEvent& event) +{ + ImageFile& rIF = GetDocument()->getImageFile(); + rIF.sqrt (rIF); + rIF.labelAdd ("Square-root Pixel Values"); + if (theApp->getSetModifyNewDocs()) + GetDocument()->Modify(TRUE); + GetDocument()->UpdateAllViews(this); +} + +void +ImageFileView::OnLog (wxCommandEvent& event) +{ + ImageFile& rIF = GetDocument()->getImageFile(); + rIF.log (rIF); + rIF.labelAdd ("Logrithm base-e Pixel Values"); + if (theApp->getSetModifyNewDocs()) + GetDocument()->Modify(TRUE); + GetDocument()->UpdateAllViews(this); +} + +void +ImageFileView::OnExp (wxCommandEvent& event) +{ + ImageFile& rIF = GetDocument()->getImageFile(); + rIF.exp (rIF); + rIF.labelAdd ("Exponent base-e Pixel Values"); + if (theApp->getSetModifyNewDocs()) + GetDocument()->Modify(TRUE); + GetDocument()->UpdateAllViews(this); +} + +void +ImageFileView::OnAdd (wxCommandEvent& event) +{ + std::vector vecIF; + theApp->getCompatibleImages (GetDocument(), vecIF); + + if (vecIF.size() == 0) { + wxMessageBox ("There are no compatible image files open for comparision", "No comparison images"); + } else { + DialogGetComparisonImage dialogGetCompare (m_frame, "Get Image to Add", vecIF, false); + + if (dialogGetCompare.ShowModal() == wxID_OK) { + ImageFile& rIF = GetDocument()->getImageFile(); + ImageFileDocument* pRHSDoc = dialogGetCompare.getImageFileDocument(); + const ImageFile& rRHSIF = pRHSDoc->getImageFile(); + ImageFileDocument* pNewDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT)); + if (! pNewDoc) { + sys_error (ERR_SEVERE, "Unable to create image file"); + return; + } + ImageFile& newImage = pNewDoc->getImageFile(); + newImage.setArraySize (rIF.nx(), rIF.ny()); + rIF.addImages (rRHSIF, newImage); + std::ostringstream os; + os << "Add image " << GetDocument()->GetFirstView()->GetFrame()->GetTitle().c_str() << " and " + << pRHSDoc->GetFirstView()->GetFrame()->GetTitle().c_str(); + wxString s = GetDocument()->GetFirstView()->GetFrame()->GetTitle() + ": "; + newImage.labelsCopy (rIF, s.c_str()); + s = pRHSDoc->GetFirstView()->GetFrame()->GetTitle() + ": "; + newImage.labelsCopy (rRHSIF, s.c_str()); + newImage.labelAdd (os.str().c_str()); + *theApp->getLog() << os.str().c_str() << "\n"; + if (theApp->getSetModifyNewDocs()) + pNewDoc->Modify(TRUE); + pNewDoc->UpdateAllViews(this); + pNewDoc->GetFirstView()->OnUpdate (this, NULL); + } + } +} + +void +ImageFileView::OnSubtract (wxCommandEvent& event) +{ + std::vector vecIF; + theApp->getCompatibleImages (GetDocument(), vecIF); + + if (vecIF.size() == 0) { + wxMessageBox ("There are no compatible image files open for comparision", "No comparison images"); + } else { + DialogGetComparisonImage dialogGetCompare (m_frame, "Get Image to Subtract", vecIF, false); + + if (dialogGetCompare.ShowModal() == wxID_OK) { + ImageFile& rIF = GetDocument()->getImageFile(); + ImageFileDocument* pRHSDoc = dialogGetCompare.getImageFileDocument(); + const ImageFile& rRHSIF = pRHSDoc->getImageFile(); + ImageFileDocument* pNewDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT)); + if (! pNewDoc) { + sys_error (ERR_SEVERE, "Unable to create image file"); + return; + } + ImageFile& newImage = pNewDoc->getImageFile(); + newImage.setArraySize (rIF.nx(), rIF.ny()); + rIF.subtractImages (rRHSIF, newImage); + std::ostringstream os; + os << "Subtract image " << GetDocument()->GetFirstView()->GetFrame()->GetTitle().c_str() << " and " + << pRHSDoc->GetFirstView()->GetFrame()->GetTitle().c_str(); + wxString s = GetDocument()->GetFirstView()->GetFrame()->GetTitle() + ": "; + newImage.labelsCopy (rIF, s.c_str()); + s = pRHSDoc->GetFirstView()->GetFrame()->GetTitle() + ": "; + newImage.labelsCopy (rRHSIF, s.c_str()); + newImage.labelAdd (os.str().c_str()); + *theApp->getLog() << os.str().c_str() << "\n"; + if (theApp->getSetModifyNewDocs()) + pNewDoc->Modify(TRUE); + pNewDoc->UpdateAllViews(this); + pNewDoc->GetFirstView()->OnUpdate (this, NULL); + } + } +} + +void +ImageFileView::OnMultiply (wxCommandEvent& event) +{ + std::vector vecIF; + theApp->getCompatibleImages (GetDocument(), vecIF); + + if (vecIF.size() == 0) { + wxMessageBox ("There are no compatible image files open for comparision", "No comparison images"); + } else { + DialogGetComparisonImage dialogGetCompare (m_frame, "Get Image to Multiply", vecIF, false); + + if (dialogGetCompare.ShowModal() == wxID_OK) { + ImageFile& rIF = GetDocument()->getImageFile(); + ImageFileDocument* pRHSDoc = dialogGetCompare.getImageFileDocument(); + const ImageFile& rRHSIF = pRHSDoc->getImageFile(); + ImageFileDocument* pNewDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT)); + if (! pNewDoc) { + sys_error (ERR_SEVERE, "Unable to create image file"); + return; + } + ImageFile& newImage = pNewDoc->getImageFile(); + newImage.setArraySize (rIF.nx(), rIF.ny()); + rIF.multiplyImages (rRHSIF, newImage); + std::ostringstream os; + os << "Multiply image " << GetDocument()->GetFirstView()->GetFrame()->GetTitle().c_str() << " and " + << pRHSDoc->GetFirstView()->GetFrame()->GetTitle().c_str(); + wxString s = GetDocument()->GetFirstView()->GetFrame()->GetTitle() + ": "; + newImage.labelsCopy (rIF, s.c_str()); + s = pRHSDoc->GetFirstView()->GetFrame()->GetTitle() + ": "; + newImage.labelsCopy (rRHSIF, s.c_str()); + newImage.labelAdd (os.str().c_str()); + *theApp->getLog() << os.str().c_str() << "\n"; + if (theApp->getSetModifyNewDocs()) + pNewDoc->Modify(TRUE); + pNewDoc->UpdateAllViews(this); + pNewDoc->GetFirstView()->OnUpdate (this, NULL); + } + } +} + +void +ImageFileView::OnDivide (wxCommandEvent& event) +{ + std::vector vecIF; + theApp->getCompatibleImages (GetDocument(), vecIF); + + if (vecIF.size() == 0) { + wxMessageBox ("There are no compatible image files open for comparision", "No comparison images"); + } else { + DialogGetComparisonImage dialogGetCompare (m_frame, "Get Image to Divide", vecIF, false); + + if (dialogGetCompare.ShowModal() == wxID_OK) { + ImageFile& rIF = GetDocument()->getImageFile(); + ImageFileDocument* pRHSDoc = dialogGetCompare.getImageFileDocument(); + const ImageFile& rRHSIF = pRHSDoc->getImageFile(); + ImageFileDocument* pNewDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT)); + if (! pNewDoc) { + sys_error (ERR_SEVERE, "Unable to create image file"); + return; + } + ImageFile& newImage = pNewDoc->getImageFile(); + newImage.setArraySize (rIF.nx(), rIF.ny()); + rIF.divideImages (rRHSIF, newImage); + std::ostringstream os; + os << "Divide image " << GetDocument()->GetFirstView()->GetFrame()->GetTitle().c_str() << " by " + << pRHSDoc->GetFirstView()->GetFrame()->GetTitle().c_str(); + wxString s = GetDocument()->GetFirstView()->GetFrame()->GetTitle() + ": "; + newImage.labelsCopy (rIF, s.c_str()); + s = pRHSDoc->GetFirstView()->GetFrame()->GetTitle() + ": "; + newImage.labelsCopy (rRHSIF, s.c_str()); + newImage.labelAdd (os.str().c_str()); + *theApp->getLog() << os.str().c_str() << "\n"; + if (theApp->getSetModifyNewDocs()) + pNewDoc->Modify(TRUE); + pNewDoc->UpdateAllViews(this); + pNewDoc->GetFirstView()->OnUpdate (this, NULL); + } + } +} + + +#ifdef HAVE_FFTW +void +ImageFileView::OnFFT (wxCommandEvent& event) +{ + ImageFile& rIF = GetDocument()->getImageFile(); + wxProgressDialog dlgProgress (wxString("FFT"), wxString("FFT Progress"), 1, m_frame, wxPD_APP_MODAL); + rIF.fft (rIF); + rIF.labelAdd ("FFT Image"); + m_bMinSpecified = false; + m_bMaxSpecified = false; + if (theApp->getSetModifyNewDocs()) + GetDocument()->Modify(TRUE); + GetDocument()->UpdateAllViews(this); +} + +void +ImageFileView::OnIFFT (wxCommandEvent& event) +{ + ImageFile& rIF = GetDocument()->getImageFile(); + wxProgressDialog dlgProgress (wxString("IFFT"), wxString("IFFT Progress"), 1, m_frame, wxPD_APP_MODAL); + rIF.ifft (rIF); + rIF.labelAdd ("IFFT Image"); + m_bMinSpecified = false; + m_bMaxSpecified = false; + if (theApp->getSetModifyNewDocs()) + GetDocument()->Modify(TRUE); + GetDocument()->UpdateAllViews(this); +} +#endif + +void +ImageFileView::OnFourier (wxCommandEvent& event) +{ + ImageFile& rIF = GetDocument()->getImageFile(); + wxProgressDialog dlgProgress (wxString("Fourier"), wxString("Fourier Progress"), 1, m_frame, wxPD_APP_MODAL); + rIF.fourier (rIF); + rIF.labelAdd ("Fourier Image"); + m_bMinSpecified = false; + m_bMaxSpecified = false; + if (theApp->getSetModifyNewDocs()) + GetDocument()->Modify(TRUE); + GetDocument()->UpdateAllViews(this); +} +void +ImageFileView::OnInverseFourier (wxCommandEvent& event) +{ + ImageFile& rIF = GetDocument()->getImageFile(); + wxProgressDialog dlgProgress (wxString("Inverse Fourier"), wxString("Inverse Fourier Progress"), 1, m_frame, wxPD_APP_MODAL); + rIF.inverseFourier (rIF); + rIF.labelAdd ("Inverse Fourier Image"); + m_bMinSpecified = false; + m_bMaxSpecified = false; + if (theApp->getSetModifyNewDocs()) + GetDocument()->Modify(TRUE); + GetDocument()->UpdateAllViews(this); +} + +void +ImageFileView::OnShuffleNaturalToFourierOrder (wxCommandEvent& event) +{ + ImageFile& rIF = GetDocument()->getImageFile(); + Fourier::shuffleNaturalToFourierOrder (rIF); + rIF.labelAdd ("Shuffle Natural To Fourier Order"); + m_bMinSpecified = false; + m_bMaxSpecified = false; + if (theApp->getSetModifyNewDocs()) + GetDocument()->Modify(TRUE); + GetDocument()->UpdateAllViews(this); +} + +void +ImageFileView::OnShuffleFourierToNaturalOrder (wxCommandEvent& event) +{ + ImageFile& rIF = GetDocument()->getImageFile(); + Fourier::shuffleFourierToNaturalOrder (rIF); + rIF.labelAdd ("Shuffle Fourier To Natural Order"); + m_bMinSpecified = false; + m_bMaxSpecified = false; + if (theApp->getSetModifyNewDocs()) + GetDocument()->Modify(TRUE); + GetDocument()->UpdateAllViews(this); +} + +void +ImageFileView::OnMagnitude (wxCommandEvent& event) +{ + ImageFile& rIF = GetDocument()->getImageFile(); + if (rIF.isComplex()) { + rIF.magnitude (rIF); + rIF.labelAdd ("Magnitude of complex-image"); + m_bMinSpecified = false; + m_bMaxSpecified = false; + if (theApp->getSetModifyNewDocs()) + GetDocument()->Modify(TRUE); + GetDocument()->UpdateAllViews(this); + } +} + +void +ImageFileView::OnPhase (wxCommandEvent& event) +{ + ImageFile& rIF = GetDocument()->getImageFile(); + if (rIF.isComplex()) { + rIF.phase (rIF); + rIF.labelAdd ("Phase of complex-image"); + m_bMinSpecified = false; + m_bMaxSpecified = false; + if (theApp->getSetModifyNewDocs()) + GetDocument()->Modify(TRUE); + GetDocument()->UpdateAllViews(this); + } +} + + ImageFileCanvas* ImageFileView::CreateCanvas (wxView *view, wxFrame *parent) { @@ -697,16 +697,16 @@ ImageFileView::CreateChildFrame(wxDocument *doc, wxView *view) wxMenu *file_menu = new wxMenu; - file_menu->Append(MAINMENU_FILE_CREATE_PHANTOM, "Cr&eate Phantom..."); - file_menu->Append(MAINMENU_FILE_CREATE_FILTER, "Create &Filter..."); + file_menu->Append(MAINMENU_FILE_CREATE_PHANTOM, "Cr&eate Phantom..."); + file_menu->Append(MAINMENU_FILE_CREATE_FILTER, "Create &Filter..."); file_menu->Append(wxID_OPEN, "&Open..."); file_menu->Append(wxID_SAVE, "&Save"); file_menu->Append(wxID_SAVEAS, "Save &As..."); file_menu->Append(wxID_CLOSE, "&Close"); - file_menu->AppendSeparator(); + file_menu->AppendSeparator(); file_menu->Append(IFMENU_FILE_PROPERTIES, "P&roperties"); - file_menu->Append(IFMENU_FILE_EXPORT, "&Export..."); + file_menu->Append(IFMENU_FILE_EXPORT, "&Export..."); file_menu->AppendSeparator(); file_menu->Append(wxID_PRINT, "&Print..."); @@ -716,56 +716,56 @@ ImageFileView::CreateChildFrame(wxDocument *doc, wxView *view) wxMenu *view_menu = new wxMenu; view_menu->Append(IFMENU_VIEW_SCALE_MINMAX, "Display Scale &Set..."); view_menu->Append(IFMENU_VIEW_SCALE_AUTO, "Display Scale &Auto..."); - - wxMenu* filter_menu = new wxMenu; - filter_menu->Append (IFMENU_FILTER_INVERTVALUES, "&Invert Values"); - filter_menu->Append (IFMENU_FILTER_SQUARE, "&Square"); - filter_menu->Append (IFMENU_FILTER_SQRT, "Square &Root"); - filter_menu->Append (IFMENU_FILTER_LOG, "&Log"); - filter_menu->Append (IFMENU_FILTER_EXP, "&Exp"); - filter_menu->AppendSeparator(); -#ifdef HAVE_FFTW - filter_menu->Append (IFMENU_FILTER_FFT, "&FFT"); - filter_menu->Append (IFMENU_FILTER_IFFT, "&IFFT"); - filter_menu->Append (IFMENU_FILTER_FOURIER, "F&ourier"); - filter_menu->Append (IFMENU_FILTER_INVERSE_FOURIER, "Inverse Fo&urier"); -#else - filter_menu->Append (IFMENU_FILTER_FOURIER, "&Fourier"); - filter_menu->Append (IFMENU_FILTER_INVERSE_FOURIER, "&Inverse Fourier"); -#endif - filter_menu->Append (IFMENU_FILTER_SHUFFLEFOURIERTONATURALORDER, "S&huffle Fourier to Natural Order"); - filter_menu->Append (IFMENU_FILTER_SHUFFLENATURALTOFOURIERORDER, "Shu&ffle Natural to Fourier Order"); - filter_menu->Append (IFMENU_FILTER_MAGNITUDE, "&Magnitude"); - filter_menu->Append (IFMENU_FILTER_PHASE, "&Phase"); - - wxMenu* image_menu = new wxMenu; - filter_menu->AppendSeparator(); - image_menu->Append (IFMENU_IMAGE_ADD, "&Add..."); - image_menu->Append (IFMENU_IMAGE_SUBTRACT, "&Subtract..."); - image_menu->Append (IFMENU_IMAGE_MULTIPLY, "&Multiply..."); - image_menu->Append (IFMENU_IMAGE_DIVIDE, "&Divide..."); - image_menu->AppendSeparator(); - image_menu->Append (IFMENU_IMAGE_SCALESIZE, "S&cale Size..."); - - wxMenu *plot_menu = new wxMenu; - plot_menu->Append (IFMENU_PLOT_ROW, "Plot &Row"); - plot_menu->Append (IFMENU_PLOT_COL, "Plot &Column"); - - wxMenu *compare_menu = new wxMenu; - compare_menu->Append (IFMENU_COMPARE_IMAGES, "Compare &Images..."); - compare_menu->Append (IFMENU_COMPARE_ROW, "Compare &Row"); - compare_menu->Append (IFMENU_COMPARE_COL, "Compare &Column"); - + + wxMenu* filter_menu = new wxMenu; + filter_menu->Append (IFMENU_FILTER_INVERTVALUES, "&Invert Values"); + filter_menu->Append (IFMENU_FILTER_SQUARE, "&Square"); + filter_menu->Append (IFMENU_FILTER_SQRT, "Square &Root"); + filter_menu->Append (IFMENU_FILTER_LOG, "&Log"); + filter_menu->Append (IFMENU_FILTER_EXP, "&Exp"); + filter_menu->AppendSeparator(); +#ifdef HAVE_FFTW + filter_menu->Append (IFMENU_FILTER_FFT, "&FFT"); + filter_menu->Append (IFMENU_FILTER_IFFT, "&IFFT"); + filter_menu->Append (IFMENU_FILTER_FOURIER, "F&ourier"); + filter_menu->Append (IFMENU_FILTER_INVERSE_FOURIER, "Inverse Fo&urier"); +#else + filter_menu->Append (IFMENU_FILTER_FOURIER, "&Fourier"); + filter_menu->Append (IFMENU_FILTER_INVERSE_FOURIER, "&Inverse Fourier"); +#endif + filter_menu->Append (IFMENU_FILTER_SHUFFLEFOURIERTONATURALORDER, "S&huffle Fourier to Natural Order"); + filter_menu->Append (IFMENU_FILTER_SHUFFLENATURALTOFOURIERORDER, "Shu&ffle Natural to Fourier Order"); + filter_menu->Append (IFMENU_FILTER_MAGNITUDE, "&Magnitude"); + filter_menu->Append (IFMENU_FILTER_PHASE, "&Phase"); + + wxMenu* image_menu = new wxMenu; + filter_menu->AppendSeparator(); + image_menu->Append (IFMENU_IMAGE_ADD, "&Add..."); + image_menu->Append (IFMENU_IMAGE_SUBTRACT, "&Subtract..."); + image_menu->Append (IFMENU_IMAGE_MULTIPLY, "&Multiply..."); + image_menu->Append (IFMENU_IMAGE_DIVIDE, "&Divide..."); + image_menu->AppendSeparator(); + image_menu->Append (IFMENU_IMAGE_SCALESIZE, "S&cale Size..."); + + wxMenu *plot_menu = new wxMenu; + plot_menu->Append (IFMENU_PLOT_ROW, "Plot &Row"); + plot_menu->Append (IFMENU_PLOT_COL, "Plot &Column"); + + wxMenu *compare_menu = new wxMenu; + compare_menu->Append (IFMENU_COMPARE_IMAGES, "Compare &Images..."); + compare_menu->Append (IFMENU_COMPARE_ROW, "Compare &Row"); + compare_menu->Append (IFMENU_COMPARE_COL, "Compare &Column"); + wxMenu *help_menu = new wxMenu; help_menu->Append(MAINMENU_HELP_ABOUT, "&About"); wxMenuBar *menu_bar = new wxMenuBar; menu_bar->Append(file_menu, "&File"); - menu_bar->Append(view_menu, "&View"); - menu_bar->Append(image_menu, "&Image"); - menu_bar->Append(filter_menu, "Fi<er"); - menu_bar->Append(plot_menu, "P&lot"); + menu_bar->Append(view_menu, "&View"); + menu_bar->Append(image_menu, "&Image"); + menu_bar->Append(filter_menu, "Fi<er"); + menu_bar->Append(plot_menu, "P&lot"); menu_bar->Append(compare_menu, "&Compare"); menu_bar->Append(help_menu, "&Help"); @@ -795,7 +795,7 @@ ImageFileView::OnCreate (wxDocument *doc, long WXUNUSED(flags) ) int x, y; // X requires a forced resize m_frame->GetSize(&x, &y); m_frame->SetSize(-1, -1, x, y); - m_frame->SetFocus(); + m_frame->SetFocus(); m_frame->Show(true); Activate(true); @@ -807,10 +807,10 @@ ImageFileView::OnDraw (wxDC* dc) { if (m_bitmap.Ok()) dc->DrawBitmap(m_bitmap, 0, 0, false); - - int xCursor, yCursor; - if (m_canvas->GetCurrentCursor (xCursor, yCursor)) - m_canvas->DrawRubberBandCursor (*dc, xCursor, yCursor); + + int xCursor, yCursor; + if (m_canvas->GetCurrentCursor (xCursor, yCursor)) + m_canvas->DrawRubberBandCursor (*dc, xCursor, yCursor); } @@ -851,7 +851,7 @@ ImageFileView::OnUpdate (wxView *WXUNUSED(sender), wxObject *WXUNUSED(hint) ) ySize = clamp (ySize, 0, 800); m_frame->SetClientSize (xSize, ySize); m_canvas->SetScrollbars(20, 20, nx/20, ny/20); - m_canvas->SetBackgroundColour(*wxWHITE); + m_canvas->SetBackgroundColour(*wxWHITE); } if (m_canvas) @@ -880,323 +880,323 @@ ImageFileView::OnClose (bool deleteWindow) } return true; } - -void -ImageFileView::OnExport (wxCommandEvent& event) -{ - ImageFile& rIF = dynamic_cast(GetDocument())->getImageFile(); - ImageFileArrayConst v = rIF.getArray(); - int nx = rIF.nx(); - int ny = rIF.ny(); - if (v != NULL && nx != 0 && ny != 0) { - if (! m_bMinSpecified || ! m_bMaxSpecified) { - double min, max; - rIF.getMinMax (min, max); - if (! m_bMinSpecified) - m_dMinPixel = min; - if (! m_bMaxSpecified) - m_dMaxPixel = max; - } - - DialogExportParameters dialogExport (m_frame, m_iDefaultExportFormatID); - if (dialogExport.ShowModal() == wxID_OK) { - wxString strFormatName (dialogExport.getFormatName ()); - m_iDefaultExportFormatID = ImageFile::convertFormatNameToID (strFormatName.c_str()); - - wxString strExt; - wxString strWildcard; - if (m_iDefaultExportFormatID == ImageFile::FORMAT_PGM || m_iDefaultExportFormatID == ImageFile::FORMAT_PGMASCII) { - strExt = ".pgm"; - strWildcard = "PGM Files (*.pgm)|*.pgm"; - } -#ifdef HAVE_PNG - else if (m_iDefaultExportFormatID == ImageFile::FORMAT_PNG || m_iDefaultExportFormatID == ImageFile::FORMAT_PNG16) { - strExt = ".png"; - strWildcard = "PNG Files (*.png)|*.png"; - } -#endif - - const wxString& strFilename = wxFileSelector (wxString("Export Filename"), wxString(""), - wxString(""), strExt, strWildcard, wxOVERWRITE_PROMPT | wxHIDE_READONLY | wxSAVE); - if (strFilename) { - rIF.exportImage (strFormatName.c_str(), strFilename.c_str(), 1, 1, m_dMinPixel, m_dMaxPixel); - *theApp->getLog() << "Exported file " << strFilename << "\n"; - } - } - } -} - -void -ImageFileView::OnScaleSize (wxCommandEvent& event) -{ - ImageFile& rIF = GetDocument()->getImageFile(); - unsigned int iOldNX = rIF.nx(); - unsigned int iOldNY = rIF.ny(); - - DialogGetXYSize dialogGetXYSize (m_frame, "Set New X & Y Dimensions", iOldNX, iOldNY); - if (dialogGetXYSize.ShowModal() == wxID_OK) { - unsigned int iNewNX = dialogGetXYSize.getXSize(); - unsigned int iNewNY = dialogGetXYSize.getYSize(); - std::ostringstream os; - os << "Scale Size from (" << iOldNX << "," << iOldNY << ") to (" << iNewNX << "," << iNewNY << ")"; - ImageFileDocument* pScaledDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT)); - if (! pScaledDoc) { - sys_error (ERR_SEVERE, "Unable to create image file"); - return; - } - ImageFile& rScaledIF = pScaledDoc->getImageFile(); - rScaledIF.setArraySize (iNewNX, iNewNY); - rScaledIF.labelAdd (os.str().c_str()); - rIF.scaleImage (rScaledIF); - *theApp->getLog() << os.str().c_str() << "\n"; - if (theApp->getSetModifyNewDocs()) - pScaledDoc->Modify(TRUE); - pScaledDoc->UpdateAllViews (this); - pScaledDoc->GetFirstView()->OnUpdate (this, NULL); - } -} - -void -ImageFileView::OnPlotRow (wxCommandEvent& event) -{ - int xCursor, yCursor; - if (! m_canvas->GetCurrentCursor (xCursor, yCursor)) { - wxMessageBox ("No row selected. Please use left mouse button on image to select column","Error"); - return; - } - - const ImageFile& rIF = dynamic_cast(GetDocument())->getImageFile(); - ImageFileArrayConst v = rIF.getArray(); - int nx = rIF.nx(); - int ny = rIF.ny(); - - if (v != NULL && yCursor < ny) { - double* pX = new double [nx]; - double* pY = new double [nx]; - for (int i = 0; i < nx; i++) { - pX[i] = i; - pY[i] = v[i][yCursor]; - } - PlotFileDocument* pPlotDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.plt", wxDOC_SILENT)); - if (! pPlotDoc) { - sys_error (ERR_SEVERE, "Internal error: unable to create Plot file"); - } else { - PlotFile& rPlotFile = pPlotDoc->getPlotFile(); - std::ostringstream os; - os << "Row " << yCursor; - std::string title("title "); - title += os.str(); - rPlotFile.addEzsetCommand (title.c_str()); - rPlotFile.addEzsetCommand ("xlabel Column"); - rPlotFile.addEzsetCommand ("ylabel Pixel Value"); - rPlotFile.addEzsetCommand ("lxfrac 0"); - rPlotFile.addEzsetCommand ("box"); - rPlotFile.addEzsetCommand ("grid"); - rPlotFile.setCurveSize (2, nx); - rPlotFile.addColumn (0, pX); - rPlotFile.addColumn (1, pY); - } - delete pX; - delete pY; - if (theApp->getSetModifyNewDocs()) - pPlotDoc->Modify(true); - pPlotDoc->UpdateAllViews(); - } -} - -void -ImageFileView::OnPlotCol (wxCommandEvent& event) -{ - int xCursor, yCursor; - if (! m_canvas->GetCurrentCursor (xCursor, yCursor)) { - wxMessageBox ("No column selected. Please use left mouse button on image to select column","Error"); - return; - } - - const ImageFile& rIF = dynamic_cast(GetDocument())->getImageFile(); - ImageFileArrayConst v = rIF.getArray(); - int nx = rIF.nx(); - int ny = rIF.ny(); - - if (v != NULL && xCursor < nx) { - double* pX = new double [ny]; - double* pY = new double [ny]; - for (int i = 0; i < ny; i++) { - pX[i] = i; - pY[i] = v[xCursor][i]; - } - PlotFileDocument* pPlotDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.plt", wxDOC_SILENT)); - if (! pPlotDoc) { - sys_error (ERR_SEVERE, "Internal error: unable to create Plot file"); - } else { - PlotFile& rPlotFile = pPlotDoc->getPlotFile(); - std::ostringstream os; - os << "Column " << xCursor; - std::string title("title "); - title += os.str(); - rPlotFile.addEzsetCommand (title.c_str()); - rPlotFile.addEzsetCommand ("xlabel Row"); - rPlotFile.addEzsetCommand ("ylabel Pixel Value"); - rPlotFile.addEzsetCommand ("lxfrac 0"); - rPlotFile.addEzsetCommand ("box"); - rPlotFile.addEzsetCommand ("grid"); - rPlotFile.setCurveSize (2, nx); - rPlotFile.addColumn (0, pX); - rPlotFile.addColumn (1, pY); - } - delete pX; - delete pY; - if (theApp->getSetModifyNewDocs()) - pPlotDoc->Modify(true); - pPlotDoc->UpdateAllViews(); - } -} - -void -ImageFileView::OnCompareCol (wxCommandEvent& event) -{ - int xCursor, yCursor; - if (! m_canvas->GetCurrentCursor (xCursor, yCursor)) { - wxMessageBox ("No column selected. Please use left mouse button on image to select column","Error"); - return; - } - - std::vector vecIFDoc; - theApp->getCompatibleImages (GetDocument(), vecIFDoc); - if (vecIFDoc.size() == 0) { - wxMessageBox ("No compatible images for Column Comparison", "Error"); - return; - } - DialogGetComparisonImage dialogGetCompare (m_frame, "Get Comparison Image", vecIFDoc, false); - - if (dialogGetCompare.ShowModal() == wxID_OK) { - ImageFileDocument* pCompareDoc = dialogGetCompare.getImageFileDocument(); - const ImageFile& rIF = GetDocument()->getImageFile(); - const ImageFile& rCompareIF = pCompareDoc->getImageFile(); - - ImageFileArrayConst v1 = rIF.getArray(); - ImageFileArrayConst v2 = rCompareIF.getArray(); - int nx = rIF.nx(); - int ny = rIF.ny(); - - if (v1 != NULL && xCursor < nx) { - double* pX = new double [ny]; - double* pY1 = new double [ny]; - double* pY2 = new double [ny]; - for (int i = 0; i < ny; i++) { - pX[i] = i; - pY1[i] = v1[xCursor][i]; - pY2[i] = v2[xCursor][i]; - } - PlotFileDocument* pPlotDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.plt", wxDOC_SILENT)); - if (! pPlotDoc) { - sys_error (ERR_SEVERE, "Internal error: unable to create Plot file"); - } else { - PlotFile& rPlotFile = pPlotDoc->getPlotFile(); - std::ostringstream os; - os << "Column " << xCursor << " Comparison"; - std::string title("title "); - title += os.str(); - rPlotFile.addEzsetCommand (title.c_str()); - rPlotFile.addEzsetCommand ("xlabel Row"); - rPlotFile.addEzsetCommand ("ylabel Pixel Value"); - rPlotFile.addEzsetCommand ("lxfrac 0"); - rPlotFile.addEzsetCommand ("curve 1"); - rPlotFile.addEzsetCommand ("color 2"); - rPlotFile.addEzsetCommand ("curve 2"); - rPlotFile.addEzsetCommand ("color 4"); - rPlotFile.addEzsetCommand ("dash 5"); - rPlotFile.addEzsetCommand ("box"); - rPlotFile.addEzsetCommand ("grid"); - rPlotFile.setCurveSize (3, nx); - rPlotFile.addColumn (0, pX); - rPlotFile.addColumn (1, pY1); - rPlotFile.addColumn (2, pY2); - } - delete pX; - delete pY1; - delete pY2; - if (theApp->getSetModifyNewDocs()) - pPlotDoc->Modify(true); - pPlotDoc->UpdateAllViews(); - } - } -} - -void -ImageFileView::OnCompareRow (wxCommandEvent& event) -{ - int xCursor, yCursor; - if (! m_canvas->GetCurrentCursor (xCursor, yCursor)) { - wxMessageBox ("No column selected. Please use left mouse button on image to select column","Error"); - return; - } - - std::vector vecIFDoc; - theApp->getCompatibleImages (GetDocument(), vecIFDoc); - - if (vecIFDoc.size() == 0) { - wxMessageBox ("No compatible images for Row Comparison", "Error"); - return; - } - - DialogGetComparisonImage dialogGetCompare (m_frame, "Get Comparison Image", vecIFDoc, false); - - if (dialogGetCompare.ShowModal() == wxID_OK) { - ImageFileDocument* pCompareDoc = dialogGetCompare.getImageFileDocument(); - const ImageFile& rIF = GetDocument()->getImageFile(); - const ImageFile& rCompareIF = pCompareDoc->getImageFile(); - - ImageFileArrayConst v1 = rIF.getArray(); - ImageFileArrayConst v2 = rCompareIF.getArray(); - int nx = rIF.nx(); - int ny = rIF.ny(); - - if (v1 != NULL && yCursor < ny) { - double* pX = new double [nx]; - double* pY1 = new double [nx]; - double* pY2 = new double [nx]; - for (int i = 0; i < nx; i++) { - pX[i] = i; - pY1[i] = v1[i][yCursor]; - pY2[i] = v2[i][yCursor]; - } - PlotFileDocument* pPlotDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.plt", wxDOC_SILENT)); - if (! pPlotDoc) { - sys_error (ERR_SEVERE, "Internal error: unable to create Plot file"); - } else { - PlotFile& rPlotFile = pPlotDoc->getPlotFile(); - std::ostringstream os; - os << "Row " << yCursor << " Comparison"; - std::string title("title "); - title += os.str(); - rPlotFile.addEzsetCommand (title.c_str()); - rPlotFile.addEzsetCommand ("xlabel Column"); - rPlotFile.addEzsetCommand ("ylabel Pixel Value"); - rPlotFile.addEzsetCommand ("lxfrac 0"); - rPlotFile.addEzsetCommand ("curve 1"); - rPlotFile.addEzsetCommand ("color 2"); - rPlotFile.addEzsetCommand ("curve 2"); - rPlotFile.addEzsetCommand ("color 4"); - rPlotFile.addEzsetCommand ("dash 5"); - rPlotFile.addEzsetCommand ("box"); - rPlotFile.addEzsetCommand ("grid"); - rPlotFile.setCurveSize (3, ny); - rPlotFile.addColumn (0, pX); - rPlotFile.addColumn (1, pY1); - rPlotFile.addColumn (2, pY2); - } - delete pX; - delete pY1; - delete pY2; - if (theApp->getSetModifyNewDocs()) - pPlotDoc->Modify(true); - pPlotDoc->UpdateAllViews(); - } - } -} - - + +void +ImageFileView::OnExport (wxCommandEvent& event) +{ + ImageFile& rIF = dynamic_cast(GetDocument())->getImageFile(); + ImageFileArrayConst v = rIF.getArray(); + int nx = rIF.nx(); + int ny = rIF.ny(); + if (v != NULL && nx != 0 && ny != 0) { + if (! m_bMinSpecified || ! m_bMaxSpecified) { + double min, max; + rIF.getMinMax (min, max); + if (! m_bMinSpecified) + m_dMinPixel = min; + if (! m_bMaxSpecified) + m_dMaxPixel = max; + } + + DialogExportParameters dialogExport (m_frame, m_iDefaultExportFormatID); + if (dialogExport.ShowModal() == wxID_OK) { + wxString strFormatName (dialogExport.getFormatName ()); + m_iDefaultExportFormatID = ImageFile::convertFormatNameToID (strFormatName.c_str()); + + wxString strExt; + wxString strWildcard; + if (m_iDefaultExportFormatID == ImageFile::FORMAT_PGM || m_iDefaultExportFormatID == ImageFile::FORMAT_PGMASCII) { + strExt = ".pgm"; + strWildcard = "PGM Files (*.pgm)|*.pgm"; + } +#ifdef HAVE_PNG + else if (m_iDefaultExportFormatID == ImageFile::FORMAT_PNG || m_iDefaultExportFormatID == ImageFile::FORMAT_PNG16) { + strExt = ".png"; + strWildcard = "PNG Files (*.png)|*.png"; + } +#endif + + const wxString& strFilename = wxFileSelector (wxString("Export Filename"), wxString(""), + wxString(""), strExt, strWildcard, wxOVERWRITE_PROMPT | wxHIDE_READONLY | wxSAVE); + if (strFilename) { + rIF.exportImage (strFormatName.c_str(), strFilename.c_str(), 1, 1, m_dMinPixel, m_dMaxPixel); + *theApp->getLog() << "Exported file " << strFilename << "\n"; + } + } + } +} + +void +ImageFileView::OnScaleSize (wxCommandEvent& event) +{ + ImageFile& rIF = GetDocument()->getImageFile(); + unsigned int iOldNX = rIF.nx(); + unsigned int iOldNY = rIF.ny(); + + DialogGetXYSize dialogGetXYSize (m_frame, "Set New X & Y Dimensions", iOldNX, iOldNY); + if (dialogGetXYSize.ShowModal() == wxID_OK) { + unsigned int iNewNX = dialogGetXYSize.getXSize(); + unsigned int iNewNY = dialogGetXYSize.getYSize(); + std::ostringstream os; + os << "Scale Size from (" << iOldNX << "," << iOldNY << ") to (" << iNewNX << "," << iNewNY << ")"; + ImageFileDocument* pScaledDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT)); + if (! pScaledDoc) { + sys_error (ERR_SEVERE, "Unable to create image file"); + return; + } + ImageFile& rScaledIF = pScaledDoc->getImageFile(); + rScaledIF.setArraySize (iNewNX, iNewNY); + rScaledIF.labelAdd (os.str().c_str()); + rIF.scaleImage (rScaledIF); + *theApp->getLog() << os.str().c_str() << "\n"; + if (theApp->getSetModifyNewDocs()) + pScaledDoc->Modify(TRUE); + pScaledDoc->UpdateAllViews (this); + pScaledDoc->GetFirstView()->OnUpdate (this, NULL); + } +} + +void +ImageFileView::OnPlotRow (wxCommandEvent& event) +{ + int xCursor, yCursor; + if (! m_canvas->GetCurrentCursor (xCursor, yCursor)) { + wxMessageBox ("No row selected. Please use left mouse button on image to select column","Error"); + return; + } + + const ImageFile& rIF = dynamic_cast(GetDocument())->getImageFile(); + ImageFileArrayConst v = rIF.getArray(); + int nx = rIF.nx(); + int ny = rIF.ny(); + + if (v != NULL && yCursor < ny) { + double* pX = new double [nx]; + double* pY = new double [nx]; + for (int i = 0; i < nx; i++) { + pX[i] = i; + pY[i] = v[i][yCursor]; + } + PlotFileDocument* pPlotDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.plt", wxDOC_SILENT)); + if (! pPlotDoc) { + sys_error (ERR_SEVERE, "Internal error: unable to create Plot file"); + } else { + PlotFile& rPlotFile = pPlotDoc->getPlotFile(); + std::ostringstream os; + os << "Row " << yCursor; + std::string title("title "); + title += os.str(); + rPlotFile.addEzsetCommand (title.c_str()); + rPlotFile.addEzsetCommand ("xlabel Column"); + rPlotFile.addEzsetCommand ("ylabel Pixel Value"); + rPlotFile.addEzsetCommand ("lxfrac 0"); + rPlotFile.addEzsetCommand ("box"); + rPlotFile.addEzsetCommand ("grid"); + rPlotFile.setCurveSize (2, nx); + rPlotFile.addColumn (0, pX); + rPlotFile.addColumn (1, pY); + } + delete pX; + delete pY; + if (theApp->getSetModifyNewDocs()) + pPlotDoc->Modify(true); + pPlotDoc->UpdateAllViews(); + } +} + +void +ImageFileView::OnPlotCol (wxCommandEvent& event) +{ + int xCursor, yCursor; + if (! m_canvas->GetCurrentCursor (xCursor, yCursor)) { + wxMessageBox ("No column selected. Please use left mouse button on image to select column","Error"); + return; + } + + const ImageFile& rIF = dynamic_cast(GetDocument())->getImageFile(); + ImageFileArrayConst v = rIF.getArray(); + int nx = rIF.nx(); + int ny = rIF.ny(); + + if (v != NULL && xCursor < nx) { + double* pX = new double [ny]; + double* pY = new double [ny]; + for (int i = 0; i < ny; i++) { + pX[i] = i; + pY[i] = v[xCursor][i]; + } + PlotFileDocument* pPlotDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.plt", wxDOC_SILENT)); + if (! pPlotDoc) { + sys_error (ERR_SEVERE, "Internal error: unable to create Plot file"); + } else { + PlotFile& rPlotFile = pPlotDoc->getPlotFile(); + std::ostringstream os; + os << "Column " << xCursor; + std::string title("title "); + title += os.str(); + rPlotFile.addEzsetCommand (title.c_str()); + rPlotFile.addEzsetCommand ("xlabel Row"); + rPlotFile.addEzsetCommand ("ylabel Pixel Value"); + rPlotFile.addEzsetCommand ("lxfrac 0"); + rPlotFile.addEzsetCommand ("box"); + rPlotFile.addEzsetCommand ("grid"); + rPlotFile.setCurveSize (2, nx); + rPlotFile.addColumn (0, pX); + rPlotFile.addColumn (1, pY); + } + delete pX; + delete pY; + if (theApp->getSetModifyNewDocs()) + pPlotDoc->Modify(true); + pPlotDoc->UpdateAllViews(); + } +} + +void +ImageFileView::OnCompareCol (wxCommandEvent& event) +{ + int xCursor, yCursor; + if (! m_canvas->GetCurrentCursor (xCursor, yCursor)) { + wxMessageBox ("No column selected. Please use left mouse button on image to select column","Error"); + return; + } + + std::vector vecIFDoc; + theApp->getCompatibleImages (GetDocument(), vecIFDoc); + if (vecIFDoc.size() == 0) { + wxMessageBox ("No compatible images for Column Comparison", "Error"); + return; + } + DialogGetComparisonImage dialogGetCompare (m_frame, "Get Comparison Image", vecIFDoc, false); + + if (dialogGetCompare.ShowModal() == wxID_OK) { + ImageFileDocument* pCompareDoc = dialogGetCompare.getImageFileDocument(); + const ImageFile& rIF = GetDocument()->getImageFile(); + const ImageFile& rCompareIF = pCompareDoc->getImageFile(); + + ImageFileArrayConst v1 = rIF.getArray(); + ImageFileArrayConst v2 = rCompareIF.getArray(); + int nx = rIF.nx(); + int ny = rIF.ny(); + + if (v1 != NULL && xCursor < nx) { + double* pX = new double [ny]; + double* pY1 = new double [ny]; + double* pY2 = new double [ny]; + for (int i = 0; i < ny; i++) { + pX[i] = i; + pY1[i] = v1[xCursor][i]; + pY2[i] = v2[xCursor][i]; + } + PlotFileDocument* pPlotDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.plt", wxDOC_SILENT)); + if (! pPlotDoc) { + sys_error (ERR_SEVERE, "Internal error: unable to create Plot file"); + } else { + PlotFile& rPlotFile = pPlotDoc->getPlotFile(); + std::ostringstream os; + os << "Column " << xCursor << " Comparison"; + std::string title("title "); + title += os.str(); + rPlotFile.addEzsetCommand (title.c_str()); + rPlotFile.addEzsetCommand ("xlabel Row"); + rPlotFile.addEzsetCommand ("ylabel Pixel Value"); + rPlotFile.addEzsetCommand ("lxfrac 0"); + rPlotFile.addEzsetCommand ("curve 1"); + rPlotFile.addEzsetCommand ("color 2"); + rPlotFile.addEzsetCommand ("curve 2"); + rPlotFile.addEzsetCommand ("color 4"); + rPlotFile.addEzsetCommand ("dash 5"); + rPlotFile.addEzsetCommand ("box"); + rPlotFile.addEzsetCommand ("grid"); + rPlotFile.setCurveSize (3, nx); + rPlotFile.addColumn (0, pX); + rPlotFile.addColumn (1, pY1); + rPlotFile.addColumn (2, pY2); + } + delete pX; + delete pY1; + delete pY2; + if (theApp->getSetModifyNewDocs()) + pPlotDoc->Modify(true); + pPlotDoc->UpdateAllViews(); + } + } +} + +void +ImageFileView::OnCompareRow (wxCommandEvent& event) +{ + int xCursor, yCursor; + if (! m_canvas->GetCurrentCursor (xCursor, yCursor)) { + wxMessageBox ("No column selected. Please use left mouse button on image to select column","Error"); + return; + } + + std::vector vecIFDoc; + theApp->getCompatibleImages (GetDocument(), vecIFDoc); + + if (vecIFDoc.size() == 0) { + wxMessageBox ("No compatible images for Row Comparison", "Error"); + return; + } + + DialogGetComparisonImage dialogGetCompare (m_frame, "Get Comparison Image", vecIFDoc, false); + + if (dialogGetCompare.ShowModal() == wxID_OK) { + ImageFileDocument* pCompareDoc = dialogGetCompare.getImageFileDocument(); + const ImageFile& rIF = GetDocument()->getImageFile(); + const ImageFile& rCompareIF = pCompareDoc->getImageFile(); + + ImageFileArrayConst v1 = rIF.getArray(); + ImageFileArrayConst v2 = rCompareIF.getArray(); + int nx = rIF.nx(); + int ny = rIF.ny(); + + if (v1 != NULL && yCursor < ny) { + double* pX = new double [nx]; + double* pY1 = new double [nx]; + double* pY2 = new double [nx]; + for (int i = 0; i < nx; i++) { + pX[i] = i; + pY1[i] = v1[i][yCursor]; + pY2[i] = v2[i][yCursor]; + } + PlotFileDocument* pPlotDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.plt", wxDOC_SILENT)); + if (! pPlotDoc) { + sys_error (ERR_SEVERE, "Internal error: unable to create Plot file"); + } else { + PlotFile& rPlotFile = pPlotDoc->getPlotFile(); + std::ostringstream os; + os << "Row " << yCursor << " Comparison"; + std::string title("title "); + title += os.str(); + rPlotFile.addEzsetCommand (title.c_str()); + rPlotFile.addEzsetCommand ("xlabel Column"); + rPlotFile.addEzsetCommand ("ylabel Pixel Value"); + rPlotFile.addEzsetCommand ("lxfrac 0"); + rPlotFile.addEzsetCommand ("curve 1"); + rPlotFile.addEzsetCommand ("color 2"); + rPlotFile.addEzsetCommand ("curve 2"); + rPlotFile.addEzsetCommand ("color 4"); + rPlotFile.addEzsetCommand ("dash 5"); + rPlotFile.addEzsetCommand ("box"); + rPlotFile.addEzsetCommand ("grid"); + rPlotFile.setCurveSize (3, ny); + rPlotFile.addColumn (0, pX); + rPlotFile.addColumn (1, pY1); + rPlotFile.addColumn (2, pY2); + } + delete pX; + delete pY1; + delete pY2; + if (theApp->getSetModifyNewDocs()) + pPlotDoc->Modify(true); + pPlotDoc->UpdateAllViews(); + } + } +} + + // PhantomCanvas PhantomCanvas::PhantomCanvas (PhantomView* v, wxFrame *frame, const wxPoint& pos, const wxSize& size, const long style) @@ -1246,14 +1246,14 @@ PhantomView::OnProperties (wxCommandEvent& event) const int idPhantom = GetDocument()->getPhantomID(); const wxString& namePhantom = GetDocument()->getPhantomName(); std::ostringstream os; - os << "Phantom " << namePhantom.c_str() << " (" << idPhantom << ")" << "\n"; - const Phantom& rPhantom = GetDocument()->getPhantom(); + os << "Phantom " << namePhantom.c_str() << " (" << idPhantom << ")" << "\n"; + const Phantom& rPhantom = GetDocument()->getPhantom(); rPhantom.printDefinitions (os); #if DEBUG rPhantom.print (os); #endif - *theApp->getLog() << os.str().c_str() << "\n"; - wxMessageBox (os.str().c_str(), "Phantom Properties"); + *theApp->getLog() << os.str().c_str() << "\n"; + wxMessageBox (os.str().c_str(), "Phantom Properties"); } @@ -1275,10 +1275,10 @@ PhantomView::OnProjections (wxCommandEvent& event) if (m_iDefaultNDet > 0 && m_iDefaultNView > 0 && sGeometry != "") { const Phantom& rPhantom = GetDocument()->getPhantom(); - ProjectionFileDocument* pProjectionDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.pj", wxDOC_SILENT)); - if (! pProjectionDoc) { - sys_error (ERR_SEVERE, "Unable to create projection document"); - return; + ProjectionFileDocument* pProjectionDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.pj", wxDOC_SILENT)); + if (! pProjectionDoc) { + sys_error (ERR_SEVERE, "Unable to create projection document"); + return; } Projections& rProj = pProjectionDoc->getProjections(); Scanner theScanner (rPhantom, sGeometry.c_str(), m_iDefaultNDet, m_iDefaultNView, m_iDefaultNSample, m_dDefaultRotation, m_dDefaultFocalLength, m_dDefaultFieldOfView); @@ -1289,7 +1289,7 @@ PhantomView::OnProjections (wxCommandEvent& event) rProj.initFromScanner (theScanner); m_dDefaultRotation /= PI; // convert back to PI units - Timer timer; + Timer timer; if (m_iDefaultTrace > Trace::TRACE_CONSOLE) { ProjectionsDialog dialogProjections (theScanner, rProj, rPhantom, m_iDefaultTrace, dynamic_cast(m_frame)); for (int iView = 0; iView < rProj.nView(); iView++) { @@ -1307,7 +1307,7 @@ PhantomView::OnProjections (wxCommandEvent& event) } } } else { - wxProgressDialog dlgProgress (wxString("Projection"), wxString("Projection Progress"), rProj.nView() + 1, m_frame, wxPD_CAN_ABORT); + wxProgressDialog dlgProgress (wxString("Projection"), wxString("Projection Progress"), rProj.nView() + 1, m_frame, wxPD_CAN_ABORT); for (int i = 0; i < rProj.nView(); i++) { theScanner.collectProjections (rProj, rPhantom, i, 1, true, m_iDefaultTrace); if (! dlgProgress.Update (i+1)) { @@ -1318,18 +1318,18 @@ PhantomView::OnProjections (wxCommandEvent& event) } std::ostringstream os; - 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(); + 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(); rProj.setCalcTime (timer.timerEnd()); rProj.setRemark (os.str()); *theApp->getLog() << os.str().c_str() << "\n"; m_frame->Lower(); ::wxYield(); - ProjectionFileView* projView = dynamic_cast(pProjectionDoc->GetFirstView()); - if (projView) { - projView->getFrame()->SetFocus(); - projView->OnUpdate (projView, NULL); - } + ProjectionFileView* projView = dynamic_cast(pProjectionDoc->GetFirstView()); + if (projView) { + projView->getFrame()->SetFocus(); + projView->OnUpdate (projView, NULL); + } if (wxView* pView = pProjectionDoc->GetFirstView()) { if (wxFrame* pFrame = pView->GetFrame()) { pFrame->SetFocus(); @@ -1338,7 +1338,7 @@ PhantomView::OnProjections (wxCommandEvent& event) theApp->getDocManager()->ActivateView (pView, true, false); } ::wxYield(); - if (theApp->getSetModifyNewDocs()) + if (theApp->getSetModifyNewDocs()) pProjectionDoc->Modify(true); pProjectionDoc->UpdateAllViews(this); } @@ -1354,20 +1354,20 @@ PhantomView::OnRasterize (wxCommandEvent& event) if (retVal == wxID_OK) { int xSize = dialogRaster.getXSize(); int ySize = dialogRaster.getYSize(); - int nSamples = dialogRaster.getNSamples(); + int nSamples = dialogRaster.getNSamples(); if (nSamples < 1) nSamples = 1; if (xSize > 0 && ySize > 0) { const Phantom& rPhantom = GetDocument()->getPhantom(); - ImageFileDocument* pRasterDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT)); - if (! pRasterDoc) { - sys_error (ERR_SEVERE, "Unable to create image file"); - return; + ImageFileDocument* pRasterDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT)); + if (! pRasterDoc) { + sys_error (ERR_SEVERE, "Unable to create image file"); + return; } ImageFile& imageFile = pRasterDoc->getImageFile(); imageFile.setArraySize (xSize, ySize); - wxProgressDialog dlgProgress (wxString("Rasterize"), wxString("Rasterization Progress"), imageFile.nx() + 1, m_frame, wxPD_CAN_ABORT); + wxProgressDialog dlgProgress (wxString("Rasterize"), wxString("Rasterization Progress"), imageFile.nx() + 1, m_frame, wxPD_CAN_ABORT); Timer timer; for (unsigned int i = 0; i < imageFile.nx(); i++) { rPhantom.convertToImagefile (imageFile, nSamples, Trace::TRACE_NONE, i, 1, true); @@ -1376,20 +1376,20 @@ PhantomView::OnRasterize (wxCommandEvent& event) return; } } - if (theApp->getSetModifyNewDocs()) + if (theApp->getSetModifyNewDocs()) pRasterDoc->Modify(true); - pRasterDoc->UpdateAllViews(this); - std::ostringstream os; - os << "Rasterize Phantom " << rPhantom.name() << ": XSize=" << xSize << ", YSize=" << ySize << ", nSamples=" << nSamples; - *theApp->getLog() << os.str().c_str() << "\n"; + pRasterDoc->UpdateAllViews(this); + std::ostringstream os; + os << "Rasterize Phantom " << rPhantom.name() << ": XSize=" << xSize << ", YSize=" << ySize << ", nSamples=" << nSamples; + *theApp->getLog() << os.str().c_str() << "\n"; imageFile.labelAdd (os.str().c_str(), timer.timerEnd()); - ImageFileView* rasterView = dynamic_cast(pRasterDoc->GetFirstView()); - if (rasterView) { - rasterView->getFrame()->SetFocus(); - rasterView->OnUpdate (rasterView, NULL); - } + ImageFileView* rasterView = dynamic_cast(pRasterDoc->GetFirstView()); + if (rasterView) { + rasterView->getFrame()->SetFocus(); + rasterView->OnUpdate (rasterView, NULL); + } - } + } } } @@ -1417,9 +1417,9 @@ PhantomView::CreateChildFrame(wxDocument *doc, wxView *view) wxMenu *file_menu = new wxMenu; file_menu->Append(MAINMENU_FILE_CREATE_PHANTOM, "Cr&eate Phantom..."); - file_menu->Append(MAINMENU_FILE_CREATE_FILTER, "Create &Filter..."); - file_menu->Append(wxID_OPEN, "&Open..."); - file_menu->Append(wxID_SAVEAS, "Save &As..."); + file_menu->Append(MAINMENU_FILE_CREATE_FILTER, "Create &Filter..."); + file_menu->Append(wxID_OPEN, "&Open..."); + file_menu->Append(wxID_SAVEAS, "Save &As..."); file_menu->Append(wxID_CLOSE, "&Close"); file_menu->AppendSeparator(); @@ -1603,10 +1603,10 @@ ProjectionFileView::OnReconstruct (wxCommandEvent& event) m_iDefaultBackprojector = Backprojector::convertBackprojectNameToID (optBackprojectName.c_str()); m_iDefaultTrace = dialogReconstruction.getTrace(); if (m_iDefaultNX > 0 && m_iDefaultNY > 0) { - ImageFileDocument* pReconDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT)); - if (! pReconDoc) { - sys_error (ERR_SEVERE, "Unable to create image file"); - return; + ImageFileDocument* pReconDoc = dynamic_cast(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT)); + if (! pReconDoc) { + sys_error (ERR_SEVERE, "Unable to create image file"); + return; } ImageFile& imageFile = pReconDoc->getImageFile(); const Projections& rProj = GetDocument()->getProjections(); @@ -1616,9 +1616,9 @@ ProjectionFileView::OnReconstruct (wxCommandEvent& event) 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); return; } - + 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); - + Timer timerRecon; if (m_iDefaultTrace > Trace::TRACE_CONSOLE) { ReconstructDialog* pDlgReconstruct = new ReconstructDialog (*pReconstruct, rProj, imageFile, m_iDefaultTrace, m_frame); @@ -1651,14 +1651,14 @@ ProjectionFileView::OnReconstruct (wxCommandEvent& event) } } delete pReconstruct; - if (theApp->getSetModifyNewDocs()) + if (theApp->getSetModifyNewDocs()) pReconDoc->Modify(true); pReconDoc->UpdateAllViews(this); - ImageFileView* rasterView = dynamic_cast(pReconDoc->GetFirstView()); - if (rasterView) { - rasterView->getFrame()->SetFocus(); - rasterView->OnUpdate (rasterView, NULL); - } + ImageFileView* rasterView = dynamic_cast(pReconDoc->GetFirstView()); + if (rasterView) { + rasterView->getFrame()->SetFocus(); + rasterView->OnUpdate (rasterView, NULL); + } std::ostringstream os; 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(); *theApp->getLog() << os.str().c_str() << "\n"; @@ -1693,7 +1693,7 @@ ProjectionFileView::CreateChildFrame(wxDocument *doc, wxView *view) wxMenu *file_menu = new wxMenu; file_menu->Append(MAINMENU_FILE_CREATE_PHANTOM, "Cr&eate Phantom..."); - file_menu->Append(MAINMENU_FILE_CREATE_FILTER, "Create &Filter..."); + file_menu->Append(MAINMENU_FILE_CREATE_FILTER, "Create &Filter..."); file_menu->Append(wxID_OPEN, "&Open..."); file_menu->Append(wxID_SAVE, "&Save"); file_menu->Append(wxID_SAVEAS, "Save &As..."); @@ -1846,7 +1846,7 @@ PlotFileCanvas::OnDraw(wxDC& dc) { if (m_pView) m_pView->OnDraw(& dc); -} +} // PlotFileView @@ -1855,8 +1855,8 @@ IMPLEMENT_DYNAMIC_CLASS(PlotFileView, wxView) BEGIN_EVENT_TABLE(PlotFileView, wxView) EVT_MENU(PJMENU_FILE_PROPERTIES, PlotFileView::OnProperties) -EVT_MENU(PLOTMENU_VIEW_SCALE_MINMAX, PlotFileView::OnScaleMinMax) -EVT_MENU(PLOTMENU_VIEW_SCALE_AUTO, PlotFileView::OnScaleAuto) +EVT_MENU(PLOTMENU_VIEW_SCALE_MINMAX, PlotFileView::OnScaleMinMax) +EVT_MENU(PLOTMENU_VIEW_SCALE_AUTO, PlotFileView::OnScaleAuto) END_EVENT_TABLE() PlotFileView::PlotFileView(void) @@ -1865,8 +1865,8 @@ PlotFileView::PlotFileView(void) } PlotFileView::~PlotFileView(void) -{ - if (m_pEZPlot) +{ + if (m_pEZPlot) delete m_pEZPlot; } @@ -1874,8 +1874,8 @@ void PlotFileView::OnProperties (wxCommandEvent& event) { const PlotFile& rPlot = GetDocument()->getPlotFile(); - std::ostringstream os; - os << "Columns: " << rPlot.getNumColumns() << ", Records: " << rPlot.getNumRecords() << "\n"; + std::ostringstream os; + os << "Columns: " << rPlot.getNumColumns() << ", Records: " << rPlot.getNumRecords() << "\n"; rPlot.printHeaders (os); *theApp->getLog() << os.str().c_str(); wxMessageDialog dialogMsg (m_frame, os.str().c_str(), "Plot File Properties", wxOK | wxICON_INFORMATION); @@ -1883,58 +1883,58 @@ PlotFileView::OnProperties (wxCommandEvent& event) } -void -PlotFileView::OnScaleAuto (wxCommandEvent& event) -{ - const PlotFile& rPlotFile = GetDocument()->getPlotFile(); - double min, max, mean, mode, median, stddev; - rPlotFile.statistics (1, min, max, mean, mode, median, stddev); - DialogAutoScaleParameters dialogAutoScale (m_frame, mean, mode, median, stddev, m_dAutoScaleFactor); - int iRetVal = dialogAutoScale.ShowModal(); - if (iRetVal == wxID_OK) { - m_bMinSpecified = true; - m_bMaxSpecified = true; - double dMin, dMax; - if (dialogAutoScale.getMinMax (&dMin, &dMax)) { - m_dMinPixel = dMin; - m_dMaxPixel = dMax; - m_dAutoScaleFactor = dialogAutoScale.getAutoScaleFactor(); - OnUpdate (this, NULL); - } - } -} - -void -PlotFileView::OnScaleMinMax (wxCommandEvent& event) -{ - const PlotFile& rPlotFile = GetDocument()->getPlotFile(); - double min; - double max; - - if (! m_bMinSpecified || ! m_bMaxSpecified) { - if (! rPlotFile.getMinMax (1, min, max)) { - *theApp->getLog() << "Error: unable to find Min/Max\n"; - return; - } - } - - if (m_bMinSpecified) - min = m_dMinPixel; - if (m_bMaxSpecified) - max = m_dMaxPixel; - - DialogGetMinMax dialogMinMax (m_frame, "Set Y-axis Minimum & Maximum", min, max); - int retVal = dialogMinMax.ShowModal(); - if (retVal == wxID_OK) { - m_bMinSpecified = true; - m_bMaxSpecified = true; - m_dMinPixel = dialogMinMax.getMinimum(); - m_dMaxPixel = dialogMinMax.getMaximum(); - OnUpdate (this, NULL); - } -} - - +void +PlotFileView::OnScaleAuto (wxCommandEvent& event) +{ + const PlotFile& rPlotFile = GetDocument()->getPlotFile(); + double min, max, mean, mode, median, stddev; + rPlotFile.statistics (1, min, max, mean, mode, median, stddev); + DialogAutoScaleParameters dialogAutoScale (m_frame, mean, mode, median, stddev, m_dAutoScaleFactor); + int iRetVal = dialogAutoScale.ShowModal(); + if (iRetVal == wxID_OK) { + m_bMinSpecified = true; + m_bMaxSpecified = true; + double dMin, dMax; + if (dialogAutoScale.getMinMax (&dMin, &dMax)) { + m_dMinPixel = dMin; + m_dMaxPixel = dMax; + m_dAutoScaleFactor = dialogAutoScale.getAutoScaleFactor(); + OnUpdate (this, NULL); + } + } +} + +void +PlotFileView::OnScaleMinMax (wxCommandEvent& event) +{ + const PlotFile& rPlotFile = GetDocument()->getPlotFile(); + double min; + double max; + + if (! m_bMinSpecified || ! m_bMaxSpecified) { + if (! rPlotFile.getMinMax (1, min, max)) { + *theApp->getLog() << "Error: unable to find Min/Max\n"; + return; + } + } + + if (m_bMinSpecified) + min = m_dMinPixel; + if (m_bMaxSpecified) + max = m_dMaxPixel; + + DialogGetMinMax dialogMinMax (m_frame, "Set Y-axis Minimum & Maximum", min, max); + int retVal = dialogMinMax.ShowModal(); + if (retVal == wxID_OK) { + m_bMinSpecified = true; + m_bMaxSpecified = true; + m_dMinPixel = dialogMinMax.getMinimum(); + m_dMaxPixel = dialogMinMax.getMaximum(); + OnUpdate (this, NULL); + } +} + + PlotFileCanvas* PlotFileView::CreateCanvas (wxView *view, wxFrame *parent) { @@ -1958,7 +1958,7 @@ PlotFileView::CreateChildFrame(wxDocument *doc, wxView *view) wxMenu *file_menu = new wxMenu; file_menu->Append(MAINMENU_FILE_CREATE_PHANTOM, "Cr&eate Phantom..."); - file_menu->Append(MAINMENU_FILE_CREATE_FILTER, "Create &Filter..."); + file_menu->Append(MAINMENU_FILE_CREATE_FILTER, "Create &Filter..."); file_menu->Append(wxID_OPEN, "&Open..."); file_menu->Append(wxID_SAVE, "&Save"); file_menu->Append(wxID_SAVEAS, "Save &As..."); @@ -1972,10 +1972,10 @@ PlotFileView::CreateChildFrame(wxDocument *doc, wxView *view) file_menu->Append(wxID_PRINT_SETUP, "Print &Setup..."); file_menu->Append(wxID_PREVIEW, "Print Pre&view"); - wxMenu *view_menu = new wxMenu; - view_menu->Append(PLOTMENU_VIEW_SCALE_MINMAX, "Display Scale &Set..."); - view_menu->Append(PLOTMENU_VIEW_SCALE_AUTO, "Display Scale &Auto..."); - + wxMenu *view_menu = new wxMenu; + view_menu->Append(PLOTMENU_VIEW_SCALE_MINMAX, "Display Scale &Set..."); + view_menu->Append(PLOTMENU_VIEW_SCALE_AUTO, "Display Scale &Auto..."); + wxMenu *help_menu = new wxMenu; help_menu->Append(MAINMENU_HELP_CONTENTS, "&Contents"); help_menu->AppendSeparator(); @@ -1983,7 +1983,7 @@ PlotFileView::CreateChildFrame(wxDocument *doc, wxView *view) wxMenuBar *menu_bar = new wxMenuBar; - menu_bar->Append(file_menu, "&File"); + menu_bar->Append(file_menu, "&File"); menu_bar->Append(view_menu, "&View"); menu_bar->Append(help_menu, "&Help"); @@ -2001,10 +2001,10 @@ PlotFileView::OnCreate (wxDocument *doc, long WXUNUSED(flags) ) m_frame = CreateChildFrame(doc, this); SetFrame(m_frame); - m_bMinSpecified = false; - m_bMaxSpecified = false; - m_dAutoScaleFactor = 1.; - + m_bMinSpecified = false; + m_bMaxSpecified = false; + m_dAutoScaleFactor = 1.; + int width, height; m_frame->GetClientSize(&width, &height); m_frame->SetTitle ("Plot File"); @@ -2025,66 +2025,66 @@ PlotFileView::OnCreate (wxDocument *doc, long WXUNUSED(flags) ) void PlotFileView::OnDraw (wxDC* dc) { - const PlotFile& rPlotFile = GetDocument()->getPlotFile(); - const int iNColumns = rPlotFile.getNumColumns(); - const int iNRecords = rPlotFile.getNumRecords(); - - if (iNColumns > 0 && iNRecords > 0) { - int xsize, ysize; - m_canvas->GetClientSize (&xsize, &ysize); - SGPDriver driver (dc, xsize, ysize); - SGP sgp (driver); - if (m_pEZPlot) - m_pEZPlot->plot (&sgp); - } + const PlotFile& rPlotFile = GetDocument()->getPlotFile(); + const int iNColumns = rPlotFile.getNumColumns(); + const int iNRecords = rPlotFile.getNumRecords(); + + if (iNColumns > 0 && iNRecords > 0) { + int xsize, ysize; + m_canvas->GetClientSize (&xsize, &ysize); + SGPDriver driver (dc, xsize, ysize); + SGP sgp (driver); + if (m_pEZPlot) + m_pEZPlot->plot (&sgp); + } } void PlotFileView::OnUpdate (wxView *WXUNUSED(sender), wxObject *WXUNUSED(hint) ) { - const PlotFile& rPlotFile = GetDocument()->getPlotFile(); - const int iNColumns = rPlotFile.getNumColumns(); - const int iNRecords = rPlotFile.getNumRecords(); - - if (iNColumns > 0 && iNRecords > 0) { - if (m_pEZPlot) - delete m_pEZPlot; - m_pEZPlot = new EZPlot; - - for (unsigned int iEzset = 0; iEzset < rPlotFile.getNumEzsetCommands(); iEzset++) - m_pEZPlot->ezset (rPlotFile.getEzsetCommand (iEzset)); - - if (m_bMinSpecified) { - std::ostringstream os; - os << "ymin " << m_dMinPixel; - m_pEZPlot->ezset (os.str()); - } - - if (m_bMaxSpecified) { - std::ostringstream os; - os << "ymax " << m_dMaxPixel; - m_pEZPlot->ezset (os.str()); - } - - m_pEZPlot->ezset("box"); - m_pEZPlot->ezset("grid"); - - double* pdXaxis = new double [iNRecords]; - rPlotFile.getColumn (0, pdXaxis); - - double* pdY = new double [iNRecords]; - for (int iCol = 1; iCol < iNColumns; iCol++) { - rPlotFile.getColumn (iCol, pdY); - m_pEZPlot->addCurve (pdXaxis, pdY, iNRecords); - } - - delete pdXaxis; - delete pdY; - } - - if (m_canvas) - m_canvas->Refresh(); + const PlotFile& rPlotFile = GetDocument()->getPlotFile(); + const int iNColumns = rPlotFile.getNumColumns(); + const int iNRecords = rPlotFile.getNumRecords(); + + if (iNColumns > 0 && iNRecords > 0) { + if (m_pEZPlot) + delete m_pEZPlot; + m_pEZPlot = new EZPlot; + + for (unsigned int iEzset = 0; iEzset < rPlotFile.getNumEzsetCommands(); iEzset++) + m_pEZPlot->ezset (rPlotFile.getEzsetCommand (iEzset)); + + if (m_bMinSpecified) { + std::ostringstream os; + os << "ymin " << m_dMinPixel; + m_pEZPlot->ezset (os.str()); + } + + if (m_bMaxSpecified) { + std::ostringstream os; + os << "ymax " << m_dMaxPixel; + m_pEZPlot->ezset (os.str()); + } + + m_pEZPlot->ezset("box"); + m_pEZPlot->ezset("grid"); + + double* pdXaxis = new double [iNRecords]; + rPlotFile.getColumn (0, pdXaxis); + + double* pdY = new double [iNRecords]; + for (int iCol = 1; iCol < iNColumns; iCol++) { + rPlotFile.getColumn (iCol, pdY); + m_pEZPlot->addCurve (pdXaxis, pdY, iNRecords); + } + + delete pdXaxis; + delete pdY; + } + + if (m_canvas) + m_canvas->Refresh(); } bool -- 2.34.1