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.3 2000/06/19 17:58:20 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 ******************************************************************************/
32 #include "projections.h"
34 // Projections are collected along an array of ndet detectors. The data
35 // for these detectors is stored in the class DetectorArray
37 typedef float DetectorValue;
42 DetectorArray (const int ndet);
43 ~DetectorArray (void);
45 const int nDet(void) const {return m_nDet;}
46 const double viewAngle(void) const {return m_viewAngle;}
47 DetectorValue* detValues(void) {return m_detValues;}
48 const DetectorValue* detValues(void) const {return m_detValues;}
50 void setViewAngle (double viewAngle)
51 { m_viewAngle = viewAngle; }
54 DetectorValue* m_detValues; /* Pointer to array of values recorded by detector */
55 int m_nDet; /* Number of detectors in array */
56 double m_viewAngle; /* View angle in radians */
58 DetectorArray& operator=(const DetectorArray& rhs);
59 DetectorArray& operator()(const DetectorArray& rhs);
72 Scanner (const Phantom& phm, const ScannerGeometry geometry, int nDet, int nView, int nSample, const double rot_anglen);
75 void collectProjections (Projections& proj, const Phantom& phm, const int start_view, const int trace);
77 void setNView (int nView);
78 const unsigned int nDet(void) const {return m_nDet;}
79 const unsigned int nView(void) const {return m_nView;}
80 const double phmLen(void) const {return m_phmLen;}
81 const double rotInc(void) const {return m_rotInc;}
82 const double detInc(void) const {return m_detInc;}
83 const double radius(void) const {return m_radius;}
87 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);
89 double projectSingleLine (const Phantom& phm, const double x1, const double y1, const double x2, const double y2);
91 double projectLineAgainstPElem (const PhantomElement& pelem, const double x1, const double y1, const double x2, const double y2);
93 void traceShowParam (const char *label, const char *fmt, int row, int color, ...);
96 ScannerGeometry m_geometry; /* Geometry of detectory */
97 unsigned int m_nDet; /* Number of detectors in array */
98 unsigned int m_nView; /* Number of rotated views */
99 unsigned int m_nSample; /* Number of rays per detector */
100 double m_detLen; /* Total length of detector array */
101 double m_rotLen; /* Rotation angle length in radians (norm 2PI) */
102 double m_detInc; /* Increment between centers of detectors */
103 double m_rotInc; /* Increment in rotation angle between views */
104 double m_radius; /* Radius of rotation. Distance from */
105 /* center of phm to center of det */
106 double m_phmLen; /* Maximum Length of phantom or area of interest */
109 double xd1,yd1,xd2,yd2; /* Coordinates of detector endpoints */
110 double xs1,ys1,xs2,ys2; /* Coordinates of source endpoints */
111 double angle; /* Starting angle */
114 static const int N_EXTRA_DETECTORS=4; /* Number of extra detectors widths when calculating detlen */
120 const static int RAYSUM_TRACE_ROW_TITLE=1;
121 const static int RAYSUM_TRACE_ROW_TITLE2=2;
122 const static int RAYSUM_TRACE_ROW_PHANT_ID=4;
123 const static int RAYSUM_TRACE_ROW_CHROMATIC=7;
124 const static int RAYSUM_TRACE_ROW_SCATTER=8;
125 const static int RAYSUM_TRACE_ROW_PHOT_STAT=9;
126 const static int RAYSUM_TRACE_ROW_NDET=12;
127 const static int RAYSUM_TRACE_ROW_NVIEW=13;
128 const static int RAYSUM_TRACE_ROW_SAMPLES=14;
129 const static int RAYSUM_TRACE_ROW_CURR_VIEW=17;
130 const static int RAYSUM_TRACE_ROW_ATTEN=18;