Incompatibility between default settings and NODDI

[AlbertoImg]

Hi there,

I would like to report an issue: Running CLI UKF (from pnlNipipeline) with the noddi option we get a terminate called without an active exception with the message:

recordFA, recordTrace, record_free_water, recordTensors parameters can only be used with tensor models

It seems that recordTensor from the default settings is not compatible with NODDI.

• DEFAULT ukf.py

ukfdefaults = ['--numTensor', 2, '--stoppingFA', 0.15, '--seedingThreshold', 0.18, '--Qm', 0.001, '--Ql', 70,
               '--Rs', 0.015, '--stepLength', 0.3, '--recordLength', 1.7, '--stoppingThreshold', 0.1,
               '--seedsPerVoxel', 10, '--recordTensors']

• CODE

ukf -i $dwi \
    --bvals $dwi_bval \
    --bvecs $dwi_bvec \
    -m $dwi_mask \
    -o $dir_dwi/$ukf_out/tracts.vtk \
  --params --noddi,--recordKappa,--recordVic,--recordViso,--Qkappa,0.01,--Qvic,0.004

• OUTPUT

Using the 2T simple model. Setting the default parameters accordingly:
"*": set by user
"-": default setting
- stoppingFA: 0.15
* seedingThreshold: 0.18
* Qm: 0.001
* Qkappa: 0.01
* Rs: 0.015
* stepLength: 0.3
* recordLength: 1.7
* Qvic = Qviso: 0.004
* stoppingThreshold: 0.1
* seedsPerVoxel: 10
Found 10 cores on your system.
Running tractography with 10 thread(s).
DWMRI_b-value 3010
Number of non-zero gradients: 202
Number of zero gradients: 14
Permuting the axis order to: 3 0 1 2
Resizing the data to: 202 140 140 92
Computing the baseline image
Dividing the signal by baseline image
Data normalization finished!

Using NODDI 2-Fiber model.
recordFA, recordTrace, record_free_water, recordTensors parameters can only be used with tensor models
terminate called without an active exception

• TEST

When I delete --recordTensors from the defaults, we are getting the following and it is still running without crushing:

Using NODDI 2-Fiber model.
Branching disabled
Using constrained filter
Seeding using Option 3, from a seeding label map.