BAYESIAN IMAGING USING GOOD'S ROUGHNESS MEASURE - IMPLEMENTATION ON A MASSIVELY PARALLEL PROCESSOR.

A constrained maximum-likelihood estimator is derived by incorporating a rotationally invariant roughness penalty proposed by I. J. Good (1981) into the likelihood functional. This leads to a set of nonlinear differential equations the solution of which is a spline-smoothing of the data. The nonlinear partial differential equations are mapped onto a grid via finite differences, and it is shown that the resulting computations possess a high degree of parallelism as well as locality in the data-passage, which allows an efficient implementation on a 48-by-48 mesh-connected array of NCR GAPP processors. The smooth reconstruction of the intensity functions of Poisson point processes is demonstrated in two dimensions.