Hello MRtrix team,
I’ve been thinking about improved ways to quantify FOD complexity, and it occurred to me that FOD entropy/redundancy, as discussed in https://cds.ismrm.org/ismrm-2005/Files/00384.pdf might be pretty optimal. Any pointers on how to compute these metrics?
There’s no command to do that in MRtrix3, mostly because we never needed it. It was sufficient for what we were doing at the time, but that was quickly swept aside with the introduction of the non-negativity constraint a year later…
That said, I think this one-liner should give you a decent approximation to the entropy. You’ll need to generate a directions files on which to evaluate the FOD first, which you can do using e.g. dirgen 300 dirs300.txt (feel free to use whatever number you feel is appropriate, but you probably want to make sure it’s sufficient to adequately sample the FOD):
Assuming I’m using the correct formula for entropy… There’s probably some other scaling factors in there that I’ve omitted, feel free to amend if you spot any errors…
Probably… It’ll reverse the contrast though, I reckon the image has a more familiar contrast with this (low entropy/high intensity in single fibre regions), but this is entirely your call…
Ha, yes, indeed, we did look at that and found some interesting correlates of complexity. We’ve just been trying to think of alternate parameterizations of “FOD complexity” to double check our intuitions about what we think might be going on there.