We introduce a parameter estimation framework for automatically and robustly personalizing aortic haemodynamic computations from four-dimensional magnetic resonance imaging data. The framework is based on a reduced-order multiscale fluid-structure interaction blood flow model, and on two calibration procedures. First, Windkessel parameters of the outlet boundary conditions are personalized by solving a system of nonlinear equations. Second, the regional mechanical wall properties of the aorta are personalized by employing a nonlinear least-squares minimization method. The two calibration procedures are run sequentially and iteratively until both procedures have converged. The parameter estimation framework was successfully evaluated on 15 datasets from patients with aortic valve disease. On average, only 1.27 ± 0.96 and 7.07 ± 1.44 iterations were required to personalize the outlet boundary conditions and the regional mechanical wall properties, respectively. Overall, the computational model was in close agreement with the clinical measurements used as objectives (pressures, flow rates, cross-sectional areas), with a maximum error of less than 1%. Given its level of automation, robustness and the short execution time (6.2 ± 1.2 min on a standard hardware configuration), the framework is potentially well suited for a clinical setting.