//Based on code from Sproxel
bool loadCSV(const QString& filename, Thermite::Volume* volume, qint32 posX, qint32 posY, qint32 posZ)
{
   PolyVox::SimpleVolume<PolyVox::Material16>* pPolyVoxVolume = volume->m_pPolyVoxVolume;

   int fscanfStatus = 0;
   FILE* fp = fopen(filename.toAscii().constData(), "rb");
   if (!fp) return false;

   // Read the dimensions
   int sizeX = 0, sizeY = 0, sizeZ = 0;
   fscanfStatus = fscanf(fp, "%d,%d,%d\n", &sizeX, &sizeY, &sizeZ);

   // Read the data
   QColor color;
   for (int y = sizeY-1; y >= 0; y--)
   {
      for (int z = 0; z < sizeZ; z++)
      {
         for (int x = 0; x < sizeX; x++)
         {
            int ir, ig, ib, ia;
            fscanfStatus = fscanf(fp, "#%02X%02X%02X%02X,", &ir, &ig, &ib, &ia);

            /*float r = ir / (float)0xff;
            float g = ig / (float)0xff;
            float b = ib / (float)0xff;
            float a = ia / (float)0xff;*/

            QColor color(ir, ig, ib, ia);

            if(ia > 0)
            {
               PolyVox::Material16 voxel;
               voxel.setMaterial(QColorToMaterial(color));

               int xToWrite = posX + x;
               int yToWrite = posY + y;
               int zToWrite = posZ + z;
               if((xToWrite >= 0) && (xToWrite < pPolyVoxVolume->getWidth()) &&
                  (yToWrite >= 0) && (yToWrite < pPolyVoxVolume->getHeight()) &&
                  (zToWrite >= 0) && (zToWrite < pPolyVoxVolume->getDepth()))
               {
                  pPolyVoxVolume->setVoxelAt(xToWrite, yToWrite, zToWrite, voxel);
               }
            }

            if (x != sizeZ-1)
               fscanfStatus = fscanf(fp, ",");
         }
         fscanfStatus = fscanf(fp, "\n");
      }
      fscanfStatus = fscanf(fp, "\n");
   }

   fclose(fp);
   return true;
}

   void DrawLine(float x1, float y1, float z1, float x2, float y2, float z2)
   {
      int i = (int)floorf(x1);
      int j = (int)floorf(y1);
      int k = (int)floorf(z1);

      int iend = (int)floorf(x2);
      int jend = (int)floorf(y2);
      int kend = (int)floorf(z2);

      int di = ((x1 < x2) ? 1 : ((x1 > x2) ? -1 : 0));
      int dj = ((y1 < y2) ? 1 : ((y1 > y2) ? -1 : 0));
      int dk = ((z1 < z2) ? 1 : ((z1 > z2) ? -1 : 0));

      float minx = floorf(x1), maxx = minx + 1.0f;
      float tx = ((x1 > x2) ? (x1 - minx) : (maxx - x1)) / abs(x2 - x1);
      float miny = floorf(y1), maxy = miny + 1.0f;
      float ty = ((y1 > y2) ? (y1 - miny) : (maxy - y1)) / abs(y2 - y1);
      float minz = floorf(z1), maxz = minz + 1.0f;
      float tz = ((z1 > z2) ? (z1 - minz) : (maxz - z1)) / abs(z2 - z1);

      float deltatx = 1.0f / abs(x2 - x1);
      float deltaty = 1.0f / abs(y2 - y1);
      float deltatz = 1.0f / abs(z2 - z1);

      SetPos(i,j,k);

      for(;;)
      {
         if(tx <= ty && tx <= tz)
         {
            if(i == iend) break;
            tx += deltatx;
            i += di;

            if(di == 1) pos.setX(pos.getX()+1);
            if(di == -1) pos.setX(pos.getX()-1);
            SetVoxel();
         } else if (ty <= tz)
         {
            if(j == jend) break;
            ty += deltaty;
            j += dj;

            if(dj == 1) pos.setY(pos.getY()+1);
            if(dj == -1) pos.setY(pos.getY()-1);
            SetVoxel();
         } else
         {
            if(k == kend) break;
            tz += deltatz;
            k += dk;

            if(dk == 1) pos.setZ(pos.getZ()+1);
            if(dk == -1) pos.setZ(pos.getZ()-1);
            SetVoxel();
         }
      }
   } 

 SetPos(i,j,k);

pos.setX(pos.getX()+1);