Eigene Dateien/uni/visualisierung/vislu/bsp1/VisLu/VObject.cpp

Go to the documentation of this file.

#include "VObject.h"

VObject::VObject(){
}

void VObject::load(string _filename){

        if (VolumeTex != 0){
                glDeleteTextures(1,&VolumeTex);
        }
        
        direction = 0;
        sliceNo = -1;
        chSlice = new int(0.0f);
        
        unsigned short *volumeData;
        
        FILE * pFile;
        pFile = fopen(_filename.c_str(),"rb");
        
        if (pFile == NULL){
                exit(1);
        }
        /*
        char charX [2], charY[2], charZ[2];
        
        fgets(charX,3,pFile);
        fgets(charY,3,pFile);
        fgets(charZ,3,pFile);
        
        unsigned short sizeX = *charX;
        unsigned short sizeY = *charY;
        unsigned short sizeZ = *charZ;
        */
        unsigned short sizeX = getc(pFile);
                                        sizeX |= getc(pFile) << 8;
        unsigned short sizeY = getc(pFile);
                                        sizeY |= getc(pFile) << 8;
        unsigned short sizeZ = getc(pFile);
                                        sizeZ |= getc(pFile) << 8;
        
        int twoSizeX = getPowerofTwo(sizeX);
        int twoSizeY = getPowerofTwo(sizeY);
        int twoSizeZ = getPowerofTwo(sizeZ);
        //      volumeData = new float[sizeX*sizeY*sizeZ*4];
        // volumeData = new float[twoSizeX*twoSizeY*twoSizeZ*4];
        volumeData = new unsigned short[twoSizeX*twoSizeY*twoSizeZ];
        histogram = new int[4096];
        
        int i = 0;
        
        for(int arrayCount = 0; arrayCount < 4096; arrayCount++){
                histogram[arrayCount] = 0;
        }
        
        for(int z=0; z < twoSizeZ; z++) 
        {                
                //iterate over the volume
                for(int y=0; y < twoSizeY; y++)
                {       
                        for(int x=0; x < twoSizeX; x++)
                        {       
                                if ((x>=sizeX) || (y>=sizeY) || (z>=sizeZ)){
                                        volumeData[i] = 0;
                                        }
                                else {                                  
                                        unsigned short integer = getc(pFile);
                                        integer |= getc(pFile) << 8;
                                        
                                        (histogram[integer])++;
                                        volumeData[i] = integer;
                                        }
                                i++;
                        }       
                }       
        }
        
        fclose(pFile);
        
        glGenTextures(1, &VolumeTex);
        glBindTexture(GL_TEXTURE_3D, VolumeTex);
        glTexParameteri(GL_TEXTURE_3D, GL_TEXTURE_WRAP_S, GL_CLAMP);
        glTexParameteri(GL_TEXTURE_3D, GL_TEXTURE_WRAP_T, GL_CLAMP);
        glTexParameteri(GL_TEXTURE_3D, GL_TEXTURE_WRAP_R, GL_CLAMP);
        glTexParameteri(GL_TEXTURE_3D, GL_TEXTURE_MAG_FILTER, GL_LINEAR);
        glTexParameteri(GL_TEXTURE_3D, GL_TEXTURE_MIN_FILTER, GL_LINEAR);
        
        glTexImage3D(GL_TEXTURE_3D, 0, GL_LUMINANCE, twoSizeX, twoSizeY, twoSizeZ, 0, GL_LUMINANCE, GL_UNSIGNED_SHORT, volumeData);
        
        //glTexImage3D(GL_TEXTURE_3D, 0, 4, twoSizeX, twoSizeY, twoSizeZ, 0, GL_RGBA, GL_FLOAT, volumeData);
                
        delete[] volumeData;
}

void VObject::draw(){

        sliceNo += *chSlice;
        glEnable(GL_BLEND);
        glBlendFunc(GL_SRC_ALPHA,GL_DST_ALPHA);
        
        glEnable(GL_TEXTURE_3D);
        glTexEnvf(GL_TEXTURE_ENV, GL_TEXTURE_ENV_MODE, GL_REPLACE);
        
        glBindTexture(GL_TEXTURE_3D, VolumeTex);
        float count = 64.0f;
        if (sliceNo >= count)
                        sliceNo = count;
        else if (sliceNo < 0)
                        sliceNo =-1;
        glRotatef(-90,1,0,0);
        /*
         glBegin(GL_QUADS);
         glTexCoord3d( 0.0, 0.0, 0.3);
         glVertex3d(0.0, 0.0,0.3);
         glTexCoord3d(0.0, 1.0,0.3);
         glVertex3d(0.0, 1.0,0.3);
         glTexCoord3d(1.0, 1.0, 0.3);
         glVertex3d(1.0, 1.0, 0.3);
         glTexCoord3d(1.0, 0.0, 0.3);
         glVertex3d(1.0, 0.0, 0.3);
         glEnd();
         */     
        if ((direction == 0) || (sliceNo == -1)){
                for (int i = 0; i < count; i++){
                        glBegin(GL_QUADS);
                        glTexCoord3d( 0.0, 0.0, (float)i/count);
                        glVertex3d(0.0, 0.0,(float)i/count);
                        glTexCoord3d(0.0, 1.0, (float)i/count);
                        glVertex3d(0.0, 1.0, (float)i/count);
                        glTexCoord3d(1.0, 1.0, (float)i/count);
                        glVertex3d(1.0, 1.0, (float)i/count);
                        glTexCoord3d(1.0, 0.0, (float)i/count);
                        glVertex3d(1.0, 0.0, (float)i/count);
                        glEnd();
                }
        }
        else if (direction == 1){
                glBegin(GL_QUADS);
                glTexCoord3d( (float)sliceNo/count,0.0, 0.0);
                glVertex3d((float)sliceNo/count,0.0, 0.0);
                glTexCoord3d((float)sliceNo/count, 1.0, 0.0);
                glVertex3d((float)sliceNo/count,1.0, 0.0);
                glTexCoord3d((float)sliceNo/count,1.0, 1.0);
                glVertex3d((float)sliceNo/count,1.0, 1.0);
                glTexCoord3d((float)sliceNo/count,0.0, 1.0);
                glVertex3d((float)sliceNo/count,0.0, 1.0);
                glEnd();                        
        }
        
        else if (direction == 2){
                glBegin(GL_QUADS);
                glTexCoord3d( 0.0,(float)sliceNo/count, 0.0);
                glVertex3d(0.0, (float)sliceNo/count, 0.0);
                glTexCoord3d(1.0, (float)sliceNo/count, 0.0);
                glVertex3d(1.0, (float)sliceNo/count, 0.0);
                glTexCoord3d(1.0, (float)sliceNo/count, 1.0);
                glVertex3d(1.0, (float)sliceNo/count, 1.0);
                glTexCoord3d(0.0, (float)sliceNo/count, 1.0);
                glVertex3d(0.0, (float)sliceNo/count, 1.0);
                glEnd();
        }
        else if (direction == 3){
                glBegin(GL_QUADS);
                glTexCoord3d( 0.0,0.0,(float)sliceNo/count);
                glVertex3d(0.0, 0.0,(float)sliceNo/count);
                glTexCoord3d(1.0, 0.0,(float)sliceNo/count);
                glVertex3d(1.0, 0.0,(float)sliceNo/count);
                glTexCoord3d(1.0, 1.0,(float)sliceNo/count);
                glVertex3d(1.0, 1.0,(float)sliceNo/count);
                glTexCoord3d(0.0, 1.0,(float)sliceNo/count);
                glVertex3d(0.0, 1.0,(float)sliceNo/count);
                glEnd();
        }
        glDisable(GL_TEXTURE_3D);
        /*
         int count1 = 0;
         int count2 = 0;
         int count3 = 0;
         
         for (int i = 0; i < points.size(); ++i){
                 //TODO: Transferfkt
                 if (points[i]->density <= 50000){
                         glColor4fv(Color1);
                         count2++;
                 }
                 else {
                         glColor4fv(Color3);
                         count3++;
                         glBegin(GL_POINTS);
                         glVertex3f((points[i]->x),(points[i]->y),(points[i]->z));
                         glEnd();
                 }
                 
         }
         */
}

int* VObject::getHistogram(void){
        return histogram;
}

void VObject::setDirection(int dir){
        direction = dir;
}

int VObject::getDirection(void){
        return direction;
}

void VObject::setsliceNo(int slice){
        sliceNo = slice;
}

int* VObject::getsliceNo(void){
        return this->chSlice;
}

void VObject::changesliceNo(int ch){
        sliceNo += ch;
}

int VObject::getPowerofTwo(int _n){
        int i = 0;
        while (pow(2.0f,i) < _n){
                i++;
        }
        return pow(2.0f,i);
}

VObject::~VObject(){
}

---