iScannerSPECTConv.hh

Go to the documentation of this file.

/*
This file is part of CASToR.

    CASToR is free software: you can redistribute it and/or modify it under the
    terms of the GNU General Public License as published by the Free Software
    Foundation, either version 3 of the License, or (at your option) any later
    version.

    CASToR is distributed in the hope that it will be useful, but WITHOUT ANY
    WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
    FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
    details.

    You should have received a copy of the GNU General Public License along with
    CASToR (in file GNU_GPL.TXT). If not, see <http://www.gnu.org/licenses/>.

Copyright 2017-2019 all CASToR contributors listed below:

    --> Didier BENOIT, Claude COMTAT, Marina FILIPOVIC, Thibaut MERLIN, Mael MILLARDET, Simon STUTE, Valentin VIELZEUF

This is CASToR version 3.0.
*/

#ifndef ISCANNERSPECTConv_HH
#define ISCANNERSPECTConv_HH 1

#include "gVariables.hh"
#include "vScanner.hh"
#include "sAddonManager.hh"

class iScannerSPECTConv : public vScanner
{
  // -------------------------------------------------------------------
  // Constructor & Destructor
  public:
    iScannerSPECTConv();
    ~iScannerSPECTConv();


  // -------------------------------------------------------------------
  // Public member functions
  public:
    // Function for automatic insertion (put the class name as the parameters and do not add semi-colon at the end of the line)
    FUNCTION_SCANNER(iScannerSPECTConv)
    void DescribeSpecific();
    int Instantiate(bool a_scannerFileIsLUT);
    int CheckParameters();
    int Initialize();
    int BuildLUT( bool a_scannerFileIsLUT );
    int GetPositionsAndOrientations( int a_index1, int a_index2,
                                     FLTNB ap_Position1[3], FLTNB ap_Position2[3],
                                     FLTNB ap_Orientation1[3], FLTNB ap_Orientation2[3],
                                     FLTNB* ap_POI1 = NULL, FLTNB* ap_POI2 = NULL );
    int GetRdmPositionsAndOrientations( int a_index1, int a_index2,
                                        FLTNB ap_Position1[3], FLTNB ap_Position2[3],
                                        FLTNB ap_Orientation1[3], FLTNB ap_Orientation2[3] );
    int GetPositionWithRandomDepth( int a_index1, int a_index2, FLTNB ap_Position1[3], FLTNB ap_Position2[3] );
    int GetTwoCorners( int a_index1, int a_index2,
                       FLTNB ap_CornerInf1[3], FLTNB ap_CornerSup1[3],
                       FLTNB ap_CornerInf2[3], FLTNB ap_CornerSup2[3] );
    int GetEdgesCenterPositions( int a_index1, int a_index2,
      FLTNB ap_pos_line_point1[3], FLTNB ap_pos_line_point2[3],
      FLTNB ap_pos_point1_x[4], FLTNB ap_pos_point1_y[4], FLTNB ap_pos_point1_z[4],
      FLTNB ap_pos_point2_x[4], FLTNB ap_pos_point2_y[4], FLTNB ap_pos_point2_z[4] );
    int GetGeometricInfoFromDataFile( string a_pathToDF );
    inline int GetSystemNbElts()
           {return m_nbCrystals;}
    void ShowHelp();
    inline FLTNB GetDetectionElementSizeTrans()
           {return m_vPixelsSizeTrans;}
    inline FLTNB GetDetectionElementSizeAxial()
           {return m_vPixelsSizeTrans;}

  // -------------------------------------------------------------------
  // Functions dedicated to Analytic Projection
  public:
    inline int IsAvailableLOR(int a_elt1, int a_elt2)
           {return 1;};
    int GetSPECTSpecificParameters( uint16_t* ap_nbOfProjections, 
                                    uint16_t* ap_nbHeads,
                                    FLTNB* ap_acquisitionZoom, 
                                    uint16_t* ap_nbOfBins,
                                    FLTNB*  ap_pixSizeXY, 
                                    FLTNB*& ap_angles, 
                                    FLTNB*& ap_CORtoDetectorDistance,
                                    int* ap_headRotDirection);
    int PROJ_SetSPECTAngles(FLTNB* ap_projectionAngles);
    int PROJ_SetSPECTCORtoDetectorDistance(FLTNB a_distance);
    inline uint16_t PROJ_GetSPECTNbProjections()
           {return m_nbOfProjections;}
    inline uint16_t PROJ_GetSPECTNbPixels()
           {return m_vNbPixelsTrans*m_vNbPixelsAxial;}
    int PROJ_GetSPECTNbBins(uint16_t* ap_nbOfBins);
    int PROJ_SetSPECTNbBins(uint16_t* ap_nbOfBins);
    int PROJ_SetSPECTNbProjections(uint32_t a_nbOfProjections);

  // -------------------------------------------------------------------
  // Private member functions
  private:
    int LoadLUT();
    int ComputeLUT();
    int ComputeFocalPositions( FLTNB a_posX, FLTNB a_posY, FLTNB a_posZ, int a_headID, int a_cryID );
  
  
  // -------------------------------------------------------------------
  // Data members
  private:
    int m_nbCrystals;                  
    int m_nbHeads;                     
    uint16_t mp_nbOfBins[2];           
    uint16_t m_nbOfProjections;        
    FLTNB* mp_projectionAngles;        
    FLTNB* mp_CORtoDetectorDistance;   
    FLTNB* mp_radius;                  
    uint32_t m_nbPixelsTrans;          
    FLTNB m_pixelsSizeTrans;           
    FLTNB m_gapSizeTrans;              
    uint32_t m_nbPixelsAxial;          
    FLTNB m_pixelsSizeAxial;           
    FLTNB m_gapSizeAxial;              
    FLTNB m_acquisitionZoom;                      
    uint32_t m_vNbPixelsTrans;         
    uint32_t m_vNbPixelsAxial;         
    FLTNB m_vPixelsSizeTrans;          
    FLTNB m_vPixelsSizeAxial;          
    FLTNB m_crystalDepth;              
    string* mp_focalModelTrans;        
    uint8_t* mp_nbCoefModelTrans;      
    FLTNB** m2p_transFocalParameters;  
    string* mp_focalModelAxial;        
    uint8_t* mp_nbCoefModelAxial;      
    FLTNB** m2p_axialFocalParameters;  
    FLTNB* mp_crystalCentralPositionX; 
    FLTNB* mp_crystalCentralPositionY; 
    FLTNB* mp_crystalCentralPositionZ; 
    FLTNB* mp_crystalOrientationX;     
    FLTNB* mp_crystalOrientationY;     
    FLTNB* mp_crystalOrientationZ;     
    FLTNB* mp_crystalFocalPositionX;   
    FLTNB* mp_crystalFocalPositionY;   
    FLTNB* mp_crystalFocalPositionZ;   
};


// Class for automatic insertion (set here the visible scanner type name, put the class name as the parameters and do not add semi-colon at the end of the line)
CLASS_SCANNER(SPECT_CONVERGENT,iScannerSPECTConv)

#endif