gradient_pulses.h

Go to the documentation of this file.

// Copyright 2025 Julien Lamy, ICube, Université de Strasbourg-CNRS.
// Part of msl, distributed under the terms of the MIT license.
 
#ifndef _d6067f26_790d_47d2_9843_9ffef6843fbb
#define _d6067f26_790d_47d2_9843_9ffef6843fbb
 
#include <algorithm>
 
#include <MrImagingFW/libMath/MrException.h>
#include <MrImagingFW/libSeqSysProp/SysProperties.h>
 
#include <MrMeasSrv/SeqIF/libRT/libRTDefines.h>
 
#include "helpers.h"
#include "GradientPulse.h"
#include "GradientSpecs.h"
#include "Vector.h"
 
namespace msl
{

namespace gradient_pulses
{

GradientPulse quickest(
    Vector3d const & k, double maxAmplitude, double minRiseTime);

GradientPulse quickest(Vector3d const & k, GradientSpecs const & specs);
    
GradientPulse plateau(
    Vector3d const & kPlateau, long plateauDuration,
    double maxAmplitude, double minRiseTime);

GradientPulse plateau(
    Vector3d const & kPlateau, long plateauDuration,
    GradientSpecs const & specs);

GradientPulse fixedDuration(
    Vector3d const & k, long rampDuration, long plateauDuration,
    double maxAmplitude, double minRiseTime);

GradientPulse fixedDuration(
    Vector3d const & k, GradientPulse const & reference,
    double maxAmplitude, double minRiseTime);

GradientPulse fixedDuration(
    Vector3d const & k, long rampDuration, long plateauDuration,
    GradientSpecs const & specs);

GradientPulse fixedDuration(
    Vector3d const & k, GradientPulse const & reference,
    GradientSpecs const & specs);

template<typename T>
GradientPulse sliceSelection(T && pulse, double maxAmplitude, double minRiseTime)
{
    auto const amplitude = 
        pulse.getBandwidth() / (pulse.getLarmorConst() * pulse.getThickness());
    if(amplitude > maxAmplitude)
    {
        throw MrException("Maximum amplitude exceeded");
    }
    
    auto rampDuration = std::lround(roundUp(
        std::max(
            minRiseTime * amplitude, double(SysProperties::getCoilCtrlLead())),
        GRAD_RASTER_TIME));
    auto plateauDuration = std::lround(
        roundUp(pulse.getDuration(), GRAD_RASTER_TIME));
    return {{0., 0., amplitude}, rampDuration, plateauDuration};
}

template<typename T>
GradientPulse sliceSelection(T && pulse, GradientSpecs const & specs)
{
    return sliceSelection(pulse, specs.maxAmplitude, specs.minRiseTime);
}
 
}
 
}
 
#endif // _d6067f26_790d_47d2_9843_9ffef6843fbb