function f_ds = cosmo_fx(ds, f, split_by, dim, check)
cosmo fmri reorient hdr
cosmo independent samples partitioner hdr