** Date Started: July 1, 1984
**
** This is part of the CTSim program
-** Copyright (C) 1983-2000 Kevin Rosenberg
+** Copyright (c) 1983-2001 Kevin Rosenberg
**
-** $Id: phantom.h,v 1.13 2000/08/25 15:59:13 kevin Exp $
+** $Id: phantom.h,v 1.22 2001/02/08 06:25:07 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
#ifndef PHANTOM_H
#define PHANTOM_H
-#include <slist>
+#include <list>
#include "ctsupport.h"
-using namespace std;
-
typedef enum {
PELEM_INVALID,
PELEM_RECTANGLE,
double* rectLimits() {return m_rectLimits;}
double* xOutline() {return m_xOutline;}
double* yOutline() {return m_yOutline;}
+ double* const xOutline() const {return m_xOutline;}
+ double* const yOutline() const {return m_yOutline;}
const double atten() const {return m_atten;}
- const double diameter() const {return m_diameter;}
const double xmin() const {return m_xmin;}
const double xmax() const {return m_xmax;}
const double ymin() const {return m_ymin;}
const double u() const {return m_u;}
const double v() const {return m_v;}
+ void printDefinition (std::ostream& os) const;
+ void printDefinition (std::ostringstream& os) const;
+
private:
PhmElemType m_type; // pelem type (box, ellipse, etc)
double m_cx, m_cy; // center of pelem
double *m_x, *m_y; // ptr to array of points in obj world coord
int m_nPoints; // number of points in outline arrays
double m_xmin, m_xmax, m_ymin, m_ymax; // pelem limits
- double m_diameter;
GRFMTX_2D m_xformPhmToObj; // map from phantom to normalized pelem coords
GRFMTX_2D m_xformObjToPhm; // map from normalized pelem coords to phantom coords
double* m_xOutline;
double* m_yOutline;
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 const int POINTS_PER_CIRCLE;
+ static const double SCALE_PELEM_EXTENT; // increase pelem limits by 0.5%
- static PhmElemType PhantomElement::convertNameToType (const char* const typeName);
+ static PhmElemType convertNameToType (const char* const typeName);
+ static const char* const convertTypeToName (PhmElemType iType);
void makeTransformMatrices ();
bool createFromFile (const char* const fname);
+ bool fileWrite (const char* const fname);
+
void addPElem (const PhantomElement& pelem);
void addPElem (const char* const composition, const double cx, const double cy, const double u, const double v, const double rot, const double atten);
- void convertToImagefile (ImageFile& im, const int in_nsample, const int trace, const int colStart, const int colCount) const;
+ void convertToImagefile (ImageFile& im, double dViewRatio, const int in_nsample, const int trace, const int colStart, const int colCount, bool bStoreAtColumnPos) const;
- void convertToImagefile (ImageFile& im, const int in_nsample, const int trace) const;
+ void convertToImagefile (ImageFile& im, double dViewRatio, const int in_nsample, const int trace) const;
+
+ void printDefinitions (std::ostream& os) const;
+ void printDefinitions (std::ostringstream& os) const;
bool fail() const {return m_fail;}
- const string& failMessage() const {return m_failMessage;}
- const string& name() const {return m_name;}
+ const std::string& failMessage() const {return m_failMessage;}
+ const std::string& name() const {return m_name;}
const int id() const {return m_id;}
#if HAVE_SGP
void addStdSheppLogan ();
void addStdSheppLoganBordered ();
- void print () const;
+ void print (std::ostream& os) const;
+ void print (std::ostringstream& os) const;
+
+ double maxAxisLength () const
+ { return maxValue<double> (m_xmax - m_xmin, m_ymax - m_ymin); }
- const double maxAxisLength () const {return (((m_xmax - m_xmin) > (m_ymax - m_ymin)) ? (m_xmax - m_xmin) : (m_ymax - m_ymin));}
+ double getDiameterBoundaryCircle() const
+ { return SQRT2 * maxAxisLength(); }
- const double diameter() const {return m_diameter;}
const double xmin() const {return m_xmin;}
const double xmax() const {return m_xmax;}
const double ymin() const {return m_ymin;}
const double ymax() const {return m_ymax;}
- slist<PhantomElement*>& listPElem() {return m_listPElem;}
- const slist<PhantomElement*>& listPElem() const {return m_listPElem;}
+ std::list<PhantomElement*>& listPElem() {return m_listPElem;}
+ const std::list<PhantomElement*>& listPElem() const {return m_listPElem;}
const int nPElem() const {return m_nPElem;}
static const int getPhantomCount() {return s_iPhantomCount;}
PhantomComposition m_composition;
int m_nPElem; // number of pelems in phantom
double m_xmin, m_xmax, m_ymin, m_ymax; // extent of pelems in pelem coordinates
- double m_diameter; // diameter of object
- mutable slist<PhantomElement*> m_listPElem; // pelem lists
- string m_name;
+ mutable std::list<PhantomElement*> m_listPElem; // pelem lists
+ std::string m_name;
int m_id;
bool m_fail;
- string m_failMessage;
-
+ std::string m_failMessage;
static const char* s_aszPhantomName[];
static const char* s_aszPhantomTitle[];
static const int s_iPhantomCount;
Phantom& operator= (const Phantom&); // assignment operator
};
-typedef slist<PhantomElement*>::iterator PElemIterator;
-typedef slist<PhantomElement*>::const_iterator PElemConstIterator;
+typedef std::list<PhantomElement*>::iterator PElemIterator;
+typedef std::list<PhantomElement*>::const_iterator PElemConstIterator;
#endif