“…These trajectories were defined on a time-domain from inner to outer surfaces with normality constraints enforced at each end. Equivolumetric depth was obtained through reparameterization of intermediate deformed surfaces to maintain constant volume along trajectories [30] . A corresponding depth image with range from 0 at the pia surface to 1 at gray-white matter boundary was then created by interpolating registration trajectories on a regular grid.…”