95 lines
3.1 KiB
C++
95 lines
3.1 KiB
C++
/*******************************************************************************
|
|
Copyright (c) 2005-2009 David Williams
|
|
|
|
This software is provided 'as-is', without any express or implied
|
|
warranty. In no event will the authors be held liable for any damages
|
|
arising from the use of this software.
|
|
|
|
Permission is granted to anyone to use this software for any purpose,
|
|
including commercial applications, and to alter it and redistribute it
|
|
freely, subject to the following restrictions:
|
|
|
|
1. The origin of this software must not be misrepresented; you must not
|
|
claim that you wrote the original software. If you use this software
|
|
in a product, an acknowledgment in the product documentation would be
|
|
appreciated but is not required.
|
|
|
|
2. Altered source versions must be plainly marked as such, and must not be
|
|
misrepresented as being the original software.
|
|
|
|
3. This notice may not be removed or altered from any source
|
|
distribution.
|
|
*******************************************************************************/
|
|
|
|
#include "OpenGLWidget.h"
|
|
|
|
#include "MaterialDensityPair.h"
|
|
#include "CubicSurfaceExtractorWithNormals.h"
|
|
#include "SurfaceMesh.h"
|
|
#include "Volume.h"
|
|
|
|
#include <QApplication>
|
|
|
|
//Use the PolyVox namespace
|
|
using namespace PolyVox;
|
|
|
|
void createSphereInVolume(Volume<MaterialDensityPair44>& volData, float fRadius)
|
|
{
|
|
//This vector hold the position of the center of the volume
|
|
Vector3DFloat v3dVolCenter(volData.getWidth() / 2, volData.getHeight() / 2, volData.getDepth() / 2);
|
|
|
|
//This three-level for loop iterates over every voxel in the volume
|
|
for (int z = 0; z < volData.getWidth(); z++)
|
|
{
|
|
for (int y = 0; y < volData.getHeight(); y++)
|
|
{
|
|
for (int x = 0; x < volData.getDepth(); x++)
|
|
{
|
|
//Store our current position as a vector...
|
|
Vector3DFloat v3dCurrentPos(x,y,z);
|
|
//And compute how far the current position is from the center of the volume
|
|
float fDistToCenter = (v3dCurrentPos - v3dVolCenter).length();
|
|
|
|
//If the current voxel is less than 'radius' units from the center then we make it solid.
|
|
if(fDistToCenter <= fRadius)
|
|
{
|
|
//Our new density value
|
|
uint8_t uDensity = MaterialDensityPair44::getMaxDensity();
|
|
|
|
//Get the old voxel
|
|
MaterialDensityPair44 voxel = volData.getVoxelAt(x,y,z);
|
|
|
|
//Modify the density
|
|
voxel.setDensity(uDensity);
|
|
|
|
//Wrte the voxel value into the volume
|
|
volData.setVoxelAt(x, y, z, voxel);
|
|
}
|
|
}
|
|
}
|
|
}
|
|
}
|
|
|
|
int main(int argc, char *argv[])
|
|
{
|
|
//Create and show the Qt OpenGL window
|
|
QApplication app(argc, argv);
|
|
OpenGLWidget openGLWidget(0);
|
|
openGLWidget.show();
|
|
|
|
//Create an empty volume and then place a sphere in it
|
|
Volume<MaterialDensityPair44> volData(64, 64, 64);
|
|
createSphereInVolume(volData, 30);
|
|
|
|
//Extract the surface
|
|
SurfaceMesh<PositionMaterialNormal> mesh;
|
|
CubicSurfaceExtractorWithNormals<MaterialDensityPair44> surfaceExtractor(&volData, volData.getEnclosingRegion(), &mesh);
|
|
surfaceExtractor.execute();
|
|
|
|
//Pass the surface to the OpenGL window
|
|
openGLWidget.setSurfaceMeshToRender(mesh);
|
|
|
|
//Run the message pump.
|
|
return app.exec();
|
|
}
|