[Camino-users] modelfit 'WARNING: Can't use normalization constant NaN'

[Marie Amelie]

Dear camino users,

I am using camino to fit the tensor model (nonlinear optimization,
constrained to be positive semi-definite) on dMRI data with the following
command:

modelfit -inputfile dwi.Bfloat -schemefile schemefile.scheme -inversion 2
-bgmask mask_Bfloat.nii.gz -maskdatatype float -outputfile dt.Bdouble

The process globally seems to work fine: the tensor fitting and then the
derived scalar maps (produced as indicated in camino's DTI tutorial) are ok
both at the visual inspection and by comparing their values to the ones
computed by a previous tensor fitting with another tool. However in the
terminal I am seeing many instances of these warnings:

Jan 17, 2018 10:50:35 AM imaging.DW_Scheme normalizeData
WARNING: Can't use normalization constant NaN

My question is: can you explain me why I am seeing these warnings? Can I
trust my tensor fit despite these warnings? Again, scalar maps and their
values look globally fine.
Could it be related to my brain/backgound mask? It is actually rather
conservative (i.e. tends to exclude some voxels of the brain rather than
include background voxels), but I can't exclude that in some instances it
might also include some background (even if perhaps not that frequently).
I am a new camino user, I apologize in advance if this is a basic question.

----------------------------------------------------------------------------
Here follows additional information.

Before the tensor fit, I performed an eddy-current correction (MRtrix'
dwipreproc) and then did with camino:
fsl2scheme -bvecfile bvec_rot.txt -bvalfile bval.txt > schemefile.scheme
image2voxel -4dimage dwi.nii.gz -outputfile dwi.Bfloat


For completeness, I report here that I also get some instances of this
error:
Jan 17, 2018 10:50:32 AM misc.LoggedException logException
WARNING: Exception in thread "main" class
optimizers.MarquardtMinimiserException:
Singular Matrix in gaussj - 2
    at optimizers.MarquardtMinimiser.gaussj(MarquardtMinimiser.java:362)
    at optimizers.MarquardtMinimiser.mrqmin(MarquardtMinimiser.java:231)
    at optimizers.MarquardtMinimiser.minimise(MarquardtMinimiser.java:156)
    at inverters.NonLinearDT_Inversion.invert(NonLinearDT_
Inversion.java:135)
    at apps.ModelFit.execute(ModelFit.java:156)
    at apps.EntryPoint.main(EntryPoint.java:290)

but from the camino FAQ webpage I think I understood their meaning.

I am running it on a Ubuntu 16.04 machine.
The git log of my camino version is:
commit af55b0acd7895b7d1dafa040df60c494369b138a
Author: Philip Cook <cookpa at mail.med.upenn.edu>
Date:   Fri Oct 6 12:00:17 2017 -0400

    ENH: Allow non-zero unweighted measurements for dtspd

---------------------------------------------------------------------------
Any help on my question is highly appreciated.
Thank you in advance.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.nitrc.org/pipermail/camino-users/attachments/20180118/f49e3ba4/attachment.html>