update07/08/2014

add the functions using stronger constraints
the output is saved in .nii format