Looks like @anege came across the solution to this issue, but the text got piped to the wrong stream deep in the depths of the Discourse server where nobody can see it
With the fixel directory format, because each fixel data file is just a 1D image file, many “workhorse” MRtrix3 commands are applicable, including mrthrehsold
. So applying a threshold of e.g. 0.5, and then feeding the result as the -mask
option of mrstats
, should provide the number of super-threshold fixels in the “count” output.
The fixelcrop
option suggested by @rgrazioplene should give the same answer, but is doing a bit more work to get there: It creates a new fixel directory, where only those fixels in the input directory that have a value of 1 in the mask (e.g. super-threshold as determined by mrthreshold
) are included in the output fixel directory; those fixels that have a value of 0 in the mask are omitted entirely. So the number of fixels present in that directory will correspond to the number of fixels that had a value of 1 in your mask; but you also have a new fixel directory containing only those fixels, which you may or may not want depending on context.