function [confusion_matrix, classes] = cosmo_confusion_matrix(ds, varargin)
cosmo config hdr
cosmo convert neighborhood hdr