Estimating high-dimensional demand systems in the presence of many binding non-negativity constraints