cosmo fmri orientation hdr

function [orient, cano_rot] = cosmo_fmri_orientation(ds)