We show how to sequentially optimize magnetic resonance imaging measurement designs over stacks of neighbouring image slices, by performing convex variational inference on a large scale non-Gaussian linear dynamical system, tracking dominating directions of posterior covariance without imposing any factorization constraints. Our approach can be scaled up to high-resolution images by reductions to numerical mathematics primitives and parallelization on several levels. In a first study, designs are found that improve significantly on others chosen independently for each slice or drawn at random.
Jean-Philippe Thiran, Friedhelm Christoph Hummel, Tobias Kober, Tom Hilbert, Erick Jorge Canales Rodriguez, Gabriel Girard, Elda Fischi Gomez, Marco Pizzolato, Gian Franco Piredda, Thomas Yu, Takuya Morishita, Elena Beanato, Alessandro Daducci, Maximilian Jonas Wessel, Chang-Hyun Park, Philipp Johannes Koch, Andéol Geoffroy Cadic-Melchior, Julia Brügger