PolyVox::CubicSurfaceExtractor<PolyVox::MaterialDensityPair44> surfaceExtractor(&volData, volData.getEnclosingRegion(), &mesh);
    surfaceExtractor.execute();

    if ( PVWorld != NULL ) {
        PVWorldNode->detachObject(PVWorld);
        scenemgr->destroyManualObject(PVWorld);
    }
    PVWorld = scenemgr->createManualObject("PVWorld");

    uint32_t noVertices = mesh.getNoOfVertices();
    uint32_t noIndices = mesh.getNoOfIndices();

    PVWorld->estimateVertexCount(noVertices);
    PVWorld->estimateIndexCount(noIndices);

    PVWorld->begin("ColouredCubicVoxel", Ogre::RenderOperation::OT_TRIANGLE_LIST);

    // vertices
    const std::vector<PolyVox::SurfaceVertex>& vVertices = mesh.getVertices();
    for (unsigned int i=0; i<noVertices; i++) {
        const PolyVox::Vector3DFloat& pos = vVertices[i].getPosition();
        PVWorld->position(pos.getX(), pos.getY(), pos.getZ());
        const PolyVox::Vector3DFloat& normal = vVertices[i].getNormal();
        PVWorld->normal(normal.getX(), normal.getY(), normal.getZ());
    }

    // indices
    const std::vector<uint32_t>& vIndices = mesh.getIndices();
    for (unsigned int i=0; i<noIndices; i++) {
        PVWorld->index( vIndices[i] );
    }

    PVWorld->end();

    if ( PVWorldNode == NULL ) {
        PVWorldNode = scenemgr->getRootSceneNode()->createChildSceneNode();
    }
    PVWorldNode->attachObject(PVWorld);
    PVWorldNode->setScale(100.0, 100.0, 100.0);
    PVWorldNode->setPosition(0.0, 0.0, 0.0);

        :m_volData(volData)
        ,m_sampVolume(volData)
        ,m_regSizeInVoxels(region)
        ,m_meshCurrent(result)

        :m_volData(volData)
        ,m_sampVolume(volData)
        ,m_meshCurrent(result)
        ,m_regSizeInVoxels(region)

    const std::vector<PolyVox::SurfaceVertex>& vVertices = mesh.getVertices();
    for (unsigned int i=0; i<noVertices; i++) {
        // position
        const PolyVox::Vector3DFloat& pos = vVertices[i].getPosition();
        PVWorld->position(pos.getX(), pos.getY(), pos.getZ());
        // normal
        const PolyVox::Vector3DFloat& normal = vVertices[i].getNormal();
        PVWorld->normal(normal.getX(), normal.getY(), normal.getZ());
        // colour/material
        Ogre::ColourValue val;
        uint8_t material = vVertices[i].getMaterial() + 0.5;
        val.r = 0.453f;
        val.g = 0.0f;
        val.b = 0.0f;
        val.a = 1.0f;
        PVWorld->colour(val);
    }

#define VOXELSCALE 100.0
#define VOXELCENTEROFFSET 0.5
// pixel size of tex atlas
#define TEXATLASSIZE 256
// pixel per single texture
#define TEXATLAS_TEXSIZE 16
// number of textures per row
#define TEXATLAS_NUMTEX (TEXATLASSIZE / TEXATLAS_TEXSIZE)
// size of a single texture in 0.0 - 1.0 values
#define TEXATLAS_TEXSIZE_NORMALIZED ((float)TEXATLAS_TEXSIZE / (float)TEXATLASSIZE)

void ColouredCubicVoxelVP(
    float4 inPosition    : POSITION,
    float4 inNormal        : NORMAL,
    float4 inColor        : COLOR,

    out float4 outClipPosition    : POSITION,
    out float4 outWorldPosition    : TEXCOORD0,
    out float4 outWorldNormal    : TEXCOORD1,
    out float4 outColor            : TEXCOORD2,

    uniform float4x4 world,
    uniform float4x4 viewProj
    )
{
    //Compute the world space position
    outWorldPosition = mul(world, inPosition);

    //Just pass though the normals without transforming them in any way. No rotation occurs.
    outWorldNormal = inNormal;

    //Compute the clip space position
    outClipPosition = mul(viewProj, outWorldPosition);

    outColor = inColor;
}

void ColouredCubicVoxelFP(
    float4 inPosition        : POSITION,
    float4 inWorldPosition    : TEXCOORD0,
    float4 inWorldNormal    : TEXCOORD1,
    float4 inColor            : TEXCOORD2,

    uniform sampler2D texAtlas    : TEXUNIT0,

    out float4 result        : COLOR)
{
    inWorldNormal = normalize(inWorldNormal);

    float3 col;

    //World position is used as texture coordinates. Choose which
    //two components of world position to use based on normal. Could
    //optionally use a different texture for each face here as well.
    float2 pos;
    if(inWorldNormal. x > 0.5)
    {
        col = tex2D(texAtlas, (inWorldPosition.yz / VOXELSCALE) + VOXELCENTEROFFSET);
    }

    if(inWorldNormal. x < -0.5)
    {
        col = tex2D(texAtlas, (inWorldPosition.yz / VOXELSCALE) + VOXELCENTEROFFSET);
    }

    // top
    if(inWorldNormal. y > 0.5)
    {
        pos = ((inWorldPosition.xz / VOXELSCALE) + VOXELCENTEROFFSET) / TEXATLAS_NUMTEX;
        if ( inColor[0] == 0.453 ) {
            col=0;
        }
        else {
            col = tex2D(texAtlas, pos);
        }
    }

    // bottom
    if(inWorldNormal. y < -0.5)
    {
        col = tex2D(texAtlas, (inWorldPosition.xz / VOXELSCALE) + VOXELCENTEROFFSET);
    }

    if(inWorldNormal. z > 0.5)
    {
        col = tex2D(texAtlas, (inWorldPosition.xy / VOXELSCALE) + VOXELCENTEROFFSET);
    }

    if(inWorldNormal. z < -0.5)
    {
        col = tex2D(texAtlas, (inWorldPosition.xy / VOXELSCALE) + VOXELCENTEROFFSET);
    }

    result = float4(col, 1.0);
}