cosmo neighborhood split hdr

function [nh_split, masks] = cosmo_neighborhood_split(nh, varargin)