** 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: projections.h,v 1.19 2001/01/03 22:00:46 kevin Exp $
+** $Id: projections.h,v 1.25 2001/02/22 18:22:40 kevin Exp $
**
**
** This program is free software; you can redistribute it and/or modify
class Scanner;
class DetectorArray;
class Array2dFileLabel;
+class fnetorderstream;
+#include "array2dfile.h"
+#include "imagefile.h"
+#include <complex>
// Projections
class Projections
{
public:
+
+ static const int POLAR_INTERP_INVALID;
+ static const int POLAR_INTERP_NEAREST;
+ static const int POLAR_INTERP_BILINEAR;
+ static const int POLAR_INTERP_BICUBIC;
+
Projections (const int nView, const int nDet);
Projections (const Scanner& scanner);
Projections ();
~Projections ();
+ static const int getInterpCount() {return s_iInterpCount;}
+ static const char* const* getInterpNameArray() {return s_aszInterpName;}
+ static const char* const* getInterpTitleArray() {return s_aszInterpTitle;}
+ static int convertInterpNameToID (const char* const interpName);
+ static const char* convertInterpIDToName (const int interpID);
+ static const char* convertInterpIDToTitle (const int interpID);
+
void initFromScanner (const Scanner& scanner);
void printProjectionData (int startView, int endView);
bool convertPolar (ImageFile& rIF, int iInterpolation);
bool convertFFTPolar (ImageFile& rIF, int iInterpolation, int iZeropad);
+ bool calcArrayPolarCoordinates (unsigned int nx, unsigned int ny, double** ppdView, double** ppdDet);
+ void interpolatePolar (ImageFileArray& v, ImageFileArray& vImag, unsigned int nx, unsigned int ny, std::complex<double>** ppcDetValue,
+ double** ppdDet, double** ppdView, unsigned int nView, unsigned int nDet, int iInterpolate);
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;}
void setDetInc (double detInc) { m_detInc = detInc;}
- void setPhmLen (double phmLen) { m_phmLen = phmLen;}
void setCalcTime (double calcTime) {m_calcTime = calcTime;}
void setRemark (const char* remark) {m_remark = remark; m_label.setLabelString(remark);}
void setRemark (const std::string& remark) {setRemark(remark.c_str());}
double detStart() const {return m_detStart;}
double rotStart() const {return m_rotStart;}
double calcTime() const {return m_calcTime;}
- double phmLen() const {return m_phmLen;}
+ double viewLen() const {return m_dViewDiameter / SQRT2;}
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;}
+ double focalLength() const {return m_dFocalLength;}
+ double phmLen() const { return m_dViewDiameter / SQRT2; }
+ void setPhmLen(double phmLen) { m_dViewDiameter = phmLen * SQRT2; }
const std::string& getFilename() const {return m_filename;}
Array2dFileLabel& getLabel() {return m_label;}
double m_rotInc; // angle between rotations
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;
+ double m_dFocalLength;
+ double m_dViewDiameter;
+ double m_dFanBeamAngle;
kuint32 m_year; // Creation date & time
kuint32 m_month;
kuint32 m_day;
const static kuint16 m_signature;
+ static const char* const s_aszInterpName[];
+ static const char* const s_aszInterpTitle[];
+ static const int s_iInterpCount;
+
bool headerWrite (fnetorderstream& fs);
bool headerRead (fnetorderstream& fs);
void newProjData ();