cosmo correlation measure hdr

function ds_sa = cosmo_correlation_measure(ds, varargin)