mitkPAVessel.h

Go to the documentation of this file.

/*============================================================================

The Medical Imaging Interaction Toolkit (MITK)

Copyright (c) German Cancer Research Center (DKFZ)
All rights reserved.

Use of this source code is governed by a 3-clause BSD license that can be
found in the LICENSE file.

============================================================================*/

#ifndef MITKVESSEL_H
#define MITKVESSEL_H

#include "mitkVector.h"
#include "mitkPAVesselMeanderStrategy.h"
#include "mitkPAInSilicoTissueVolume.h"
#include "mitkPAVector.h"
#include "mitkPAVesselProperties.h"
#include "mitkPAVesselDrawer.h"

#include <MitkPhotoacousticsLibExports.h>

//Includes for smart pointer usage
#include "mitkCommon.h"
#include "itkLightObject.h"

namespace mitk {
  namespace pa {
    class MITKPHOTOACOUSTICSLIB_EXPORT Vessel : public itk::LightObject
    {
    public:
      mitkClassMacroItkParent(Vessel, itk::LightObject);
      mitkNewMacro1Param(Self, VesselProperties::Pointer);

      typedef void (VesselMeanderStrategy::*CalculateNewVesselPositionCallback)
      (Vector::Pointer, double, std::mt19937*);

      void ExpandVessel(mitk::pa::InSilicoTissueVolume::Pointer volume,
        CalculateNewVesselPositionCallback calculateNewPosition, double bendingFactor, std::mt19937* rng);

      bool CanBifurcate();

      Vessel::Pointer Bifurcate(std::mt19937* rng);

      bool IsFinished();

      itkGetConstMacro(VesselProperties, VesselProperties::Pointer);

    protected:
      Vessel(VesselProperties::Pointer parameters);
      ~Vessel() override;

    private:

      const double MINIMUM_VESSEL_RADIUS = 0.1;
      const double NEW_RADIUS_MINIMUM_RELATIVE_SIZE = 0.6;
      const double NEW_RADIUS_MAXIMUM_RELATIVE_SIZE = 0.8;

      VesselMeanderStrategy::Pointer m_VesselMeanderStrategy;
      bool m_Finished;
      double m_WalkedDistance;

      std::uniform_real_distribution<> m_RangeDistribution;
      std::uniform_real_distribution<> m_SignDistribution;
      std::uniform_real_distribution<> m_RadiusRangeDistribution;

      int GetSign(std::mt19937* rng);

      VesselProperties::Pointer m_VesselProperties;

      VesselDrawer::Pointer m_VesselDrawer;
    };

    MITKPHOTOACOUSTICSLIB_EXPORT bool Equal(const Vessel::Pointer leftHandSide, const Vessel::Pointer rightHandSide, double eps, bool verbose);
  }
}
#endif // MITKVESSEL_H