Inter-subject non-rigid deformations are modeled as elastic waves. The underlying transformation is governed by elastodynamics wave equation. The proposed dynamic elasticity model recovers both global shape differences and large local deformations. The proposed registration scheme pursues in a step-by-step approach thus avoiding local optimal transformations.