A method for segmenting an object of interest from an image of a patient having such object. Each one of a plurality of training shapes is distorted to overlay a reference shape with a parameter Θi being a measure of the amount of distortion required to effect the overlay. A vector of the parameters Θi is obtained for every one of the training shapes through the minimization of a cost function along with an estimate of uncertainty for every one of the obtained vectors of parameters Θi, such uncertainty being quantified as a covariance matrix Σi. A statistical model represented as {circumflex over (f)}H (Θ,Σ) is generated with the sum of kernels having a mean Θi and covariance Σi. The desired object of interest in the image of the patient is identified by positioning of the reference shape on the image and distorting the reference shape to overlay the obtained image with a parameter Θ being a measure of the amount of distortion required to effect the overlay. An uncertainty is quantified as a covariance matrix Σ and an energy function E=Eshape+Eimage is computed to obtain the probability of the current shape in the statistical shape model Eshape(Θ,Σ)=−log({circumflex over (f)}H) and the fit in the image Eimage.