Abraham,
The accuracy of handling analytical multiplicity is likely in the eye of the beholder. In my opinion, three major concerns arise from the voxel-level p-threshold of 0.001, which can result in substantial reproducibility challenges. First, the conventional cluster-based approach implicitly assumes a uniform brain distribution across the brain, a premise that is substantially violated. Second, the resulting penalty is overly harsh and lacks adaptability. Lastly, there is an exclusive emphasis on identifying "false positives," while "false negatives" are totally overlooked.
As a recommended alternative, consider adopting a "highlight, but don’t hide" strategy, as exemplified here.
Gang Chen