function c = cosmo_corr(x, y, corr_type)
cosmo convert neighborhood hdr
cosmo correlation measure hdr