X-Git-Url: http://git.kpe.io/?a=blobdiff_plain;f=include%2Fphantom.h;h=a82381f78fdfc014e8ddebe70ed604a173ac33fd;hb=27a474e0622ebb7229fd5705552021f63d8f932d;hp=32eedf6e9f933ee8e617e19eeef89b50266c5ba9;hpb=595e63c804284d460ce4d032c3848b75bc57186d;p=ctsim.git

diff --git a/include/phantom.h b/include/phantom.h
index 32eedf6..a82381f 100644
--- a/include/phantom.h
+++ b/include/phantom.h
@@ -9,7 +9,7 @@
 **  This is part of the CTSim program
 **  Copyright (C) 1983-2000 Kevin Rosenberg
 **
-**  $Id: phantom.h,v 1.2 2000/06/19 20:08:09 kevin Exp $
+**  $Id: phantom.h,v 1.9 2000/07/20 11:17:31 kevin Exp $
 **
 **  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,7 +56,7 @@ class PhantomElement
 
     ~PhantomElement (void);
 
-    bool isPointInside (double x, double y, const CoordType coord_type);
+    bool isPointInside (double x, double y, const CoordType coord_type) const;
 
     bool clipLineNormalizedCoords (double& x1, double& y1, double& x2, double& y2) const;
  
@@ -80,9 +80,9 @@ class PhantomElement
 
  private:
     PhmElemType m_type;	     // pelem type (box, ellipse, etc)
-    double m_atten;	     // X-ray attenuation coefficient
     double m_cx, m_cy;	     // center of pelem 
     double m_u, m_v;		     // size of pelem 
+    double m_atten;	     // X-ray attenuation coefficient
     double m_rot;		     // pelem rotation angle (in radians) 
     double *m_x, *m_y;	     // ptr to array of points in obj world coord 
     int m_nPoints;		     // number of points in outline arrays 
@@ -95,6 +95,9 @@ class PhantomElement
     double  m_rectLimits[4];
 
     static const int POINTS_PER_CIRCLE = 360;
+    static const double SCALE_PELEM_EXTENT=0.005;  // increase pelem limits by 0.5% 
+
+    static PhmElemType PhantomElement::convertNameToType (const char* const typeName);
 
     void makeTransformMatrices (void);
 
@@ -104,8 +107,10 @@ class PhantomElement
 
     void calcEllipsePoints (double x[], double y[], const int pts, const double u, const double v);
 
-    int numCirclePoints (double theta) const;
+    static int numCirclePoints (double theta);
 
+    PhantomElement (const PhantomElement& rhs);        // copy constructor
+    PhantomElement& operator= (const PhantomElement&); // assignment operator
 };
 
 
@@ -120,7 +125,30 @@ typedef enum {
 class Phantom
 {
  public:
+    typedef enum {
+      PHM_INVALID,
+      PHM_HERMAN,               /* Herman head phantom */
+      PHM_BHERMAN,              /* Bordered Herman head phantom */
+      PHM_ROWLAND,              /* Rowland head phantom */
+      PHM_BROWLAND,             /* Bordered Rowland head phantom */
+      PHM_UNITPULSE             /* Unit pulse phantom */
+    } PhantomID;
+
+    static const char PHM_HERMAN_STR[];
+    static const char PHM_BHERMAN_STR[];
+    static const char PHM_ROWLAND_STR[];
+    static const char PHM_BROWLAND_STR[];
+    static const char PHM_UNITPULSE_STR[];
+
+    static const char PHM_HERMAN_TITLE_STR[];
+    static const char PHM_BHERMAN_TITLE_STR[];
+    static const char PHM_ROWLAND_TITLE_STR[];
+    static const char PHM_BROWLAND_TITLE_STR[];
+    static const char PHM_UNITPULSE_TITLE_STR[];
+
     Phantom (void);
+    Phantom (const char* const phmName);
+
     ~Phantom (void);
 
     void setComposition (PhantomComposition composition)
@@ -129,7 +157,9 @@ class Phantom
     const PhantomComposition getComposition (void) const
 	{ return m_composition; }
 
-    void create (const int phmid);
+    bool createFromPhantom (const char* const phmName);
+
+    bool createFromPhantom (const PhantomID phmid);
 
     bool createFromFile (const char* const fname);
 
@@ -141,12 +171,18 @@ class Phantom
 
     void convertToImagefile (ImageFile& im, const int in_nsample, const int trace) const;
 
+    bool fail(void) const             {return m_fail;}
+    const string& failMessage(void) const {return m_failMessage;}
+    const string& name(void) const     {return m_name;}
+    const PhantomID id(void) const     {return m_id;}
+
 #if HAVE_SGP
     void show (void) const;
     void draw (void) const;
 #endif
     
     void addStdHerman (void);
+    void addStdHermanBordered (void);
     void addStdRowland (void);
     void addStdRowlandBordered (void);
 
@@ -170,6 +206,18 @@ class Phantom
     double m_diameter;                        // diameter of object
     mutable slist<PhantomElement*> m_listPElem;      // pelem lists
     string m_name;
+    PhantomID m_id;
+    bool m_fail;
+    string m_failMessage;
+
+    // Standard phantomsa
+    static PhantomID convertNameToPhantomID (const char* const phmName);
+    static const char* convertPhantomIDToName (const PhantomID phmID);
+
+    void init(void);
+
+    Phantom (const Phantom& rhs);        // copy constructor
+    Phantom& operator= (const Phantom&); // assignment operator
 };
 
 typedef slist<PhantomElement*>::iterator PElemIterator;