This note presents a framework for generating tissue probability maps that represent the average shape of a number of subjects' brain images. The procedure is formulated as finding maximum a posteriori estimates within a probabilistic generative model. Estimating the parameters involves alternating between estimating the deformations that match tissue class images of individual subjects to template, and updating the template according to the latest estimates of the deformations. A multinomial matching criterion is used, such that multiple tissue class images (e.g. grey and white matter) are registered simultaneously with the current template estimate. In order to generalise the resulting template to a broader range of subjects, a template blurriness prior is included within the model.