Accurate spatial precipitation estimates are important in various fields. An approach to simulate spatial rainfall fields conditioned on radar and rain gauge data is proposed. Unlike the commonly used Kriging methods, which provide a Kriged mean field, the output of the proposed approach is an ensemble of estimates that represents the estimation uncertainty. The approach is robust to nonlinear error in radar estimates and is shown to have some advantages, especially when estimating the extremes.