D:/Timo/Visualisierung/LU/Aufgabe1/src/Volume.cpp

Go to the documentation of this file.

#include "Volume.h"
#include <QtDebug>
#include <QFileDialog>
#include <QMessageBox>

Volume::Volume() : width(1), height(1), depth(1), voxels(1), histogramm(1), heightHist(0)
{
        
};

Volume::~Volume(void)
{
};

const int Volume::getWidth() const
{
        return width;
};

const int Volume::getHeight() const
{
        return height;
};

const int Volume::getDepth() const
{
        return depth;
};

const int Volume::getSize() const
{
        return int(voxels.size());
};

const unsigned short* Volume::getData() const
{
        return &(voxels.front());
}

const std::vector<unsigned int>* Volume::getHistogramm() const
{
        return &histogramm;
}

const unsigned int Volume::getHistHeight() const 
{
        return heightHist;
}

const int Volume::getDataRange() const
{
        return 4096;
}

void Volume::loadDatFile() 
{
         QString fileName = QFileDialog::getOpenFileName(0,
         "Load Data", curDir,
         "Vis Data (*.dat);;All Files (*)");
         if(fileName!=0) {
                 curDir = fileName;
                 qDebug() << "Pfad vom File: " << fileName.toStdString().c_str();
                 loadDatFile(fileName.toStdString().c_str());
         }
}

void Volume::loadDatFile(const std::string &fileName)
{
        qDebug() << "- Loading file '" << fileName.c_str() << "' ... ";
        FILE *file = NULL;
        
        fopen_s(&file, fileName.c_str(), "rb");

        if(!file)
        {
                QMessageBox::information(0, "Unable to load file", "");
        }
        else
        {
                unsigned short tempShort;

                fread(&tempShort, sizeof(unsigned short), 1, file);
                width = int(tempShort);
                fread(&tempShort, sizeof(unsigned short), 1, file);
                height = int(tempShort);
                fread(&tempShort, sizeof(unsigned short), 1, file);
                depth = int(tempShort);

                int size;
                // histogramm initialisieren
                histogramm.resize(4096);
                for(unsigned int j=0;j<histogramm.size();j++) {
                        histogramm[j] = 0;
                }

                if(width % 2)
                {
                        int realWidth = width;
                        width = realWidth + realWidth % 2;
                        size = width * height * depth;
                        voxels.resize(size);

                        unsigned short *voxelPointer = &(voxels.front());
                        for(int i = 0;i < height * depth;++i) {
                                fread(voxelPointer + i * width, sizeof(unsigned short), realWidth, file);
                                voxelPointer[i * width + realWidth] = 0;
                        }
                }
                else
                {
                        size = width * height * depth;
                        voxels.resize(size);
                        fread((void*)&(voxels.front()), sizeof(unsigned short), size, file);
                }
                fclose(file);

                for(int i = 0;i < size;++i)
                {
                        histogramm[voxels[i]]++;
                        voxels[i] <<= 4;
                }

                heightHist = 0;
                for(unsigned int k = 0;k <histogramm.size();k++)
                {       
                        if(heightHist<histogramm[k])
                        {
                                heightHist = histogramm[k];
                        }
                }
        
                qDebug() << "- File loaded. Size: " << width << height << depth;

                emit loadedDatFile(this);
                QString temp = QString(fileName.c_str());
                QStringList tempList = temp.split("/");
                emit loadedDatFileName(tempList.last());
        }
}

---