We propose an approach for boundary finding where the correspondence of a subset of boundary points to a model is simultaneously determined. Global shape parameters derived from the statistical variation of object boundary points in a training set are used to model the abject. A Bayesian formulation, based on this prior knowledge and the edge information of the input image, is employed to find the object boundary with its subset paints in correspondence with boundaries in the training set or the mean boundary. We compared the use of a generic smoothness prior and a uniform independent prior with the training set prior in order to demonstrate the power of this statistical information. A number of experiments were performed on both synthetic and real medical images of the brain, and heart to evaluate the approach, including the validation of the dependence of the method on image quality, different initialization and prior information.