r186: *** empty log message ***
[ctsim.git] / include / projections.h
index 8d5acacc9c3d0adc0a46a4fd5b1cca175fdcbb73..8a916166170327ee92adcde5dffe93d6b7993c21 100644 (file)
@@ -9,7 +9,7 @@
 **  This is part of the CTSim program
 **  Copyright (C) 1983-2000 Kevin Rosenberg
 **
-**  $Id: projections.h,v 1.6 2000/06/27 10:48:11 kevin Exp $
+**  $Id: projections.h,v 1.13 2000/08/31 08:38:58 kevin Exp $
 **
 **
 **  This program is free software; you can redistribute it and/or modify
@@ -32,7 +32,7 @@
 class Scanner;
 class DetectorArray;
 class Array2dFileLabel;
-
+class ostringstream;
 
 // Projections
 class Projections
@@ -40,20 +40,22 @@ class Projections
  public:
   Projections (const int nView, const int nDet);
   Projections (const Scanner& scanner);
-  Projections (void);
-  ~Projections (void);
+  Projections ();
+  ~Projections ();
 
   void initFromScanner (const Scanner& scanner);
 
-  void printProjectionData (void);
-  void printScanInfo (void) const;
+  void printProjectionData ();
+  void printScanInfo (ostringstream& os) const;
 
+  bool read (const string& fname);
   bool read (const char* fname);
   bool write (const char* fname);
+  bool write (const string& fname);
   bool detarrayRead (fnetorderstream& fs, DetectorArray& darray, const int view_num);
   bool detarrayWrite (fnetorderstream& fs, const DetectorArray& darray, const int view_num);
 
-  bool reconstruct (ImageFile& im, const char* const filterName, double filt_param, const char* const interpName, int interp_param, const char* const backprojName, const int trace);
+  bool reconstruct (ImageFile& im, const char* const filterName, double filt_param, const char* const filterMethodName, const int zeropad, const char* frequencyFilterName, const char* const interpName, int interp_param, const char* const backprojName, const int trace) const;
 
   void setNView (int nView);  // used in MPI to restrict # of views
   void setRotInc (double rotInc) { m_rotInc = rotInc;}
@@ -61,24 +63,31 @@ class Projections
   void setPhmLen (double phmLen) { m_phmLen = phmLen;}
   void setCalcTime (double calcTime) {m_calcTime = calcTime;}
   void setRemark (const char* remark) {m_remark = remark;}
-  void setRemark (const string remark) {m_remark = remark;}
-
-  const double detStart(void) const {return m_detStart;}
-  const double rotStart(void) const {return m_rotStart;}
-  const double calcTime(void) const {return m_calcTime;}
-  const double phmLen(void) const {return m_phmLen;}
-  const char*  remark(void) const {return m_remark.c_str();}
-  const double detInc(void) const {return m_detInc;}
-  const double rotInc(void) const {return m_rotInc;}
-  const int nDet(void) const {return m_nDet;}
-  const int nView(void) const {return m_nView;}
-  const string& getFilename(void) const {return m_filename;}
-  Array2dFileLabel& getLabel(void) {return m_label;}
-  const Array2dFileLabel& getLabel(void) const {return m_label;}
+  void setRemark (const string& remark) {m_remark = remark;}
+
+  double detStart() const {return m_detStart;}
+  double rotStart() const {return m_rotStart;}
+  double calcTime() const {return m_calcTime;}
+  double phmLen() const {return m_phmLen;}
+  const char*  remark() const {return m_remark.c_str();}
+  double detInc() const {return m_detInc;}
+  double rotInc() const {return m_rotInc;}
+  int nDet() const {return m_nDet;}
+  int nView() const {return m_nView;}
+  int geometry() const {return m_geometry;}
+  double focalLength() const {return m_focalLength;}
+  double fieldOfView() const {return m_fieldOfView;}
+
+  const string& getFilename() const {return m_filename;}
+  Array2dFileLabel& getLabel() {return m_label;}
+  const Array2dFileLabel& getLabel() const {return m_label;}
 
   DetectorArray& getDetectorArray (const int iview)
       { return (*m_projData[iview]); }
   
+  const DetectorArray& getDetectorArray (const int iview) const
+      { return (*m_projData[iview]); }
+  
  private:
   int m_headerSize;            // Size of disk file header 
   int m_geometry;              // Geometry of scanner 
@@ -92,6 +101,8 @@ class Projections
   double m_detStart;           // distance of beginning detector to center phantom
   double m_detInc;             // increment between detectors 
   double m_phmLen;             // Length of phantom edge (phm is square) 
+  double m_focalLength;
+  double m_fieldOfView;
   kuint32 m_year;                   // Creation date & time
   kuint32 m_month;
   kuint32 m_day;
@@ -103,12 +114,12 @@ class Projections
 
   const static kuint16 m_signature = ('P'*256 + 'J');
 
-  bool headerRead (void);
-  bool headerWrite (void);
+  bool headerRead ();
+  bool headerWrite ();
   bool headerRead (fnetorderstream& fs);
   bool headerWrite (fnetorderstream& fs);
-  void newProjData (void);
-  void deleteProjData (void);
+  void newProjData ();
+  void deleteProjData ();
 
   void init (const int nView, const int nDet);