Robust and scalable Bayesian analysis of spatial neural tuning function data