X-Git-Url: http://git.kpe.io/?p=ctsim.git;a=blobdiff_plain;f=include%2Fscanner.h;h=7f8f95e7e2aafcb29a4e8323167b0c66bf91bfb6;hp=550573a68e6ecdaccfd662d4ac4a36777e889407;hb=a05f3cb550877e94aa118cc04b361c0c8fdb3dc3;hpb=93eb50659a64db3e37753d36c9923e82c7f7ca8a

diff --git a/include/scanner.h b/include/scanner.h
index 550573a..7f8f95e 100644
--- a/include/scanner.h
+++ b/include/scanner.h
@@ -9,7 +9,7 @@
 **  This is part of the CTSim program
 **  Copyright (C) 1983-2000 Kevin Rosenberg
 **
-**  $Id: scanner.h,v 1.11 2000/08/02 18:06:00 kevin Exp $
+**  $Id: scanner.h,v 1.14 2000/08/31 08:38:58 kevin Exp $
 **
 **
 **  This program is free software; you can redistribute it and/or modify
@@ -45,12 +45,12 @@ class DetectorArray
 {
  public:
   DetectorArray (const int ndet);
-  ~DetectorArray (void);
+  ~DetectorArray ();
 
-  const int nDet(void) const {return m_nDet;}
-  const double viewAngle(void) const {return m_viewAngle;}
-  DetectorValue* detValues(void) {return m_detValues;}
-  const DetectorValue* detValues(void) const {return m_detValues;}
+  const int nDet() const {return m_nDet;}
+  const double viewAngle() const {return m_viewAngle;}
+  DetectorValue* detValues() {return m_detValues;}
+  const DetectorValue* detValues() const {return m_detValues;}
 
   void setViewAngle (double viewAngle)
       { m_viewAngle = viewAngle; }
@@ -74,23 +74,28 @@ class Scanner
   static const int Scanner::GEOMETRY_EQUIANGULAR;
 
   
-  Scanner (const Phantom& phm, const char* const geometryName, int nDet, int nView, int nSample, const double rot_anglen);
+  Scanner (const Phantom& phm, const char* const geometryName, int nDet, int nView, int nSample, const double rot_anglen, double dFieldOfView, double dFocalLength);
   ~Scanner();
   
-  void collectProjections (Projections& proj, const Phantom& phm, const int start_view = 0, const int trace = TRACE_NONE, SGP* pSGP = NULL);
+  void collectProjections (Projections& proj, const Phantom& phm, const int trace = Trace::TRACE_NONE, SGP* pSGP = NULL);
 
-  void setNView (int nView);
+  void collectProjections (Projections& proj, const Phantom& phm, const int iStartView, const int iNumViews, bool bStoreAtViewPosition, const int trace = Trace::TRACE_NONE, SGP* pSGP = NULL);
 
-  const bool fail(void) const {return m_fail;}
-  const string& failMessage(void) const {return m_failMessage;}
-  const unsigned int nDet(void) const {return m_nDet;}
-  const unsigned int nView(void) const {return m_nView;}
-  const double phmLen(void) const {return m_phmLen;}
-  const double rotInc(void) const {return m_rotInc;}
-  const double detInc(void) const {return m_detInc;}
-  const double radius(void) const {return m_radius;}
+  void setNView (int nView);
 
-  static const int getGeometryCount() {return s_iGeometryCount;}
+  bool fail() const {return m_fail;}
+  const string& failMessage() const {return m_failMessage;}
+  unsigned int nDet() const {return m_nDet;}
+  unsigned int nView() const {return m_nView;}
+  double phmLen() const {return m_phmLen;}
+  double rotInc() const {return m_rotInc;}
+  double detInc() const {return m_detInc;}
+  double detLen() const {return m_detLen;}
+  double focalLength() const {return m_dFocalLength;}
+  double fieldOfView() const {return m_dFieldOfView;}
+  int geometry() const {return m_idGeometry;}
+
+  static int getGeometryCount() {return s_iGeometryCount;}
   static const char** getGeometryNameArray() {return s_aszGeometryName;}
   static const char** getGeometryTitleArray() {return s_aszGeometryTitle;}
   static int convertGeometryNameToID (const char* const geometryName);
@@ -104,46 +109,67 @@ class Scanner
   unsigned int m_nDet;		/* Number of detectors in array */
   unsigned int m_nView;		/* Number of rotated views */
   unsigned int m_nSample;	/* Number of rays per detector */
-  double m_detLen;		/* Total length of detector array */
-  double m_rotLen;		/* Rotation angle length in radians (norm 2PI) */
-  double m_detInc;		/* Increment between centers of detectors */
-  double m_rotInc;		/* Increment in rotation angle between views */
-  double m_radius;		/* Radius of rotation.  Distance from center of phm to center of det */
-  double m_phmLen;              /* Maximum Length of phantom or area of interest */
+  double m_dFieldOfView;        // Field of View
+  double m_dFocalLength;        // Focal Length
+  double m_dFieldOfViewRatio;   // Field of View Ratio to diameter phmLen
+  double m_dFocalLengthRatio;   // Focal Length ratio to radius phmLen
+  double m_detLen;		// Total length of detector array 
+  double m_rotLen;		// Rotation angle length in radians (norm 2PI)
+  double m_detInc;		// Increment between centers of detectors 
+  double m_rotInc;		// Increment in rotation angle between views 
+  double m_phmLen;              // Maximum Length of phantom or area of interest
+  double m_dXCenter;            // Center of Phantom
+  double m_dYCenter;            
+  double m_dAngularDetIncrement;
+  double m_dAngularDetLen;
+
   int m_trace;
   struct {
     double xd1,yd1,xd2,yd2;	/* Coordinates of detector endpoints */
     double xs1,ys1,xs2,ys2;	/* Coordinates of source endpoints */
     double angle;		/* Starting angle */
+    double dAngularDet;
   } m_initPos;
 
+  GRFMTX_2D m_rotmtxIncrement;
+
+#if HAVE_SGP
+  SGP* m_pSGP;                  // Pointer to graphics device
+  double m_dXMinWin;            // Extent of graphics window
+  double m_dXMaxWin;
+  double m_dYMinWin;
+  double m_dYMaxWin;
+  double m_dTextHeight;
+#endif
+
   static const char* s_aszGeometryName[];
   static const char* s_aszGeometryTitle[];
   static const int s_iGeometryCount;
-  static const int N_EXTRA_DETECTORS=0;    // Number of extra detectors widths when calculating detlen 
 
-  void projectSingleView (const Phantom& phm, DetectorArray& darray, const double xd1, const double yd1, const double xd2, const double yd2, const double xs1, const double ys1, const double xs2, const double ys2, SGP* pSGP);
+  void projectSingleView (const Phantom& phm, DetectorArray& darray, const double xd1, const double yd1, const double xd2, const double yd2, const double xs1, const double ys1, const double xs2, const double ys2, const double dDetAngle);
 
-  double projectSingleLine (const Phantom& phm, const double x1, const double y1, const double x2, const double y2, SGP* pSGP);
+  double projectSingleLine (const Phantom& phm, const double x1, const double y1, const double x2, const double y2);
 
-  double projectLineAgainstPElem (const PhantomElement& pelem, const double x1, const double y1, const double x2, const double y2, SGP* pSGP);
+  double projectLineAgainstPElem (const PhantomElement& pelem, const double x1, const double y1, const double x2, const double y2);
 
-  void traceShowParam (const char *label, const char *fmt, int row, int color, ...);
+  void traceShowParam (const char* szLabel, const char *fmt, int row, int color, ...);
+  void traceShowParamXOR (const char* szLabel, const char *fmt, int row, int color, ...);
+  void traceShowParamRasterOp (int iRasterOp, const char* szLabel, const char* fmt, int row, int color, va_list va);
 
 
 };
 
-const static int RAYSUM_TRACE_ROW_TITLE=1;
-const static int RAYSUM_TRACE_ROW_TITLE2=2;
-const static int RAYSUM_TRACE_ROW_PHANT_ID=4;
-const static int RAYSUM_TRACE_ROW_CHROMATIC=7;
-const static int RAYSUM_TRACE_ROW_SCATTER=8;
-const static int RAYSUM_TRACE_ROW_PHOT_STAT=9;
-const static int RAYSUM_TRACE_ROW_NDET=12;
-const static int RAYSUM_TRACE_ROW_NVIEW=13;
-const static int RAYSUM_TRACE_ROW_SAMPLES=14;
-const static int RAYSUM_TRACE_ROW_CURR_VIEW=17;
-const static int RAYSUM_TRACE_ROW_ATTEN=18;
+const static int PROJECTION_TRACE_ROW_TITLE=0;
+const static int PROJECTION_TRACE_ROW_TITLE2=1;
+const static int PROJECTION_TRACE_ROW_PHANT_ID=3;
+const static int PROJECTION_TRACE_ROW_GEOMETRY=4;
+const static int PROJECTION_TRACE_ROW_FOCAL_LENGTH=5;
+const static int PROJECTION_TRACE_ROW_FIELD_OF_VIEW=6;
+const static int PROJECTION_TRACE_ROW_NDET=7;
+const static int PROJECTION_TRACE_ROW_NVIEW=8;
+const static int PROJECTION_TRACE_ROW_SAMPLES=9;
+const static int PROJECTION_TRACE_ROW_CURR_VIEW=10;
+const static int PROJECTION_TRACE_ROW_ATTEN=11;