function res = cosmo_distatis(ds, varargin)
cosmo dissimilarity matrix measure hdr
cosmo find local extrema hdr