X-Git-Url: http://git.kpe.io/?p=ctsim.git;a=blobdiff_plain;f=include%2Fscanner.h;h=d302b67b657ce1dcbf660d2df85075c2e3dc2924;hp=162e89686744e2a34a9f04374e3333ea0fe7f007;hb=f13a8c004b8f182b42d9e4df2bcd7c7f030bf1ad;hpb=c5e7140bd08b8c8f527713e8dc861bcb7ee5f633 diff --git a/include/scanner.h b/include/scanner.h index 162e896..d302b67 100644 --- a/include/scanner.h +++ b/include/scanner.h @@ -7,10 +7,7 @@ ** Date Started: July 1, 1984 ** ** This is part of the CTSim program -** Copyright (C) 1983-2000 Kevin Rosenberg -** -** $Id: scanner.h,v 1.17 2001/02/08 06:25:07 kevin Exp $ -** +** Copyright (C) 1983-2009 Kevin Rosenberg ** ** This program is free software; you can redistribute it and/or modify ** it under the terms of the GNU General Public License (version 2) as @@ -56,9 +53,9 @@ class DetectorArray { 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 @@ -68,32 +65,42 @@ class DetectorArray class Scanner { public: - 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 dFocalLengthRatio, double dViewRatio, double dScanRatio); + static const int GEOMETRY_INVALID; + static const int GEOMETRY_PARALLEL; + static const int GEOMETRY_EQUILINEAR; + static const int GEOMETRY_EQUIANGULAR; + static const int GEOMETRY_LINOGRAM; + + + Scanner (const Phantom& phm, const char* const geometryName, int nDet, + int nView, int iOffsetView, int nSample, const double rot_anglen, + double dFocalLengthRatio, double dCenterDetectorRatio, double dViewRatio, double dScanRatio); ~Scanner(); - - 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 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, const int iOffsetView, bool bStoreAtViewPosition, const int trace = Trace::TRACE_NONE, SGP* pSGP = NULL); + + void collectProjections (Projections& proj, const Phantom& phm, const int iStartView, const int iNumViews, const int iOffsetView, int iStorageOffset, const int trace = Trace::TRACE_NONE, SGP* pSGP = NULL); void setNView (int nView); + void setOffsetView (int iOffsetView); 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;} - + unsigned int offsetView() const {return m_iOffsetView;} + unsigned int startView() const {return m_startView;} double rotInc() const {return m_rotInc;} double detInc() const {return m_detInc;} double detLen() const {return m_detLen;} + double detStart() const {return m_detStart;} double focalLength() const {return m_dFocalLength;} + double sourceDetectorLength() const {return m_dSourceDetectorLength;} + double centerDetectorLength() const {return m_dCenterDetectorLength;} + double viewDiameter() const {return m_dViewDiameter;} double scanDiameter() const {return m_dScanDiameter;} double fanBeamAngle() const {return m_dFanBeamAngle;} @@ -101,46 +108,52 @@ class Scanner 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 const char* const* getGeometryNameArray() {return s_aszGeometryName;} + static const char* const* 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_dFocalLength; // Focal Length + unsigned int m_nDet; /* Number of detectors in array */ + unsigned int m_nView; /* Number of rotated views */ + unsigned int m_iOffsetView; + unsigned int m_startView; + unsigned int m_nSample; /* Number of rays per detector */ + double m_dFocalLength; // Focal Length, distance from source to center + double m_dSourceDetectorLength; // Distance from source to detectors + double m_dCenterDetectorLength; // Distance from center to detectors double m_dViewDiameter; // Diameter of area being processed double m_dScanDiameter; // Diamer of area being scanned double m_dViewRatio; // View Ratio to diameter phantom - double m_dFocalLengthRatio; // Focal Length ratio to radius phantom + double m_dFocalLengthRatio; // Source to Center Length as ratio to viewDiameter radius + double m_dCenterDetectorRatio; // Center to Detector Length as ratio of viewDiameter radius double m_dScanRatio; // Scan length to view length ratio double m_dFanBeamAngle; - 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_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_detStart; double m_dXCenter; // Center of Phantom - double m_dYCenter; + 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 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; -#if HAVE_SGP +#ifdef HAVE_SGP SGP* m_pSGP; // Pointer to graphics device double m_dXMinWin; // Extent of graphics window double m_dXMaxWin; @@ -149,8 +162,8 @@ class Scanner double m_dTextHeight; #endif - static const char* s_aszGeometryName[]; - static const char* s_aszGeometryTitle[]; + static const char* const s_aszGeometryName[]; + static const char* const s_aszGeometryTitle[]; static const int s_iGeometryCount; 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);