~ubuntu-branches/ubuntu/precise/openwalnut/precise

//---------------------------------------------------------------------------

//

// Project: OpenWalnut ( http://www.openwalnut.org )

//

// Copyright 2009 OpenWalnut Community, BSV@Uni-Leipzig and CNCF@MPI-CBS

// For more information see http://www.openwalnut.org/copying

//

// This file is part of OpenWalnut.

//

// OpenWalnut is free software: you can redistribute it and/or modify

// it under the terms of the GNU Lesser General Public License as published by

// the Free Software Foundation, either version 3 of the License, or

// (at your option) any later version.

//

// OpenWalnut is distributed in the hope that it will be useful,

// but WITHOUT ANY WARRANTY; without even the implied warranty of

// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

// GNU Lesser General Public License for more details.

//

// You should have received a copy of the GNU Lesser General Public License

// along with OpenWalnut. If not, see <http://www.gnu.org/licenses/>.

//

//---------------------------------------------------------------------------

 

#ifndef WINTEGRATIONPARAMETERIZATION_H

#define WINTEGRATIONPARAMETERIZATION_H

 

#include <vector>

 

#include "WRasterParameterization.h"

 

/**

 * Stores the direction if a line in a separate dataset for each voxel.

 */

class WIntegrationParameterization: public WRasterParameterization

{

public:

 

    /**

     * Default constructor.

     *

     * \param grid the grid used for the new dataset.

     */

    explicit WIntegrationParameterization( boost::shared_ptr< WGridRegular3D > grid );

 

    /**

     * Destructor.

     */

    virtual ~WIntegrationParameterization();

 

    /**

     * This method allows this parameterization to update. It gets called for every voxel which is rasterized by the WRasterAlgorithm.

     *

     * \param voxel the voxel to parameterize

     * \param voxelIdx the voxel index in the common grid calculated using "voxel" (this is for convenience)

     * \param axis  Along which axis the traversal takes place. Since when walking in e.g. X-direction there are not supporting voxels in the

     * \param value the new voxel value

     * \param start Start point of the line segment (used to computed the distance)

     * \param end   End point of the line segment (used to computed the distance)

     */

    virtual void parameterizeVoxel( const WValue< int >& voxel, size_t voxelIdx, const int axis, const double value,

                                    const WPosition& start,

                                    const WPosition& end );

 

    /**

     * Gets the dataset representing the parameterization.

     *

     * \return the dataset.

     */

    virtual boost::shared_ptr< WDataSetScalar > getDataSet();

 

    /**

     * Gets called for each new line getting rasterized. Here, it is used to reset the internal length integrator

     *

     * \param line the new line.

     */

    virtual void newLine( const WLine& line );

 

    /**

     * Gets called for each new line segment getting rasterized, as one segment can have multiple voxels.

     *

     * \param start start point of the new line segment

     * \param end end point of the new line segment

     */

    virtual void newSegment( const WPosition& start, const WPosition& end );

 

protected:

 

    /**

     * Stores the current length of the fiber at each voxel.

     */

    std::vector< double > m_lengthValues;

 

    /**

     * The current length of a line.

     */

    double m_curLength;

 

private:

};

 

#endif  // WINTEGRATIONPARAMETERIZATION_H