mitk2DBackProjectionUsingGPU.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 __mitk2DBackProjectionUsingGPU_h
#define __mitk2DBackProjectionUsingGPU_h

#include "mitkBackProjectionUsingGPU.h"
#include <GL/glew.h>
#include <GL/glut.h>
#include <Cg/cgGL.h>

#define BLOCKNUM 7
#define FANWINDOWSIZE 7






class MITK_RECONSTRUCTION_CG_API mitk2DBackProjectionUsingGPU : public mitkBackProjectionUsingGPU
{
public:
    MITK_TYPE(mitk2DBackProjectionUsingGPU,mitkBackProjectionUsingGPU)

    mitk2DBackProjectionUsingGPU();  

    virtual void PrintSelf(ostream& os);

    virtual void InputProjectionData(float * data);

    virtual void Draw (int mBottom);

    virtual void Rotate(float angle);

    virtual void ReadResult(int mSubSliceNum,float * result);

    virtual void SetParameter(int xrayNum,int angleNum,int sliceNum,double begenAngle,double endAngle,double SOD,double SDD,double lengthDetector,double *fov,int interpolationModel,int width,int height);

    virtual bool CreateTexture();

     
private:
    void initWindow();

    void initFBO();

    void initCG();

    bool _createTexture();

    void AttachFBO();
     
    void EnableFshader();

    void EnableVshader();

    void SetDtheta(float m) {dthetaDegree = m;}

    void SetFrustum(float zn,float zf,float zc, float detectorLength);

    void DestroyAll();

    void CheckForCgError(const char *situation);

    void _bindCurrentFBO();
    void _unBind();
    void _scanTexParams(CGprogram m_prog);
    void _setTexture(char *name, GLuint tex);
    void _bindTextures();



    void SetXrayNum(int m){xrayNum = m;}
    void SetAngleNum(int m){angleNum = m;}
    void SetWidth(int m){width = m;}
    void SetHeight(int m){height = m;}
    void SetImgslice(int m){imgslice = m;}
    void SetFov(double *m)
    {
        fov[0] = m[0]; 
        fov[1] = m[1]; 
        fov[2] = m[2]; 
        fov[3] = m[3];
    }
    void SetFov(double xmin,double xmax,double ymin,double ymax)
    {
        fov[0] = xmin; 
        fov[1] = xmax; 
        fov[2] = ymin; 
        fov[3] = ymax;
    }
    int xrayNum;
    int angleNum;
    int width;
    int height;
    int imgslice;
    double fov[4];



    struct TexParam 
    {
        const char *name;
        int unit;
        int target;
        int tex;
    };
    std::vector<TexParam> m_texParams;

    float dthetaDegree;
    float zNear;
    float zFar;     //L
    float zCenter;  //D
    float nearDetectorLengthX;
    GLfloat frustum[16];

    GLuint windowHandle;    
    CGcontext cgContext;
    CGprofile fProfile ,vProfile;
    CGprogram fProgram ,vProgram;
    CGparameter mvpMatrixParam,preResultTexParam,
        texMatrixParam[BLOCKNUM];

    GLuint *projectDataID,  fbotex, *ResultID;
    GLuint fb;
    GLint m_preFboId;

    
protected:
    virtual ~mitk2DBackProjectionUsingGPU();
      
private:

    mitk2DBackProjectionUsingGPU(const mitk2DBackProjectionUsingGPU&);
    void operator=(const mitk2DBackProjectionUsingGPU&);
};


//#define DEFINED_mitk2DBackProjectionUsingGPU

#endif

---