X-Git-Url: http://git.kpe.io/?p=ctsim.git;a=blobdiff_plain;f=include%2Fscanner.h;h=79eee68d0969287709637f61dc2f931e4d1d5d17;hp=1a01c89312cc2670dc16b3c3195680d08073f314;hb=1e88cf0f7fa4f690ea9f110e8ed3f2b5338d0a10;hpb=16d093f9366acfccee0f701f736840bcd47c0709

diff --git a/include/scanner.h b/include/scanner.h
index 1a01c89..79eee68 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.6 2000/07/18 14:51:06 kevin Exp $
+**  $Id: scanner.h,v 1.12 2000/08/25 15:59:13 kevin Exp $
 **
 **
 **  This program is free software; you can redistribute it and/or modify
@@ -29,9 +29,12 @@
 #ifndef SCANNER_H
 #define SCANNER_H
 
+#include "trace.h"
+
 class Projections;
 class Phantom;
 class PhantomElement;
+class SGP;
 
 // Projections are collected along an array of ndet detectors.  The data
 // for these detectors is stored in the class DetectorArray
@@ -61,24 +64,22 @@ class DetectorArray
   DetectorArray (const DetectorArray& rhs);             // copy constructor
 };
 
+
 class Scanner
 {
  public:
-    typedef enum {
-	GEOMETRY_INVALID,
-	GEOMETRY_PARALLEL,
-	GEOMETRY_EQUILINEAR,
-	GEOMETRY_EQUIANGLE
-    } GeometryID;
-  
-    static const char GEOMETRY_PARALLEL_STR[] = "parallel";
-    static const char GEOMETRY_EQUILINEAR_STR[] = "equilinear";
-    static const char GEOMETRY_EQUIANGLE_STR[] = "equiangle";
+  static const int Scanner::GEOMETRY_INVALID;
+  static const int Scanner::GEOMETRY_PARALLEL;
+  static const int Scanner::GEOMETRY_EQUILINEAR;
+  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, const int trace);
+  void collectProjections (Projections& proj, const Phantom& phm, const int trace = TRACE_NONE, SGP* pSGP = NULL);
+
+  void collectProjections (Projections& proj, const Phantom& phm, const int iStartView, const int iNumViews, bool bStoreAtViewPosition, const int trace = TRACE_NONE, SGP* pSGP = NULL);
 
   void setNView (int nView);
 
@@ -89,22 +90,29 @@ class Scanner
   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;}
-
-
+  const double detLen(void) const {return m_detLen;}
+
+  static const 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);
+  static const char* convertGeometryIDToName (const int idGeometry);
+  static const char* convertGeometryIDToTitle (const int idGeometry);
+  
  private:
   bool m_fail;
   string m_failMessage;
-  GeometryID m_idGeometry;
+  int m_idGeometry;
   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_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
   int m_trace;
   struct {
     double xd1,yd1,xd2,yd2;	/* Coordinates of detector endpoints */
@@ -112,17 +120,17 @@ class Scanner
     double angle;		/* Starting angle */
   } m_initPos;
 
-  static const int N_EXTRA_DETECTORS=4;           /* Number of extra detectors widths when calculating detlen */
+  static const char* s_aszGeometryName[];
+  static const char* s_aszGeometryTitle[];
+  static const int s_iGeometryCount;
 
-  static GeometryID convertGeometryNameToID (const char* const geometryName);
-  
-  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);
+  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);
 
-  double projectSingleLine (const Phantom& phm, const double x1, const double y1, const double x2, const double y2);
+  double projectSingleLine (const Phantom& phm, 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);
+  double projectLineAgainstPElem (const PhantomElement& pelem, const double x1, const double y1, const double x2, const double y2, SGP* pSGP);
 
-  void traceShowParam (const char *label, const char *fmt, int row, int color, ...);
+  void traceShowParam (SGP* pSGP, const char *label, const char *fmt, int row, int color, ...);
 
 
 };