Add option to take top <n> voxels in make_masks.py
Fix top voxel identification in maskfactory