We propose a method for performing simulation-based inference (SBI) of function-valued parameters, and apply it to spatial inference problems in the geosciences.