include/scanner.h

   1 /*****************************************************************************
   2 ** FILE IDENTIFICATION
   3 **
   4 **   Name:          scanner.h
   5 **   Purpose:       Scanner header file
   6 **   Programmer:    Kevin Rosenberg
   7 **   Date Started:  July 1, 1984
   8 **
   9 **  This is part of the CTSim program
  10 **  Copyright (C) 1983-2000 Kevin Rosenberg
  11 **
  12 **  $Id: scanner.h,v 1.14 2000/08/31 08:38:58 kevin Exp $
  13 **
  14 **
  15 **  This program is free software; you can redistribute it and/or modify
  16 **  it under the terms of the GNU General Public License (version 2) as
  17 **  published by the Free Software Foundation.
  18 **
  19 **  This program is distributed in the hope that it will be useful,
  20 **  but WITHOUT ANY WARRANTY; without even the implied warranty of
  21 **  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  22 **  GNU General Public License for more details.
  23 **
  24 **  You should have received a copy of the GNU General Public License
  25 **  along with this program; if not, write to the Free Software
  26 **  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
  27 ******************************************************************************/
  28
  29 #ifndef SCANNER_H
  30 #define SCANNER_H
  31
  32 #include "trace.h"
  33
  34 class Projections;
  35 class Phantom;
  36 class PhantomElement;
  37 class SGP;
  38
  39 // Projections are collected along an array of ndet detectors.  The data
  40 // for these detectors is stored in the class DetectorArray
  41
  42 typedef float DetectorValue;
  43
  44 class DetectorArray
  45 {
  46  public:
  47   DetectorArray (const int ndet);
  48   ~DetectorArray ();
  49
  50   const int nDet() const {return m_nDet;}
  51   const double viewAngle() const {return m_viewAngle;}
  52   DetectorValue* detValues() {return m_detValues;}
  53   const DetectorValue* detValues() const {return m_detValues;}
  54
  55   void setViewAngle (double viewAngle)
  56       { m_viewAngle = viewAngle; }
  57
  58  private:
  59   DetectorValue* m_detValues;   // Pointer to array of values recorded by detector
  60   int m_nDet;                   // Number of detectors in array */
  61   double m_viewAngle;           // View angle in radians
  62
  63   DetectorArray& operator=(const DetectorArray& rhs);   // assignment
  64   DetectorArray (const DetectorArray& rhs);             // copy constructor
  65 };
  66
  67
  68 class Scanner
  69 {
  70  public:
  71   static const int Scanner::GEOMETRY_INVALID;
  72   static const int Scanner::GEOMETRY_PARALLEL;
  73   static const int Scanner::GEOMETRY_EQUILINEAR;
  74   static const int Scanner::GEOMETRY_EQUIANGULAR;
  75
  76
  77   Scanner (const Phantom& phm, const char* const geometryName, int nDet, int nView, int nSample, const double rot_anglen, double dFieldOfView, double dFocalLength);
  78   ~Scanner();
  79
  80   void collectProjections (Projections& proj, const Phantom& phm, const int trace = Trace::TRACE_NONE, SGP* pSGP = NULL);
  81
  82   void collectProjections (Projections& proj, const Phantom& phm, const int iStartView, const int iNumViews, bool bStoreAtViewPosition, const int trace = Trace::TRACE_NONE, SGP* pSGP = NULL);
  83
  84   void setNView (int nView);
  85
  86   bool fail() const {return m_fail;}
  87   const string& failMessage() const {return m_failMessage;}
  88   unsigned int nDet() const {return m_nDet;}
  89   unsigned int nView() const {return m_nView;}
  90   double phmLen() const {return m_phmLen;}
  91   double rotInc() const {return m_rotInc;}
  92   double detInc() const {return m_detInc;}
  93   double detLen() const {return m_detLen;}
  94   double focalLength() const {return m_dFocalLength;}
  95   double fieldOfView() const {return m_dFieldOfView;}
  96   int geometry() const {return m_idGeometry;}
  97
  98   static int getGeometryCount() {return s_iGeometryCount;}
  99   static const char** getGeometryNameArray() {return s_aszGeometryName;}
 100   static const char** getGeometryTitleArray() {return s_aszGeometryTitle;}
 101   static int convertGeometryNameToID (const char* const geometryName);
 102   static const char* convertGeometryIDToName (const int idGeometry);
 103   static const char* convertGeometryIDToTitle (const int idGeometry);
 104
 105  private:
 106   bool m_fail;
 107   string m_failMessage;
 108   int m_idGeometry;
 109   unsigned int m_nDet;          /* Number of detectors in array */
 110   unsigned int m_nView;         /* Number of rotated views */
 111   unsigned int m_nSample;       /* Number of rays per detector */
 112   double m_dFieldOfView;        // Field of View
 113   double m_dFocalLength;        // Focal Length
 114   double m_dFieldOfViewRatio;   // Field of View Ratio to diameter phmLen
 115   double m_dFocalLengthRatio;   // Focal Length ratio to radius phmLen
 116   double m_detLen;              // Total length of detector array
 117   double m_rotLen;              // Rotation angle length in radians (norm 2PI)
 118   double m_detInc;              // Increment between centers of detectors
 119   double m_rotInc;              // Increment in rotation angle between views
 120   double m_phmLen;              // Maximum Length of phantom or area of interest
 121   double m_dXCenter;            // Center of Phantom
 122   double m_dYCenter;
 123   double m_dAngularDetIncrement;
 124   double m_dAngularDetLen;
 125
 126   int m_trace;
 127   struct {
 128     double xd1,yd1,xd2,yd2;     /* Coordinates of detector endpoints */
 129     double xs1,ys1,xs2,ys2;     /* Coordinates of source endpoints */
 130     double angle;               /* Starting angle */
 131     double dAngularDet;
 132   } m_initPos;
 133
 134   GRFMTX_2D m_rotmtxIncrement;
 135
 136 #if HAVE_SGP
 137   SGP* m_pSGP;                  // Pointer to graphics device
 138   double m_dXMinWin;            // Extent of graphics window
 139   double m_dXMaxWin;
 140   double m_dYMinWin;
 141   double m_dYMaxWin;
 142   double m_dTextHeight;
 143 #endif
 144
 145   static const char* s_aszGeometryName[];
 146   static const char* s_aszGeometryTitle[];
 147   static const int s_iGeometryCount;
 148
 149   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);
 150
 151   double projectSingleLine (const Phantom& phm, const double x1, const double y1, const double x2, const double y2);
 152
 153   double projectLineAgainstPElem (const PhantomElement& pelem, const double x1, const double y1, const double x2, const double y2);
 154
 155   void traceShowParam (const char* szLabel, const char *fmt, int row, int color, ...);
 156   void traceShowParamXOR (const char* szLabel, const char *fmt, int row, int color, ...);
 157   void traceShowParamRasterOp (int iRasterOp, const char* szLabel, const char* fmt, int row, int color, va_list va);
 158
 159
 160 };
 161
 162 const static int PROJECTION_TRACE_ROW_TITLE=0;
 163 const static int PROJECTION_TRACE_ROW_TITLE2=1;
 164 const static int PROJECTION_TRACE_ROW_PHANT_ID=3;
 165 const static int PROJECTION_TRACE_ROW_GEOMETRY=4;
 166 const static int PROJECTION_TRACE_ROW_FOCAL_LENGTH=5;
 167 const static int PROJECTION_TRACE_ROW_FIELD_OF_VIEW=6;
 168 const static int PROJECTION_TRACE_ROW_NDET=7;
 169 const static int PROJECTION_TRACE_ROW_NVIEW=8;
 170 const static int PROJECTION_TRACE_ROW_SAMPLES=9;
 171 const static int PROJECTION_TRACE_ROW_CURR_VIEW=10;
 172 const static int PROJECTION_TRACE_ROW_ATTEN=11;
 173
 174
 175
 176 #endif