X-Git-Url: http://git.kpe.io/?p=ctsim.git;a=blobdiff_plain;f=include%2Fprojections.h;fp=include%2Fprojections.h;h=39cb09a7a67446a3723a7cb06a6921bc6ddc68a7;hp=610cae76d3dc9b63189fba45a15e65739e3bfc9b;hb=b361677a2f7d5b443641faec70b057f2a84dc77e;hpb=af74e42547cbc35de784e30ebdbba746ae0fc9e5 diff --git a/include/projections.h b/include/projections.h index 610cae7..39cb09a 100644 --- a/include/projections.h +++ b/include/projections.h @@ -9,7 +9,7 @@ ** This is part of the CTSim program ** Copyright (c) 1983-2001 Kevin Rosenberg ** -** $Id: projections.h,v 1.27 2001/03/02 02:08:14 kevin Exp $ +** $Id: projections.h,v 1.28 2001/03/10 23:14:15 kevin Exp $ ** ** ** This program is free software; you can redistribute it and/or modify @@ -38,6 +38,36 @@ class fnetorderstream; #include "imagefile.h" #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 + int m_iViewNum; // View number for this raysum + int m_iDetNum; // Detector number for this raysum + + bool lessThanT (ParallelRaysumCoordinate& rCompare) + {return m_dT < rCompare.m_dT; } + bool lessThanTheta (ParallelRaysumCoordinate& rCompare) + {return m_dTheta < rCompare.m_dTheta; } +}; + +class ParallelRaysums { +private: + ParallelRaysumCoordinate** m_ppCoordinates; + int m_iNumCoordinates; + +public: + ParallelRaysums (Projections* pProjections); + ~ParallelRaysums (); + + ParallelRaysumCoordinate** getCoordinates() {return m_ppCoordinates;} + int getNumCoordinates() const {return m_iNumCoordinates;} + void getLimits (double* dMinT, double* dMaxT, double* dMinTheta, double* dMaxTheta) const; +}; + + + // Projections class Projections { @@ -74,6 +104,8 @@ class Projections bool detarrayRead (fnetorderstream& fs, DetectorArray& darray, const int view_num); bool detarrayWrite (fnetorderstream& fs, const DetectorArray& darray, const int view_num); + bool interpolateToParallel(); + 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);