Given spatially located observed random variables (x, z = {(xi, zi)}i, we propose a new method for non-parametric estimation of the potential functions of a Markov random field p(x|z), based on a roughness penalty approach. The new estimator maximizes the penalized log-pseudolikelihood function...