X-Git-Url: http://git.kpe.io/?a=blobdiff_plain;f=libctsim%2Fprojections.cpp;h=35d8418a070edd3029684246f8e8eea645d1ffff;hb=711cae0ee02e046370fdb4d6c6f440596ff71980;hp=a29d7c5a6ca967ab5e569349a46acc0ef2082797;hpb=e4c1f7f8eb87558c3abf3bf1d20732361f425351;p=ctsim.git

diff --git a/libctsim/projections.cpp b/libctsim/projections.cpp
index a29d7c5..35d8418 100644
--- a/libctsim/projections.cpp
+++ b/libctsim/projections.cpp
@@ -8,7 +8,7 @@
 **  This is part of the CTSim program
 **  Copyright (C) 1983-2000 Kevin Rosenberg
 **
-**  $Id: projections.cpp,v 1.15 2000/07/15 08:36:13 kevin Exp $
+**  $Id: projections.cpp,v 1.16 2000/07/28 08:28:08 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
@@ -515,7 +515,6 @@ Projections::reconstruct (ImageFile& im, const char* const filterName, double fi
     cout << "Reconstruct: filter="<<filterName<< ", interp="<<interpName<<", backproject="<<backprojectName<<endl;
 
 #if HAVE_SGP
-  SGP_ID gid;
   int nVecFilter = filter.getNFilterPoints();
   double plot_xaxis [nVecFilter]; 			// array for plotting 
 
@@ -523,14 +522,17 @@ Projections::reconstruct (ImageFile& im, const char* const filterName, double fi
     int i;
     double f;
     double filterInc = filter.getFilterIncrement();
-    for (i = 0, f = filter.getFilterMin(); i < nVecFilter; i++, f += filterInc) 
+    for (i = 0, f = filter.getFilterMin(); i < nVecFilter; i++, f += filterInc)
       plot_xaxis[i] = f;
 
     if (filter.getFilter()) {
-      gid = ezplot (plot_xaxis, filter.getFilter(), nVecFilter);
+      SGPDriver sgpDriver ("Filter Function");
+      SGP sgp (sgpDriver);
+      EZPlot ezplot (sgp);
+
+      ezplot.addCurve (plot_xaxis, filter.getFilter(), nVecFilter);
       cio_put_str ("Press any key to continue");
       cio_kb_getc ();
-      sgp2_close (gid);
     }
   }
   if (trace >= TRACE_TEXT) {
@@ -555,26 +557,30 @@ Projections::reconstruct (ImageFile& im, const char* const filterName, double fi
 
 #ifdef HAVE_SGP
     if (trace >= TRACE_PLOT)  {
-      ezset  ("clear.");
-      ezset  ("xticks major 5.");
-      ezset  ("xlabel ");
-      ezset  ("ylabel ");
-      ezset  ("xlength .5.");
-      ezset  ("box.");
-      ezset  ("grid.");
-      ezset  ("ufinish yes.");
-      ezplot (detval, plot_xaxis, m_nDet);
-      ezset  ("clear.");
-      ezset  ("xticks major 5.");
-      ezset  ("xlabel ");
-      ezset  ("ylabel ");
-      ezset  ("ustart yes.");
-      ezset  ("xporigin .5.");
-      ezset  ("xlength .5.");
-      ezset ("box");
-
-      ezset ("grid");
-      gid = ezplot (filteredProj, plot_xaxis, n_filteredProj);
+      SGPDriver sgpDriverProj ("Projection");
+      SGP sgpProj (sgpDriverProj);
+      EZPlot ezplotProj (sgpProj);
+
+      ezplotProj.ezset  ("clear");
+      ezplotProj.ezset  ("xticks major 5.");
+      ezplotProj.ezset  ("xlabel ");
+      ezplotProj.ezset  ("ylabel ");
+      ezplotProj.ezset  ("xlength .5.");
+      ezplotProj.ezset  ("box.");
+      ezplotProj.ezset  ("grid.");
+      ezplotProj.addCurve (detval, plot_xaxis, m_nDet);
+      ezplotProj.ezset  ("clear.");
+      ezplotProj.ezset  ("xticks major 5.");
+      ezplotProj.ezset  ("xlabel ");
+      ezplotProj.ezset  ("ylabel ");
+      ezplotProj.ezset  ("xporigin .5.");
+      ezplotProj.ezset  ("xlength .5.");
+      ezplotProj.ezset ("box");
+      ezplotProj.ezset ("grid");
+      ezplotProj.addCurve (filteredProj, plot_xaxis, n_filteredProj);
+      ezplotProj.plot ();
+      cout << "Press enter to continue\n";
+      cio_kb_getc();
     }
 #endif  //HAVE_SGP
 
@@ -597,7 +603,6 @@ Projections::reconstruct (ImageFile& im, const char* const filterName, double fi
       char str[256];
       printf ("Do you want to exit with current pic (y/n) -- ");
       fgets(str, sizeof(str), stdin);
-      sgp2_close (sgp2_get_active_win());
       if (tolower(str[0]) == 'y') {
 	break;
       }