1 /*****************************************************************************
5 ** Purpose: Scanner header file
6 ** Programmer: Kevin Rosenberg
7 ** Date Started: July 1, 1984
9 ** This is part of the CTSim program
10 ** Copyright (C) 1983-2000 Kevin Rosenberg
12 ** $Id: scanner.h,v 1.16 2000/12/16 03:29:02 kevin Exp $
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.
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.
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 ******************************************************************************/
39 // Projections are collected along an array of ndet detectors. The data
40 // for these detectors is stored in the class DetectorArray
42 typedef float DetectorValue;
47 DetectorArray (const int ndet);
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;}
55 void setViewAngle (double viewAngle)
56 { m_viewAngle = viewAngle; }
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
63 DetectorArray& operator=(const DetectorArray& rhs); // assignment
64 DetectorArray (const DetectorArray& rhs); // copy constructor
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;
77 Scanner (const Phantom& phm, const char* const geometryName, int nDet, int nView, int nSample, const double rot_anglen, double dFieldOfView, double dFocalLength);
80 void collectProjections (Projections& proj, const Phantom& phm, const int trace = Trace::TRACE_NONE, SGP* pSGP = NULL);
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);
84 void setNView (int nView);
86 bool fail() const {return m_fail;}
87 const std::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;}
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);
107 std::string m_failMessage;
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
123 double m_dAngularDetIncrement;
124 double m_dAngularDetLen;
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 */
134 GRFMTX_2D m_rotmtxIncrement;
137 SGP* m_pSGP; // Pointer to graphics device
138 double m_dXMinWin; // Extent of graphics window
142 double m_dTextHeight;
145 static const char* s_aszGeometryName[];
146 static const char* s_aszGeometryTitle[];
147 static const int s_iGeometryCount;
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);
151 double projectSingleLine (const Phantom& phm, const double x1, const double y1, const double x2, const double y2);
153 double projectLineAgainstPElem (const PhantomElement& pelem, const double x1, const double y1, const double x2, const double y2);
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);
162 const static int PROJECTION_TRACE_ROW_PHANT_ID=0;
163 const static int PROJECTION_TRACE_ROW_GEOMETRY=1;
164 const static int PROJECTION_TRACE_ROW_FOCAL_LENGTH=2;
165 const static int PROJECTION_TRACE_ROW_FIELD_OF_VIEW=3;
166 const static int PROJECTION_TRACE_ROW_NDET=4;
167 const static int PROJECTION_TRACE_ROW_NVIEW=5;
168 const static int PROJECTION_TRACE_ROW_SAMPLES=6;
169 const static int PROJECTION_TRACE_ROW_CURR_VIEW=7;
170 const static int PROJECTION_TRACE_ROW_ATTEN=8;