NSM constrained approximation of Lipschitz functions from data