X-Git-Url: http://git.kpe.io/?p=ctsim.git;a=blobdiff_plain;f=src%2Fviews.cpp;h=04adab9fba869a643254e697eae5e55609690544;hp=d40fd5edd988240a2f38ff43d1fda38b4a95f770;hb=fd1d136a94a6d20013f38d6a997bdfefad0f5e98;hpb=efc398029330f4d9170257db29e4420e50f25625

diff --git a/src/views.cpp b/src/views.cpp
index d40fd5e..04adab9 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.33 2000/12/20 14:52:30 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
@@ -90,7 +90,7 @@ ImageFileCanvas::DrawRubberBandCursor (wxDC& dc, int x, int y)
 	ImageFileArrayConst v = rIF.getArray();
 	int nx = rIF.nx();
 	int ny = rIF.ny();
-
+	
 	dc.SetLogicalFunction (wxINVERT);
 	dc.SetPen (*wxGREEN_PEN);
 	dc.DrawLine (0, y, nx, y);
@@ -103,7 +103,7 @@ ImageFileCanvas::GetCurrentCursor (int& x, int& y)
 {
 	x = m_xCursor;
 	y = m_yCursor;
-
+	
 	if (m_xCursor >= 0 && m_yCursor >= 0)
 		return true;
 	else
@@ -120,7 +120,7 @@ ImageFileCanvas::OnMouseEvent(wxMouseEvent& event)
     PrepareDC(dc);
     
     wxPoint pt(event.GetLogicalPosition(dc));
-
+	
     if (event.RightIsDown()) {
 		const ImageFile& rIF = m_pView->GetDocument()->getImageFile();
 		ImageFileArrayConst v = rIF.getArray();
@@ -134,7 +134,7 @@ ImageFileCanvas::OnMouseEvent(wxMouseEvent& event)
 		} else
 			*theApp->getLog() << "Mouse out of image range (" << pt.x << "," << pt.y << ")\n";
     }
-    else if (event.LeftIsDown()) {
+    else if (event.LeftIsDown() || event.LeftUp()) {
 		const ImageFile& rIF = m_pView->GetDocument()->getImageFile();
 		ImageFileArrayConst v = rIF.getArray();
 		int nx = rIF.nx();
@@ -147,11 +147,13 @@ ImageFileCanvas::OnMouseEvent(wxMouseEvent& event)
 			DrawRubberBandCursor (dc, pt.x, pt.y);
 			m_xCursor = pt.x;
 			m_yCursor = pt.y;
+		} else
+			*theApp->getLog() << "Mouse out of image range (" << pt.x << "," << pt.y << ")\n";
+	}
+    if (event.LeftUp()) {
 			std::ostringstream os;
 			os << "Selected column" << pt.x << " and row" << pt.y << "\n";
 			*theApp->getLog() << os.str().c_str();
-		} else
-			*theApp->getLog() << "Mouse out of image range (" << pt.x << "," << pt.y << ")\n";
     }
 }
 
@@ -198,17 +200,20 @@ void
 ImageFileView::OnScaleAuto (wxCommandEvent& event)
 {
     const ImageFile& rIF = GetDocument()->getImageFile();
-    DialogAutoScaleParameters dialogAutoScale (m_frame, rIF, m_dAutoScaleFactor);
+    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;
-		dialogAutoScale.getMinMax (&dMin, &dMax);
-		m_dMinPixel = dMin;
-		m_dMaxPixel = dMax;
-		m_dAutoScaleFactor = dialogAutoScale.getAutoScaleFactor();
-		OnUpdate (this, NULL);
+		double dMin, dMax;
+		if (dialogAutoScale.getMinMax (&dMin, &dMax)) {
+			m_dMinPixel = dMin;
+			m_dMaxPixel = dMax;
+			m_dAutoScaleFactor = dialogAutoScale.getAutoScaleFactor();
+			OnUpdate (this, NULL);
+		}
     }
 }
 
@@ -225,7 +230,7 @@ ImageFileView::OnScaleMinMax (wxCommandEvent& event)
     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;
@@ -281,7 +286,7 @@ ImageFileView::CreateChildFrame(wxDocument *doc, wxView *view)
 	wxMenu *plot_menu = new wxMenu;
 	plot_menu->Append (IFMENU_PLOT_ROW, "Plot &Row");
 	plot_menu->Append (IFMENU_PLOT_COL, "Plot &Column");
-
+	
     wxMenu *help_menu = new wxMenu;
     help_menu->Append(MAINMENU_HELP_ABOUT, "&About");
     
@@ -330,7 +335,7 @@ 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);
@@ -414,12 +419,12 @@ ImageFileView::OnPlotRow (wxCommandEvent& event)
 		*theApp->getLog() << "No row selected. Please use left mouse button on image to select row\n";
 		return;
 	}
-
+	
     const ImageFile& rIF = dynamic_cast<ImageFileDocument*>(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];
@@ -431,20 +436,21 @@ ImageFileView::OnPlotRow (wxCommandEvent& event)
 		if (! pPlotDoc) {
 			sys_error (ERR_SEVERE, "Internal error: unable to create Plot file");
 		} else {
-			PlotFile& rPlot = pPlotDoc->getPlotFile();
+			PlotFile& rPlotFile = pPlotDoc->getPlotFile();
 			std::ostringstream title;
 			title << "Row " << yCursor;
-			rPlot.setTitle(title.str());
-			rPlot.setXLabel("Column");
-			rPlot.setYLabel("Pixel Value");
-			rPlot.setCurveSize (2, nx);
-			rPlot.addColumn (0, pX);
-			rPlot.addColumn (1, pY);
+			rPlotFile.setTitle(title.str());
+			rPlotFile.setXLabel("Column");
+			rPlotFile.setYLabel("Pixel Value");
+			rPlotFile.setCurveSize (2, nx);
+			rPlotFile.addColumn (0, pX);
+			rPlotFile.addColumn (1, pY);
 		}
 		delete pX;
 		delete pY;
+		pPlotDoc->Modify(true);
 	}
-
+	
 }
 
 void
@@ -455,12 +461,12 @@ ImageFileView::OnPlotCol (wxCommandEvent& event)
 		// os << "No column selected. Please use left mouse button on image to select column\n";
 		return;
 	}
-
+	
     const ImageFile& rIF = dynamic_cast<ImageFileDocument*>(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];
@@ -479,18 +485,19 @@ ImageFileView::OnPlotCol (wxCommandEvent& event)
 		if (! pPlotDoc) {
 			sys_error (ERR_SEVERE, "Internal error: unable to create Plot file");
 		} else {
-			PlotFile& rPlot = pPlotDoc->getPlotFile();
+			PlotFile& rPlotFile = pPlotDoc->getPlotFile();
 			std::ostringstream title;
 			title << "Column " << xCursor;
-			rPlot.setTitle(title.str());
-			rPlot.setXLabel("Row");
-			rPlot.setYLabel("Pixel Value");
-			rPlot.setCurveSize (2, nx);
-			rPlot.addColumn (0, pX);
-			rPlot.addColumn (1, pY);
+			rPlotFile.setTitle(title.str());
+			rPlotFile.setXLabel("Row");
+			rPlotFile.setYLabel("Pixel Value");
+			rPlotFile.setCurveSize (2, nx);
+			rPlotFile.addColumn (0, pX);
+			rPlotFile.addColumn (1, pY);
 		}
 		delete pX;
 		delete pY;
+		pPlotDoc->Modify(true);
 	}
 }
 
@@ -571,7 +578,11 @@ PhantomView::OnProjections (wxCommandEvent& event)
 		
 		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));
+			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);
 			if (theScanner.fail()) {
@@ -648,7 +659,11 @@ PhantomView::OnRasterize (wxCommandEvent& event)
 			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));
+			if (! pRasterDoc) {
+				sys_error (ERR_SEVERE, "Unable to create image file");
+				return;
+			}
 			ImageFile& imageFile = pRasterDoc->getImageFile();
 			
 			imageFile.setArraySize (xSize, ySize);
@@ -883,7 +898,11 @@ 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<ImageFileDocument*>(theApp->getDocManager()->CreateDocument("untitled.if", wxDOC_SILENT));
+			ImageFileDocument* pReconDoc = dynamic_cast<ImageFileDocument*>(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();
 			imageFile.setArraySize (m_iDefaultNX, m_iDefaultNY);
@@ -1135,7 +1154,8 @@ PlotFileCanvas::OnDraw(wxDC& dc)
 {
     if (m_pView)
         m_pView->OnDraw(& dc);
-}
+}
+
 
 // PlotFileView
 
@@ -1143,10 +1163,12 @@ 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)
 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)
 {
 }
 
@@ -1166,6 +1188,57 @@ 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, 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)
 {
@@ -1202,6 +1275,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 *help_menu = new wxMenu;
     help_menu->Append(MAINMENU_HELP_CONTENTS, "&Contents");
     help_menu->AppendSeparator();
@@ -1209,7 +1286,8 @@ 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");
     
     subframe->SetMenuBar(menu_bar);
@@ -1246,10 +1324,10 @@ PlotFileView::OnCreate(wxDocument *doc, long WXUNUSED(flags) )
 void 
 PlotFileView::OnDraw (wxDC* dc)
 {
-    const PlotFile& rPlot = GetDocument()->getPlotFile();
-    const int iNColumns = rPlot.getNumColumns();
-	const int iNRecords = rPlot.getNumRecords();
-
+    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);
@@ -1257,38 +1335,50 @@ PlotFileView::OnDraw (wxDC* dc)
 		SGP sgp (driver);
 		const PlotFile& rPhantom = GetDocument()->getPlotFile();
 		EZPlot plot (sgp);
-
-		if (! rPlot.getTitle().empty()) {
+		
+		if (! rPlotFile.getTitle().empty()) {
 			std::string s("title ");
-			s += rPlot.getTitle();
+			s += rPlotFile.getTitle();
 			plot.ezset (s);
 		}
-		if (! rPlot.getXLabel().empty()) {
+		if (! rPlotFile.getXLabel().empty()) {
 			std::string s("xlabel ");
-			s += rPlot.getXLabel();
+			s += rPlotFile.getXLabel();
 			plot.ezset (s);
 		}
-		if (! rPlot.getYLabel().empty()) {
+		if (! rPlotFile.getYLabel().empty()) {
 			std::string s("ylabel ");
-			s += rPlot.getYLabel();
+			s += rPlotFile.getYLabel();
 			plot.ezset (s);
 		}
+		
+		if (m_bMinSpecified) {
+			std::ostringstream os;
+			os << "ymin " << m_dMinPixel;
+			plot.ezset (os.str());
+		}
 
-	plot.ezset("box");
-	plot.ezset("grid");
+		if (m_bMaxSpecified) {
+			std::ostringstream os;
+			os << "ymax " << m_dMaxPixel;
+			plot.ezset (os.str());
+		}
 
+		plot.ezset("box");
+		plot.ezset("grid");
+		
 		double* pdXaxis = new double [iNRecords];
-		rPlot.getColumn (0, pdXaxis);
-
+		rPlotFile.getColumn (0, pdXaxis);
+		
 		double* pdY = new double [iNRecords];
 		for (int iCol = 1; iCol < iNColumns; iCol++) {
-			rPlot.getColumn (iCol, pdY);
+			rPlotFile.getColumn (iCol, pdY);
 			plot.addCurve (pdXaxis, pdY, iNRecords);
 		}
-
+		
 		delete pdXaxis;
 		delete pdY;
-
+		
 		plot.plot();
 	}
 }
@@ -1298,7 +1388,7 @@ void
 PlotFileView::OnUpdate(wxView *WXUNUSED(sender), wxObject *WXUNUSED(hint) )
 {
     if (m_canvas)
-      m_canvas->Refresh();
+		m_canvas->Refresh();
 }
 
 bool 
@@ -1312,14 +1402,14 @@ PlotFileView::OnClose (bool deleteWindow)
     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;
 }