Intensity correction was implemented using an in-house procedure that performs an edge-aware local normalization by the mean of the signal within each image for each subject.  The heuristic signal likelihood at each voxel and the global mean of the signal