src/Volume.cpp

Go to the documentation of this file.

#include "Volume.h"


Volume::Volume()
{
  release();
}

Volume::Volume( const std::string &strFilename ) 
{
  release();
        load( strFilename );
}

Volume::~Volume(void)
{
  release();
}


void Volume::load( const std::string &strFilename )
{
        std::cout << "- Loading file '" << strFilename << "' ... " << std::endl;
        
  FILE *fp = NULL;
        fopen_s( &fp, strFilename.c_str(), "rb" );

        if (!fp)
                std::cerr << "- Error loading file!" << std::endl << std::endl;
        else
        {

                /* who needs this??
    char vcPath[1024];
                char *pFileName = NULL;
                GetFullPathName(strFilename.c_str(),1024,vcPath,&pFileName);
                char vcDrive[1024], vcDirectory[1024], vcFilename[1024], vcExtension[1024];
                _splitpath_s(vcPath,vcDrive,vcDirectory,vcFilename,vcExtension);
                const std::string strAdditionalFilename =  std::string(vcDrive)+std::string(vcDirectory)+std::string(vcFilename)+std::string(".ini");

                char vpSpacingX[1024],vpSpacingY[1024],vpSpacingZ[1024];
                GetPrivateProfileString("DatFile","oldDat Spacing X","1.0",vpSpacingX,256,strAdditionalFilename.c_str());
                GetPrivateProfileString("DatFile","oldDat Spacing Y","1.0",vpSpacingY,256,strAdditionalFilename.c_str());
                GetPrivateProfileString("DatFile","oldDat Spacing Z","1.0",vpSpacingZ,256,strAdditionalFilename.c_str());
    */

                //unsigned short uWidth, uHeight, uDepth;

                unsigned short sizeX, sizeY, sizeZ;
                fread( &sizeX, sizeof(unsigned short), 1, fp );
                fread( &sizeY, sizeof(unsigned short), 1, fp );
                fread( &sizeZ, sizeof(unsigned short), 1, fp );

                m_sizeX = int(sizeX);
                m_sizeY = int(sizeY);
                m_sizeZ = int(sizeZ);

    std::cout << " - Volume-dimensions: (" << 
      m_sizeX << ", " <<
      m_sizeY << ", " <<
      m_sizeZ << ")" << std::endl;

                const int slice      = sizeX * sizeY;
                const int volumeSize = slice * sizeZ;
                m_voxels.resize( volumeSize );

                std::vector<unsigned short> data;
                data.resize( volumeSize );

                // read file into data-vector and close it
    fread( (void*) &(data.front()), sizeof(unsigned short), volumeSize, fp );
                fclose( fp );

                std::cout << "- File loaded." << std::endl;

    //** create density data for histogram
    memset(m_densityArray, 0, sizeof(m_densityArray));
    int densityIndex = 0;

                // process data and store it in m_voxels
    for ( int z=0; z<m_sizeZ; z++ )
                {
                        for ( int y=0; y<m_sizeY; y++ )
                        {
                                for ( int x=0; x<m_sizeX; x++ )
                                {
                                        //we convert the data to float values in an interval of [0..1]
          const float fValue = std::min( 1.0f, float( data[ x + y*m_sizeX + z*slice ] ) / 4095.0f );
          m_voxels[ x + y*m_sizeX + z*slice ] = Voxel( fValue );

          // Creating Histogram
          densityIndex = (int) ( Voxel(fValue).getValue() * DENSITY_COUNT );
          m_densityArray[densityIndex]++; //increment arrayIndexValue

                            if( m_densityArray[densityIndex] > m_maxDensityIndex )
                            {
                                    m_maxDensityIndex = m_densityArray[densityIndex];
                            }
        }
                        }
      std::cout << "\r- Preparing data (" << (z*100) / (m_sizeZ - 1) << "%) ...";
                }
                std::cout << std::endl << "- Data prepared." << std::endl;
        }
}

void Volume::release()
{
  printf("Unload Volume data\n");

  memset(m_densityArray, 0, sizeof(m_densityArray));
  m_voxels.clear();
  m_sizeX = m_sizeY = m_sizeZ = 1;
}

---