cosmo interval neighborhood hdr

function nbrhood = cosmo_interval_neighborhood(ds, label, varargin)