afni.nimh.nih.gov

3dAllineate does have mask for the source dataset, but I'm not sure if that's what you want. Are you trying to fill holes? Is this for ROIs with a particular integer value for each region type, or are you trying to fill holes of function with floating point values.