OpenWalnut  1.4.0
WFiber.h
00001 //---------------------------------------------------------------------------
00002 //
00003 // Project: OpenWalnut ( http://www.openwalnut.org )
00004 //
00005 // Copyright 2009 OpenWalnut Community, BSV@Uni-Leipzig and CNCF@MPI-CBS
00006 // For more information see http://www.openwalnut.org/copying
00007 //
00008 // This file is part of OpenWalnut.
00009 //
00010 // OpenWalnut is free software: you can redistribute it and/or modify
00011 // it under the terms of the GNU Lesser General Public License as published by
00012 // the Free Software Foundation, either version 3 of the License, or
00013 // (at your option) any later version.
00014 //
00015 // OpenWalnut is distributed in the hope that it will be useful,
00016 // but WITHOUT ANY WARRANTY; without even the implied warranty of
00017 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
00018 // GNU Lesser General Public License for more details.
00019 //
00020 // You should have received a copy of the GNU Lesser General Public License
00021 // along with OpenWalnut. If not, see <http://www.gnu.org/licenses/>.
00022 //
00023 //---------------------------------------------------------------------------
00024 
00025 #ifndef WFIBER_H
00026 #define WFIBER_H
00027 
00028 #include <vector>
00029 
00030 #include "../math/WLine.h"
00031 #include "../math/linearAlgebra/WPosition.h"
00032 
00033 
00034 class WFiberTest;
00035 
00036 /**
00037  * Represents a neural pathway.
00038  */
00039 class WFiber : public WLine
00040 {
00041 friend class WFiberTest;
00042 public:
00043     /**
00044      * Constructs a new fiber out of WPositions.
00045      *
00046      * \param points Reference to the points which belong to this fiber
00047      */
00048     explicit WFiber( const std::vector< WPosition > &points );
00049 
00050     /**
00051      * Creates an empty fiber.
00052      */
00053     WFiber();
00054 
00055     /**
00056      * This is the Smaller thresholded distance as described by Zhang: http://dx.doi.org/10.1109/TVCG.2008.52 .
00057      *
00058      * \param thresholdSquare Threshold upto which the distances should be ignored given as square for reasons of performance.
00059      * \param q First fiber
00060      * \param r Second fiber
00061      *
00062      * \return The minimum of dt(Q,R) and dt(R,Q)
00063      */
00064     static double distDST( double thresholdSquare, const WFiber &q, const WFiber &r );
00065 
00066     /**
00067      * This is the Larger thresholded distance as described by Zhang: http://dx.doi.org/10.1109/TVCG.2008.52 .
00068      *
00069      * \param thresholdSquare Threshold upto which the distances should be ignored given as square for reasons of performance.
00070      * \param q First fiber
00071      * \param r Second fiber
00072      *
00073      * \return The maximum of dt(Q,R) and dt(R,Q)
00074      */
00075     static double distDLT( double thresholdSquare, const WFiber &q, const WFiber &r );
00076 
00077 protected:
00078 private:
00079 };
00080 #endif  // WFIBER_H