We consider the Bayesian inference of nonlinear, large-scale regression problems in which the parameters model the spatial distribution of some property. A homoscedastic Gaussian sampling distribution is supposed as well as certain assumptions about the regression function. Propriety of the posterior and the existence of its moments are explored when using improper prior distributions expressing different levels of prior knowledge, ranging from a purely noninformative prior over intrinsic Gaussian Markov random field priors to a partition prior. The considered class of problems includes magnetic resonance fingerprinting (MRF). We apply an approximate Bayesian inference to this particular application and demonstrate its practicability in dimensions up to $$10^5$$ or larger. The benefit of incorporating substantial prior knowledge is illustrated. By analyzing simulated realistic MRF data, it is shown that MAP estimates can significantly improve the results achieved with maximum likelihood estimation.