cosmo spherical neighborhood hdr

function nbrhood = cosmo_spherical_neighborhood(ds, varargin)