Efficient Multivariate Skellam Shrinkage For Denoising Photon-limited Image Data: An Empirical Bayes Approach