The present invention provides a method for obtaining an estimation of the shape, position and/or orientation of one or more existing teeth of a patient or of one or more teeth to be included in a dental restoration destined to replace one or more missing teeth in a partially edentulous patient. The method involves adapting a virtual teeth setup to the intra-oral anatomical situation of the patient, wherein said virtual teeth setup comprises separated surface meshes of individual teeth positioned in a dental arch or segment thereof. The virtual teeth setup is adapted by optimizing an energy function, which represents a quality measure for said virtual teeth setup, and using a statistical model, which describes for a given dentition or segment thereof a probability distribution for at least the shapes of individual teeth, the relations between shapes of neighbouring teeth and/or relations between positions and/or orientations of neighbouring teeth. The adapted virtual teeth setup resulting from said optimized energy function is subsequently used to estimate said sought for shape, position and/or orientation.
