About MCPC3DS

[youyiyiyi]

Do you know what is the possible reason of MCP3DS with an output of all NAN?
Also when running: mask3 = robustmask(romeovoxelquality(phase; mag)); # Using magnitude and phase
, it complains about the loss of key word TEs.
Thanks!

[korbinian90]

For the NAN problem
Can you check if your inputs to mcpc3ds contain NaNs?
for example: sum(.!isfinite.(phase))
phase[.!isfinite.(phase)] .= 0 might help in that case before using mcpc3ds.
If that doesn't fix your problem, can you share your dataset for me to have a look?

romeovoxelquality
When your data is 4D, for example with multiple echoes, then you need to specify the echo time:

TEs = [1,2,3]
mask3 = robustmask(romeovoxelquality(phase; mag, TEs));