Model Calibration via Deformation

Dynamical computer models often exhibit space-time features that are partially misaligned or misshapen when compared to observational data. Whether due to approximate numerical schemes, incomplete physics, or estimated boundary conditions, the goal of calibrating these models to field data involves optimally aligning model output with observed features. The traditional approach to correcting model discrepancy is to introduce an additive and/or multiplicative bias. Especially for dynamical models, systematic bias may alternatively be viewed as deformation bias. We introduce an expanded approach to model calibration in the presence of space-time feature discrepancy. Borrowing ideas from the image warping literature, we propose a nonlinear deformation of the computer model that optimally aligns with observed images; probabilistically this manifests as a transformation of model coordinate space with a variational penalty on the likelihood function. We apply the approach to a dynamical magnetosphere-ionosphere...

[1]  A. O'Hagan,et al.  Bayesian calibration of computer models , 2001 .

[2]  Richard L. Smith,et al.  Spatial Forecast Verification: Image Warping , 2010 .

[3]  Finn Lindgren,et al.  Analyzing the Image Warp Forecast Verification Method on Precipitation Fields from the ICP , 2010 .

[4]  Stephan R. Sain,et al.  Fast Sequential Computer Model Calibration of Large Nonstationary Spatial-Temporal Processes , 2013, Technometrics.

[5]  G. Shaddick,et al.  Modeling Nonstationary Processes Through Dimension Expansion , 2010, 1011.2553.

[6]  A. O'Hagan,et al.  Bayesian emulation of complex multi-output and dynamic computer models , 2010 .

[7]  C. Goodall.  Procrustes methods in the statistical analysis of shape , 1991 .

[8]  Paul Dupuis,et al.  Variational problems on ows of di eomorphisms for image matching , 1998 .

[9]  P. Guttorp,et al.  Studies in the history of probability and statistics XLIX On the Matérn correlation family , 2006 .

[10]  Y. Amit,et al.  Towards a coherent statistical framework for dense deformable template estimation , 2007 .

[11]  J. Fedder,et al.  Modeling seasonal variations of auroral particle precipitation in a global-scale magnetosphere-ionosphere simulation , 2009 .

[12]  B. Silverman,et al.  Nonparametric Regression and Generalized Linear Models: A roughness penalty approach , 1993 .

[13]  Finn Lindgren,et al.  An image warping approach to spatio-temporal modelling , 2005 .

[14]  Karl Rohr,et al.  Radial basis functions with compact support for elastic registration of medical images , 2001, Image Vis. Comput..

[15]  J. Rougier.  Efficient Emulators for Multivariate Deterministic Functions , 2008 .

[16]  James O. Berger,et al.  A Framework for Validation of Computer Models , 2007, Technometrics.

[17]  A. Gelfand,et al.  Gaussian predictive process models for large spatial data sets , 2008, Journal of the Royal Statistical Society. Series B, Statistical methodology.

[18]  Mevin B. Hooten,et al.  Assessing First-Order Emulator Inference for Physical Parameters in Nonlinear Mechanistic Models , 2011 .

[19]  Yali Amit,et al.  A Nonlinear Variational Problem for Image Matching , 1994, SIAM J. Sci. Comput..

[20]  D. Cox.  Nonparametric Regression and Generalized Linear Models: A roughness penalty approach , 1993 .

[21]  A. O'Hagan,et al.  Gaussian process emulation of dynamic computer codes , 2009 .

[22]  Grace Wahba,et al.  Spline Models for Observational Data , 1990 .

[23]  K. Mardia,et al.  A penalized likelihood approach to image warping , 2001 .

[24]  Timothy F. Cootes,et al.  On representing edge structure for model matching , 2001, Proceedings of the 2001 IEEE Computer Society Conference on Computer Vision and Pattern Recognition. CVPR 2001.

[25]  John Lyon,et al.  The Lyon-Fedder-Mobarry (LFM) global MHD magnetospheric simulation code , 2004 .

[26]  Hans-Peter Meinzer,et al.  Statistical shape models for 3D medical image segmentation: A review , 2009, Medical Image Anal..

[27]  Stephan R. Sain,et al.  Emulating and calibrating the multiple‐fidelity Lyon–Fedder–Mobarry magnetosphere–ionosphere coupled computer model , 2015 .

[28]  Michael A. West,et al.  A dynamic modelling strategy for Bayesian computer model emulation , 2009 .

[29]  D. Higdon,et al.  Computer Model Calibration Using High-Dimensional Output , 2008 .

[30]  Peter Z. G. Qian,et al.  Bayesian Hierarchical Modeling for Integrating Low-Accuracy and High-Accuracy Experiments , 2008, Technometrics.

[31]  J. Lyon,et al.  Effects of the low-latitude ionospheric boundary condition on the global magnetosphere , 2010 .

[32]  Fred L. Bookstein,et al.  Principal Warps: Thin-Plate Splines and the Decomposition of Deformations , 1989, IEEE Trans. Pattern Anal. Mach. Intell..

[33]  Matthew J. Heaton,et al.  Parameter tuning for a multi-fidelity dynamical model of the magnetosphere , 2013, 1303.6992.

[34]  Hermann Ney,et al.  Deformation Models for Image Recognition , 2007, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[35]  Nur Arad,et al.  Image Warping Using Few Anchor Points and Radial Functions , 1995, Comput. Graph. Forum.

[36]  P. Guttorp,et al.  Nonparametric Estimation of Nonstationary Spatial Covariance Structure , 1992 .

[37]  B. Brown,et al.  Forecasts of Spatial Fields , 2012 .