** This is part of the CTSim program
** Copyright (C) 1983-2000 Kevin Rosenberg
**
-** $Id: scanner.h,v 1.1 2000/06/17 20:12:14 kevin Exp $
+** $Id: scanner.h,v 1.16 2000/12/16 03:29:02 kevin Exp $
**
**
** This program is free software; you can redistribute it and/or modify
#ifndef SCANNER_H
#define SCANNER_H
-#include "projections.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
{
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; }
private:
- DetectorValue* m_detValues; /* Pointer to array of values recorded by detector */
- int m_nDet; /* Number of detectors in array */
- double m_viewAngle; /* View angle in radians */
+ DetectorValue* m_detValues; // Pointer to array of values recorded by detector
+ int m_nDet; // Number of detectors in array */
+ double m_viewAngle; // View angle in radians
+
+ DetectorArray& operator=(const DetectorArray& rhs); // assignment
+ DetectorArray (const DetectorArray& rhs); // copy constructor
};
-typedef enum {
- DETECTOR_PARALLEL,
- DETECTOR_EQUIANGLE,
- DETECTOR_EQUILINEAR
-} ScannerGeometry;
-
class Scanner
{
public:
- Scanner (const Phantom& phm, const ScannerGeometry geometry, int nDet, int nView, int nSample, const double rot_anglen);
+ 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, 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::TRACE_NONE, SGP* pSGP = NULL);
+
+ void collectProjections (Projections& proj, const Phantom& phm, const int iStartView, const int iNumViews, bool bStoreAtViewPosition, const int trace = Trace::TRACE_NONE, SGP* pSGP = NULL);
+
+ void setNView (int nView);
+
+ bool fail() const {return m_fail;}
+ const std::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);
+ static const char* convertGeometryIDToName (const int idGeometry);
+ static const char* convertGeometryIDToTitle (const int idGeometry);
+
+ private:
+ bool m_fail;
+ std::string m_failMessage;
+ 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_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;
- 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;}
+#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;
- private:
- 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, const double dDetAngle);
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);
- 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);
- ScannerGeometry m_geometry; /* Geometry of detectory */
- 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 */
-
- struct {
- double xd1,yd1,xd2,yd2; /* Coordinates of detector endpoints */
- double xs1,ys1,xs2,ys2; /* Coordinates of source endpoints */
- double angle; /* Starting angle */
- } m_initPos;
+};
- static const int N_EXTRA_DETECTORS=4; /* Number of extra detectors widths when calculating detlen */
+const static int PROJECTION_TRACE_ROW_PHANT_ID=0;
+const static int PROJECTION_TRACE_ROW_GEOMETRY=1;
+const static int PROJECTION_TRACE_ROW_FOCAL_LENGTH=2;
+const static int PROJECTION_TRACE_ROW_FIELD_OF_VIEW=3;
+const static int PROJECTION_TRACE_ROW_NDET=4;
+const static int PROJECTION_TRACE_ROW_NVIEW=5;
+const static int PROJECTION_TRACE_ROW_SAMPLES=6;
+const static int PROJECTION_TRACE_ROW_CURR_VIEW=7;
+const static int PROJECTION_TRACE_ROW_ATTEN=8;
- int m_trace;
-};
#endif
projects / ctsim.git / blobdiff