function msk = cosmo_match(haystack, needle, varargin)
cosmo mask dim intersect hdr
cosmo measure clusters hdr