Function two_shells_voxels() is referred here but I don't see where it is used. You mean read_isbi2013_2shell()?