animaMCMConstants.h

Go to the documentation of this file.

#pragma once
#include <cmath>

namespace anima
{

// The following units are in place to keep things right: distances are in millimeters, time in seconds, magnetic field in Tesla

const double DiffusionGyromagneticRatio = 267.513e6;

inline double GetBValueFromAcquisitionParameters(double smallDelta, double bigDelta, double gradientStrength)
{
    double alpha = DiffusionGyromagneticRatio * smallDelta * gradientStrength;
    return alpha * alpha * (bigDelta - smallDelta / 3.0);
}

inline double GetGradientStrengthFromBValue(double bValue, double smallDelta, double bigDelta)
{
    double alpha = std::sqrt(bValue / (bigDelta - smallDelta / 3.0));
    return alpha / (DiffusionGyromagneticRatio * smallDelta);
}

const double DiffusionSmallDelta = 10.0e-3;

const double DiffusionBigDelta = 40.0e-3;

const double MCMZeroLowerBound = 0.0;

const double MCMEpsilon = 1.0e-2;

const double MCMCompartmentsFractionUpperBound = 1.0e7;

const double MCMAxialDiffusivityAddonLowerBound = 5.0e-4;

const double MCMDiffusivityLowerBound = 1.0e-6;

const double MCMDiffusivityUpperBound = 3.0e-3;

const double MCMRadialDiffusivityUpperBound = 1.0e-3;

const double MCMFreeWaterDiffusivityLowerBound = 2.0e-3;

const double MCMFreeWaterDiffusivityUpperBound = 4.0e-3;

const double MCMIsotropicRestrictedWaterDiffusivityLowerBound = 0.25e-4;

const double MCMIsotropicRestrictedWaterDiffusivityUpperBound = 2.0e-3;

const double MCMPolarAngleUpperBound = M_PI;

const double MCMAzimuthAngleUpperBound = 2.0 * M_PI;

const double MCMConcentrationUpperBound = 128.0;

const double MCMFractionUpperBound = 1.0;

const double MCMTissueRadiusLowerBound = 1.0e-4;

const double MCMTissueRadiusUpperBound = 4.01e-2;

}