mitkCTReconstructionFilter.h

/*=========================================================================

  Program:   3DMed
  Date:      $Date: 2014-02-25 18:30:00 +0800 $
  Version:   $Version: 4.6.0 $
  Copyright: MIPG, Institute of Automation, Chinese Academy of Sciences

=========================================================================*/


#ifndef __mitkCTReconstructionFilter_h
#define __mitkCTReconstructionFilter_h

#include "mitkFilter.h"
#include "mitkVolume.h"
#include "mitkCTProjectionData.h"
#include "mitkRCPtr.h"
#include "mitkRawDataToVolumeFilter.h"
#define PI 3.14159265358979


#ifndef USE_SIMD
#define USE_SIMD
#endif

class MITK_RECONSTRUCTION_API mitkCTReconstructionFilter : public mitkRawDataToVolumeFilter  
{
public:
    MITK_TYPE(mitkCTReconstructionFilter,mitkRawDataToVolumeFilter)

    virtual void PrintSelf(ostream& os);

    //enum InterpolationModel
    //{
    //  NEAREST,
    //  BILINEAR
    //};

    void SetBeginAngle(double angle)  { beginAngle = angle; }
    
    double GetBeginAngle() const { return beginAngle; }
    
    void SetEndAngle(double angle)  { endAngle = angle; }
    
    double GetEndAngle() const { return endAngle; }

    void SetInterpolationModel(int m){mInterp = m;}

    int GetInterpolationModel(){return mInterp;}

    void SetInput(mitkCTProjectionData *inData) {m_InData = inData;}
    
    mitkCTProjectionData* GetInput() {return m_InData;}
    
    mitkVolume* GetOutput();

    void SetOutput(mitkVolume *m) {m_OutData = m;}

    void SetIsLogOperation(bool m) {isLogOperation = m;}

    bool GetIsLogOperation() {return isLogOperation;}

    void SetEmptyScanImage(mitkVolume *m){emptyImage = m;}

    mitkVolume* GetEmptyScanImage(){return emptyImage;}
    

    float GetFilterFactor(){return m_filterFactor;}
    int GetFilterType(){return m_filterType;}
    void SetFilterType(int type=4, float factor=0.46){m_filterType=type; m_filterFactor=factor;}

protected:

    mitkCTReconstructionFilter();

    virtual ~mitkCTReconstructionFilter();

protected:  
    mitkRCPtr<mitkCTProjectionData> m_InData;

    mitkRCPtr<mitkVolume> emptyImage;
 // mitkVolume *m_OutData;
    bool isLogOperation;
    int mInterp;
    double beginAngle;
    double endAngle;
    float m_filterFactor;//set the factor to change filter window in ct filter
    int m_filterType;// m_filterType:1 for ramp, 2 for Blackman, 3 for Hanning, 4 for Hamming, 5 for self define (1-m_filterFactor+m_filterFactor*cos(2*PI*ii/L))

private:
    mitkCTReconstructionFilter(const mitkCTReconstructionFilter&);
    void operator=(const mitkCTReconstructionFilter&);

};


//#define DEFINED_mitkCTReconstructionFilter


#endif 

---