Yes, you will need to install it. Nipype provides a standardised way to call functions from various packages, but the packages still need to be present in your local environment as it won’t download them for you. It looks like mindboggle distributes a Docker container, which might be of use !
Thank you for your interest in Mindlogger! We have disabled thickinthehead in the automated shape analysis pipeline because we found a possible dependency on isotropic resolution. While it is available as a command-line argument and in Mindlogger’s python library, we recommend using this shape measure only in special circumstances.
@binarybottle: Thank you!
I have noticed that from here: https://github.com/nipy/mindboggle/issues/149
So by special circumstances, do you mean when the images are isotropic 1 mm^3?
I have my T1 images in isotropic 1 mm^3 space, so can I assume that the revised algorithm will work similarly to what you have published?
Thank you, @aghayoor. This is precisely why I caution against use of thickinthehead for general-purpose brain image processing. It assumes not only 1mm^3 isometric voxels but also a minimum cortical thickness. In all the cases I have run, the cortex was several voxels thick, so I didn’t run into this problem, but if you are running these on brain images with very thin cortices, thickinthehead could fail.
Re: line 379, only “label_volume_thickness” and “output_table” are returned by the function, not “thickness”. if the “if statement” is not passed, then -1 values are saved in label_volume_thickness and nothing is written to output_table for that label.
@binarybottle: Thank you for your review. Do you think can you offer a more generalized version of the algorithm? Do you have any suggestions to fix the above bug? I will be glad to help making the code more reliable. Thank you.
I tried a quick fix by upsampling the resolution to simulate thicknesses of multiple voxels, but this required too much memory, so I removed that function. I haven’t had time to revisit this problem, but welcome any pull requests that tackle this issue!
@binarybottle: Is it possible if you put back that function and make its usage optional by a flag?
Using this option can be recommended only in cases that the input subject scan might include very thin cortical thicknesses in some cortical regions. For example in processing of the brain scans in some neurodegenerative disease studies like late AD.
It’s easy enough to reset the upsample factor, but this raises a new issue – the aliasing will be so great on curved regions that the estimate of volume/area will be increasingly inaccurate. You are welcome to play with this function with greater upsampling but I don’t want to alter the code if it encourages users to run into memory and accuracy issues…