Implement automatic extraction of BSup pulses from the Bremen sequence:

Here is the formula used for the calulcation:

BSInvTime[0]               =  ((TI + 2.0 * log(0.25 + 0.75 * exp(-0.5 * TI * R1opt)) / R1opt));

BSInvTime[1]               =  ((TI + 2.0 * log(0.75 + 0.25 * exp(-0.5 * TI * R1opt)) / R1opt));

The free parameters can be calculated based the Siemens WIP-parameters as found in the phoenix protocol (starting with “### ASCCONV BEGIN ###” in the Dicom file):

R1opt = sWiPMemBlock.adFree[1] ( 800ms in the data you sent)
He noted that R1opt is actually not BSparam1, but actually 1/BSparam1.

TI = TI as provided in sequence (will differ for multi-TI data) - sWiPMemBlock.adFree[2] (-10ms in the data you sent)

You can get the start-TI and increment from the WIP parameter for multi-TI/PLD as well as some settings for the FOCI pulses (which might allow for inversion efficiency estimation). However, these might vary with different...