** This is part of the CTSim program
** Copyright (C) 1983-2000 Kevin Rosenberg
**
-** $Id: views.cpp,v 1.32 2000/12/20 14:39:09 kevin Exp $
+** $Id: views.cpp,v 1.34 2000/12/20 20:08:48 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
ImageFileArrayConst v = rIF.getArray();\r
int nx = rIF.nx();\r
int ny = rIF.ny();\r
-\r
+ \r
dc.SetLogicalFunction (wxINVERT);\r
dc.SetPen (*wxGREEN_PEN);\r
dc.DrawLine (0, y, nx, y);\r
{\r
x = m_xCursor;\r
y = m_yCursor;\r
-\r
+ \r
if (m_xCursor >= 0 && m_yCursor >= 0)\r
return true;\r
else\r
PrepareDC(dc);
wxPoint pt(event.GetLogicalPosition(dc));
-\r
+ \r
if (event.RightIsDown()) {\r
const ImageFile& rIF = m_pView->GetDocument()->getImageFile();\r
ImageFileArrayConst v = rIF.getArray();\r
} else\r
*theApp->getLog() << "Mouse out of image range (" << pt.x << "," << pt.y << ")\n";\r
}\r
- else if (event.LeftIsDown()) {\r
+ else if (event.LeftIsDown() || event.LeftUp()) {\r
const ImageFile& rIF = m_pView->GetDocument()->getImageFile();\r
ImageFileArrayConst v = rIF.getArray();\r
int nx = rIF.nx();\r
DrawRubberBandCursor (dc, pt.x, pt.y);\r
m_xCursor = pt.x;\r
m_yCursor = pt.y;\r
+ } else\r
+ *theApp->getLog() << "Mouse out of image range (" << pt.x << "," << pt.y << ")\n";\r
+ }\r
+ if (event.LeftUp()) {\r
std::ostringstream os;\r
os << "Selected column" << pt.x << " and row" << pt.y << "\n";\r
*theApp->getLog() << os.str().c_str();\r
- } else\r
- *theApp->getLog() << "Mouse out of image range (" << pt.x << "," << pt.y << ")\n";\r
}\r
}
ImageFileView::OnScaleAuto (wxCommandEvent& event)
{
const ImageFile& rIF = GetDocument()->getImageFile();
- DialogAutoScaleParameters dialogAutoScale (m_frame, rIF, m_dAutoScaleFactor);
+ double min, max, mean, mode, median, stddev;\r
+ rIF.statistics(min, max, mean, mode, median, stddev);\r
+ 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;
- dialogAutoScale.getMinMax (&dMin, &dMax);
- m_dMinPixel = dMin;
- m_dMaxPixel = dMax;
- m_dAutoScaleFactor = dialogAutoScale.getAutoScaleFactor();
- OnUpdate (this, NULL);
+ double dMin, dMax;\r
+ if (dialogAutoScale.getMinMax (&dMin, &dMax)) {
+ m_dMinPixel = dMin;
+ m_dMaxPixel = dMax;
+ m_dAutoScaleFactor = dialogAutoScale.getAutoScaleFactor();
+ OnUpdate (this, NULL);\r
+ }
}
}
if (m_bMaxSpecified)
max = m_dMaxPixel;
- DialogGetImageMinMax dialogMinMax (m_frame, rIF, min, max);
+ DialogGetMinMax dialogMinMax (m_frame, "Set Image Minimum & Maximum", min, max);
int retVal = dialogMinMax.ShowModal();
if (retVal == wxID_OK) {
m_bMinSpecified = true;
wxMenu *plot_menu = new wxMenu;\r
plot_menu->Append (IFMENU_PLOT_ROW, "Plot &Row");\r
plot_menu->Append (IFMENU_PLOT_COL, "Plot &Column");\r
-
+
wxMenu *help_menu = new wxMenu;
help_menu->Append(MAINMENU_HELP_ABOUT, "&About");
{
if (m_bitmap.Ok())
dc->DrawBitmap(m_bitmap, 0, 0, false);
-\r
+ \r
int xCursor, yCursor;\r
if (m_canvas->GetCurrentCursor (xCursor, yCursor))\r
m_canvas->DrawRubberBandCursor (*dc, xCursor, yCursor);\r
*theApp->getLog() << "No row selected. Please use left mouse button on image to select row\n";\r
return;\r
}\r
-\r
+ \r
const ImageFile& rIF = dynamic_cast<ImageFileDocument*>(GetDocument())->getImageFile();\r
ImageFileArrayConst v = rIF.getArray();\r
int nx = rIF.nx();\r
int ny = rIF.ny();\r
-\r
+ \r
if (v != NULL && yCursor < ny) {\r
double* pX = new double [nx];\r
double* pY = new double [nx];\r
if (! pPlotDoc) {\r
sys_error (ERR_SEVERE, "Internal error: unable to create Plot file");\r
} else {\r
- PlotFile& rPlot = pPlotDoc->getPlotFile();\r
- rPlot.setTitle("Row Plot");\r
- rPlot.setXLabel("Column");\r
- rPlot.setYLabel("Pixel Value");\r
- rPlot.setCurveSize (2, nx);\r
- rPlot.addColumn (0, pX);\r
- rPlot.addColumn (1, pY);\r
- delete pX;\r
- delete pY;\r
+ PlotFile& rPlotFile = pPlotDoc->getPlotFile();\r
+ std::ostringstream title;\r
+ title << "Row " << yCursor;\r
+ rPlotFile.setTitle(title.str());\r
+ rPlotFile.setXLabel("Column");\r
+ rPlotFile.setYLabel("Pixel Value");\r
+ rPlotFile.setCurveSize (2, nx);\r
+ rPlotFile.addColumn (0, pX);\r
+ rPlotFile.addColumn (1, pY);\r
}\r
+ delete pX;\r
+ delete pY;\r
+ pPlotDoc->Modify(true);\r
}\r
-\r
+ \r
}\r
void\r
// os << "No column selected. Please use left mouse button on image to select column\n";\r
return;\r
}\r
-\r
+ \r
const ImageFile& rIF = dynamic_cast<ImageFileDocument*>(GetDocument())->getImageFile();\r
ImageFileArrayConst v = rIF.getArray();\r
int nx = rIF.nx();\r
int ny = rIF.ny();\r
-\r
+ \r
if (v != NULL && xCursor < nx) {\r
- double* pVec = new double [ny];\r
+ double* pX = new double [ny];\r
+ double* pY = new double [ny];\r
double minVal = v[xCursor][0];\r
double maxVal = minVal;\r
for (int i = 0; i < ny; i++) {\r
- if (minVal < v[xCursor][i])\r
- minVal = v[xCursor][i];\r
- else if (maxVal > v[xCursor][i])\r
- maxVal = v[xCursor][i];\r
- pVec[i] = v[xCursor][i];\r
+ double y = v[xCursor][i];\r
+ if (minVal < y)\r
+ minVal = y;\r
+ else if (maxVal > y)\r
+ maxVal = y;\r
+ pX[i] = i;\r
+ pY[i] = y;\r
+ }\r
+ PlotFileDocument* pPlotDoc = dynamic_cast<PlotFileDocument*>(theApp->getDocManager()->CreateDocument("untitled.plt", wxDOC_SILENT));\r
+ if (! pPlotDoc) {\r
+ sys_error (ERR_SEVERE, "Internal error: unable to create Plot file");\r
+ } else {\r
+ PlotFile& rPlotFile = pPlotDoc->getPlotFile();\r
+ std::ostringstream title;\r
+ title << "Column " << xCursor;\r
+ rPlotFile.setTitle(title.str());\r
+ rPlotFile.setXLabel("Row");\r
+ rPlotFile.setYLabel("Pixel Value");\r
+ rPlotFile.setCurveSize (2, nx);\r
+ rPlotFile.addColumn (0, pX);\r
+ rPlotFile.addColumn (1, pY);\r
}\r
+ delete pX;\r
+ delete pY;\r
+ pPlotDoc->Modify(true);\r
}\r
}\r
\r
if (m_iDefaultNDet > 0 && m_iDefaultNView > 0 && sGeometry != "") {
const Phantom& rPhantom = GetDocument()->getPhantom();
- ProjectionFileDocument* pProjectionDoc = dynamic_cast<ProjectionFileDocument*>(theApp->getDocManager()->CreateDocument("untitled.pj", wxDOC_SILENT));
+ ProjectionFileDocument* pProjectionDoc = dynamic_cast<ProjectionFileDocument*>(theApp->getDocManager()->CreateDocument("untitled.pj", wxDOC_SILENT));\r
+ if (! pProjectionDoc) {\r
+ sys_error (ERR_SEVERE, "Unable to create projection document");\r
+ return;\r
+ }
Projections& rProj = pProjectionDoc->getProjections();
Scanner theScanner (rPhantom, sGeometry.c_str(), m_iDefaultNDet, m_iDefaultNView, m_iDefaultNSample, m_dDefaultRotation, m_dDefaultFocalLength, m_dDefaultFieldOfView);
if (theScanner.fail()) {
nSamples = 1;
if (xSize > 0 && ySize > 0) {
const Phantom& rPhantom = GetDocument()->getPhantom();
- ImageFileDocument* pRasterDoc = dynamic_cast<ImageFileDocument*>(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT));
+ ImageFileDocument* pRasterDoc = dynamic_cast<ImageFileDocument*>(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT));\r
+ if (! pRasterDoc) {\r
+ sys_error (ERR_SEVERE, "Unable to create image file");\r
+ return;\r
+ }
ImageFile& imageFile = pRasterDoc->getImageFile();
imageFile.setArraySize (xSize, ySize);
m_iDefaultBackprojector = Backprojector::convertBackprojectNameToID (optBackprojectName.c_str());
m_iDefaultTrace = dialogReconstruction.getTrace();
if (m_iDefaultNX > 0 && m_iDefaultNY > 0) {
- ImageFileDocument* pReconDoc = dynamic_cast<ImageFileDocument*>(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT));
+ ImageFileDocument* pReconDoc = dynamic_cast<ImageFileDocument*>(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT));\r
+ if (pReconDoc) {\r
+ sys_error (ERR_SEVERE, "Unable to create image file");\r
+ return;\r
+ }
ImageFile& imageFile = pReconDoc->getImageFile();
const Projections& rProj = GetDocument()->getProjections();
imageFile.setArraySize (m_iDefaultNX, m_iDefaultNY);
{
if (m_pView)
m_pView->OnDraw(& dc);
-}
+}\r
+
// PlotFileView
BEGIN_EVENT_TABLE(PlotFileView, wxView)
EVT_MENU(PJMENU_FILE_PROPERTIES, PlotFileView::OnProperties)
+EVT_MENU(PLOTMENU_VIEW_SCALE_MINMAX, PlotFileView::OnScaleMinMax)\r
+EVT_MENU(PLOTMENU_VIEW_SCALE_AUTO, PlotFileView::OnScaleAuto)\r
END_EVENT_TABLE()
PlotFileView::PlotFileView(void)
-: wxView(), m_canvas(NULL), m_frame(NULL)
+: wxView(), m_canvas(NULL), m_frame(NULL), m_bMinSpecified(false), m_bMaxSpecified(false)
{
}
}
+void \r
+PlotFileView::OnScaleAuto (wxCommandEvent& event)\r
+{\r
+ const PlotFile& rPlotFile = GetDocument()->getPlotFile();\r
+ double min, max, mean, mode, median, stddev;\r
+ rPlotFile.statistics (1, min, max, mean, mode, median, stddev);\r
+ DialogAutoScaleParameters dialogAutoScale (m_frame, mean, mode, median, stddev, m_dAutoScaleFactor);\r
+ int iRetVal = dialogAutoScale.ShowModal();\r
+ if (iRetVal == wxID_OK) {\r
+ m_bMinSpecified = true;\r
+ m_bMaxSpecified = true;\r
+ double dMin, dMax;\r
+ if (dialogAutoScale.getMinMax (&dMin, &dMax)) {\r
+ m_dMinPixel = dMin;\r
+ m_dMaxPixel = dMax;\r
+ m_dAutoScaleFactor = dialogAutoScale.getAutoScaleFactor();\r
+ OnUpdate (this, NULL);\r
+ }\r
+ }\r
+}\r
+\r
+void \r
+PlotFileView::OnScaleMinMax (wxCommandEvent& event)\r
+{\r
+ const PlotFile& rPlotFile = GetDocument()->getPlotFile();\r
+ double min, max;\r
+\r
+ if (! m_bMinSpecified && ! m_bMaxSpecified) {\r
+ if (! rPlotFile.getMinMax (1, min, max)) {\r
+ *theApp->getLog() << "Error: unable to find Min/Max\n";\r
+ return;\r
+ }\r
+ }\r
+ \r
+ if (m_bMinSpecified)\r
+ min = m_dMinPixel;\r
+ if (m_bMaxSpecified)\r
+ max = m_dMaxPixel;\r
+ \r
+ DialogGetMinMax dialogMinMax (m_frame, "Set Y-axis Minimum & Maximum", min, max);\r
+ int retVal = dialogMinMax.ShowModal();\r
+ if (retVal == wxID_OK) {\r
+ m_bMinSpecified = true;\r
+ m_bMaxSpecified = true;\r
+ m_dMinPixel = dialogMinMax.getMinimum();\r
+ m_dMaxPixel = dialogMinMax.getMaximum();\r
+ OnUpdate (this, NULL);\r
+ }\r
+}\r
+\r
+\r
PlotFileCanvas*
PlotFileView::CreateCanvas (wxView *view, wxFrame *parent)
{
file_menu->Append(wxID_PRINT_SETUP, "Print &Setup...");
file_menu->Append(wxID_PREVIEW, "Print Pre&view");
+ wxMenu *view_menu = new wxMenu;\r
+ view_menu->Append(PLOTMENU_VIEW_SCALE_MINMAX, "Display Scale &Set...");\r
+ view_menu->Append(PLOTMENU_VIEW_SCALE_AUTO, "Display Scale &Auto...");\r
+ \r
wxMenu *help_menu = new wxMenu;
help_menu->Append(MAINMENU_HELP_CONTENTS, "&Contents");
help_menu->AppendSeparator();
wxMenuBar *menu_bar = new wxMenuBar;
- menu_bar->Append(file_menu, "&File");
+ menu_bar->Append(file_menu, "&File");\r
+ menu_bar->Append(view_menu, "&View");
menu_bar->Append(help_menu, "&Help");
subframe->SetMenuBar(menu_bar);
void
PlotFileView::OnDraw (wxDC* dc)
{
- const PlotFile& rPlot = GetDocument()->getPlotFile();\r
- const int iNColumns = rPlot.getNumColumns();\r
- const int iNRecords = rPlot.getNumRecords();\r
-\r
+ const PlotFile& rPlotFile = GetDocument()->getPlotFile();\r
+ const int iNColumns = rPlotFile.getNumColumns();\r
+ const int iNRecords = rPlotFile.getNumRecords();\r
+ \r
if (iNColumns > 0 && iNRecords > 0) {\r
int xsize, ysize;\r
m_canvas->GetClientSize (&xsize, &ysize);\r
SGP sgp (driver);\r
const PlotFile& rPhantom = GetDocument()->getPlotFile();\r
EZPlot plot (sgp);\r
-\r
- if (! rPlot.getTitle().empty()) {\r
+ \r
+ if (! rPlotFile.getTitle().empty()) {\r
std::string s("title ");\r
- s += rPlot.getTitle();\r
+ s += rPlotFile.getTitle();\r
plot.ezset (s);\r
}\r
- if (! rPlot.getXLabel().empty()) {\r
+ if (! rPlotFile.getXLabel().empty()) {\r
std::string s("xlabel ");\r
- s += rPlot.getXLabel();\r
+ s += rPlotFile.getXLabel();\r
plot.ezset (s);\r
}\r
- if (! rPlot.getYLabel().empty()) {\r
+ if (! rPlotFile.getYLabel().empty()) {\r
std::string s("ylabel ");\r
- s += rPlot.getYLabel();\r
+ s += rPlotFile.getYLabel();\r
plot.ezset (s);\r
}\r
+ \r
+ if (m_bMinSpecified) {\r
+ std::ostringstream os;\r
+ os << "ymin " << m_dMinPixel;\r
+ plot.ezset (os.str());\r
+ }\r
\r
- plot.ezset("box");\r
- plot.ezset("grid");\r
+ if (m_bMaxSpecified) {\r
+ std::ostringstream os;\r
+ os << "ymax " << m_dMaxPixel;\r
+ plot.ezset (os.str());\r
+ }\r
\r
+ plot.ezset("box");\r
+ plot.ezset("grid");\r
+ \r
double* pdXaxis = new double [iNRecords];\r
- rPlot.getColumn (0, pdXaxis);\r
-\r
+ rPlotFile.getColumn (0, pdXaxis);\r
+ \r
double* pdY = new double [iNRecords];\r
for (int iCol = 1; iCol < iNColumns; iCol++) {\r
- rPlot.getColumn (iCol, pdY);\r
+ rPlotFile.getColumn (iCol, pdY);\r
plot.addCurve (pdXaxis, pdY, iNRecords);\r
}\r
-\r
+ \r
delete pdXaxis;\r
delete pdY;\r
-\r
+ \r
plot.plot();\r
}\r
}
PlotFileView::OnUpdate(wxView *WXUNUSED(sender), wxObject *WXUNUSED(hint) )
{
if (m_canvas)
- m_canvas->Refresh();
+ m_canvas->Refresh();
}
bool
m_canvas = NULL;
wxString s(wxTheApp->GetAppName());
if (m_frame)
- m_frame->SetTitle(s);
+ m_frame->SetTitle(s);
SetFrame(NULL);
Activate(false);
if (deleteWindow) {
- delete m_frame;
- return true;
+ delete m_frame;
+ return true;
}
return true;
}