Curl-constrained Gradient Estimation for Image Recovery from Highly Incomplete Spectral Data