Which dwi2response algorithm will work best

dieloz · October 27, 2021, 6:15am

Hi Philip,

The default behavior of dwi2fod csd when you run dwi2response tournier on multishell data is indeed automatically selecting largest b-value

github.com

MRtrix3/mrtrix3/blob/430465d64ad7cedb71ece60cb81428e415190029/core/dwi/shells.cpp#L250


      
              throw Exception ("User did not select b=0 shell, but the command requires the presence of b=0 data");
          
            if (force_without_bzero && bzero_selected)
              throw Exception ("User selected b=0 shell, but the command is not compatible with b=0 data");
          
          } else {
          
            if (force_singleshell && !is_single_shell()) {
              if (count() == 1 && has_bzero())
                throw Exception ("No non b=0 data found, but the command requires a non b=0 shell");
              WARN ("Multiple non-zero b-value shells detected, automatically selecting largest b-value: b=" + str(largest().get_mean()) + " with " + str(largest().count()) + " volumes");
              to_retain[count()-1] = true;
              if (has_bzero())
                to_retain[0] = true;
            } else {
              // default: keep everything
              to_retain.clear (true);
            }
          
            if (force_with_bzero && !has_bzero())
              throw Exception ("No b=0 data found, but the command requires the presence of b=0 data");

I am also curious to know what could be the rationalle behind the focusing on the largest b-value since there are also views against this approach vs the msmt_csd.
Best,
Diego