X-Git-Url: http://git.kpe.io/?p=ctsim.git;a=blobdiff_plain;f=include%2Fprojections.h;h=3efbde920f64aab989976457266eafedfd082ac1;hp=35325dfe370ca34e1f0b5227d4819e2aae17680a;hb=f13a8c004b8f182b42d9e4df2bcd7c7f030bf1ad;hpb=11ea72090da18aa290e2c3ab217f1d27736bca50 diff --git a/include/projections.h b/include/projections.h index 35325df..3efbde9 100644 --- a/include/projections.h +++ b/include/projections.h @@ -7,10 +7,7 @@ ** Date Started: July 1, 1984 ** ** This is part of the CTSim program -** Copyright (c) 1983-2001 Kevin Rosenberg -** -** $Id: projections.h,v 1.24 2001/02/20 17:44:14 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 @@ -37,6 +34,61 @@ class fnetorderstream; #include "array2dfile.h" #include "imagefile.h" #include +#include + + +//used for rebinning divergent beam projections to parallel +class ParallelRaysumCoordinate { +public: + double m_dT; // Distance from center of origin + double m_dTheta; // perpendicular angle to origin + double m_dRaysum; + + static bool compareByTheta (ParallelRaysumCoordinate* a, ParallelRaysumCoordinate* b) + { return a->m_dTheta == b->m_dTheta ? b->m_dT > a->m_dT : b->m_dTheta > a->m_dTheta; } + + // sort first by T, then Theta + static bool compareByT (ParallelRaysumCoordinate* a, ParallelRaysumCoordinate* b) + { return a->m_dT == b->m_dT ? b->m_dTheta > a->m_dTheta : b->m_dT > a->m_dT; } + +}; + +class ParallelRaysums { +public: + + enum { + THETA_RANGE_UNCONSTRAINED = 0, + THETA_RANGE_NORMALIZE_TO_TWOPI, + THETA_RANGE_FOLD_TO_PI, + }; + + ParallelRaysums (const Projections* pProjections, int iThetaRange); + ~ParallelRaysums (); + + typedef std::vector CoordinateContainer; + typedef CoordinateContainer::iterator CoordinateIterator; + + CoordinateContainer& getCoordinates() {return m_vecpCoordinates;} + int getNumCoordinates() const {return m_iNumCoordinates;} + void getLimits (double* dMinT, double* dMaxT, double* dMinTheta, double* dMaxTheta) const; + CoordinateContainer& getSortedByT(); + CoordinateContainer& getSortedByTheta(); + double interpolate (double* pdX, double* pdY, int n, double dXValue, int* iLastFloor = NULL); + void getThetaAndRaysumsForT (int iT, double* pdTheta, double* pdRaysum); + void getDetPositions (double* pdDetPos); + +private: + CoordinateContainer m_vecpCoordinates; + CoordinateContainer m_vecpSortedByT; + CoordinateContainer m_vecpSortedByTheta; + ParallelRaysumCoordinate* m_pCoordinates; + int m_iNumCoordinates; + int m_iNumView; + int m_iNumDet; + int m_iThetaRange; +}; + + // Projections class Projections @@ -54,18 +106,23 @@ class Projections ~Projections (); static const int getInterpCount() {return s_iInterpCount;} - static const char** getInterpNameArray() {return s_aszInterpName;} - static const char** getInterpTitleArray() {return s_aszInterpTitle;} + 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); + bool initFromSomatomAR_STAR (int iNViews, int iNDets, unsigned char* pData, unsigned long lDataLength); void printProjectionData (int startView, int endView); void printProjectionData (); void printScanInfo (std::ostringstream& os) const; + int Helical180LI(int interpView); + int Helical180LI_Equiangular(int interpView); + int HalfScanFeather(void); + bool read (const std::string& fname); bool read (const char* fname); bool write (const char* fname); @@ -73,11 +130,15 @@ class Projections bool detarrayRead (fnetorderstream& fs, DetectorArray& darray, const int view_num); bool detarrayWrite (fnetorderstream& fs, const DetectorArray& darray, const int view_num); + Projections* interpolateToParallel() const; + 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 calcArrayPolarCoordinates (unsigned int nx, unsigned int ny, double** ppdView, double** ppdDet, + int iNumDetWithZeros, double dZeropadRatio, double dDetInc); void interpolatePolar (ImageFileArray& v, ImageFileArray& vImag, unsigned int nx, unsigned int ny, std::complex** ppcDetValue, - double** ppdDet, double** ppdView, unsigned int nView, unsigned int nDet, int iInterpolate); + double** ppdDet, double** ppdView, unsigned int nView, unsigned int nDet, unsigned int nDetWithZeros, + 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; @@ -99,6 +160,8 @@ class Projections int nView() const {return m_nView;} int geometry() const {return m_geometry;} double focalLength() const {return m_dFocalLength;} + double sourceDetectorLength() const { return m_dSourceDetectorLength;} + double viewDiameter() const {return m_dViewDiameter; } double phmLen() const { return m_dViewDiameter / SQRT2; } void setPhmLen(double phmLen) { m_dViewDiameter = phmLen * SQRT2; } @@ -108,7 +171,7 @@ class Projections DetectorArray& getDetectorArray (const int iview) { return (*m_projData[iview]); } - + const DetectorArray& getDetectorArray (const int iview) const { return (*m_projData[iview]); } @@ -117,20 +180,21 @@ class Projections static bool copyViewData (const char* const filename, std::ostream& os, int startView, int endView); static bool copyViewData (const std::string& filename, std::ostream& os, int startView, int endView); - + private: - int m_headerSize; // Size of disk file header - int m_geometry; // Geometry of scanner - class DetectorArray **m_projData; // Pointer to array of detarray_st pointers - std::string m_remark; // description of raysum data - int m_nDet; // number of detectors in array - int m_nView; // number of rotated views - double m_calcTime; // time required to calculate raysums - double m_rotStart; // starting view rotation - double m_rotInc; // angle between rotations - double m_detStart; // distance of beginning detector to center phantom - double m_detInc; // increment between detectors + int m_headerSize; // Size of disk file header + int m_geometry; // Geometry of scanner + class DetectorArray **m_projData; // Pointer to array of detarray_st pointers + std::string m_remark; // description of raysum data + int m_nDet; // number of detectors in array + int m_nView; // number of rotated views + double m_calcTime; // time required to calculate raysums + double m_rotStart; // starting view rotation + double m_rotInc; // angle between rotations + double m_detStart; // distance of beginning detector to center phantom + double m_detInc; // increment between detectors double m_dFocalLength; + double m_dSourceDetectorLength; double m_dViewDiameter; double m_dFanBeamAngle; kuint32 m_year; // Creation date & time @@ -144,8 +208,8 @@ class Projections const static kuint16 m_signature; - static const char* s_aszInterpName[]; - static const char* s_aszInterpTitle[]; + static const char* const s_aszInterpName[]; + static const char* const s_aszInterpTitle[]; static const int s_iInterpCount; bool headerWrite (fnetorderstream& fs);