I'm not saying the erosion method is best, but it is pretty good. It is definitely faster than the BB method. The ball and box method measures thickness in a very different way, so nodules or bumps would have different measures. Both are somewhat sensitive to the sampling resolution. The BB method will work about twice as fast if "boxes" aren't included with "-balls_only". For macaque resolutions at 0.5 or 0.25 mm resolution data, upsampling isn't really required, but the thickness precision is more limited with lower thickness values. The best way is to try and see. Thickness precision is limited by voxel size for all methods really.
We are working on a distance transform method that is just about ready for release too. The distance transform part was started in an AFNI hackathon (Code Convergence) with Chris Rorden. He wrote up a very nice version as Depth3D, that is both fast and accurate. The method computes the shortest distance to an edge of the mask at every voxel in a very efficient way. Transforming depth to thickness can be done as in the other methods by projecting 2*maximum depth either volumetrically or by volume to surface projection along normal vectors. Then smoothing the data in the volume or surface gives the final map.
[
github.com]
[
afni.nimh.nih.gov]